自相关滤波方法及其在山东齐河—禹城铁矿勘查应用
Application of autocorrelation filtering to iron ore prospecting in Qihe-Yucheng area, Shandong Province
通讯作者: 于长春(1964-),教授级高级工程师,长期从事航空物探方法技术研究和资料解释工作 。Email:bjycc@126.com
责任编辑: 王萌
收稿日期: 2021-08-31 修回日期: 2021-09-26
基金资助: |
|
Received: 2021-08-31 Revised: 2021-09-26
作者简介 About authors
吴成平(1982-),教授级高级工程师,主要从事航空物探数据处理和综合解释工作。Email:
山东齐河—禹城地区属于典型的厚覆盖区,深部隐伏铁矿在地表引起的地球物理异常相对较弱,因此提取铁矿引起的局部弱异常是该地区找矿预测的关键。自相关滤波方法采用一定波长的高通滤波压制区域场的影响,达到提取局部弱异常的目的。首先,通过建立不同规模和深度的铁矿模型和产生背景场的倾斜岩体模型,在此基础上正演获得叠加磁场,对该磁场进行自相关滤波方法处理,有效地提取了铁矿体产生的弱异常。而后,对齐河—禹城地区的实测数据进行处理,处理结果与钻孔资料进行对比,验证了方法的有效性。自相关滤波方法提取弱异常缩小了找矿范围,对圈定找矿有利部位和钻探验证具有重要指导意义。
关键词:
Deeply concealed iron ores cause relatively weak geophysical anomalies on the ground surface in the Qihe-Yucheng area, Shandong Province, a typical thickly covered area. Therefore, it is the key to the ore prospecting in this area to extract local weak anomalies caused by iron ores. Autocorrelation filtering can extract local weak anomalies by reducing the influence of a regional field using a high-pass filter of a certain wavelenth. The steps of the method are as follows. First, establish a model of incline plutons generating background field and models of iron ores of different scales and depths. Based on this, obtain the superimposed magnetic field through forward modeling. Then carry out autocorrelation filtering processing of the superimposed magnetic field to effectively extract the weak anomalies of iron ores. The autocorrelation filtering method was used to process the measured data in the Qihe-Yucheng area and the processing results were compared to the borehole data, verifying that the autocorrelation filtering method is effective. The extraction of weak anomalies using the autocorrelation filtering method allows for the ore prospecting scope to be reduced, which can provide guidance on the determination of prospecting favorable locations and drilling verification.
Keywords:
本文引用格式
吴成平, 于长春, 熊盛青, 徐剑春, 徐东礼, 乔春贵.
WU Cheng-Ping, YU Chang-Chun, XIONG Sheng-Qing, XU Jian-Chun, QIAO Chun-Gui.
0 引言
自相关滤波方法是数字信号处理中非常重要的处理方法,多用于识别有意义的弱异常,如识别矿体、含油气构造体、断裂构造等。该方法在重磁领域的应用仍在不断发展过程中,最早G.D.Garland[1]基于泊松定理,利用相关系数分析研究重磁异常。姜枚等[2]阐述了统计分析方法在区域重磁资料解释中的应用,指出计算相关函数时延迟距离或窗口选择的重要性。秦葆瑚[3]应用自相关滤波方法降低噪声、突出弱异常,对实测磁数据用插值切割法与自相关滤波方法进行对比试验,发现两者获取的剩余异常相似,但后者异常形态受滤波半径影响更小。熊盛青[4]提出自协方差、互协方差为基础的结构逻辑法提取非平稳场中的弱异常,应用于铜多金属矿和金矿勘查。刘天佑[5]采用剖面相关滤波对扬子地台断裂不同延拓高度数据进行处理,取得较好的效果。韩兆红等[6]在测线上自动寻找剩余局部异常最大点处的观测值作为参考模型,改造自相关滤波因子,并处理实测资料,圈定新的找矿远景区。自相关滤波方法在重磁领域应用还不是很多,在深部铁矿勘查中未见应用。山东齐河—禹城厚覆盖区铁矿在地表引起的异常微弱,利用自相关滤波方法对研究区内数据进行局部弱异常提取,取得了较好的效果,为厚覆盖区找矿提供了支撑,拓展了自相关滤波方法的应用领域。
1 原理及公式
其中:Z0为窗口中心点的异常值,即计算点Z(i,j); l=1,2,…,R; k=0,1,…,l; i=1,2,…,M; j=1,2,…,N; M、N表示研究区大小(测线数为M,每条测线的测点数为N)。
计算点Z(i,j)处经过自相关滤波后的异常值Zloc(i,j)为:
2 模型试验
图1
图1
地质模型在xy、xz、yz平面投影
a—模型在xy平面投影;b—模型在xz平面投影;c—模型在yz平面投影
Fig.1
xy,xz and yz plane profiles of geological model
a—xy plane profile of theoretical geological model; b—xz plane profile of theoretical geological model; c—yz plane profile of theoretical geological model;
表1 矿体模型参数
Table 1
编号 | 质心坐标/ (x/m,y/m,z/m) | 长/m | 宽/m | 高/m | 磁化强度/ (10-3A·m-1) |
---|---|---|---|---|---|
① | (1050,1050,-650) | 100 | 100 | 100 | 33600 |
② | (2025,2425,-475) | 50 | 50 | 50 | 33600 |
③ | (2950,2450,-1050) | 100 | 100 | 100 | 33600 |
图2
图2
自相关滤波方法模型试验
a—地面正演ΔT异常;b—自相关滤波方法提取局部弱异常
Fig.2
The theoretical model test of self-correlation method
a—ground forward ΔT field; b—extracted local weak anomaly map by using self auto correlation filting method
3 实际应用
3.1 数据处理
山东齐河—禹城地区位于山东省西北部,地处黄河以北。2011年在山东齐河—禹城地区完成了1∶5万航磁测量,测线方向为SN向;2017年,又进行了1∶5万航磁测量,测线方向与2011年相同,但测线位置穿插于2011年测线,相当于进行了加密测量。以2017年航磁数据为基准,对上述航磁数据进行调平融合处理[10],得到了1∶2.5万航磁数据。数据网格化采用的网格距为50 m×50 m。
在自相关滤波方法处理之前对数据进行了化极处理。根据钻探资料,研究区内铁矿体埋深在700~1500 m,厚度几米至上百米不等。通过分析已知铁矿产生的航磁异常以及对铁矿体产生的异常进行正演模拟,该地区铁矿引起的航磁异常宽度一般为1~2 km,因此设定自相关半径为20个数据点。自相关滤波方法使用的软件为基于中国自然资源航空物探遥感中心GeoProbe Ver 4.0软件平台[11]自主开发的插件。
3.2 应用效果
图3
经自相关滤波处理后,上述3个异常区由岩体引起的背景场基本消除,认为提取的规模和强度相对较大的弱异常是地下磁性体局部凸起或磁性不均的反映,是找矿的重要部位。需要注意的是:①这里指的弱异常是相对于原叠加异常强度而言;②提取的局部弱异常中也有部分异常是地面干扰引起,其波长短、相对杂乱,在解释中应根据异常特征,并结合遥感、实地踏勘等资料加以区分。
本区已施工多个钻孔[16⇓-18],可以看到见矿钻孔均位于提取的局部弱异常范围内(图4),而且局部异常强度相对较大。例如:李屯地区ZK1、ZK0701、ZK5、ZK6、ZK7等见矿钻孔所处的局部异常;潘店地区PZK01见矿钻孔所处的局部异常;以及大张地区ZK001、ZK002、ZK003、DZK01所处的局部异常。而未见矿钻孔多位于强度较弱的局部异常,甚至不在局部异常范围内,如图4中ZK3、ZK4在提取的局部异常边部,ZK2801则不在提取的局部异常范围内。钻孔资料表明,自相关滤波方法提取的局部弱异常,尤其那些强度相对较大的局部弱异常与深部铁矿密切相关,对该地区寻找深部铁矿具有重要意义,可指导后续钻探工作。
图4
图4
自相关滤波方法处理结果及钻孔分布
Fig.4
The result of self-correlation autocorrection filting and borehole distribution
4 结论
1)应用自相关滤波方法对磁异常进行处理,压制了背景场的影响,明显突出了局部弱异常。不仅从叠加磁场中提取出有意义的异常,而且自相关滤波方法提取的局部弱异常,其异常范围比原异常更接近磁性地质体的实际边界。
2)在齐河—禹城厚覆盖区,铁矿埋深大导致在地表观测的铁矿体引起的地球物理信号相对弱,而且铁矿引起的异常湮没在大面积的磁性岩体引起的异常中,铁矿预测的难度大。通过自相关滤波方法提取弱异常显著缩小了找矿范围,为该地区找矿靶区圈定和钻探验证等工作提供了重要依据。
3)自相关滤波方法提取的局部弱异常突出矿致异常等有用信息的同时,也突出了部分地表人文干扰异常,在地质解释过程中建议使用遥感、实地踏勘等资料进行排除。
致谢
山东省煤田地质规划勘察研究院、中化地质矿山总局山东地质勘察院、山东省地调查院等单位及相关技术人员提供了研究区钻孔资料,匿名审稿专家和编辑部对本文提出了宝贵意见和建议,在此表示由衷感谢。
参考文献
Combined analysis of gravity and magnetic anomalies
[J]. ,DOI:10.1190/1.1437650 URL [本文引用: 1]
试谈统计分析方法在区域重磁资料解释中的某些应用
[J]. ,
Some applications of statistical analysis method in regional gravity and magnetic data interpretation
[J]. ,
用自相关滤波法提取弱异常
[J]. ,
Extracting weak anomalies by self-correlation filtering
[J]. ,
航空物探勘查金属矿产时提取弱信息的解释方法
[J]. ,
Interpretation method of extracting weak signals of aero geophysical data for prospecting metallic deposits
[J]. ,DOI:10.3390/geosciences11020067 URL [本文引用: 1]
相关滤波方法及其在扬子地台断裂分析中的应用
[J]. ,
The relative filtering and its application to fault analysis of the yangtze platform
[J]. ,
自相关滤波法提取重磁场中弱异常
[J]. ,
Extracting weak anomaly in gravity and magnetic field with self-correlation filtering method
[J]. ,
位场数据处理中的最小曲率扩边和补空方法研究
[J]. ,
The research to the extending edge and interpolation based on the minimum curvature method inpotential field data processing
[J]. ,
磁(重力)异常局部场与区域场分离处理中的扩边方法研究
[J]. ,
A study of the scheme of extending edge in the processing of separating local field from regional field for magnetic/gravity anomaly
[J]. ,
航磁弱异常区磁场水平调整——以黔东地区为例
[J]. ,
Leveling aeromagnetic data for the weak magnetic field:A case study in eastern Guizhou
[J]. ,
插件技术在GeoProbe地球物理软件平台中的应用
[J]. ,
The application of the plugin technology to geophysical software platform(Geoprobe)
[J]. ,
广东莲花山断裂带南西段整装勘查区矿床预测模型
[J]. ,
Mineral prediction model of the integrative exploration area in the southwest section of the Lianhuashan fault zone, Guangdong Province
[J]. ,DOI:10.1086/624061 URL [本文引用: 1]
山东省齐河——禹城地区矽卡岩型铁矿成矿预测
[J]. ,
Metallogenic prognosis of skarn-type iron ore deposits in Qihe-Yucheng area, Shandong province
[J]. ,
鲁西齐河地区岩(矿)石物性特征及应用
[J]. ,DOI:10.11867/j.issn.1001-8166.2019.10.1099 [本文引用: 1]
为系统掌握鲁西齐河地区岩(矿)石物性特征,服务于基础地质和矿产资源调查,在该地区及周边实地采集了磁性、密度等物性参数,并根据不同岩(矿)石分类进行统计分析。结果表明,磁铁矿的磁性最大,基性侵入岩的磁性大于中酸性侵入岩,而沉积岩的磁性最弱。新生界与下伏地层和岩体的密度差异明显,其界面起伏变化可以引起明显的重力异常。磁铁矿的密度最大,黄铁矿次之。太古界泰山岩群变质岩密度较大,也可引起一定强度的重力异常。通过物性特征与地球物理场之间的联系,综合其他资料进行了岩性填图,圈定变质岩和侵入岩范围,确定了强磁性地质体范围作为找矿靶区,并基于物性特征进行了重磁联合反演,推断矿体的埋深、厚度等信息,与钻探验证结果吻合。研究结果提高了对该地区岩(矿)石物性特征的认识程度,为地球物理资料的反演解释和地质构造、矿产研究等提供更可靠的依据。
Physical characteristics of rocks and ores and their application in Qihe area,Western Shandong
[J]. ,
山东齐河—禹城地区重磁场特征及找矿预测
[J]. ,
Magnetic field characteristics and ore prediction in Qihe-Yucheng area of Shandong Province
[J]. ,
山东齐河—禹城地区李屯富铁矿床地球化学特征及地质意义
[J]. ,
Geochemical characteristics and geological significance of Litun iron ore deposit in Qihe-Yucheng area, Shandong Province
[J]. ,
/
〈 | 〉 |