山东省归来庄金矿床异常分布及其时空演化.docxVIP

山东省归来庄金矿床异常分布及其时空演化 归庄金矿床是芦溪地区首次发现的大型新金矿床，其成矿过程基本清晰。能否扩大矿床规模,如何进一步指导矿区周边今后的找矿工作是摆在研究者面前的难题。通过在时空域内对归来庄金矿床的金异常进行地质统计学研究,以及总结该区金异常分布及其时空演化规律,分析金的赋存部位,试图科学地指导该区今后的找矿工作,为扩大矿床规模提供理论与方法依据。研究区为山东省平邑县归来庄金矿区及其周边地区,面积约1 000 km2。 1 准平稳假设 在时空域内,区域化变量Au在研究邻域内满足平稳或准平稳假设,其数学期望值未知,即E[Z(x)]=m未知,故可用普通克立格法进行研究。研究中用到的Au的数据来自预测区内Au的265个已知取样点。 1.1 数据分布特征 由于所有资料都是经过处理的,不是区域化变量的原始值(不带任何单位),所以这些资料只有指证意义。 对Au的数据作直方图并对其拟合(图1),可以看出Au呈单峰分布,变化范围在-1.88～52.29,变化幅度大,数据主要集中在小值区;大值区的数据虽然很少,但对数据的整体分布影响显著;数据分布明显右偏,需对特异值进行处理。经逐步估计邻域法处理原始数据后Au数据的分布轻微正偏,但基本可以认为服从正态分布。 1.2 最佳结构拟合过程 区域化变量的时空分布特征分析主要是通过计算区域化变量在不同方向的变异函数,并对反映其时空结构的变异函数实验曲线进行拟合来实现。 经估计邻域法处理特异值后Au的数据分布仍轻微正偏。两两相关半变异函数(pairwise relative semivariogram)对稀少的、有外来值(特异值)且服从正偏分布的数据很有用。将两两相关半变异函数作为时空结构度量方法,对5个不同方向(0°,22.5°,45°,90°,135°)Au的数据分布进行结构分析,结果表明Au在不同的方向基本显示各向同性,可以当作各向同性处理。 对Au进行套合结构分析的过程其实就是一个结合金矿成矿过程对Au的时空变异实验曲线进行拟合的过程。由于归来庄金矿床的成因主要与次火山杂岩有关,次火山杂岩体对金的成矿起主导作用,所以在研究中只考虑次火山岩这一主要成矿地质作用,且只用一个套合结构对归来庄金矿床Au的时空结构进行拟合(图2)。从图2可以看出,现时相Au的变异函数实验曲线可以用一个块金效应为5个单位、基台值为200个单位、变程为12 000 m的球状模型的曲线拟合。块金效应的大小取决于Au的区域背景值;基台值远大于块金效应说明该区的Au数据在某些地方明显高于背景值,指示有Au异常存在,说明该区有找矿前景;变程主要与该区铜石次火山杂岩体有关。该杂岩体出露面积超过30 km2,最大长度大于12 000 m。 归来庄金矿床现时相Au的时空变异模型可简化为: γ(h)=C0+C(1.5h/a?0.5h3/a3)=5+97.5(h/4000?h3/120003)0＜h≤aγ(h)=C0+C(1.5h/a-0.5h3/a3)=5+97.5(h/4000-h3/120003)0＜h≤a 该时空结构参数与该区地质情况相对应,可以用来进行普通克立格法估计。 1.3 时空变异结构模型验证 根据Au的时空结构变异特征,结合归来庄金矿床具体地质情况,同时考虑到原始数据间距较大等因素,制定适合本区的普通克立格法估计方案。将研究区分为90×80个待估区,待估区大小为400 m×325 m,搜索直径为12 000 m,各向异性比为1。 为了验证Au的时空变异结构模型在该研究区内的正确性,根据普通克立格法的设计方案,用265个原始数据值作为检验标准进行交叉验证。验证结果表明其估计精度较高,因而用该模型和方案对该区现时相Au进行估计,其结果是可信的。 对研究区内的7 200个待估区Au的平均品位进行普通克立格法估计,估计结果见图3。图3显示,该区区域化变量Au的普通克立格法估计高值区集中在归来庄金矿区,归来庄金矿田南部与次火山杂岩体南部边缘接触部位及铜石镇西与次火山岩杂岩体的西部边缘接触部位也有较高的Au品位出现。 1.4 金矿床au异常区地质异常及定位预测 获取研究区内Au估计量的最终目的是了解Au的时空结构,分析Au异常的分布规律,为扩大矿床规模、进一步指导今后的找矿工作服务。根据Au品位的时空结构所确定的异常在地质异常理论中属于“专属致矿地质异常”,这类异常的分布区可以称为“找矿可行地段”。根据地质异常理论,山东省平邑县归来庄地区的Au异常(Au数值大于4个单位的区域)属于局部地质异常,与归来庄金矿田相一致。该区Au异常具有地质异常的基本特性: (1)不连续性和突变性。从现时相Au的普通克立格法估计的等值线图上可以看出,Au异常的分布具有突变性,Au的数值在小于4个单位时等值线显得平缓、稀薄,而大于4个单位后等值线突