CN111708094B - 一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统 - Google Patents
一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统 Download PDFInfo
- Publication number
- CN111708094B CN111708094B CN202010696996.XA CN202010696996A CN111708094B CN 111708094 B CN111708094 B CN 111708094B CN 202010696996 A CN202010696996 A CN 202010696996A CN 111708094 B CN111708094 B CN 111708094B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- data
- depth
- point
- wide
- 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
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000004576 sand Substances 0.000 title claims abstract description 40
- 229910052770 Uranium Inorganic materials 0.000 title claims abstract description 35
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 title claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 238000012937 correction Methods 0.000 claims abstract description 23
- 238000002386 leaching Methods 0.000 claims abstract description 9
- 238000005070 sampling Methods 0.000 claims description 20
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 11
- 229910052802 copper Inorganic materials 0.000 claims description 11
- 239000010949 copper Substances 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 230000003068 static effect Effects 0.000 claims description 7
- 230000005684 electric field Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 4
- 230000007480 spreading Effects 0.000 claims description 4
- 230000002087 whitening effect Effects 0.000 claims description 4
- 230000005288 electromagnetic effect Effects 0.000 claims description 3
- 238000011065 in-situ storage Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 4
- 239000002689 soil Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 3
- 229910052782 aluminium Inorganic materials 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000009933 burial Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000010287 polarization Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000033558 biomineral tissue development Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000010291 electrical method Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统，所述识别方法包括：采集地浸砂岩型铀矿沉积盆地中沿剖面不同测点处电压值的时间序列；将所述不同测点处电压值的时间序列转换为视电阻率；对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面；对所述分布断面进行深度的标定和校正；基于标定和校正后的分布断面对砂体进行划分。本发明中的上述方法能够提高适用性和探测精度。
Description
技术领域
本发明涉及地球物理勘探领域，特别是涉及一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统。
背景技术
传统的地浸砂岩型铀矿勘探常采用音频大地电磁法、大地电磁法或可控源音频大地电磁法探测深部成矿环境，识别与成矿相关的有利砂体。该类方法通过接收来自天然场或人工发射的不同频率的电磁场，从而反映地下岩石、地层电性结构及特征。
在实际方法应用中，由于上述几种传统电磁测深方法在信噪比、频点密度以及视电阻率计算上的限值，导致在识别电性差异小、厚度较薄的泥砂岩时效果难以达到勘探所需的精度。为了进一步提高深部勘探效率，使地球物理测深方法在砂岩型铀矿中发挥更大作用，一种广域电磁法探测技术被提出。
发明内容
本发明的目的是提供一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统，提高适用性和探测精度。
为实现上述目的，本发明提供了如下方案：
一种基于广域电磁法的砂岩型铀矿砂体识别方法，所述识别方法包括：
采集地浸砂岩型铀矿沉积盆地中沿剖面不同测点处电压值的时间序列；
将所述不同测点处电压值的时间序列转换为视电阻率；
对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面；
对所述分布断面进行深度的标定和校正；
基于标定和校正后的分布断面对砂体进行划分。
可选的，所述采集地浸砂岩型铀矿沉积盆地中沿剖面不同测点处电压值的时间序列具体包括：
确定采样点；
确定广域电磁仪接收端频率；
布设广域电磁仪发射极；
在所述采样点上采用广域电磁仪的铜棒或不极化电极采集电压值的时间序列。
可选的，所述将所述不同测点处电压值的时间序列转换为视电阻率具体包括：
确定单点不同频率的电阻率和相位数据；
删除所述单点不同频率的电阻率和相位数据中的各频点叠加前产生的“飞点”数据；
获取测点点位数据、发射电流数据；
基于所述测点点位数据、发射电流数据以及删除“飞点”数据后的单点不同频率的电阻率和相位数据确定视电阻率。
可选的，所述确定单点不同频率的电阻率和相位数据采用以下公式：
其中，ρa为测点处视电阻率，
为装置系数；
为测点处电压差，由测点处电场强度Ex和接收端电极距MN相乘得到，I为供电电流；
为电磁效应函数，r为收发距，-ik＝(-1-i)/δ，其中δ为均匀半空间的波数。
可选的，所述对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面具体包括：
删除或平移视电阻率中不连续的频点、频段；
对整条剖面上的各测点开展静位移校正，去除浅部不均匀体对单点整支数据的干扰；
对静位移校正后的视电阻率进行反演计算，得到以深度为纵轴、水平距离为横轴的电阻率剖面数据；
将电阻率剖面数据中的深度转换为标高，得到以标高为纵轴，剖面距离为横轴的电阻率剖面数据；
将所述以标高为纵轴，剖面距离为横轴的电阻率剖面数据进行插值；
将插值后的电阻率剖面数据进行白化，得到电阻率随深度的分布断面图。
可选的，对所述分布断面进行深度的标定和校正具体包括：
根据岩心编录或电阻率测井曲线结果中对应砂岩、砂砾岩、砾岩层的顶底界面深度，采用线性拟合方法对反演结果中电阻率剖面数据中的深度进行标定和校准。
可选的，基于标定和校正后的分布断面对砂体进行划分具体包括：
根据层位的横向连续性和相邻剖面的一致性等特征，划分整个剖面上砂体的展布形态。
本发明另外提供一种基于广域电磁法的砂岩型铀矿砂体识别系统，所述系统包括：
数据采集模块，用于采集地浸砂岩型铀矿沉积盆地中沿剖面不同测点处电压值的时间序列；
转换模块，用于将所述不同测点处电压值的时间序列转换为视电阻率；
反演计算模块，用于对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面；
标定和校正模块，用于对所述分布断面进行深度的标定和校正；
划分模块，用于基于标定和校正后的分布断面对砂体进行划分。
可选的，所述数据采集模块具体包括：
采样点确定单元，用于确定采样点；
接收端频率调节单元，用于确定广域电磁仪接收端频率；
发射极布设单元，用于布设广域电磁仪发射极；
电压值的时间序列采集单元，用于在所述采样点上采用广域电磁仪的铜棒或不极化电极采集电压值的时间序列。
可选的，所述转换模块具体包括：
电阻率和相位数据确定单元，用于确定单点不同频率的电阻率和相位数据；
“飞点”数据删除单元，用于删除所述单点不同频率的电阻率和相位数据中的各频点叠加前产生的“飞点”数据；
点位数据、发射电流数据获取单元，用于获取测点点位数据、发射电流数据；
视电阻率确定单元，用于基于所述测点点位数据、发射电流数据以及删除“飞点”数据后的单点不同频率的电阻率和相位数据确定视电阻率。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明的一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统，采用广域电磁数据采集系统在砂岩型铀矿盆地上进行数据采集、处理、反演工作，基于广域电磁法在接收数据频率密度、信噪比以及视电阻率计算公式上所具有的优势，以及利用已知钻孔进行先验信息校正技术，圈定与成矿相关的有利砂体，从而降低勘探风险，提高铀矿找矿的经济效益。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例一种基于广域电磁法的砂岩型铀矿砂体识别方法流程图；
图2为本发明实施例利用钻孔揭露结果对广域电磁结果记性深度校正示意图；
图3为本发明实施例一种基于广域电磁法的砂岩型铀矿砂体识别系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统，提高适用性和探测精度。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例一种基于广域电磁法的砂岩型铀矿砂体识别方法流程图，如图1所示，所述方法包括：
步骤101：采集地浸砂岩型铀矿沉积盆地中沿剖面不同测点处电压值的时间序列。
具体的，步骤101中共包含以下几点：
步骤1011：确定采样点。
其中，电场采集间隔在同一剖面上为100米连续采样，即每隔100米作为一个采样点。各采样点上采用铜棒或不极化电极与地面土壤耦合的方式采集信号。各采样点处采用铜棒采样时直接将铜棒尽量深的插入地下，当用电阻表测量两个铜棒之间接地电阻大于5000Ω时可以采用盐水浇注方式降低接地电阻，保持其与土壤较好的耦合状态。采用不极化电极时先挖深10cm以上的小坑，保持坑内底部有松软湿润的土壤，并去除根、石子等杂物，将不极化电极移入坑内，使底部与土壤充分接触，设置好后用土壤覆盖极坑。
步骤1012：确定广域电磁仪接收端频率。
即，设置接收端频率组合。由于砂体展布与垂向分辨能力相关，因此在接收时采用尽量密的频率密度。采集时首先根据勘探任务选取最低频率数值。其原则与常规大地电磁法相同，具体操作步骤为：
a)明确目标体最大埋深，选择比该最大埋设略大的数值作为勘探深度d；
b)了解测区平均电阻率
计算方法可根据经验或收集测区周边地区钻孔电阻率测井曲线、电法测量断面结果估计；
c)根据频率域电磁测深法勘探深度估算公式
计算得到最低频率f。
然后选择所需最低频率至设备能够发送到的最高频率范围内的所有频组，从而获取测量频率密度最高的单点数据。
步骤1013：布设广域电磁仪发射极。
发射偶极位置至最近接收点的距离要求大于勘探深度的3倍以上，实际布设时采用的5倍以上收发距。发射极布设时，为保证良好的供电效果，一般要求尽量降低接地电阻。采用铝板布设，单个铝板面积约1㎡，阵列式排布，两偶极直接间距约1～2km。
步骤1014：在所述采样点上采用广域电磁仪的铜棒或不极化电极采集电压值的时间序列。
选择JSDY-180广域电磁系统，该发射系统具有大功率、发射伪随机波的特征。本次实例中发射端电流一般大于150A。选用的接收机应具备解码伪随机信号的高灵敏度系统，本次实例中所采用的接收机为JSGY-2广域电磁接收系统。
在此，为了保证记录的各接收道电场值稳定，通过采集软件实时关注各道各频点数据，保证多次叠加时偏差不大于10％。
步骤102：将所述不同测点处电压值的时间序列转换为视电阻率。
具体包括：
步骤1021：确定单点不同频率的电阻率和相位数据。
即，合成单点数据，同一测点采集的不同频组数据文件加载后在软件中操作进行合并，从而获得单点不同频率的电阻率、相位数据，具体采用以下公式：
其中，ρa为测点处视电阻率，
为装置系数；
为测点处电压差，由测点处电场强度Ex和接收端电极距MN相乘得到，I为供电电流；
为电磁效应函数，r为收发距，-ik＝(-1-i)/δ，其中δ为均匀半空间的波数。
步骤1022：删除所述单点不同频率的电阻率和相位数据中的各频点叠加前产生的“飞点”数据。
即，将步骤1021合成的数据在“点数据处理”中加载，查看、删除各频点叠加前产生的“飞点”数据。
步骤1023：获取测点点位数据、发射电流数据。
步骤1024：基于所述测点点位数据、发射电流数据以及删除“飞点”数据后的单点不同频率的电阻率和相位数据确定视电阻率。
具体的，是在“视电阻率计算”中分别导入RTK采集的测点点位数据、发射电流数据以及②采集的电场数据，然后点击“开始计算”，将计算结果导出。
步骤103：对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面。
利用步骤:102计算的视电阻率数据，载入到WFEM2D_Pro反演软件中进行反演模拟计算，从而获得剖面下方电阻率随深度的分布断面。采用该系统对广域电磁二维数据进行删除“飞点”、校正、反演、成图，其处理流程与原则常用的大地电磁结果相似，具体步骤为：
步骤1031：删除或平移单点曲线数据中不连续的频点、频段。
步骤1032：然后利用软件的校正功能，对整条剖面上的各测点开展静位移校正，去除浅部不均匀体对单点整支数据的干扰。
步骤1033：反演计算是在参数设置完成后由计算机自动运行完成的。其中反演过程中采用的主要参数为：水平分辨指数为0.3～0.5；垂向分辨指数上由于沉积地层中期望得到更多分层信息，因此一般选用不大于1的值；反演迭代次数一般大于5次；反演终止拟合差不易大于10％。本次实例采用的分辨率系数分别为水平方向0.5，垂直方向0.85。将上述反演参数在WFEM2D_Pro软件中设置好，点击运行二维反演即可自动进行迭代运算，完成反演过程。
步骤1034：经过迭代反演计算后将得到以深度为纵轴、水平距离为横轴的电阻率剖面数据，该数据格式为“水平距离/米深度/米电阻率/欧姆米”。
步骤1035：利用野外RTK记录的高程数据，将步骤1034中数据结果的深度转换为标高。
步骤1036：以纵轴为标高，横轴为剖面距离，对反演结果模型文件用Surfer软件进行插值。然后再次利用RTK记录的高程数据对插值后数据进行白化，获得带地形的电阻率分布特征断面图。
步骤104：对所述分布断面进行深度的标定和校正。
通过资料收集或自测方式获取一口以上的、经过剖面或剖面附近的钻孔，根据其岩心编录或电阻率测井曲线结果中对应砂岩、砂砾岩、砾岩层的位置，对反演结果中的中高阻电性层位进行标定和校准。本次实例中主要的砂岩层顶底界面位置与步骤103中的1036步骤获得电阻率分布特征断面图中10～30欧姆米的电性层位范围对应。记录砂岩层顶底界面深度[d]＝(d1、d2、d3、d4…dn)，记录断面反演电阻率层位顶底界面深度[D]＝(D1、D2、D3、D4…Dn)(如图2所示)。当断面层数与实际编录结果不相符时，可将编录结果中的若干薄层看做一个合并层，并记录合并层顶底面信息)。将上述两组数据[d]、[D]进行最小二乘拟合。该步骤通过编写MATLAB代码实现，具体代码如下：
d＝[d1,d2,d3,d4,…,dn]；％定义钻孔揭露砂岩顶底界面
D＝[D1,D2,D3,D4,…,Dn]；％定义反演断面推断砂岩顶底界面
[a,b]＝polyfit(d,D,1)；％获得反演结果深度方程系数a,b
对获得的拟合方程y＝ax+b，将步骤103中的1034步骤获得反演结果数据中“深度”一列数据代入x中进行转换，用获得转换后的y数据替换原有深度数据。再重复步骤1035、1036过程，重新绘制剖面电阻率断面图。
步骤105：基于标定和校正后的分布断面对砂体进行划分。
利用步骤104获取的电阻率断面结果和解译原则，通过沿等值线变化较快区域追索的方式划分主要砂体顶底界面，并根据沉积地区层位的横向连续性和相邻剖面的一致性等特征，划分整个测区所有剖面上中高阻层位展布形态，并将其标记为砂岩型铀矿成矿有利砂体，最终从剖面上获取所需的找矿信息，主要包括砂体厚度、砂体顶底板埋深等。
图3为本发明实施例基于广域电磁法的砂岩型铀矿砂体识别系统结构示意图，如图3所示，所述系统包括：
数据采集模块201，用于采集地浸砂岩型铀矿沉积盆地中沿剖面不同测点处电压值的时间序列；
转换模块202，用于将所述不同测点处电压值的时间序列转换为视电阻率；
反演计算模块203，用于对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面；
标定和校正模块204，用于对所述分布断面进行深度的标定和校正；
划分模块205，用于基于标定和校正后的分布断面对砂体进行划分。
其中，所述数据采集模块201具体包括：
采样点确定单元2011，用于确定采样点；
接收端频率调节单元2012，用于确定广域电磁仪接收端频率；
发射极布设单元2013，用于布设广域电磁仪发射极；
电压值的时间序列采集单元2014，用于在所述采样点上采用广域电磁仪的铜棒或不极化电极采集电压值的时间序列。
所述转换模块202具体包括：
电阻率和相位数据确定单元2021，用于确定单点不同频率的电阻率和相位数据；
“飞点”数据删除单元2022，用于删除所述单点不同频率的电阻率和相位数据中的各频点叠加前产生的“飞点”数据；
点位数据、发射电流数据获取单元2023，用于获取测点点位数据、发射电流数据；
视电阻率确定单元2024，用于基于所述测点点位数据、发射电流数据以及删除“飞点”数据后的单点不同频率的电阻率和相位数据确定视电阻率。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种基于广域电磁法的砂岩型铀矿砂体识别方法，其特征在于，所述识别方法包括：
采集地浸砂岩型铀矿沉积盆地中沿剖面不同测点处电压值的时间序列；
将所述不同测点处电压值的时间序列转换为视电阻率；
对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面；
对所述分布断面进行深度的标定和校正；
基于标定和校正后的分布断面对砂体进行划分；
对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面具体包括：
删除或平移视电阻率中不连续的频点、频段；
对整条剖面上的各测点开展静位移校正，去除浅部不均匀体对单点整支数据的干扰；
对静位移校正后的视电阻率进行反演计算，得到以深度为纵轴、水平距离为横轴的电阻率剖面数据；
将电阻率剖面数据中的深度转换为标高，得到以标高为纵轴，剖面距离为横轴的电阻率剖面数据；
将所述以标高为纵轴，剖面距离为横轴的电阻率剖面数据进行插值；
将插值后的电阻率剖面数据进行白化，得到电阻率随深度的分布断面图；
对所述分布断面进行深度的标定和校正具体包括：
根据岩心编录或电阻率测井曲线结果中对应砂岩、砂砾岩、砾岩层的顶底界面深度，采用线性拟合方法对反演结果中电阻率剖面数据中的深度进行标定和校准。
2.根据权利要求1所述的基于广域电磁法的砂岩型铀矿砂体识别方法，其特征在于，所述采集地浸砂岩型铀矿沉积盆地中沿剖面不同测点处电压值的时间序列具体包括：
确定采样点；
确定广域电磁仪接收端频率；
布设广域电磁仪发射极；
在所述采样点上采用广域电磁仪的铜棒或不极化电极采集电压值的时间序列。
3.根据权利要求1所述的基于广域电磁法的砂岩型铀矿砂体识别方法，其特征在于，所述将所述不同测点处电压值的时间序列转换为视电阻率具体包括：
确定单点不同频率的电阻率和相位数据；
删除所述单点不同频率的电阻率和相位数据中的各频点叠加前产生的“飞点”数据；
获取测点点位数据、发射电流数据；
基于所述测点点位数据、发射电流数据以及删除“飞点”数据后的单点不同频率的电阻率和相位数据确定视电阻率。
4.根据权利要求3所述的基于广域电磁法的砂岩型铀矿砂体识别方法，其特征在于，所述确定单点不同频率的电阻率和相位数据采用以下公式：
其中，ρa为测点处视电阻率，
为装置系数；
为测点处电压差，由测点处电场强度Ex和接收端电极距MN相乘得到，I为供电电流；
为电磁效应函数，r为收发距，-ik＝(-1-i)/δ，其中，k为均匀半空间的波数。
5.根据权利要求1所述的基于广域电磁法的砂岩型铀矿砂体识别方法，其特征在于，基于标定和校正后的分布断面对砂体进行划分具体包括：
根据层位的横向连续性和相邻剖面的一致性特征，划分整个剖面上砂体的展布形态。
6.一种基于广域电磁法的砂岩型铀矿砂体识别系统，其特征在于，所述系统包括：
数据采集模块，用于采集地浸砂岩型铀矿沉积盆地中沿剖面不同测点处电压值的时间序列；
转换模块，用于将所述不同测点处电压值的时间序列转换为视电阻率；
反演计算模块，用于对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面；
标定和校正模块，用于对所述分布断面进行深度的标定和校正；
划分模块，用于基于标定和校正后的分布断面对砂体进行划分；
对所述视电阻率进行反演计算，得到剖面下方电阻率随深度的分布断面具体包括：
删除或平移视电阻率中不连续的频点、频段；
对整条剖面上的各测点开展静位移校正，去除浅部不均匀体对单点整支数据的干扰；
对静位移校正后的视电阻率进行反演计算，得到以深度为纵轴、水平距离为横轴的电阻率剖面数据；
将电阻率剖面数据中的深度转换为标高，得到以标高为纵轴，剖面距离为横轴的电阻率剖面数据；
将所述以标高为纵轴，剖面距离为横轴的电阻率剖面数据进行插值；
将插值后的电阻率剖面数据进行白化，得到电阻率随深度的分布断面图；
对所述分布断面进行深度的标定和校正具体包括：
根据岩心编录或电阻率测井曲线结果中对应砂岩、砂砾岩、砾岩层的顶底界面深度，采用线性拟合方法对反演结果中电阻率剖面数据中的深度进行标定和校准。
7.根据权利要求6所述的基于广域电磁法的砂岩型铀矿砂体识别系统，其特征在于，所述数据采集模块具体包括：
采样点确定单元，用于确定采样点；
接收端频率调节单元，用于确定广域电磁仪接收端频率；
发射极布设单元，用于布设广域电磁仪发射极；
电压值的时间序列采集单元，用于在所述采样点上采用广域电磁仪的铜棒或不极化电极采集电压值的时间序列。
8.根据权利要求6所述的基于广域电磁法的砂岩型铀矿砂体识别系统，其特征在于，所述转换模块具体包括：
电阻率和相位数据确定单元，用于确定单点不同频率的电阻率和相位数据；
“飞点”数据删除单元，用于删除所述单点不同频率的电阻率和相位数据中的各频点叠加前产生的“飞点”数据；
点位数据、发射电流数据获取单元，用于获取测点点位数据、发射电流数据；
视电阻率确定单元，用于基于所述测点点位数据、发射电流数据以及删除“飞点”数据后的单点不同频率的电阻率和相位数据确定视电阻率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010696996.XA CN111708094B (zh) | 2020-07-20 | 2020-07-20 | 一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010696996.XA CN111708094B (zh) | 2020-07-20 | 2020-07-20 | 一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111708094A CN111708094A (zh) | 2020-09-25 |
| CN111708094B true CN111708094B (zh) | 2023-04-18 |
Family
ID=72546748
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010696996.XA Active CN111708094B (zh) | 2020-07-20 | 2020-07-20 | 一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111708094B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112859181B (zh) * | 2021-01-08 | 2024-04-05 | 四川省自然资源投资集团物探勘查院有限公司 | 一种广域电磁法局部异常分离方法 |
| CN113960674B (zh) * | 2021-10-14 | 2023-11-21 | 湖北省水文地质工程地质勘察院有限公司 | 一种广域电磁法二维反演方法 |
| CN113835126B (zh) * | 2021-11-02 | 2024-10-22 | 新疆维吾尔自治区煤田灭火工程局 | 一种温度自校正自然电位测量系统及方法 |
| CN114415245A (zh) * | 2022-01-21 | 2022-04-29 | 山东省地质矿产勘查开发局第六地质大队(山东省第六地质矿产勘查院) | 一种适用于中深部蚀变岩型金矿找矿方法 |
| CN114609683B (zh) * | 2022-03-15 | 2024-07-12 | 核工业北京地质研究院 | 一种砂岩型铀矿砂体有利区确定方法及系统 |
| CN114966864B (zh) * | 2022-06-24 | 2024-03-08 | 山东省地质矿产勘查开发局第三地质大队(山东省第三地质矿产勘查院、山东省海洋地质勘查院) | 一种变质岩区晶质石墨矿找矿系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109782360A (zh) * | 2019-02-19 | 2019-05-21 | 核工业航测遥感中心 | 一种低阻覆盖区深部砂体探测方法 |
| CN111045114A (zh) * | 2019-12-31 | 2020-04-21 | 核工业北京地质研究院 | 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7262602B2 (en) * | 2005-07-08 | 2007-08-28 | David Gary Meyer | Active geophysical prospecting system |
-
2020
- 2020-07-20 CN CN202010696996.XA patent/CN111708094B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109782360A (zh) * | 2019-02-19 | 2019-05-21 | 核工业航测遥感中心 | 一种低阻覆盖区深部砂体探测方法 |
| CN111045114A (zh) * | 2019-12-31 | 2020-04-21 | 核工业北京地质研究院 | 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法 |
Non-Patent Citations (7)
| Title |
|---|
| I. C. ROACH, S. JAIRETH AND M. T. COSTELLOE.Applying regional airborne electromagnetic (AEM) surveying to understand the architecture of sandstonehosted uranium mineral systems in the Callabonna Sub-basin, Lake Frome region, South Australia.Australian Journal of Earth Sciences.2014,第61卷全文. * |
| 何继善 ; .广域电磁测深法研究.中南大学学报(自然科学版).2010,(第03期),全文. * |
| 何继善.广域电磁法和伪随机信号电法.北京：高等教育出版社,2010,(第ISBN 978-7-04-028704-2版),全文. * |
| 张叶鹏 ; 谢亮 ; 李艳 ; 严家斌 ; 王红 ; 孔辉 ; .基于全区视电阻率的CSAMT应用研究――以西藏帮浦―笛给铅锌矿区为例.矿产与地质.2020,(第03期),全文. * |
| 柳建新,赵然,郭振威.电磁法在金属矿勘查中的研究进展.地球物理学进展.2019,第34卷(第1期),全文. * |
| 王会波 ; 宋振涛 ; 李新站 ; 王均 ; 王生云 ; .电磁法物探在北方中新生代盆地产铀性评价中试验研究.地质找矿论丛.2018,(第03期),全文. * |
| 王延忠 ; 沈金松 ; 周正武 ; 符超 ; 路肖尧 ; .井震约束的广域电磁数据反演及其在保靖页岩气储层评价中的应用.工程地球物理学报.2016,(第03期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111708094A (zh) | 2020-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111708094B (zh) | 一种基于广域电磁法的砂岩型铀矿砂体识别方法及系统 | |
| Sørensen et al. | An integrated approach for hydrogeophysical investigations: new technologies and a case history | |
| Dlugosch et al. | Assessment of the potential of a new generation of surface nuclear magnetic resonance instruments | |
| CN110187398B (zh) | 一种寻找井间目标体的多电极系探测方法 | |
| Aizebeokhai et al. | The use of the multiple-gradient array for geoelectrical resistivity and induced polarization imaging | |
| CN103913776A (zh) | 盾构施工中球状孤石的探测方法 | |
| Legault et al. | Synthetic model testing and distributed acquisition dc resistivity results over an unconformity uranium target from the Athabasca Basin, northern Saskatchewan | |
| Saey et al. | Electrical conductivity depth modelling with a multireceiver EMI sensor for prospecting archaeological features | |
| CN102419452A (zh) | 电性源瞬变电磁法全场区探测方法 | |
| Cubbage et al. | A modified Wenner array for efficient use of eight-channel resistivity meters | |
| Kuras et al. | Capacitive resistivity imaging with towed arrays | |
| Zhe et al. | Multichannel, full waveform and flexible electrode combination resistivity-imaging system | |
| Khalil et al. | 2D and 3D resistivity inversion of Schlumberger vertical electrical soundings in Wadi El Natrun, Egypt: A case study | |
| CN110187394A (zh) | 双场源电磁测深法获取地层电阻率各向异性的方法及装置 | |
| Tao et al. | Evaluating the joint use of GPR and ERT on mapping shallow subsurface features of karst critical zone in southwest China | |
| Asry et al. | Groundwater exploration using 2-D geoelectrical resistivity imaging technique at Sungai. Udang, Melaka | |
| CN114764151A (zh) | 一种大地电磁分频层析反演方法 | |
| CN112485829A (zh) | 一种电阻率法计算动态含水量方法 | |
| Kirsch et al. | Geoelectrical methods | |
| CN114280673B (zh) | 基于静态效应的物探仪及数据采集方法 | |
| Kelly et al. | Continuous electrical imaging for mapping aquifer recharge along reaches of the Namoi River in Australia | |
| Stummer | New developments in electrical resistivity imaging | |
| Oldenburg et al. | Direct current resistivity methods | |
| Sørensen et al. | Pulled Array Continuous Electrical Sounding Pa-Cves | |
| CN112882124A (zh) | 一种三维铀成矿有利部位的圈定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |