摘要
钾盐是农作物增产的主要肥料成分,也是国民经济诸多行业不可替代的重要化工原料。全球钾矿储量很不均衡,主要分布在北美、欧洲、中东和亚洲的泰国等15个国家。中国钾资源总体贫乏,已探明钾矿的工业储量仅占世界的0.47%,2010年以后钾盐缺口已增加到1000万吨以上。中国是世界钾盐进口的主要国家之一,占全球进口总量的18.6%,占国内消费量的80~90%。建立起我国强有力的钾肥资源保障体系,以保证农业可持续发展具有深远的战略意义。
由于钾盐资源储量相对不足和钾盐消费的急速增长,我国不得不进入国际钾盐市场,但是世界较有潜力的钾盐矿区基本被几个寡头垄断,所能进入的主要是前苏联和东南亚的个别国家。老挝有丰富的矿藏,其政府制定了一系列地矿勘测、矿藏开发和矿业发展计划及招商项目,寄希望于国际社会的支持,并已引起高度关注。2006年开始国内已经陆续有多家公司与老挝政府签订了钾盐矿勘探合同,并取得了令人瞩目的成就。2008年老挝国家主席朱马里访华期间与我国达成了“资源换资产、全面推进双边经贸合作”的共识,为中国企业率先“登陆”国际舞台提供了政策依据。因此,老挝矿业发展前景看好。
2006~2012年吉林大学地球探测科学与技术学院承担了四项老挝钾盐矿勘探任务。作者有幸三次奔赴现场,全程参加了各个勘探环节,积累了一定的工作经验,选择该项研究具有一定的基础。作者使用2006年在老挝甘蒙他曲实测的数据,并着眼于重力数据处理方法特别是小波分析在寻找钾盐矿中的应用研究,希望能为缓解我国钾盐矿的对外依赖程度做出自己应有的贡献。具体工作与成果如下:
(1)建立了三个不同的重力异常模型,利用二维离散小波/小波包变换分别对它们的区域场和局部场进行分离。试验结果表明:用相同的模型,选择不同的小波进行重力场分离具有良好的效果,不同小波分析之间差异较小;小波包分解比小波分解更精确,但小波分解能够直接确定低频成分,而小波包分解需要建立与实际地质情况更相符的模型来确定低频系数。
(2)对多种规则地质体引起的重力异常进行了区域场和局部场分解,作出了各种小波系在所给模型重力场分解中的精度评价,通过对比找出了最佳小波及其分解尺度。试验结果表明:在小波分析中Coiflets小波系的精度较高,而在小波包分解中Daubechies小波系的精度较高。
(3)基于对理论模型的研究,选择相应的二维小波包对老挝钾盐勘探区的实测重力数据进行了区域场和局部异常的分离,并与测井资料做了对比解释。
(4)对研究区重力资料作了解析延拓、方向导数、视密度反演、欧拉反褶积等常规数据处理和三维分布特征分析,并与小波分析的结果作了对比,总结了各方法的优势和局限性,且在数据的二次处理过程中,发现了一个新的找矿靶区。
(5)将钻探结果与地球物理、地球化学、地质等资料进行对比,进一步验证了重力勘探在寻找钾盐矿中的效果最好。
(6)对钾盐矿勘探工作中的重力数据处理方法进行了总结,指出了存在的问题和不足,提出了继续研究的建议。
Potash is the crop yield nutrient, is also the national economy industries irreplaceableimportant chemical raw materials. Global potash reserves is not balanced, mainly in NorthAmerica, Europe, the Middle East and Asia, Thailand and other15countries. China'sPotassium Resources indigence, has proven reserves of potassium mineral industry onlyaccounts for0.47%of the world, since2010, has increased to more than10000000tons ofgap. China is one of the major countries of world potash imports accounted for about18.6%oftotal imports, domestic consumption, accounting for80~90%. Establish our country strongpotash resources guarantee system, in order to ensure the sustainable development ofagriculture has far-reaching strategic significance.
The potash resources reserves is relatively insufficient and potassium consumption ofrapid growth, China has to enter the international potassium market, but the world has morepotential sylvite area basically is a few oligarch monopoly, can enter the main is the formerSoviet Union and the individual countries in Southeast asia. Laos has rich mineral resources,the government has formulated a series of mineral exploration, mineral development andmining industry development planning and investment projects, in the hope that the support ofthe international community, and has attracted a lot of attention. Beginning in2006, has beena domestic companies and the government of Laos signed potash exploration contracts, andobtained the achievement attract people's attention. In2008the Lao President Jumari duringhis visit to China and China reached a "resource for assets, promote bilateral economic andtrade cooperation" consensus, for Chinese enterprises to take the lead in" landing"international stage provides policy basis. Therefore, Lao mining development prospects.
2006-2012Jilin University assumed four Laos sylvite mine exploration mission. Theauthor has favour three went to the scene, participate in the whole of the exploration process,accumulated certain experience, the study has certain base. The authors use in Laos in2006Gammon not measured data, and focus on gravity data processing method, especiallywavelet analysis in search for potash application study, hope to be able to reduce ourdependence on external potassium salt mine to make their due contribution. Specific work andresults are as follows:
(1) The establishment of three different gravity anomaly model, the use oftwo-dimensional discrete wavelet/wavelet packet transform respectively to their regionalfield and local field separation. The test results show that: with the same model, selecting different wavelet gravity field separation with good effect, little difference between differentwavelet analysis; wavelet packet decomposition than wavelet decomposition is more precise,but can directly determine the low frequency component of wavelet decomposition andwavelet packet decomposition, need to build more accord with actual geological conditionsthe model to determine the frequency coefficient of.
(2) To a variety of rules geological body caused by gravity abnormality by regional fieldand local field decomposition, make a variety of wavelet coefficients in the model of gravityfield in the decomposition of accuracy evaluation, through the comparison of the optimalwavelet and scale. The test results show that: the wavelet analysis Coiflets wavelet serieshigh precision, and the wavelet packet decomposition of Daubechies wavelet series highprecision.
(3) Based on the model, corresponding to the selected wavelet packet to Laos sylviteexploration of the measured gravity data of regional field and local anomaly separation, andwith well logging data to do comparative interpretation.
(4) In the study area and the analytic continuation of gravity data, direction derivative,apparent density inversion, Euler deconvolution and conventional data processing andthree-dimensional distribution characteristic analysis and wavelet analysis, and the resultswere compared, summed up the advantages and limitations of each method, and the data ofthe two process, found a new prospecting target area.
(5) The results of drilling and geophysical, geochemical, geological data are compared,further validation of gravity prospecting in searching for the best effect of potash mine.
(6) To potash in the exploration of gravity data processing method were summarized,pointed out the problem of existence and inadequacy, put forward suggestions for continuedresearch.
引文
[1]童鸿频.老挝钾盐勘探开发投资项目的风险分析[J].大众科技,2008,105(5):104-105
[2]吴江,张滨,何永等.中国农资:对话老挝钾盐开发[J].中国农资,2010,02:24
[3]邓尔新.钾盐成矿理论的若干问题[J].地质地球化学,1980,03:P37-50
[4]袁见齐,霍承禹,蔡克勤.高山深盆的成盐环境-一种新的成盐模式的剖析[J].地质评论,1983,29(2):159-164
[5]马万栋.塔里木盆地西部盐岩的地球化学特征及钾盐远景区预测研究,[D].青海:中国科学院研究生院(青海盐湖研究所),2006
[6]张建林,中国-老挝古近纪钾盐矿对比,[J].云南地质,2006,25(3):309-316
[7]续培信,老挝万象钾盐矿开采方法的探讨,[J].中国井矿盐,2005,36:28-31
[8]王建明.重磁异常变换中的延拓问题[J].新疆有色金属,2006,2:11-16
[9]赵学智,陈统坚,叶邦彦等.基于卷积型小波包变换的信号消噪算法[J].数据采集与处理,2003,18(3):292-295
[10]张花花.基于最优小波包基的信号增强算法研究及应用[J].电子测试,2008,8:15-20
[11]方东红,曾昭发,陈家林.基于小波分析的重磁数据求导方法及应用[J].吉林大学学报(地球科学版).2008,38(6):1049-1054
[12]赵学智,陈统坚,叶邦彦等.基于卷积型小波包变换的信号消噪算法[J].数据采集与处理,2003,18(3):292-295
[13]张花花.基于最优小波包基的信号增强算法研究及应用[J].电子测试,2008,8:15-20
[14]汪炳柱,王硕儒.二维位场向上延拓和向下延拓的样条函数法[J].物探化探计算技术,1998,20(2):125-129
[15]栾文贵,熊光楚.重、磁场解析延拓中的一些问题[J].桂林冶金地质学院学报,1986,1:83-93
[16]范美宁.欧拉反褶积方法的研究及应用[D].吉林:吉林大学地球探测科学与技术学院,2006
[17]史辉,刘天佑,Dawi Muna Ghaboush.利用欧拉反褶积法估计二度磁性体深度与位置[J].2005,29(3):230-233
[18]范美宁,孙运生,田庆君.关于欧拉反褶积方法计算中的一点改进[J].物探化探计算技术,2003,8(3):42-46
[19]赵玉岩.区域地球物理地球化学资料在浅覆盖区地质调查中的应用方法研究[D].吉林大学博士论文,2009F
[20]陈高,李仲远,尹正武等.复杂地区高精度综合地球物理勘探技术方法及应用[P].中国石化勘探南方分公司,2006-01~2008-12
[21]管志宁,安玉林,陈维雄.曲线与曲面上磁场向上延拓和分量转换[J].地球物理学报,1985,28(4):420-428
[22]孟军海,陈文义,霍成胜.用窗口平均法进行重力异常划分的可能性[J].青海地质(2001增刊),:48-52
[23]熊光楚.自调节趋势分析法[J].物探与化探,2000,24(4):268-277
[24]文百红,程方道.用于划分磁异常的新方法—插值切割法[J].中南矿冶学院学报,1990,21(3):229-236.
[25]丁科.插值切割法用于低应变反射波法波形曲线的解释[J].勘察科学技术,2001,1:59-63.
[26]张瑞钢,钱家忠,赵卫东,马雷.对应分析法在地下水化学特征分析中的应用[J].合肥工业大学学报(自然科学版),2008,31(10):1552-1555.
[27]邱勇.卡尔曼滤波在重力场数据处理中的应用[D].吉林大学硕士论文,2009.
[28]吴振慈,刘兴德.利用区域布格重力异常解释深部构造问题的讨论[J].物探与化探1985,9(6):426-430
[29]王继伦.物探数据处理中线性滤波的基本矛盾和作用[J].地质与勘探,1979,3:41-48
[30]石磊,郭良辉,孟小红.优化滤波在青藏高原东北缘重力异常处理上的应用[C].中国地球物理学会第二十七届年会论文集,2011,695-696
[31]韩道范.BP神经网络在反演多层密度分界面中的应用研究[J].物探化探计算技术,1994,16(3):249~252.
[32]宋清华.基于小波神经网络的机器宠物之语音情感识别研究[D].东南大学,2009.
[33]金英雄.肾交感神经活动分析方法的改进及应用[D].南开大学,2006.
[34]朱自强.同时反演两个三维密度界面的拟神经网络BP算法[J].石油物探,1995,34(1):76~85.
[35]朱自强.重磁反演的拟人工神经网络BP算法及其应用[J].中南矿冶学院学报,1994,25(3):288-293.
[36]耿喜哲.小波神经网络在重磁资料反演中的应用前景[J].物探与化探,2001,25(2):102~108.
[37]陈东敬.张新兵.带模拟退火算法的BP神经网络在伊朗某地区进行了重力资料的反演[J].勘探地球物理进展,2005,28(3):215-218.
[38]裘士忠.拟神经网络反演在煤田中应用[J].中国煤田地质,20051(5):104~106.
[39]曾政祥,吴良才,周玲,危志明.基于MATLAB的遗传神经网络在GPS高程转换中的应用[J].矿山测量,2006年,3:(37-39).
[40]闻新,周露,王丹力,熊晓英等,《MATLAB神经网络应用设计》[M].科学出版社,2001年.
[41]闻新,周露,李翔,张宝伟等,《MATLAB神经网络仿真与应用》[M].科学出版社,2003年.
[42]罗成汉,基于MATLAB神经网络工具箱的BP网络实现[J].计算机仿真,2004,21(5):109~111.
[43]杜志敏.双窗口模型结合小波多尺度分析识别人类基因组Isochore的边界[D].河北工业大学,2007.
[44]邓承志.基于小波变换图像去噪研究[D].江西师范大学,2005.
[45]董晓马.基于小波分析的框架结构故障诊断方法研究[D].东南大学,2003.
[46]王长军.基于小波分析的结构模态识别和破损诊断方法研究[D].武汉理工大学,2002.
[47]邢海英.利用小波算法建立棉纤维表面分形数学模型分析研究棉纤维的含糖量.[D].天津工业大学,2003.
[48]温省云.小波变换在isochors边界识别中的应用[D]天津大学,2004.
[49]叶振凯.基于小波分析的交流电力测功机测试系统研究[D]湖南大学,2004.
[50]金鑫.虚拟仪器及小波分析技术在电容设备绝缘介损在线检测中的应用研究[D]广西大学,2005.
[51]党晓强.小波变换技术在电容设备在线监测中的应用研究[D]四川大学,2003.
[52]黄晶,基于小波和改进相关分析的介损值算法研究及其应用[D].湖南大学,2006.
[53]袁永南.基于小波变换的地基云图云状分割[D].京信息工程大学,2011.
[54]崔雪梅.以复小波提取局部放电信号特征的原理及方法研究[D].重庆大学,2004.
[55]胡婕等.基于小波变换的污水处理过程测量数据去噪研究[C].第17届中国过程控制会议,2006.
[56]侯遵泽,杨文采.重力异常多尺度反演研究及应用[C].中国地球物理,2011
[57]刘金兰,李庆春,何秦娥.小波多分辨分析分离山西断陷盆地重磁场的研究[C].中国地球物理,2008
[58]王明,孟小红,郭良辉,赵博.多尺度分析法在南海磁异常分离中的应用[C].中国地球物理,2008
[59]杨文采,施志群,侯遵泽,程振炎离散小波变换与重力异常多重分解[J].地球物理学报,2001,44(4)534-542
[60]李宗杰,杨林,王勤聪.二维小波变换在位场数据处理中的应用试验研究[J].石油物探,1997,36(3)70-79
[61]李宗杰,杨林,王勤聪.小波变换在位场数据处理中的应用.[J].石油物探,1997,36(2)86-93
[62]飞思科技产品研发中心.小波分析理论与MATLAB7实现[M].北京:电子工业出版社,2005
[63]杨长保,聂兰仕,孙鹏远.改进的模糊神经网络模型在储层产能预测中的应用研究.[J].吉林大学学报(地球科学版)2003,33(1)48-55
[64]肖锋.重磁数据处理系统的开发及其在新疆地壳构造研究中的应用[D].吉林大学硕士毕业论文,2004.
[65]董焕成.《重磁勘探教程》北京:地质出版社,1995
[66]郭志宏,一种实用的等值线型数据网格化方法[J].物探与化探,2001,25(3):203-208
[67]罗亦泳,张立亭,陈竹安等.基于sufer的数据网格化与体积计算精度分析[J].测绘科学,2009,34(5):97-100
[68]李振海,汪海洪.重力数据网格化方法比较[J].大地测量与地球动力学,2010,30(1):140-143
[69]汪隆六,三角形线性插值在区域重力测量数据网格化中的应用[J].物化探计算技术,1986,8(2):141-146
[70]尹春蕾,刘四新. Matlab环境下地球物理波动方程的求解.[C].中国地球物理,2009.
[71]张凤旭,张凤琴,刘财,孟令顺.基于余弦变换的匹配滤波方法分离重磁异常[J].石油地球物理勘探,2006,41(2)216-227
[72]王家林.石油重磁解释[M].北京:石油工业出版社,1991.
[73]越南推进老挝钾盐项目[J].化工矿物与加工,2008(11):42.
[74]韩元红,马海州,张西营,袁小龙,李建森.沙空那空盆地龙湖钾盐矿区盐岩上覆泥岩矿物学及地球化学研究[J].2011增刊:787-788.
[75]连静,王珂,李光鑫.基于边缘的小波图像融合算法[J].通信学报,2007,28(4):17-23.