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

河南嵩县地区

金矿资源潜力估算龚庆杰2009.11.4中国地质大学（北京）地球化学教研室Tel-mail:qjiegong@163.comGeoMDIS&GeoExpl培训河南嵩县地区金矿资源潜力估算一、地质概况二、化探数据概况三、化探数据处理四、制作元素地球化学图五、制作元素地球化学异常图六、元素资源量估算七、元素地球化学组合异常图八、剥蚀系数计算与资源量估算一、地质概况面积：110*70km2地层：Ar太华群地层：Pt1宽坪群地层：Pt21熊耳群地层：Pt22汝阳群地层：Pt22管道口群地层：Pt3栾川群地层：寒武系地层：KENQ岩体：五丈山花岗岩岩体：花山花岗岩岩体：合峪花岗岩岩体：小型斑岩体构造：近EW向断裂构造：近NW向断裂构造：近NE向断裂构造：近SN向断裂矿床：几十处金矿床矿床：几处钼矿床矿床：几处铅锌矿床二、化探数据概况1:20化探数据 本实例：原始数据 数据点：1774点 数据项：12项 元素项：AgAsSbHg CuPbZnAu WSnMoBi三、化探数据处理离散数据点位图离散数据分布检验数据校正处理（滤波）离散散数数据据点点位位图图取样样空空间间范范围围，，存存在在空空白白区区规则则网网格格取取样样，，取取样样间间距距2km分布检验非对对数数正正态态分分布布存在在奇奇异异数数据据点点（（单单点点异异常常））表现现在在两两端端（（极极低低值值与与极极高高值值））分布检验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倍倍网网格格间间距距））采用用中中位位数数插插值值方方法法获获得得网网格格数数据据将网网格格数数据据转转换换为为表表数数据据获得得无无奇奇异异点点的的滤滤波波数数据据本实实例例：：滤滤波波中中位位三种种数数据据的的比比较较原始数据删低中位minmax<min>maxmin<minminmaxAg1026800104244292400As2.11335152.152.518.8Sb0.122.3690.1740.203.18Hg109400110010600三种种数数据据的的比比较较原始数据删低中位minmax<min>maxmin<minminmaxCu5525111181010101Pb8.613279358.6313.51242Zn39.161007839.1748.3615Au0.025541110.300.342.7三种种数数据据的的比比较较原始数据删低中位minmax<min>maxmin<minminmaxW0.3236870.470.6537.6Sn0.736251.2-20.911Mo0.0961.28110.3370.425.4Bi0.0228.1170.06-10.055.55结论论：：中中位位数数滤滤波波数数据据稳稳健健利用用中中位位数数滤滤波波数数据据滤波分布检验非对对数数正正态态分分布布无低低端端奇奇异异数数据据点点AuAgAs滤波分布检验非对对数数正正态态分分布布存在在高高端端奇奇异异点点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.9Cu313649Hg405070MoPb65.994149.7Sb0.931.011.4SnW4.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.60