河南嵩县地区金矿资源潜力估算

河南嵩县地区金矿资源潜力估算龚庆杰2009.11.4中国地质大学（北京）地球化学教研室Tel:010-82320895E-mail:[e-mail]&GeoExpl培训河南嵩县地区金矿资源潜力估算一、地质概况二、化探数据概况三、化探数据处理四、制作元素地球化学图五、制作元素地球化学异常图六、元素资源量估算七、元素地球化学组合异常图八、剥蚀系数计算与资源量估算一、地质概况面积：110*70km2地层：Ar太华群地层：Pt1宽坪群地层：Pt21熊耳群地层：Pt22汝阳群地层：Pt22管道口群地层：Pt3栾川群地层：寒武系地层：KENQ岩体：五丈山花岗岩岩体：花山花岗岩岩体：合峪花岗岩岩体：小型斑岩体构造：近EW向断裂构造：近NW向断裂构造：近NE向断裂构造：近SN向断裂矿床：几十处金矿床矿床：几处钼矿床矿床：几处铅锌矿床二、化探数据概况1:20化探数据本实例：原始数据数据点：1774点数据项：12项元素项：AgAsSbHgCuPbZnAuWSnMoBi三、化探数据处理离散数据点位图离散数据分布检验数据校正处理（滤波）离散数据点位图1.取样空间范围，存在空白区2.规则网格取样，取样间距2km分布检验1.非对数正态分布2.存在奇异数据点（单点异常）表现在两端（极低值与极高值）分布检验AuAgAs分布检验SbHgCu分布检验PbZnW分布检验SnMoBi结论：存在奇异数据点（单点异常）办法：滤波滤波——删除低端数据直接删除低端奇异数据（置为空值）删除原则：①低于或等于检出限②孤立数据点（从小到大排序时明显偏小）按此原则本实例中删除数据个数为：Au-1；Ag-6；Sb-2；Cu-1；As-0；Hg-0；W-1；Sn-5；Mo-1；Bi-6；Pb-0；Zn-0；获得无低端奇异点的滤波数据本实例：滤波删低滤波——中位数中位数滤波范围：5点即可（1.1倍网格间距）采用中位数插值方法获得网格数据将网格数据转换为表数据获得无奇异点的滤波数据本实例：滤波中位三种数据的比较原始数据删低中位minmaxminmaxminminminmaxAg1026800104244292400As2.11335152.152.518.8Sb0.122.3690.1740.203.18Hg109400110010600三种数据的比较原始数据删低中位minmaxminmaxminminminmaxCu5525111181010101Pb8.613279358.6313.51242Zn39.161007839.1748.3615Au0.025541110.300.342.7三种数据的比较原始数据删低中位minmaxminmaxminminminmaxW0.3236870.470.6537.6Sn0.736251.2-20.911Mo0.0961.28110.3370.425.4Bi0.0228.1170.06-10.055.55结论：中位数滤波数据稳健利用中位数滤波数据滤波分布检验1.非对数正态分布2.无低端奇异数据点AuAgAs滤波分布检验1.非对数正态分布2.存在高端奇异点SbHgCu滤波分布检验PbZnW滤波分布检验SnMoBi滤波分布检验结论：仍存在奇异数据仅表现在高端（固有性质）四、制作元素地球化学图数据：中位数滤波网格数据方法：累频，19级特点：与数据绝对大小无关优点：每一元素均存在异常；背景与异常对比清晰缺点：出现频率高的相同数据可能划归不同级元素地球化学图AuAgAsSb元素地球化学图HgCuPbZn元素地球化学图WSnMoBi五、制作元素地球化学异常图异常与背景的区分：①定值（异常下限）；②变值（趋势面）本实例采用定值（异常下限）的原因：地球化学景观相同；数据范围不大确定异常下限的方法：①累频方式；②均值标准差方式；③分形方式本实例选用累频方式（其它方式专门讨论）累频方式确定异常下限及浓度分带异常下限：85%背景：0%~85，颜色索引9异常三级分带：外带：85%~92%，颜色索引272中带：92%~98%，颜色索引274内带：98%~100%，颜色索引278计算异常下限异常下限元素异常下限中带下限内带下限Ag140180350As10.611.614.1Au2.94.812.5Bi0.550.771.9Cu313649Hg405070Mo1.92.88.1Pb65.994149.7Sb0.931.011.4Sn3.94.87.4W4.986.712.4Zn131148197浓度分带元素地化异常图AuAgAsSb元素地化异常图HgCuPbZn元素地化异常图WSnMoBi六、元素资源量估算1243567Au异常图1243567异常序号异常名称异常面积数据个数内带均值丰度估计面金属量1祈雨沟13.15320.82501002窑沟36.83917.26191693前河、店房39.141023.79484484牛头沟、萑香洼31.22918.86772275上宫、虎沟11.06316192426康山7.72225.52041047松里沟——316.419747Au异常6个内带丰度估计法资源量预测类比法思想计算公式研究区元素质量=异常面积*深度*密度*异常均值当深度与密度一致时仅取决于异常面积和异常均值预测区控制区元素质量预测资源量成矿系数元素质量已知总储量预测区控制区均值异常面积预测资源量成矿系数均值异常面积已知总储量**面金属量法资源量预测类比法思想计算公式成矿元素质量=异常面积*单位面积面金属量*深度*密度面金属量=异常面积*异常超量=异常面积*(均值-下限)当深度与密度一致时仅取决于异常面积和异常超量预测区控制区超量异常面积预测资源量成矿系数超量异常面积已知总储量**异常序号异常名称丰度估计面金属量丰度估计资源量面金属量资源量1祈雨沟25010030*30*2窑沟61916974513前河、店房9484481141344牛头沟、萑香洼67722781685上宫、虎沟1924223136康山20410424317松里沟197472414Au异常6个内带的资源量预测剥蚀系数随剥蚀深度不同，异常规模也不同。因此不需引入剥蚀系数，如丰度估计法、面金属量法。但异常在垂向上若具有对称性，即矿体异常规模强，头晕和尾晕弱，难以区分头晕和尾晕。直接导致类比时垂向深度相同，异常均值向下是增强或减弱，即垂向不均匀。因此引入：剥蚀系数七、元素地球化学组合异常图选定成矿指示元素成矿元素：Au待选指示元素：AgAsSbHgCuPbZnWSnMoBi元素地球化学组合异常图异常空间套合AuAgAsHg组合异常图AuSbCu组合异常图AuPbZn组合异常图AuWSn组合异常图AuMoBi组合异常图八、剥蚀系数计算与资源量估算头晕：AsSbHg矿体：AuAgCuPbZn尾晕：WSnMoBi数据处理：①均匀化（除以均值），消除量纲；②头晕、矿体、尾晕各自累加；③剥蚀系数=(体/体头+尾/尾体+2尾/尾头)/4剥蚀系数计算示例矿体异常值晕比值剥蚀分带尾晕矿体头晕尾/尾头体/体头尾/尾体系数无矿01200.330.000.08矿头1240.20.330.330.27矿体2420.50.670.330.50矿尾4210.80.670.670.73无矿21011.000.670.92GeoExpl软件实现数据均匀化多元素叠加分析剥蚀系数分布检验剥蚀系数网格化（最近点、1.1倍网格间距）等值区与等值线剥蚀系数与金异常空间对比分析金矿床与剥蚀系数考虑剥蚀系数的资源量预测预测区控制区剥蚀系数元素质量预测资源量剥蚀系数元素质量已知总储量//控制区控制区剥蚀系数元素质量已知总储量控制区剥蚀系数元素质量已知总储量成矿保存系数/异常序号异常名称丰度估计面金属量剥蚀系数丰度资源量面金属量资源量1祈雨沟2501000.60230*30*2窑沟6191690.75559403前河等9484480.5821181394牛头沟等6772270.57286725上宫等192420.45530176康山2041040.38039497松里沟197470.5442616Au异常6个内带的资源量预测预测资源量丰度估计面金属量已探明序号异常名称无剥蚀剥蚀无剥蚀剥蚀储量1祈雨沟30*30*30*30*30*2窑沟745951407+3前河等11411813413943+4牛头沟等8186687244+5上宫等23301317316康山24393149167松里沟2426141626资源量预测结果对比方法选择：考虑剥蚀系数的丰度估计法类似体积品位法