多项式趋势面分析在磷块岩深...用：以瓮福磷矿英坪矿段为例_罗俊峰.pdf

第 43 卷 第 4 期 矿矿 物物 学学 报报 Vol.43,No.4 2023 年 8 月 ACTA MINERALOGICA SINICA Aug.,2023 文章编号：文章编号：1000-4734(2023)04-0533-14 多项式趋势面分析在磷块岩深部勘探中应用：以瓮福磷矿英坪矿段为例 罗俊峰1，周光红2*，苏之良1，文瑞成3，蔡健龙3，欧义3，罗泰义4（1.贵州省有色金属和核工业地质勘查局 地质矿产勘查院，贵州 贵阳 550005；2.贵州师范学院 地理与资源学院，贵州 贵阳 550018；3.瓮福（集团）有限责任公司，贵州 福泉 550509；4.中国科学院地球化学研究所 矿床地球化学国家重点实验室，贵州 贵阳 550081）摘要摘要：磷块岩是重要的磷矿资源，主要发育于大陆边缘的台地相，在一定空间范围内磷块岩厚度具有稳定而缓慢变化的特点，多项式趋势面分析方法非常适于刻画这种空间演化特征。本次研究收集了瓮福磷矿田高坪矿区英坪矿段 B 矿层未经历后期断层破坏或抬升剥蚀的 88 项勘探工程资料，以 B 矿层厚度为因变量，以工程的地理坐标为自变量，开展 15 次多项式趋势面分析方法在沉积磷矿床深部勘探的示范性研究。获得如下主要认识：1）坐标选择：根据勘探范围的空间尺度，通过坐标平移进行坐标简化，能更加突出因变量的空间演化特征；2）数据变换：因变量的数据变换对变量之间的相关程度没有影响，自变量的等比例缩放不影响变量之间的相关程度，但极差正规化和极差标准化对变量之间的相关程度有微弱影响；3）最优趋势面判定及异常特征：以 B 矿层厚度等值线图为基础，结合趋势面的拟合度和F检验值的转折变化，判定3次多项式趋势面为最优趋势面（拟合度为77.39%，F 检验值=39.11F0.001(7,87)3.94），相应 B 矿层厚度的正异常主要集中在英坪矿段的中部东侧；4）古地貌特征：最优趋势面趋势值的负值表征了晚陡山沱期英坪矿段为从西缘瓮福古岛向东逐渐变深的潟湖，北侧和南侧分别为谭家院水下高地和湾子半岛，远东方向可能受上塘青平障壁岛限制；5）找矿预测：根据已有的岩相古地理研究成果和勘探工程资料，结合最优多项式趋势面特征，综合预测英坪矿段东侧具有寻找厚大 B 矿层磷块岩的较大潜力。研究表明，基于多项式趋势面分析在磷块岩矿床深部勘探中结果可靠，具有重要应用价值。关键词关键词：多项式趋势面分析；深部勘探；磷块岩；英坪矿段；瓮福磷矿 中图分类号中图分类号：P612 文献标识码文献标识码：A doi:10.16461/ki.1000-4734.2023.43.041 第一第一作者：作者：罗俊峰，男，1985 年生，高级工程师，从事地质找矿工作.E-mail: Application of polynomial trend surface analysis in the exploration of phosphorite in depth:Taking the Yingping oreblock of Wengfu deposit as an example LUO Jun-feng1,ZHOU Guang-hong2*,SU Zhi-liang1,WEN Rui-cheng3,CAI Jian-long3,OU Yi3,LUO Tai-yi4(1.Institute of Geology and Mineral Exploration,Guizhou Bureau of Non-ferrous Metals and Nuclear Industry,Guiyang 550005,China;2.School of Geography and Resources,Guizhou Education University,Guiyang 550018,China;3.Wengfu(Group)Co.,Ltd.,Fuquan 550509,China;4.The State Key Laboratory of Geochemistry,Institute of Geochemistry,Chinese Academy of Sciences,Guiyang 550081,China)Abstract:Phosphorite is an important phosphate ore resource,and is mainly developed in the platform facies of continental margin.The thickness of phosphorite is stable and slowly changing in a certain spatial range.Therefore,the polynomial trend surface analysis method is very suitable to describe this spatial evolution characteristics.In this study,收稿日期：收稿日期：2023-03-28 基金项目基金项目：贵州师范学院国家科技部和国家自然科学基金项目奖励补助资金项目（编号：黔科合平台人才20185778-04）；国家自然科学基金项目（编号：41802027;U1812402）；贵州省高等学校数字乡村创新团队（编号：黔教技2023076 号）*通讯作者，E-mail: 534 矿 物 学 报 2023 年 data of 88 items associated with the B ore layer,which were not obviously disturbed by faults or uplift denudation,were collected from exploration projects of the Yingping oreblock Gaoping Mining area,Wengfu Phosphate Mine.Taking the thickness of the B ore layer as a dependent variable and the geographical coordinates of the exploration projects as independent variables,a demonstration study on the 1-5 orders polynomial trend surface analysis method applied in the exploration of sedimentary phosphate deposit in depth has been carried out.The main following results and understandings are obtained.1)Coordinate selection:Based om the spatial scale of exploration area,the coordinate simplification through the coordinate translation can highlight the spatial evolution characteristics of dependent variables;2)Data transversion:The data transversion of dependent variable and the proportional scaling of independent variable have no effect on the degree of correlation between variables,but the range normalization and range standardization have weak effect on the correlation between variables;3)Optimal trend surface discrimination and anomaly characteristics:Based on the contour map of B ore layer thickness,combined with the collaborative changes of goodness-of-fit and F-test value of the trend surface,the 3 order polynomial trend surface is determined to be the optimal trend surface(with the goodness-of-fit of 77.39%and F-test value(F)of 39.11 which is larger than F0.001(7,87)(about 3.94),and the positive anomalies of the B ore layer thickness are mainly concentrated in the eastern central part of the Yingping oreblock.4)Palaeogeomorphic characteristics:The contour map of negative value of the optimal trend surface indicates that the Yingping oreblock of the Late Doushantuo period is a lagoon gradually deepening eastward from the Wengfu Island in the western margin,the Tanjiayuan underwater highland and Wanzi Peninsula in the north and south sides,respectively,and may be restricted by the Shangtang-Qingping barrier island in the far-east side;5)Exploration prediction:Based on the previous research results of lithofacies palaeogeography and exploration projects,combined with characteristics of the optimal polynomial trend surface,it is comprehensively predicted that the east side of the Yingping oreblock has great potential to find thick and large B phosphorite ore layer.The results show that the polynomial trend surface analysis is reliable in the exploration of the deep part of phosphorite deposit and is of important application value.Keywords:polynomial trend surface analysis;exploration in depth;phosphorite;Yingping oreblock;Wengfu phosphorite deposit 磷是地球生命的基本