z0001-2013万象地区钾盐矿首期勘查地球物理方法研究工作

项目名称：老挝万象地区钾盐矿首勘查地球物理方法研承担单位：吉林大单位行政：单位技术：赵报告编写人：主要研究人员：吴燕冈肖锋杨国东杨长保管彦武路鹏郇恒飞李秋辰李红雨李霖常畅陈曦杨刘强项目起止时间：20132月-20163《老挝万象地区钾盐矿首期勘查地球物理方法研究》课题来源于OrientMiningLimited(东方矿业)《老挝万象地区钾盐矿勘探》项目中的一部分。课题的任务是，在老挝人民 中部的万象地区东经 .15"～.86175250.38～180044.95范围内(勘探区边界为不规则形状，高件的编绘。计坐标的范围是横坐标Y=18255～18265.5km，纵坐标X=1979～1994km，勘探总面了野外任务，完成的具体工作量和精度统计见表0-1和表0-2。0-1完成的野外工作量统计表50-2实达勘探精度统计表重力勘探精度(105ms2201331320136月20最终的平面找矿靶区，其中I类靶区5个、II类靶区7个，为钾盐矿钻探提供较可靠的依 TOC\o"1-2"\h\z\u前 工作目的与任 地质任务和勘探部 工区自然地理条 管理措 人员编制和主要仪器装 工作时间安 完成的工作量与精 GPS测量工作方法与质量评 GPS测量施工前的准 GPS测量野外工作方法技 GPS测量资料整理与质量评 GPS测量结论与建 重力勘探工作方法与质量评 重力勘探施工前的准 重力勘探野外工作方 重力地形改正工作方 重力数据的整理计 重力资料的质量评 重力异常处理解 原始资料的预处 解析延拓及初步解 剩余异常求取和沉积盆地推 视密度计算及其三维分布特 断裂或构造边界的识 结论与建 结 存在的问 建 提交资 提交原始资 提交归档资 6.3提交成果资 参考文 附件1老挝万象地区钾盐矿首期勘查地球物理方法研究GPS归档资附件2老挝万象地区钾盐矿首期勘查地球物理方法研究地形改正归档资附件3老挝万象地区钾盐矿首期勘查地球物理方法研究重力归档资附112233前《老挝万象地区钾盐矿首期勘查地球物理方法研究》课题(合同编号：2013220101000225)是OrientMiningLimited(东方矿业)《老挝万象地区钾正常勘探区开展1:50000(局部加密成1:25000)的高精度GPS、重力共设计测线22579160152.8km2。完成测量和重力测点各778个，质量优于设计要求。工作目的与任《老挝万象地区钾盐矿首期勘查地球物理方法研究》课题来源于OrientMiningLimited(东方矿业)《老挝万象地区钾盐矿勘探》项目中的重力勘探部ORIENTALMININGLIMITED（東方礦業）是亿泰资本旗下的一个略开发，在2013年该公司经老挝批准，获得了在老挝万象市118平方公里的钾盐吉林大学是教育部直属重点综合性大学，1995年首批通过国家教委“211工程”，2001年被列入“985工程”国家重点建设的大学之一。原长春地质学院，是一学院拥有重点学科，是我国固体地球科学领域科技实力最强的单位之一，具有甲级研究与勘查资质，在寻找隐伏矿产资源方面有成技术，此课题次已经是第五次地质任务和勘探老挝地处大湄公河次区域(GMS)的中心位置，是次区域经济合作机制的成员国铁、锡、铝、铅、锌、锰、煤、盐和宝石等矿藏。进入21世纪后，老挝将矿藏的全面开发利用列入10年发展规划和20年发展的主要开发和发展领域，并制定了一系解释结果提供钾盐矿目标靶区的坐标和范围；④结合实测高程数据进行1:25000的地质、万象钾盐勘查118km2，距离万象市东约工区位于老挝人民中部的万象盆地，在首都万象万象钾盐勘查118km2，距离万象市东约1-1勘探区交通位置118km21-2，1-1。1-2勘探区范围和角点位置1-1工区角点坐标(ITRF1996UTMZONEABCDE制面积约向四周扩大1km，同时尽量保持测区边界整齐规则。研究区设计坐标的范根据OrientMiningLimited(东方矿业)的意见，按照1:50000的比×500m的正方形。共设计测线22条，线号由西向东从10～52双号编排，给加密线预留100～160方法、质检率和勘探工作量见表1-2：1-2设计野外工作量一览表57929608579296085792960816081681608168160816873937776ABABCDE 图-3设计的理论点编号及分布图工区自然地理条 屋 6老挝自北往南分为上寮、中寮、下寮3部分。行政上划分14个。中寮包括万象和甘蒙(他曲)两省，占面积的19.8％，人口的21.1％，扼最狭的腰部、当13地带，工农业生产居首位。1200与越南两国边界的自然分界线。长山山脉的是高原和山区。甘蒙省的石灰石峡谷是国和邻国提供了各种用途的木材。山脚下成片的乔木林有效地了山上的急流，保护老挝地处北纬14～23度，东经100～108度，属热带、带季风气候，全年分为雨季和旱季两季。雨季从5～6月份开始直到10月份，经常有雷雨，这个季节的气温在25℃～30℃，空气湿度接近100％，一般是北部省份比南部雨量要小。旱季从11月份开始到次年4月份，这个季节雨完全停了，有季风，平均气温为15℃，万象最凉爽的时候是12月份到2月份，气温在25℃度左右。4到6月份很热，平均气温达到35℃，而这时河1865段它的水流宽达几十公里。湄公河是老挝交通的大动脉和经济交流的，尤其是与中国和的物资交流。老挝境内主要河流还有北部的南塔河、南乌河、中部的南俄河和田区很难施工，见图1-4。1-4老挝万象附近地貌图锑、锡、锰、铝土矿、钾盐、石膏、煤、高岭土、、宝石、石墨、石棉、滑石、20管理措制度；专题组严格按照甲方的工作安排开始施工；外业结束后，经专题组初步验收工过程中各个环节进行全面的质量与评价，严格按照各项规章制度执行，发现问题意识，上岗前进行安全；仪器专人，并做好日常；专题的各个班组、障体系由专题统畴安排，确保经费。组建了一支能够协调的精干技术队人员编制和主要仪器装根据《合同》及《设计书》的要求，为保证项目的顺利实施和确保任务目标实现，及按时优质完成该项目外业任务，本专题组织了一支高素质的技术队伍来完成任入GPS1GPS量3探3正1-3主要技术人员素质一览表15生生刘生生肖4生李生杨生12常生12路肖作的顺利进行提供了保障。投入的主要仪器设备详见表1-4：1-4投入的主要仪器设备一览表1CG-3241MAGELLAN导航453521736151272工作时间安1-5工作时间安排表据；随时进行数据整理提供资料给甲方使用，提交初步研究成果。地完成的工作量与吉林大学的师生还是克服了重重，按时保质保量地完成了野外任务，共计完成测量和重力测点各778个，质量优于设计要求。野外各项具体工作量和精度统计如表1-6、表1-7和表1-8，详细情况见本归档资料的附件1《2013老挝GPS归档资料附1-6野外工作量统计表37464615561556461-7野外工作观测精度统计表近区地改中误差(105s2中区地改中误差(105s2远区地改中误差(105s2测点重力值观测精度(105ms2有效异常等值线距(105ms2)(105ms2(105ms2(105ms2(105ms2(105ms2GPS测量工作方法与质量幅1:10万地形图作为工作参考用图。中民GB/T18314-2009(代替GB/T18314-2001)《全球定GB/T18341-2001GB12898-91DZ/T0153-95中民地质矿产行业标准DZ/T0171-1997《大比例尺重力勘查规范中民地质矿产行业标准DZ0004-91《重力技术规定(1:50000)8618#和8623)，仪器外观见图2-1。 a2(ba2(b南方灵锐S82v(8617#，8623#，8618#和8359#

