刘邦定,陈新跃,罗小亚,谭仕敏,马小双
(1.湖南科技大学,湖南湘潭 411201;2.湖南省地质调查院,湖南长沙 410116)
地球化学找矿是隐伏、半隐伏矿、难识别矿找矿获取成矿直接信息最主要的手段,是传统宏观矿化露头向微观找矿的延续(薛顺荣等,2008;张再天等,2013)。针对当前地表找矿难度越来越大的现状,运用水系沉积物地球化学测量方法寻找有色金属矿产发挥了越来越重要的作用(夏祥标等,2009;李玉芹等,2011)。本文在分析湘南坪宝地区1∶5万水系沉积物地球化学特征的基础上,结合该区内铅锌成矿地质条件,进行了找矿靶区预测,为进一步指导找矿提供有利依据。
湘南坪宝地区位于南岭东西向构造带中段北缘,处于郴州-蓝山NE向基底构造岩浆岩带中段,亦是耒阳-临武SN向构造带中段,是南岭多金属成矿带的重要组成部分,具有优越的成矿地质条件,矿产资源丰富(胡志科,2002;柏道远等,2008;乔玉生等,2010)。
区内出露地层为晚古生代泥盆系上统、石炭系、二叠系为主的浅海相碳酸盐岩夹海陆交互相砂页岩(许以明等,2007),以及少量零星的古近系、第四系覆盖(图1)。
区内构造发育,其构造骨架由一系列SN-NNENE向的复式褶皱和逆冲断裂构成向北撒开,向南收敛,北宽南窄,向北西微凹的弧形褶断带,同时伴有一系列NW-NNW向的横断层。这些大致等距的横断层与SN-NNE向的逆冲断层相互交织呈网格状构造格架(童潜明,1990;李建中等,2005;付旭,2012)。
区内岩浆活动频繁,岩体成群成带展布,共出露约百余处。一般呈岩脉产出,亦有岩豆、岩盘、岩瘤、岩株等。多为浅成-超浅成,以燕山早期的中酸性、酸性英安斑岩、石英斑岩、花岗斑岩、花岗闪长斑岩为主。这些岩体的展布,NNE向成岩带,NWW向成岩列,纵横交汇部位成岩群(图1)。矿床一般位于NNE向岩带与NWW向岩列的交汇部位(龚述清,2006;刘悟辉,2007;田旭峰等,2007)。
坪宝地区是南岭中段多金属成矿带的重要成矿区,矿产资源丰富,并以内生矿产为主。现有矿床15处,矿(化)点41处。矿种多达13种,其中以有色金属铅、锌、铜、钼、铋、钨、锡为主,其次是贵金属金、银等(潘卓,2011)。内生矿床在空间分布上具有明显的集群性,主要分布在黄沙坪、宝山、大坊等地。矿床的分布和岩浆岩关系密切,多产于岩体接触带及其附近。
图1 坪宝地区地质略图Fig.1 Geological sketch map of Huangshaping-Baoshan area
湖南宝山地区矿产地质调查项目完成了1∶5万桂阳幅460km2水系沉积物的采集,共采集样品1920件,每个采样点控制的汇水面积在0.125~0.25km2之间,平均采样密度4.17个/km2,样点分布于一级水系或二级水系上游、分枝口及沟口,采样物质以细砂为主,每个样品在采样点附近30~50m范围内进行多点(不少于3点)采集,合并为一个样品,样品重量过筛后不少于300g。分析测试方法采用 P-GF-AAS、ES、ICP-MS、离子选择性电极、HGAFS、CV-AFS等,分析地点在华北地矿局514地质队检测中心。
对研究区1920个水系沉积物样品分析了W、Sn、Mo、Bi、Cu、Pb、Zn、Au、Ag、As、Sb、Hg、Cr、Ni、Co、F、La、V、Mn等十九种元素,其异常下限的确定采用原始数据进行对数变换后连续剔除均值±3倍标准方差外数据,计算出对数数据算术平均值(X)、标准方差(S),由X+KS(K为置信度常数,取值范围1~2)后反对数求真数得到研究区各元素的异常下限值(T),结果如表1。
随着多元正态分布理论在地球化学数据处理中的应用,综合化探异常在地球化学找矿和成矿预测中的优势得到充分的发挥,已逐渐成为一种常规的化探异常处理方法。基于多元正态分布理论的多指标综合化探异常处理方法,实际上是多变量统计方法的直接推广与应用,根据其数学原理的不同,可以分为多种的具体方法,目前通常使用的是相关分析方法、聚类分析方法和因子分析方法(朴寿成等,2003;宋运红等,2008)。运用不同的分析方法确定的元素组合及组合化探异常也不尽相同。
2.3.1 相关分析
作为地质作用的微观结果,地球化学信息必然与地质信息相关连。相关分析是一种简单而直接的研究元素亲合性的方法。利用DGSInfo软件对研究区的19种元素进行相关分析,计算了各元素之间的相关系数,计算之前首先对原始数据进行标准化,计算结果见表2。
由表2可知:W、Cu、Bi元素呈显著正相关,Pb、Zn、Cu、W、Sb、Bi元素之间呈正强相关,Au、Ag、Sb元素之间呈明显正相关,Mn、Co、Ni之间相关性也较好。
2.3.2 R聚类分析
聚类分析以变量之间的相似程度为基础,将变量分成不同级别的类或点群,直观地对变量进行分类处理。R聚类分析可以了解个别变量之间以及各变量组合之间的亲疏程度;对研究区19种元素进行R聚类分析,得出R聚类分析图如图2所示。
据图2可见R=0.6为界可分四类。结合异常所处地质矿产条件分析其组合分布规律。
第一类 Pb、Zn、W、Cu、Au、Ag、Bi、Sb、As:其异常分布与岩体有关,Pb、Zn、W密切相关,反映Pb、Zn、W、Cu、Au、Ag、Bi、Sb、As 元素的富集主要与岩浆热液成矿作用有关,组合异常的出现是测区寻找Pb、Zn多金属矿床的重要地球化学找矿标志。
第二类 Mo、Sn、Cr、V、Hg:该类反映了中低温热液作用有关的矿化信息。
表1 异常下限统计一览表Table 1 Statistics of lower limits of anomalies
表2 相关系数矩阵Table 2 Correlation coefficient matrix
图2 R聚类分析图Fig.2 R-clustering analysis of stream sediment samples
第三类 F 与 La:与 Mo、Sn、Cr、V、Hg 和 Mn、Co、Ni均有一定的相关性,反映其成因复杂,既与岩浆活动有关,又受地层控制,为双控元素。
第四类Mn、Co、Ni:由谱系图上看,三者之间相关性好,与其它元素相关性较差,其富集主要与二叠系地层有关,与岩浆热液活动关系不密切,为层控元素。
2.3.3 因子分析
因子分析实际上是一种降维分析,是将存在复杂关系的较多变量依据某种内在联系生成几个综合因子的一种多元统计方法,每个因子指示出某种地质体上的元素共生组合及成因关系,从而对地球化学现象进行内在的成因分析,不同的元素组合反映不同的地质地球化学信息。本文应用DGSInfo软件对湘南坪宝地区1∶5万水系沉积物地球化学数据进行因子分析,得出表3,根据表3中所反映的不同元素组合来确定地球化学类型。
由表3可知,19种元素总体可以分为两种元素组合,其特征如下:第一组F1因子代表的是Au、Ag、Sn、Cu、Pb、Zn、W、Sb、Bi等元素组合,该组合包含了高、中、低温元素,其反映了与岩浆热液活动相关,从元素的亲和性看,Pb、Zn、Cu、Ag等表现为亲硫性,而W表现为亲氧性,常常易于在中酸性岩体中富集。因此,F1因子应与研究区的中酸性岩体有关。第二组 F2 因子代表的是 Co、Cr、Ni、V、Mn 等元素组合,其富集与地层关系密切,根据野外实地考察,二叠系地层多见堆积型铁锰矿富集,经化学分析同时伴生 Co、Cr、Ni、V 等元素,因此,F2 因子应与地层关系密切。
利用DGSInfo软件对水系沉积物样品分析结果进行处理,首先对原始数据进行展点(带属性),然后对元素逐个进行数字特征计算,再对数据进一步网格数据等值线处理,设置等值线套区、绘制色阶、光滑等值线、等值线层值等,得到单元素异常图;根据不同的算法进行不同处理,特别是对异常下限值的算法不同得出的异常图也是有所不同的。
表3 初始因子矩阵结果Table 3 Matrix results from initial factors
通过对19种元素进行相关分析、R聚类分析和因子分析结合研究区内地质特征可以基本得出这样一个结论,以 Pb、Zn、Cu、W、Bi、Sb 等 6 种元素为主做为一组,以 Au、Ag、Sb等 3种元素为主做为一组,以Mn、Co、Ni等3种元素为主做为一组;该研究区以寻找Pb、Zn、Cu多金属矿为主,本文选取Pb、Zn、Cu、W、Bi、Sb 为主的 6 种元素进行分析研究。图 3为 Pb、Zn、Cu、W、Bi、Sb 的异常剖析图,由剖析图既可以直观看出单元素异常分布空间,又可以看出多种元素的分布规律以及共生关系或伴生关系。
对图3异常剖析图取得如下四点认识:①从剖析图上可以清楚看到 Pb、Zn、Cu、W、Bi、Sb 六种元素异常区基本相吻合,Pb、Zn异常区面积略大,并且具有明显的浓度集中区,特别是图幅东边麦子坪-宝山-黄沙坪一带异常吻合良好,说明各元素之间共生或伴生关系密切。②元素异常主要集中在麦子坪-宝山-黄沙坪一带呈南北向展布,并且高度吻合,与研究区内南北向构造相吻合,因此异常应受南北向构造控制。③麦子坪-宝山-黄沙坪一带岩浆活动强烈,多呈岩脉产出,多为浅成-超浅成的中酸性、酸性英安斑岩、石英斑岩、花岗斑岩、花岗闪长斑岩等,异常应与岩浆活动有关。④区内已知矿床有宝山、黄沙坪、大坊等,异常区与区内已知矿床吻合良好。
通过对元素进行相关分析、R聚类分析、因子分析以及元素异常剖析,基本可以确立元素的组合特征,生成元素组合异常图,圈定元素组合异常;对区内寻找Pb、Zn、Cu多金属矿提供指导。在图4中共圈定了8个组合异常区,结合区内地层、构造、岩体特征,运用GIS叠加功能对区内组合异常区进行分析,并对异常区进行了野外实地调查,区内Ⅳ、Ⅵ、Ⅷ号组合异常区分别与大坊、宝山、黄沙坪等已知矿床相吻合;Ⅰ、Ⅱ、Ⅴ号组合异常区出露地层以二叠系为主,区内近南北向断裂构造以及南北向褶皱构造发育,未见岩体出露,但异常区曾经民采堆积型铁锰矿,故该三号组合异常可能受矿山开采影响。Ⅲ、Ⅶ号组合异常区成矿地质条件与Ⅵ号组合异常区宝山矿床极为相似,组合异常区内地层以石炭系为主,区内NE向褶皱构造以及近SN向、NW向断裂构造发育,区内岩浆岩较为发育,主要为燕山早期超浅成中酸性小岩体;经对Ⅲ、Ⅶ号异常区内异常高值点采样位置实地调查,认为矿山开采对其影响不大,路线地质调查中地表见褐铁矿化、方解石化、大理岩化、灰岩重结晶等围岩蚀变,具有良好的找矿前景,下一步工作可选取Ⅲ、Ⅶ号组合异常区开展土壤剖面测量工作,进行异常检查。
在水系沉积物异常的基础上,经过野外实地调查研究,结合研究区地层、构造、岩浆岩等地质特征,通过多元信息叠加对区内异常进行综合分析研究, 并与研究区内宝山大型铅锌多金属矿床进行对比分析,圈定了麦子坪-金子岗-梧桐铅锌多金属矿找矿靶区和禾茂村-马鞍岭铅锌多金属矿找矿靶区。
通过对宝山铅锌多金属矿床分析,认为研究区内寻找铅锌多金属矿床的几个标志为:①矿床多处于断裂交叉部位或褶皱背斜部位,其上有覆盖层。②有岩体存在,并与中酸性岩体关系密切。③ 地层主要为石炭系灰岩,主要为石磴子灰岩、梓门桥灰岩以及壶天群灰岩。④围岩蚀变较强,主要有矽卡岩化、钾长石化、绿帘石化、萤石化、硅化、绿泥石化、方解石化、绢云母化、高岭石化、大理岩化、灰岩重结晶等。⑤ 为多元素组合异常,异常值高,且异常分带清楚,浓集中心明显。
麦子坪-金子岗-梧桐铅锌多金属找矿靶区位于黄沙坪-宝山成矿带的北段,出露地层以石炭系为主,分别为天鹅坪组砂岩、粉砂岩;石磴子组灰岩;测水组砂页岩;梓门桥组灰岩;壶天群灰岩、白云岩等。区内NE向褶皱构造以及近南北向、北西向断裂构造发育。区内岩浆活动较为强烈,地表多处出露燕山早期超浅成中酸性花岗闪长斑岩呈岩脉产出。围岩蚀变较强,主要有硅化、绿泥石化、方解石化、大理岩化、灰岩重结晶等。区内重力异常属炎陵-蓝山重力梯度带南侧的舒缓部位,剩余重力异常图上反映出区内存在金子岗-昭金祠南北向重力低异常呈带状展布。航磁异常有C-86-25号金子岗磁异常,地磁测量在本区为正负异常紧密相伴的线状异常(秦葆瑚,1984;曾钦旺等,2002)。水系沉积物异常为Ⅲ号 Pb、Zn、Cu、W、Bi、Sb 元素组合异常(图4),组合元素间套合良好,异常值高,且异常分带清楚,浓集中心明显,呈近南北向展布。
本区地层岩性、构造、岩浆岩对成矿较为有利,成矿地质条件与宝山矿床极为相似,是区内寻找“宝山式”矿床的极为有利地段,根据区内地球物理重力、磁性异常特征及水系沉积物组合异常和区内成矿地质条件,推测该区可能隐伏有矽卡岩型铅锌铜矿(化)体存在。
禾茂村-马鞍岭铅锌多金属矿找矿靶区位于坪宝地区中部,出露地层有泥盆系、石炭系、二叠系以及第四系,其岩性主要为锡矿山组灰岩;马栏边组灰岩;天鹅坪组砂岩、粉砂岩;石磴子组灰岩;测水组砂页岩;梓门桥组灰岩;壶天群灰岩、白云岩;栖霞组灰岩等。区内构造为近SN紧闭型褶皱及SN向逆冲断层,次级构造为NW向的横断层与近SN向逆冲断层相切,形成“井”字型构造格架。区内岩浆活动强烈,出露岩体为燕山早期闪斜煌斑岩、花岗斑岩。区内围岩蚀变强烈,主要有硅化、碳酸盐化、绿泥石化、灰岩重结晶等。区内重力异常属炎陵-蓝山重力梯度带南侧的舒缓部位,剩余重力异常图上为正负异常边缘的零值区和负值区。航磁有C-85-21号异常,存在两个封闭的正异常和一个负异常。地磁测量在本区中部为低缓的局部异常,南部为走向SN作EW向排列的正负交替异常群。水系沉积物异常为Ⅶ号 Pb、Zn、Cu、W、Bi、Sb 元素组合异常(图4),Cu、Pb、Zn异常面积较大,异常值高,且异常分带清楚,浓集中心明显,呈NE向带状分布。
图4 Pb、Zn、Cu、W、Bi、Sb 组合异常图Fig.4 Composite anomalies of elements Pb,Zn,Cu,W,Bi and Sb
本区具有与黄沙坪矿床较为相似的地层岩性、构造、岩浆岩条件,区内物化探异常明显,套合性较好,其特征与黄沙坪矿床较为相似,故该区是寻找隐伏高-中温热液交代型和裂隙充填交代型铅锌多金属矿床的有利地段。
湘南坪宝地区是寻找铅锌多金属矿资源潜力较大地区,成矿地质条件优越。本文通过DGSInfo软件对研究区内水系沉积物样品分析数据进行相关分析、R聚类分析和因子分析确定元素组合特征,圈定了 8 个 Pb、Zn、Cu、W、Bi、Sb 元素组合异常区,异常区内元素套合性较好,异常强度高,浓集中心明显。结合研究区成矿地质条件,运用GIS多元信息叠加功能对区内组合异常区进行综合分析研究,圈定了麦子坪-金子岗-梧桐铅锌多金属找矿靶区与禾茂村-马鞍岭铅锌多金属矿找矿靶区,为区内开展下一步工作寻找铅锌多金属矿提供有意义的帮助。
Bai Dao-yuan,Ma Tie-qiu,Wang Xian-hui,Zhang Xiao-yang,Chen Bi-he.2008.Progress in the study of Mesozoic tectono-magmatism and mineralization in the central segment of the Nanling Mountains-Summary of major achievements of the 1 ∶250,000 geological survey in southeastern Hunan[J].Geology in China,35(3):436-455(in Chinese with English abstract)
Fu Xu.2012.A study on diagenetic and metallogenic depths of the Huangshaping lead-zinc deposit in hunan province[J].Beijing:China U-niversity Geosciences:1-78(in Chinese with English abstract)
Gong Shu-qing.2006.Discussion on the deep lead-zinc resource potentials in Huangshaping mine[J].Mining research and development,27(4):83-86(in Chinese with English abstract)
Hu Zhi-ke.2002.Discussion on marks of prospecting information and effective prospecting methods in Huangshaping Baoshan oref ield[J].Geology and Mineral Resources of South China,(4):60-64(in Chinese with English abstract)
Li Jian-zhong,Zhang Yi-jun,Cai Xin-hua,Xu Hui-chang.2005.Analysis on the prospecting ore of potentiality in Huangshaping Pb-Zn deposit,Hunan province[J].Geology and Mineral Resources of South China,84(4):23-26(in Chinese with English abstract)
Li Yu-qin,Shen Heng-li,Wang Xue-zhen.2011.Geochemical characteristic of stream sediment survey and metallogenic prognosis of Dulan area Qinghai province china[J].Acta Minalogica Sinica,31(3):615-620(in Chinese with English abstract)
Liu Wu-hui.2007.Mineralization mechanism and prediction of Huangshaping Pb-Zn polymetallic deposit[D].changsha:Central South University:1-149(in Chinese with English abstract)
Nanjing University,Department of Geology.1981.Granites of different geological periods and their relationship with mineralization,South China[M].Beijing:Science Press:214-216(in Chinese)
Pan Zhuo.2011.The geological characteristics,ore-controlling factors and deposit genesis of the W-Mo polymetallic deposit in the south of Huangshaping mine,Hunan Province[D].Changsha:Central South University:1-52(in Chinese with English abstract)
Piao Shou-cheng,Liu Shu-tian.2003.Features of geochemical anomalies and ore prediction in Xiaoshiren gold deposit,Jilin province[J].Geology and prospecting,(2):26-29(in Chinese with English abstract)
Qiao Yu-sheng,Gong Shu-qing.2010.The deep ministry of huangshaping-Baoshan distrut looks for the mine latent capacity to analyse[J].China Science And Technoloy,(24):275-276(in Chinese with English abstract)
Qin Bao-hu.1984.The interpretations of regional geavity and magnetic anomal ies and their applications of minerogenetic prognoses in Southern Hunan[J].Hunan Geology,(2):28-33(in Chinese with English abstract)
Song Yun-hong,Jia Da-cheng,Li Jian,Bai Xue-shan,Che Yan-an,Li Yong-chan.2008.Characteristics of the anomalies and geochem ical exploration data processing of the Sangou area in chengde,Hebei province[J].Journa l of Jilin Un iversity(Earth Sc ience Edition),38:167-169(in Chinese with English abstract)
Tian Xu-feng,Gong Shu-qing.2007.Future prospect of tungsten and molybdenum resources in Huangshaping lead and zinc mine[J].Journal of University of South China(Science and Technology),21(2):110-114(in Chinese with English abstract)
Tong Qian-ming.1990.Analysis of metallogenic regularity and concealed deposits prediction of Baoshan-Xianghualing lead zinc polymetallic deposits[R].changsha:Geological research institute of Hunan Province:1-183(in Chinese)
Xia Xiang-biao,Zheng Lai-lin.2009.Geochemical character istics and ore search prospects of Puchiya area in Tibet based on stream sediment survey[J].Geophysical and Geochemical Exploration,33(6):626-629(in Chinese with English abstract)
Xu Yi-ming,Gong Shu-qing ,Jiang Yuan-cheng,Lei Ze-heng,Li Yusheng.2007.Hunan Huangshaping lead zinc mine of the deep analysis of the prospects for lead zinc mine resources[J].Geology and prospecting,43(1):38-43(in Chinese with English abstract)
Xue Shun-rong,Xiao Ke-yan,Ding Jian-hua.2008.The application of geochemical information to ore exploration in Shangri-laarea yunnan province[J].Geophysical and Geochemical Exploration,32(5):537-540(in Chinese with English abstract).
Yang Shi-yi,Liu Hou-qun,Zhang Xiu-lan,Chen chang-jiang.1986.A study of two type porphyry plutons and related mineralization in Pingbao district,South Hunan[J].Theses memoir of Nanling Geology and mineral resources,(2):157-167(in Chinese)
Zeng Qin-wang,Rao Jia-rong,Luo Jian-lan.2002.The effects of applying integra ted methods to the discovery of the Liutang concealed leadzinc deposit[J].Geophysical and Geochemical Exploration,26(3):179-191(in Chinese with English abstract)
Zhang Zai-tian,Pan Jia-yong,Xia Fei,Liu Guo-qi.2013.Geochemical data processing of stream sediment and target areas prediction for mentallogenic of Jiangxi yuhuashan.[J].Computing Techniques for Geophysical and Geochemical Exploration,35(1):86-92(in Chinese with English abstract)
[附中文参考文献]
柏道远,马铁球,王先辉,张晓阳,陈必河.2008.南岭中段中生代构造-岩浆活动与成矿作用研究进展[J].中国地质,35(3):436-455
付 旭.2012.湖南黄沙坪铅锌矿成岩成矿深度估算[D].北京:中国地质大学:1-78
龚述清.2006.黄沙坪铅锌矿深部铅锌资源潜力探讨[J].矿业研究与开发,27(4):83-86
胡志科.2002.坪宝矿田找矿信息标志及有效找矿方法探析[J].华南地质与矿产,(4):60-64
李建中,张怡军,蔡新华,徐惠长.2005.湖南省黄沙坪铅锌矿区找矿潜力分析[J].华南地质与矿产,84(4):23-26
李玉芹,沈恒丽,王学贞.2011.都兰地区水系沉积物测量地球化学特征及找矿预测[J].矿物学报,31(3):615-620
刘悟辉.2007.黄沙坪铅锌多金属矿床成矿机理及其预测研究[D].长沙:中南大学:1-149
南京大学.1981.华南不同时代花岗岩类及其与成矿的关系[M].北京:科学出版社:1-395
潘 卓.2011.湖南黄沙坪矿区南部钨钼多金属矿床地质特征、控矿因素及矿床成因[D].长沙:中南大学:1-52
朴寿成,刘树田.2003.吉林小石人金矿地球化学异常特征及成矿预测[J].地质与勘探,39(2):26-29
乔玉生,龚述清.2010.坪宝地区深部找矿潜力分析[J].中国科技纵横,(24):275-276
秦葆瑚.1984.湘南区域重磁异常的地质解释及其在成矿预测中的应用[J].湖南地质,(2):28-33
宋运红,贾大成,李 健,白雪山,车延安,李永占.2008.河北承德三沟地区化探数据处理与异常特征[J].吉林大学学报(地球科学版),38:167-169
田旭峰,龚述清.2007.黄沙坪铅锌矿深部钨钼资源远景浅析[J].南华大学学报(自然科学版),21(2):110-114
童潜明.1990.宝山-香花岭铅锌为主的多金属矿成矿规律及隐伏矿床预测研究[R].长沙:湖南省地质研究所:1-183
夏祥标,郑来林.2009.普迟亚地区水系沉积物测量地球化学特征及找矿方向[J].物体与化探,33(6):626-629
许以明,龚述清,江元成,雷泽恒,李玉生.2007.湖南黄沙坪铅锌矿深边部找矿前景分析[J].地质与勘探,43(1):38-43
薛顺荣,肖克炎,丁建华.2008.香格里拉地区勘查地球化学信息找矿应用[J].物探与化探,32(5):537
杨世仪,刘姤群,张秀兰,陈长江.1986.湘南坪宝地区两类斑岩及其成矿作用研究[J].南岭地质矿产文集,(2):157-167
曾钦旺,饶家荣,骆检兰.2002.综合方法在柳塘隐伏铅锌矿床的应用效果[J].物探与化探,26(3):179-191
张再天,潘家永,夏 菲,张 勇,刘国奇.2013.江西玉华山地区水系沉积物化探数据处理及找矿靶区划分[J].物探化探计算技术,35(1):86-92