6.2mm中、长距离测量的误差d应符合下列要求：观测距离D5km时，d1D5km时，d2a2(ba2(b(XX)2(YY在GPS校准场上进行，按仪器操作要求正确安置和操作仪器，记录测量数据，经解算GPS所得点位坐标与标准点的坐标在同(XX)2(YY (2-网误差为：7.4mm。基线长度误差的不确定度分量 32(0.01c106)2 (2-式中：c为相邻点距离，mm估计基线长度误差在c50000mm时，10%，故自由度v150GPS安置误差的不确定度分量uGPS安置采用强制归心孔，其安置误差可控制在0.1mm范围内，安置误差u20.1mm/20%，则自由度v2GPS分辨力为1mm，故u30.29x0.291mm0.29mm,v3uc(d)

(2-veff (2-取置信概率为99%k99t9957)2.68,U99k99uc(d)2.681.04的内噪声，主要是由于硬件的不完善(如钟差、信号通道的时别输入两台，根据两台的观测数据，利用相对定位的原理解算相应的基线 0通过测试仪器的噪声水平为：0.001在无功分器的情况下，GPS的内部噪声水平，也可利用长度精确已知的超短GPS测量野外工作方法投影方式：六度带投长半轴a= e2：0.00669437999013高程系统：Lao97高程系统测区位于老挝人民中部的万象盆地，在首都万象市的东部约20km处。测区范围为东经.15"～.86175250.38～180044.95。测量的投影方式采用投影，投影的子午线为105°E。隔500m流动作业，每点观测5～10分钟，流动站相对于基准站的基线中误差为PDOP周期和足够的数目。为了获得垂直方向上好的结果，同时检查VDOP值。按照设计GPS控制网采用D级网，从2013313日至18日，历时6天，完成了GPS控制网的收料、踏勘、埋桩、静态观测、平差计算、绘图等任务。使用已知老挝国家控制点2036作为测区平面控制测量的起算数据，使用已知老挝国家控制点2036362012-1GPS10(OML01、OML02、OML03、OML04、OML05、OML06、OML07、OML08、OML09OML10)，建临时转点3个(ZD1、ZD2和ZD3)，经过平差计算，解算出各控制点的最终结果坐标见表2-2，2-2，GPS2-3。2-12-2GPSABCDE 3实测GPS图2-3中浅蓝色折线为平面坐标联测路线，浅蓝色和绿色结合为高程联测路线。黑色据，蓝色空心三角符号表示测得平面和高程数据的自建GPS控制点，绿色空心圆圈符号表根据设计要求开展1:5万面积性高精度GPS，重力测点的平面和高程定位采用快速静态测量模式，工作中使用4台国产南方灵锐S82双频GPS，仪器编号分别为8359#，8617#，8618#和8623#，其中8623#仪器作为基准站使用。从2013年319日～47日，历时20天，完成的野外重力测点和加密点的图2-4GPS野外工作2-2GPS8ABABCD河公E湄 5GPSGPS测量资料整理与质量评GPS静态观测依据中民GB/T18314-2009《全球(GPS)测量规范》中的相关规定，所有各项限差要求如下。GPS静态测量的所有观测资料，由GPS随机软件TBC进行相位差分计算。资料整理工作与野外施工同步进行，使用便携nn式中：v代表较差；n代表点检查点的

(2-等单项指标均符合设计和规范要求。精度统计见表2-3。2-3GPS2-6实测高程等值线图根据《老挝万象地区钾盐矿首期勘查地球物理方法研究高精度重力勘探工作设计GPSPSGPS/GPS/水准点的网格间距、地区重力异常的分辨率和DTM的分辨率提出要求；算，这就需要收集的老挝测区内的控制点。重力勘探工作方法与质量施工所用基础资料：收集到2011年编绘的1:100000地形图1幅，编号中民地质矿产行业标准DZ0004-91《重力技术规定(1:50000)SystemOPERATIONMANUAL重力勘探施工前的准本次重力野外工作中使用了三台生产的CG-5型高精度石英弹簧重力仪(编号用过程中应经常进行必要保养和，对仪器做了各项性能测试如下：CG-5型重力器必须充电48h后，使仪器掉格速度趋于稳定才能开始各种试验。静态试本次工作静态试验是在老挝野外工区驻地进行的，2013年3月12～13日对三台重力仪进行充电稳定，14～18日分别对143#、553#和763#仪器做了静态试验。三台仪器连续240.1210-5ms2，符合设计和规范的要求。三台重力仪的静态试验数据详见附件3《2013老挝重力归档资料》第七部分的《重力仪的静态测试资统计结果优于设计见表3-1，零点掉格曲线分别见图3-2、图3-3和图3-4。3-1重力仪静态试验统计表(10-5ms23060-60-60-3-2143#3-3553

3-4763#验。根据中民地质矿产行业标准DZ/T0171-1997《大比例尺重力勘查规372.8408～2.962110-5ms22010动

mm测边(段)数(当只在两个点上观测时，n1)。3-2重力仪动态试验结果统计表(10-(10-10333-5143#3-6553#3-7763#n2i1jMn2i1jM一致

(3-式中：Vij为某点(i)上某台仪器(j)之观测值与多台仪器在该点之观测值的平均值的差m为参加一致性观测的仪器台n为一致性试验的观测点M为差值Vij的总个数Mmn。3-3s2 3-8143#、553763#本次工作中使用143#、553#和763#重力仪，比例因子测定是2012年6月10日采程往返重复观测法在高崖口国家1#和4#标准点上完成的，两点间的重力差值为105.265×105ms2。3-9各台仪器都取得合格的独立读数差8个。达到了规范要求，比例因子测试原3-9比例因子计算公式为C

(3-式中：G为格值标定场两标准点间的重力差值，105ms2；g为重力仪在两点间读出iigi1 (3-i/[n(i/[n(n

(3-式中：i为第igi)与平均重力差g的差值；n重力仪比例因子相对变化的计算公式为

CC

(3-式中：C原为重力仪原标定的比例因子；C新为重力仪新标定的比差的2.5倍，比例因子相对变化都很小，比例因子是稳定的，仍可使用原比例因子，三台仪器的原比例因子均为1。计算过程详见附件3《2013老挝重力归档资料》第七部3-4原比例因子(105ms2每单位新测比例因子(105ms2每单位≥6≤2808080重力勘探野外工作方横坐标为Y=18257.361895kmX=1990.557555kmH=168.303mN0105ms23-10工区的重力基点位置重力以检查仪器是否处于正常工作状态。工区完成的工作量统计见表3-5，普通测点与检查点，以及加密点的分布与GPS测量相同，见图2-5，工作场景见图3-11。3-5重力野外工作量统计表152.857929401608152.8739373-11重力野外工作场景性零点漂移改正。由于每个工区只设一个重力基点，其相对重力值设为0105ms2gSS2S1TT (3- T 式中：S1、T1为起始基点上的重力仪观测值、观测时间；S2、T2为闭合基点上的重力仪观测值、观测时间；Si、Ti为某测点上的重力仪观测值、观测时间。使用重力仪观测得到的相对重力值，经上述线性零点漂移改正后(所使用的CG-5型重17.2679～5.3557105ms222.6236105ms2，平均值为-4.9823105ms2。由图可见，异常大致呈北东，西重力地形改正工作方地形改正工作均选择方域模式，其精度分配按照中民地质矿产行业标准余各区都用平板公式计算，且远区一般限制在20km之内，图3-13是用方域各区分割测点ababd(,cA远区型，由操作员在野外实地目估测点周围20m远处东、南、西和北4个方位相对于测点的高差，而南东、南西、北西和北东4个方位相对于测点的高差是由东、南、西和北4个方位g

Rd

d

hhh

232

(R

hhh

(3-GRln1 1 Rh2 2 R2hhhR2hhRh2 2 2节点的高程；h2为阴影部分斜边节点的高程。由上式可知，当改正半径R、划分的方位数n和中间层给定以后，只要测到了R距离上的地形与测点的高差h1和h2，就可以很方便的算出这个斜顶面三角棱柱体的野外近区地形改正数据采用操作员目估的方法。总计完成近区地形改正测abcdi,i)，方柱Ahi,jA1 1 gDGi d

j

d0

2)3/

(3-i

j

( 由于22r2

1(r2h2)10(222)3 0(r22)3 1 h21 12 r r1

212 11h

(3- r K 212 ij1ij (3- 1 r i1j ij

ij Kij为各方柱体所应用的系数，它包含了引力系数G及数值积分所用的通节点的常数为1以便保证各区域之间无缝隙和)；rij表示柱体到地改点的距离。源于全球遥感影像和DEMsrtm(航天飞机地形测绘计划)数据，为及其周围20km内的点节点高程数据，共获取节点高程数据353686个。取数节点位置详见3-14中浅蓝色圆圈符号所示，图中蓝色圆点符号为测点实际位置。 图-14中区地形改正高程取数位置图20m×20m区地改的边界相匹配，网格化后的工区及其周围2km的高程影像见图3-15；然后按公式中区地形改正质量检查采用移动方格网节点的方法，由于网格距是20m，所以纵向100m×100m20km3-16。方格网节点的方法，由于网格距是100m，所以纵向和横向移动均为50m。远区地形改正计算点数与移动网格后检查点数相同，质检率100%，全部参与精度统计。形改正值(20～2000m)及远区地形改正值(2000～20000m)2《20130.0018～0.1408×10-5ms20.0146×10-5ms2，总体来看研究3-6(105ms2(105ms2(0～205555重力数据的整理纬度校正值g的计算公式(单位105ms2g5203.3258sin20sin0 (3-式中：0为重力总基点的纬度，本区为北纬 为重力测点的纬度B布格校正值g的计算公式(单位105ms2Bg

0.0007cos20.72107H0.0419103

R20.3086

1

h

(3- 式中：为重力测点的纬度；H为测点的海拔高程，单位m；h为测点相对于总基点的高差，单位m；为中间层密度，单位kg/m3；R为中间层校正半径，定位20km，也是地改布格重力异常gB的计算公式(单位105ms2gBgCggB (3-图3-18。由布格重力异常图可见，工区的布格重力异常变化范围是-12.443～3.073105ms2，变化量为15.516105ms2，大致呈北东，西北高，东南转南低。3-18实测布格重力异常等值线图重力资料的质量n2n22g 2约为23105ms2，实测重力观测值的各项技术指标都优于设计。3-7测点观测精度统计表10ms10ms10ms10ms10ms--g纬度校正误差 gg0.000814sin20d0.000038105ms2 (3-式中：0为重力总基点的纬度，17.Nd为重力测点沿南北方向的定位均方误差，按点位中误差0.08m计算，详见本项目报告的附件1《2013老挝GPS归档布格校正误差 的计算公式为 g

0.308610.0007cos21.44107H0.0419103 0.010810-5m

(3-式中：17.～18.，取测点最小纬度17.；H160.380～188.160160.380；为中间层密度，本工区选取2.2103kg/m3；h为测点相对于总基点的高差，-7.923～19.857，取19.857m；R为中间层校正半径，也是地形改正的最大半径，老挝GPS归档资料》。值。近区地形改正均误差采用下式计算：t t式中：t为近区地形改正均方误差，105ms2；i为第i差，105ms2n2n2t

= = 105ms2 (3-2式中：t为中区地形改正均方误差，105ms2；i为第2差，105ms2n2in2i22t远 = = (3-式中：t为远区地形改正均方误差，105ms2；i为第i差，105ms2 t 0.000852+0.000112 =0.000000734616=0.0008610-5ms2(3- 3-8总地形校正精度统计表

t远为远区地形改(105ms2(105ms2近区地形校正(0～20布格重力异常误差gB的计算公式(单位105ms22222 B0.036620.036620.010820.0008620.038105ms2

(3- 式中：C为测点重力值均方误差；g为布格校正误差；g为地形改正均方误差，由于远区未做地形改正，故 等于近区和中区的总地形改正误差； 预期的效果。工区实测异常质量与设计的分项精度要求对比见表3-9。(105ms2(105ms2(105ms2(105ms2(105ms2重力异常处理解处理、转换的目的与作用如下：但是，有些解释方法需要Vxz、Vzz。为此，必须由实际异常计算它们。处理、转换的主要内容如下：原始资料的预处它是一种较灵活和常用的法，对任何类型的数据几乎都可以使用。对于大多数数由于野外测量时点线距是按500m500m和局部中点加密来布设格化以后的点线精确的处理解释结果，我们选择100m100m网格化点线距，这已经是勉为其难了。解析延拓及初步延 e-延 e-

(4-波器，它起到以低频成分为主的区域异常，相对突出浅层局部场的作用。为了了解本区基底情况，对重力异常作了100m、200m、300m、500m、800m和1000m等六个不同高度的向上延拓，各个不同延拓高度的重力异常示于图4-1。向上延拓某一高度的异常，大致反映该高度1/2深度以下密度不均匀的响应，各延拓图件大致反映50m、100m、150m、250m、400m和500m深度以下密度不均匀的影响。对比上述图4-1中各图可知，地表0～100m之间异常形态最为复杂，研究区基底的主体大致为北东向的斜坡状，呈西北低的梯度带，在研究区的中部又有多个局部500m 4-1向上延拓重力异常图300～500m；中部的北东转东西向的负异常中心强度次之，其异常中心变化不大，上延300m以后基本，说明该处可能存在低密度矿体，150m对更弱，异常中心变化也不大，上延300m以后基本，说明其可能存在低密度矿150m为了了解本区浅部地情况，对重力异常作了50m和100m等两个不同深度的向下延拓，各个不同延拓深度的重力异常示于图4-2。 4-2向下延拓重力异常图图4-2中的两幅图可见，前面讨论的负异常重力低中心仍然存在，且更加突出，更加低圈闭的顶深在50-100m左右。除此之外，在图4-2中还可见另外几个负异常圈闭，其异常强度较弱，埋藏较浅，可剩余异常求取和沉积盆地推常规的趋势面分析的实质是最小二乘意义下的多项式拟合。假设区域重力异常egGGregAxBy (4-每一个重力测点上区域重力异常Greg可以由测点坐标x，y以及多项式系数A，BGresGGregGAxByC (4-当G2res取得最小值时， 与G在最小二乘意义下拟合得最好。G2res取得

GresGres/A0；GresGres/B0；GresGres/C0 (4- /AAx2Bx2Cx2Gx /BAxyBy2CyGyGresGres/CAxByCG (4-利用迭代算法可以计算出(4-5)中的未知系数A，B和C，从而确定了最小二乘意义2、3 4-3(a)；(b)域重力异常，就可得到剩余重力异常，见图4-4(b)。 4-4对比图4-3和图4-4可见，区域场都过于圆滑，且局部场规模较大。负异常中心在ZK4和ZK5井附近，可能是有利的找矿区域，ZK9井附近也是有利的找矿区域。比较而GregAAxAx2AyAxyAy2Ax2yA (4- 有限元素法引入形状函数Nix,y，上式记作边界上八个离散节点值Gi之和Gregx,yNix, (4-测区边界上Gi的分布在理坐标下的分布如图4-5所示，中心点坐标XC,YC,长和宽分别为2a和2b。从左下角节点起，节点序号为1，沿逆时针方向节点 点的最大序号为8。图4-5真实坐标系下的二阶多项式单 图4-6参考坐标系下的二阶多项式单xxC，yyC (4- Ni,1i1iii1/4 (4-Ni,121i/2 (4-Ni,1i12/2，i=4,8 (4-将形状函数Ni,代入公式(4-12)中，内插出测区内每一个点的区域重Greg,Ni, (4-x,Ni,xi；y,Ni, (4-研究区利用二阶有限元法和三阶有限元法计算的区域重力异常见图4-7(a)和图图4-8(b)。 4-7二阶有限元分解异常图 4-8三阶有限元分解异常图给定正交尺度函数t和小波函数t，其二尺度关系tt

2h0k2tkk2h1k2tkk

(4-(4-式中h0k、h1k是多分辨率分析中的滤波w2ntw2n1t

2h0kwn2tk2h1kwn2tk

(4-(4-式中h0k、h1k仍然是式(4-14)和式(4-15)中的滤波器系数。当n=0w0ttw1tt以上定义的函数集合wnt 为由w0tt所确定的小波包，由此，wnt 是包括尺度函数w0t和小波母函数w1t在内的一个具有一定联系的HhejkHhejk i1

制数表示为：ni ，其中i=0或者1 21H021令

H02H' H2 wnt的频域表达 'WnHk2k (k0或者 (4-k 4-9小波四阶多分辨分析异常图(a)区域异常；(b)部异 4-10五阶小波包分解异常图(a)区域异常；(b)部异测区D内重力异常gx,y由区Rx,y和局部Lx,y组成，局部场范围较小，其所在区域为DzDzD，则有gx,yRx,yLx,gmingx,yLx,y

x,yx,yDzx,yDDz

(4-(4-(4-用切割算子A作用于gx,y，得第一次切R1x,yAgx, x,y (4-切割算A应使切割区域R1x,ygminR1x,yR1x,ygx,

x,yDzx,yDDz

(4-(4-R1x,y重复使用切割算子A得第二次切割R2x,y依次迭代有limRn1x,yRnx,y (4-nRnxy作为所求区域场的近似RxyRnxyx,yD；从观测LxygxyRnxyxyD掉区域重力异常，就可得到剩余重力异常，见图4-11(b)。 插值切割4-11ZK10很理想，或者地质情况比较复杂。但是在ZK4井附近负异常圈闭明显，只是圈闭范围Naudy和Dreyer（1968）提出了一种非线性滤波方法。它利用异常宽度作为唯一的滤ady滤波通过在空间域进行滑动窗口平均，弥补目标常频谱较宽情况下线性滤波方的不足为了更适于重力异常分离，对Naudy滤波的算法进行了改进。通过单拱、多拱形态的的2D重力模型及3D组合模型进行验证，试验表明，改进的Naudy滤波分离的区域异常与理论不太可靠。计算结果见图4-12。 4-12非线性分解异常图(a)区域异常；(b)部异对比分析图12可见，ZK4井至ZK5井一线负异常圈闭明显，至zk7井附近结束，这个负异常圈闭与这三口的钻探结果吻ZK9井至ZK14井一线的北东转南北向的负异常圈闭也是找矿的有利ZK8井附近异常圈闭较好。ZK10井以北为串珠状负异常圈ZK10500-1000m至负异常圈闭中心，效果可能要好得多。coefficentofcoefficentof upwardcontinuationheight/point4-13掉区域重力异常，就可得到剩余重力异常，见图4-14(b)。 向上延拓利用剩余重力异常推测钾盐沉积盆地的位置和范围：4-14(b)中围，见图4-15。 利用重力垂向导数异常推测钾盐沉积盆地的位置和范围：f(xyz)的nnf(xyzznnf(xyhhn，因此由深度因子hexp(h可以求出异常n~ ()n(22)n ~ 2(22

(4-(4-300m、500m800m的垂向二阶导数，它们的零值线圈闭见图4-16。 4-16(a)原始异常；(b)向上延拓100m；(c)向上延拓200m；(d)向上延拓300m；(e)向上延拓500m；(f)向上延拓800m用剩余重力异常预测的低密度钾盐靶区范围(图4-15)和布格重力异常向上延拓800m的垂向二阶导数(图4-16(f))吻合得非常好，从而增强了计算结果的可信度。特别是视密度计算及其三维分布特V2Gtg

x 1x 1xb h

2

(4-的重力值g就应该是所有棱柱体重力值的叠加，即Ng(x,y,z)gi(x,y,

(4-zx2y2x2y2x2y2zx2y2x2y2x2y2giGixlny

ylnx

y设一共有MN个棱柱体，则实测重力场值可写为：g(x,y)gi(x,

(4-g(u,v)gi(u,

(4- bvehr hri(uxvygi(u,v)2Gi sin .sin . (4- Si(u,v)4sinau.sinbv.11ehr 2 H(u,v) D(u,v)ei(uxvy u2柱体中心坐标。u2(,v(uv(uvSi(u,v)为尺寸因子，i(u,v)，Di(u,v)(u)g(u,v) 2GiSi(u,v).Hi(u,v)Di(u,

(4-深，即每个棱柱体的g(u,v)2GSi(u,v)Hi(u,v)iDi(u,

(4-g*(u,v)2GiDi(u,

(4- g*(x,y)2Gidi(x,

(4-由于在频域中已经进行了向下延拓和求导。这时的场值g*(x,y)已经是棱柱体顶g*(x,y)2 g*(x, (4- 2i即为这一点小柱体的视密度值。对每个直立棱柱体均计算i，就得到了地密无关；反演值与理论值最大误差小于15％。为了了解不同深度的密度不均匀分布情况，对研究区做了0-100m、100-200m、的重力异常反演该层位的视密度，各个层位异常及区域异常分布情况见图4-17。 0-100m100-200m视密度变化幅度最大，是寻找钾盐的最有利层位，主要低密度圈闭负异常圈闭，深度为300m左右。400-600m伸较深，可作为寻找钾岩的靶区。X=1988.5km、X=1987.5km、X=1985.5km、X=1981.5km、X=1981km，见图4-18。4-18上述9条剖面的立置见图4-19研究区其它位置的视密度剖面图，可根据需要图4-19部分剖面立置和视密度分布断裂或构造边界的识若sx方向的夹角为T°cosT°sinT°i(ucosvsin)

(4-(4- 可见i(ucosvsin)是计算任意方向的方向导数的波数响出沿x方向的方向导数对于x方向的高频成分有放大作用，因此实践中常用水平方向导数度的异常，分别做了水平0、45、90、135四个不同方向的导数，分别反映南北方 (a)0；(b)45 (a)0；(b)45针对利用重磁场资料解决地质构造的精细划分和解释问题,提出向小子域滤多的识别断裂构造的信息。这就是向小子域滤波的基本原理[22]。 4-22位场数f(i,j0)xyz3个方向n阶差分算n阶归一化差DN2可表x方向 fni,jfn1ir,jrfn1ir,jr) fn(i,j)fn1(ir,j,r)fn1(ir,j,

(4-(4- fn(i,j)fn1(i,j,0)fn1(i,j, fn/ (4- 其中n阶总差

f2nf2nf2nnA (4- f2nf2nf2n f,f式中fnfnfnfn进行90°相移f,ffnF

n

nF1

yn x y 4-234-20～4-214-22～4-23步识别出不同方向展布的断裂或构造线多条，示于图4-24。250m100mABABCDE 图-24综合识别的断裂或构造线结论与建根据OrientMiningLimited(东方矿业)就实施《老挝万象地区钾盐结局部重力低圈闭，其范围都不大，目的层顶深在50～100m左右。1000m研究区中心ZK4井附近存在一个北东转东西向的负异常中心，强度次之，其异常中心变化不大，上延500m以后，说明该处可能存在低密度矿体，由浅至深产状较陡，目的层底深约250m以内，也是相对有利的找矿靶区；外，中心上下层位置变化也不大，上延300m以后基本，说明其可能存在低密度矿体，由浅至深产状较陡，目的层底深约150m以内，也应对其有所关注。600m600m靶区范围，见图4-15。可开展进一步的勘探，从而可以减少工作的盲目性。垂向二阶导数的零值线可以大致圈定地的边界。分析对比不同延拓高度的垂向400m用剩余重力异常预测的低密度钾盐靶区范围(图4-15)和布格重力异常向上延拓800m的垂4-16(f))吻合得非常好，从而增强了计算结果的可信度。特别是前者还视密度计算得到其三维分布特征可知：0-100m视密度分布比较凌乱，低密度圈闭很多且规模较小，可能对应浅部的小型沉积盆地。100-200m视密度变化幅ZK2ZK16研究区西北部找矿的有利部位；在ZK10的东北存在三个低异常圈闭呈北西向串珠状展布，它们是研究区南部找矿的有利部位。200-400m视密度分布逐渐变得相对平视密度分布更加平缓，局部低密度圈闭还有显示，说明这些异常延伸较深，可作为寻找钾岩的靶区。600-800m视密度变化相当平缓；再向下800-1000m视密度分布特别平缓，寻找钾岩比较。望的7条东西向的剖面，各剖面视密度图中深蓝色指示了低密度体的位置和深度，为寻找时提取，见图4-18和图4-19。250m100m700mI5-1。 推断的A1自ZK9井起向ZK14井方向延伸，由多个呈串珠状的负异常圈闭中心组成，呈北西向转北北西向展布。在二次曲面趋势分析、五阶小波包分析、非线性分解