云南省富宁县大湾金矿土壤地球化学异常评价及找矿方向
|
摘要
微细浸染型金矿是目前国内外重要的金矿床类型之一,目前,我国微细浸染型金矿累计资源量已超过1000t,该类型矿床也己经成为中国重要的黄金资源。滇黔桂接壤地区是20多年来发现的一大金矿成矿集中区,其金矿类型主要为微细粒浸染型。研究区大湾金矿位于富宁县,地处滇黔桂金三角区,由于富宁地区金矿已开采多年,金矿资源已日渐枯竭,急需在矿床周边继续寻找可持续开采的金矿床。
通过收集大湾地区基础地质资料,对矿区土壤进行1:10000地球化学土壤测量,应用聚类分析、主因子分析等传统数学地质方法研究各元素的共生组合规律,并利用主因子得分异常来表征主要成矿元素的组合异常。对比传统方法和求和法的多重分形方法来研究区域内主要成矿元素地球化学异常场的地球化学异常下限,并对研究区域内地球化学异常做出评价,结合地质特征探讨其找矿前景。研究结果表明:
(1)系统聚类分析结果与因子分析一致。Cu、Co、Ni同属一类,与因子分析的F1因子相对应;As、Au、Sb属一类,与因子分析的F2因子相对应;Hg单独为一类,与因子分析的F3因子相对应。
(2)本区的成矿元素组合主要为低温元素As、Au、Sb, Au是主成矿元素,Sb、Hg可能为伴生元素。
(3)研究区内7种元素应用多重分形方法所确定的下限值均高于用传统方法计算的异常下限,Au元素不具有多重分形特征。具有多重分形特征的元素,可以利用分形法求该元素的异常下限。
(4)对该地区异常图的研究表明,分形法确定的异常面积集中,且与构造、矿化点联系紧密。
(5)研究区金异常区域主要为三个区域:龙三盘以西、龙三盘以南及那桑,异常面积分别为1.0、0.38、0.5平方公里。沿着区内北西向的断层分布,说明金矿化的形成可能与这条贯穿该区域的断层有相当密切的关系。
(6)Au、As、Sb综合异常区域与金异常区的分布相一致,异常面积分别为2.5、0.68、0.87平方公里。这与大湾地区构造发育,基性岩活动强烈,具有较好的成矿条件有关;主成矿元素的异常与矿化点联系紧密、金矿化带往北西、南东均有延伸趋势,远景较好。
Carlin type gold deposit is one important type of gold deposit at home and abroad. At present, Carlin type gold deposit resources have accumulated more than 1000t, the type of deposit also has become an important gold resource. Yunnan-Guizhou-Guangxi border area is a major gold deposits which was found in concentrated areas in 20 years, the main types of gold deposit is Carlin type. Dawan area is located in Funing gold deposit in Yunnan, Guizhou Golden Triangle area, the Fu Ning gold deposit has been mined area for many years,the resources of gold deposit has been an increasing depletion, the urgent need to continue to look for sustainable exploitation of deposits around the gold deposit.
Through the collection of basic geological data of Dawan Area, surveying the 1:10000 geochemical soil around the mine, applying cluster analysis, principal factor analysis and other traditional methods of mathematical geology of the symbiotic combination of the elements regularity, and using the main factor scores to represent the exception forming elements of the combination of abnormalities. Compared to the traditional method and multifractal statistical methods of summation to the study area forming elements within the field of geochemical anomalies in geochemical anomaly threshold, and geochemical anomalies within the study area to make assessment of the geological features of its exploration Prospects. The results show that:
(1) The results of cluster analysis and factor analysis is consistent, Cu, Co, Ni is part of the same class, and factor analysis of the F1 factor corresponds; As, Au, Sb, Hg is part of the same class, and factor analysis of the F2 factor corresponds.
(2) The combination of this area is mainly forming elements As, Au, Sb, Hg, As-Au is temperature element, and the main ore-forming elements, Sb, Hg may be the associated elements.
(3) 7 elements in the study area by using multifractal methods is higher than anomaly threshold of the traditional method of anomaly threshold, Au element does not have multifractal characteristics. We can use fractal Method anomaly threshold of the element of multifractal characteristics .
(4) The study of anomaly map in the area shows that the abnormal area which was determined by the fractal method is small, and associated closely with the structure, mineralization point.
(5) Major regional gold anomaly in the study area is divided into three regions: the west of the Long Sanpan, the south of Long Sanpan and Nasang that abnormal area are 1.0,0.38,0.5 square kilometers. Along northwest the fault distribution, indicating that the formation of gold mineralization may be related to this fault runs through the region have very close relationship.
(6) The distribution of Au, As, Hg, Sb integrated abnormal area and gold anomalies abnormal area is consistent. The abnormal area are 2.5,0.68,0.87 square kilometers. Related to good mineralization conditions,Dawan construction fully developed, basic rock strongly actived,; The anomalies of the main mineralization elements and mineralization points are closely connected, gold mineralization belt are an extension of the trend in northwest,southeast . The Prospect is good
通过收集大湾地区基础地质资料,对矿区土壤进行1:10000地球化学土壤测量,应用聚类分析、主因子分析等传统数学地质方法研究各元素的共生组合规律,并利用主因子得分异常来表征主要成矿元素的组合异常。对比传统方法和求和法的多重分形方法来研究区域内主要成矿元素地球化学异常场的地球化学异常下限,并对研究区域内地球化学异常做出评价,结合地质特征探讨其找矿前景。研究结果表明:
(1)系统聚类分析结果与因子分析一致。Cu、Co、Ni同属一类,与因子分析的F1因子相对应;As、Au、Sb属一类,与因子分析的F2因子相对应;Hg单独为一类,与因子分析的F3因子相对应。
(2)本区的成矿元素组合主要为低温元素As、Au、Sb, Au是主成矿元素,Sb、Hg可能为伴生元素。
(3)研究区内7种元素应用多重分形方法所确定的下限值均高于用传统方法计算的异常下限,Au元素不具有多重分形特征。具有多重分形特征的元素,可以利用分形法求该元素的异常下限。
(4)对该地区异常图的研究表明,分形法确定的异常面积集中,且与构造、矿化点联系紧密。
(5)研究区金异常区域主要为三个区域:龙三盘以西、龙三盘以南及那桑,异常面积分别为1.0、0.38、0.5平方公里。沿着区内北西向的断层分布,说明金矿化的形成可能与这条贯穿该区域的断层有相当密切的关系。
(6)Au、As、Sb综合异常区域与金异常区的分布相一致,异常面积分别为2.5、0.68、0.87平方公里。这与大湾地区构造发育,基性岩活动强烈,具有较好的成矿条件有关;主成矿元素的异常与矿化点联系紧密、金矿化带往北西、南东均有延伸趋势,远景较好。
Carlin type gold deposit is one important type of gold deposit at home and abroad. At present, Carlin type gold deposit resources have accumulated more than 1000t, the type of deposit also has become an important gold resource. Yunnan-Guizhou-Guangxi border area is a major gold deposits which was found in concentrated areas in 20 years, the main types of gold deposit is Carlin type. Dawan area is located in Funing gold deposit in Yunnan, Guizhou Golden Triangle area, the Fu Ning gold deposit has been mined area for many years,the resources of gold deposit has been an increasing depletion, the urgent need to continue to look for sustainable exploitation of deposits around the gold deposit.
Through the collection of basic geological data of Dawan Area, surveying the 1:10000 geochemical soil around the mine, applying cluster analysis, principal factor analysis and other traditional methods of mathematical geology of the symbiotic combination of the elements regularity, and using the main factor scores to represent the exception forming elements of the combination of abnormalities. Compared to the traditional method and multifractal statistical methods of summation to the study area forming elements within the field of geochemical anomalies in geochemical anomaly threshold, and geochemical anomalies within the study area to make assessment of the geological features of its exploration Prospects. The results show that:
(1) The results of cluster analysis and factor analysis is consistent, Cu, Co, Ni is part of the same class, and factor analysis of the F1 factor corresponds; As, Au, Sb, Hg is part of the same class, and factor analysis of the F2 factor corresponds.
(2) The combination of this area is mainly forming elements As, Au, Sb, Hg, As-Au is temperature element, and the main ore-forming elements, Sb, Hg may be the associated elements.
(3) 7 elements in the study area by using multifractal methods is higher than anomaly threshold of the traditional method of anomaly threshold, Au element does not have multifractal characteristics. We can use fractal Method anomaly threshold of the element of multifractal characteristics .
(4) The study of anomaly map in the area shows that the abnormal area which was determined by the fractal method is small, and associated closely with the structure, mineralization point.
(5) Major regional gold anomaly in the study area is divided into three regions: the west of the Long Sanpan, the south of Long Sanpan and Nasang that abnormal area are 1.0,0.38,0.5 square kilometers. Along northwest the fault distribution, indicating that the formation of gold mineralization may be related to this fault runs through the region have very close relationship.
(6) The distribution of Au, As, Hg, Sb integrated abnormal area and gold anomalies abnormal area is consistent. The abnormal area are 2.5,0.68,0.87 square kilometers. Related to good mineralization conditions,Dawan construction fully developed, basic rock strongly actived,; The anomalies of the main mineralization elements and mineralization points are closely connected, gold mineralization belt are an extension of the trend in northwest,southeast . The Prospect is good
引文
[1]吴锁平.黄金:金属之王[J].中国金属通报,2010,41,14-17.
[2]张宏宾,黄映聪,赵存法,葛俊逸.云南省革档金矿流体包裹体地球化学特征[J].有色金属,2010,3(62):25-30.
[3]杨昌华,刘星,罗荣生.广南堂上金矿多种金矿化的发现[J].云南地质,2006,25(3):286-290.
[4]罗刚,杨小峰.云南广南地区老寨湾微细粒浸染型金矿床地质特征与成矿规律[J].2010,29(9),1362-1370.
[5]吕世琨,李爱平,杨绍文,雷跃增,赵禹.滇东南那能金矿带之断裂控矿[J].云南地质,2003,23(1):47-55.
[6]朱江.卡林型金矿找矿案例分析及思考[J].2010,24(4),354-358.
[7]Teal L.and Jackson M.,Geologic overview of the Carlin Trend gold deposits and descriptions of recent deep discoveries.STGNewsletter,1997,(31):1-25.
[8]周余国,刘继顺等.从滇黔桂“金三角”区域地层地球化学演化特征探讨卡林型金矿的物质来源[J].2009,16(2),199-208.
[9]朱赖民,胡瑞忠,刘显凡,刘家军,金景福,何明友.关于黔西南微细浸染型金矿床成因的一些初步认识[J」.矿产与地质,1997(a),11(5).
[10]刘显凡,苏文超,朱赖民.滇黔桂微细浸染型金矿深源流体成矿机理探讨[J].地质与勘探,1999,(01).
[11]杜琦.以我国一些斑岩铜矿为例试论斑岩矿床成矿与地层(围岩)的关系[J].矿床地质,1984,3(2):21-26.
[12]钟文旭,计政科,何永平.新疆区域铀矿地球物理场特征分析[J].铀矿地质,1985,1(5):58-61.
[13]陈衍景,邓健,富士谷胡受奚.上宫金矿地质地球化学特征和成矿模式[J].黄金地质科技,1993,4,1-9.
[14]郭保健,李永峰,王志光,叶会寿.熊耳山Au-Ag-Pb-Mo矿集区成矿模式与找矿方向[J].地质与勘探,2005,41(5):43-47.
[15]王学求.矿产勘查地球化学:过去的成就与未来的挑战[J].地学前缘,2003,10(1),239-247.
[16]罗正传.土壤地球化学测量在智利阿尔法矿区的应用效果[J].物探与化探.2010,34(4),473-475.
[17]熊亮,朱杰勇,朱林生,李侃等.土壤地球化学测量在扬子地台西缘会泽铅锌矿带筇竹寺组中的应用[J].2010,34(4),467-471.
[18]胡长征.地球化学土壤测量在找矿中的应用[J].科技创新导报,2010,24,88.
[19]中华人民共和国地质矿产行业标准(DZ/T 0011-91),地球化学普查规范.
[20]张彪.内蒙古鲁北地区地球化学异常识别与综合评价[D].中国地质大学,地球探测与信息技术系,2008.
[21]Garrett,R.G..A cry from the heart[J].Explore,1989,66:18-20.
[22]Govett,G.J.S.,Goodfellow W.D.,Chapman R.P.and Chork C.Y..Exploration Geochemistry distribution of elements and recognition of anomalies [J]. Math, 1975, 7:415-446.
[23]Miesch A.T..Estimation of the geochemical threshold and its statistical significance[J] .Geochemistry Exploration,1981,16:49-76.
[24] Sinclair A.J..A foundmental approach to threshold estimation in exploration geochemistry[J].Geochemistry Exploration,1991,41:1-22. Geochemistry Exploration,1989,32:355-357.
[25]Stainly C.R.and Sinclair A.J..Comparison of probability plts and gap statistics in the selection of threshold for exploration geochemistry data[J].1989,32:355-357.
[26]张济忠.分形[M].北京:清华大学出版社,1999.
[27]刘玲玲,孔正义,张超.分形理论在地震波表征中的应用[J].2010,36(5),157-158.
[28]Mandelbrot,B.B.1982,the Fratal Geometry of Nature.Freeman[M],San Francisco.
[29]Shaw,H.R.and Gatiier,A.E..1984.Empirical laws of order among rivers,faults,and earthquake.U.S. Geological Surgery Open File Report,84-356:51.
[30]Shaw,H.R..1986.On the graphical interpretation of Paleaseismic data U.S. Geological Surgery Open File Report,86-394:92.
[31]Chiles,J.P..1988.Fractal and statistical method for modeling of a fracture network.Math.Geol.,20:631-654.
[32]陈颙.分形与分维在地球科学中的应用[M].北京:学术期刊出版设.1989.
[33]李长江,徐有浪,蒋叙良.论矿床的分形性质[J].浙江地质,1994,10(2):25-32.
[34]李长江,麻土华.化探数据处理新技术[J].地质找矿论丛,1997,12(4):57-64.
[35]A.Ford,T.G.Bleinkinsop. Combining fractal analysis of mineral deposit clustering with weights of evidence to evaluate patterns of mineralization: Application to copper deposits of the Mount Isa Inlier, NW Queensland, Australia [J]. Ore Geology Reviews,2008,6(33):435-450.
[36]Tureotte D.L.,Fractal in geology and geophysies. pure Appl Geophys,1989,131:171-196.
[37]孟宪国,赵鹏大.地质数据的分形结构[J].地球科学—中国地质大学学报,1991,16(2),207-211.
[38]杜佑仁.分形(Fractal)理论及在地学中的应用简介[J].地质科技情报,1991,10(1),69-73.
[39]谢淑云,鲍征宇.多重分形与地球化学元素的分布规律[J].地质地球化学,2003,31(3),97-101.
[40]Qiuming Cheng,Agterberg F.P..The separation of geochemical anomalies from backgroundby fractal methods[J].Geochem Explor.1994,51(2):109-130.
[41]Allegre C.J.and Lewin E.Scaling laws and geochemical distribution.Earth and Planetary Science Letter[J].1995,132:1-13.
[42]李长江,麻土华.矿产勘查中的分形、混沌与ANN[M].北京:地质出版社.1999.
[43]申维.分形混沌与矿产预测[M].北京:地质出版社.2002.
[44]刘方成,秦德先,燕永锋,李保珠,秦来勇,洪托.地质统计学及其在云南金宝山铂钯矿的应用[J].地质与资源.2005,14(4):293-295.
[45]肖斌,赵鹏大,侯景儒.现代地质统计学的新进展[J].世界地质,1999,18(9):81-91.
[46]夏邦栋.地质学概论[M].北京:高等教育出版社,1992.
[47]姚娟.云南老寨湾金矿床成矿物质来源分析及矿床成因探讨[D],成都:成都理工大学.2008.
[48]吕贻中,李保国.土壤学[M].北京:中国农业出版社,2006.
[49]夏邦栋.普通地质学[M].北京:地质出版社,1995.
[50]庞绪贵,姜相洪,战金成,朱继托,曹秀华,王增辉,崔元俊.山东黄河流域下游土壤粒度和深度试验成果[J].地质通报,2005,24(8):767-771.
[51]AllegreC.J.,LewinE..Scaling laws and geochemical distribution[J].Earth and Planetary science Letters,1995,132:1-13.
[52]宋廷山、葛金田.统计学—以excel为分析工具[J].北京:北京大学出版社,2009.
[53]邱永华,金爱武,杜亮亮,金晓春,宋艳冬.影响毛竹林冬春笋产量的因子分析——基于对遂昌县三仁乡的调查[J].竹子研究汇刊,2010,29(3),26-31.
[54]包赞.基于因子分析法评价安徽省各城市经济综合实力[J].经济论坛,2010,10,88-90.
[55]王华义,谭卫红,李春燕.钢铁企业高温作业工人血压影响因素的因子分析[J].职业与健康,2010,26(19),2171-2173.
[56]许平,雷芬丽,文晓涛,贺振.碳酸盐岩储层流体识别因子分析[J].油气地球物理,2010,8(4),1-4.
[57]谭荣波,梅晓仁.SPSS统计分析实用教程[M].北京:科学出版社,2007.
[2]张宏宾,黄映聪,赵存法,葛俊逸.云南省革档金矿流体包裹体地球化学特征[J].有色金属,2010,3(62):25-30.
[3]杨昌华,刘星,罗荣生.广南堂上金矿多种金矿化的发现[J].云南地质,2006,25(3):286-290.
[4]罗刚,杨小峰.云南广南地区老寨湾微细粒浸染型金矿床地质特征与成矿规律[J].2010,29(9),1362-1370.
[5]吕世琨,李爱平,杨绍文,雷跃增,赵禹.滇东南那能金矿带之断裂控矿[J].云南地质,2003,23(1):47-55.
[6]朱江.卡林型金矿找矿案例分析及思考[J].2010,24(4),354-358.
[7]Teal L.and Jackson M.,Geologic overview of the Carlin Trend gold deposits and descriptions of recent deep discoveries.STGNewsletter,1997,(31):1-25.
[8]周余国,刘继顺等.从滇黔桂“金三角”区域地层地球化学演化特征探讨卡林型金矿的物质来源[J].2009,16(2),199-208.
[9]朱赖民,胡瑞忠,刘显凡,刘家军,金景福,何明友.关于黔西南微细浸染型金矿床成因的一些初步认识[J」.矿产与地质,1997(a),11(5).
[10]刘显凡,苏文超,朱赖民.滇黔桂微细浸染型金矿深源流体成矿机理探讨[J].地质与勘探,1999,(01).
[11]杜琦.以我国一些斑岩铜矿为例试论斑岩矿床成矿与地层(围岩)的关系[J].矿床地质,1984,3(2):21-26.
[12]钟文旭,计政科,何永平.新疆区域铀矿地球物理场特征分析[J].铀矿地质,1985,1(5):58-61.
[13]陈衍景,邓健,富士谷胡受奚.上宫金矿地质地球化学特征和成矿模式[J].黄金地质科技,1993,4,1-9.
[14]郭保健,李永峰,王志光,叶会寿.熊耳山Au-Ag-Pb-Mo矿集区成矿模式与找矿方向[J].地质与勘探,2005,41(5):43-47.
[15]王学求.矿产勘查地球化学:过去的成就与未来的挑战[J].地学前缘,2003,10(1),239-247.
[16]罗正传.土壤地球化学测量在智利阿尔法矿区的应用效果[J].物探与化探.2010,34(4),473-475.
[17]熊亮,朱杰勇,朱林生,李侃等.土壤地球化学测量在扬子地台西缘会泽铅锌矿带筇竹寺组中的应用[J].2010,34(4),467-471.
[18]胡长征.地球化学土壤测量在找矿中的应用[J].科技创新导报,2010,24,88.
[19]中华人民共和国地质矿产行业标准(DZ/T 0011-91),地球化学普查规范.
[20]张彪.内蒙古鲁北地区地球化学异常识别与综合评价[D].中国地质大学,地球探测与信息技术系,2008.
[21]Garrett,R.G..A cry from the heart[J].Explore,1989,66:18-20.
[22]Govett,G.J.S.,Goodfellow W.D.,Chapman R.P.and Chork C.Y..Exploration Geochemistry distribution of elements and recognition of anomalies [J]. Math, 1975, 7:415-446.
[23]Miesch A.T..Estimation of the geochemical threshold and its statistical significance[J] .Geochemistry Exploration,1981,16:49-76.
[24] Sinclair A.J..A foundmental approach to threshold estimation in exploration geochemistry[J].Geochemistry Exploration,1991,41:1-22. Geochemistry Exploration,1989,32:355-357.
[25]Stainly C.R.and Sinclair A.J..Comparison of probability plts and gap statistics in the selection of threshold for exploration geochemistry data[J].1989,32:355-357.
[26]张济忠.分形[M].北京:清华大学出版社,1999.
[27]刘玲玲,孔正义,张超.分形理论在地震波表征中的应用[J].2010,36(5),157-158.
[28]Mandelbrot,B.B.1982,the Fratal Geometry of Nature.Freeman[M],San Francisco.
[29]Shaw,H.R.and Gatiier,A.E..1984.Empirical laws of order among rivers,faults,and earthquake.U.S. Geological Surgery Open File Report,84-356:51.
[30]Shaw,H.R..1986.On the graphical interpretation of Paleaseismic data U.S. Geological Surgery Open File Report,86-394:92.
[31]Chiles,J.P..1988.Fractal and statistical method for modeling of a fracture network.Math.Geol.,20:631-654.
[32]陈颙.分形与分维在地球科学中的应用[M].北京:学术期刊出版设.1989.
[33]李长江,徐有浪,蒋叙良.论矿床的分形性质[J].浙江地质,1994,10(2):25-32.
[34]李长江,麻土华.化探数据处理新技术[J].地质找矿论丛,1997,12(4):57-64.
[35]A.Ford,T.G.Bleinkinsop. Combining fractal analysis of mineral deposit clustering with weights of evidence to evaluate patterns of mineralization: Application to copper deposits of the Mount Isa Inlier, NW Queensland, Australia [J]. Ore Geology Reviews,2008,6(33):435-450.
[36]Tureotte D.L.,Fractal in geology and geophysies. pure Appl Geophys,1989,131:171-196.
[37]孟宪国,赵鹏大.地质数据的分形结构[J].地球科学—中国地质大学学报,1991,16(2),207-211.
[38]杜佑仁.分形(Fractal)理论及在地学中的应用简介[J].地质科技情报,1991,10(1),69-73.
[39]谢淑云,鲍征宇.多重分形与地球化学元素的分布规律[J].地质地球化学,2003,31(3),97-101.
[40]Qiuming Cheng,Agterberg F.P..The separation of geochemical anomalies from backgroundby fractal methods[J].Geochem Explor.1994,51(2):109-130.
[41]Allegre C.J.and Lewin E.Scaling laws and geochemical distribution.Earth and Planetary Science Letter[J].1995,132:1-13.
[42]李长江,麻土华.矿产勘查中的分形、混沌与ANN[M].北京:地质出版社.1999.
[43]申维.分形混沌与矿产预测[M].北京:地质出版社.2002.
[44]刘方成,秦德先,燕永锋,李保珠,秦来勇,洪托.地质统计学及其在云南金宝山铂钯矿的应用[J].地质与资源.2005,14(4):293-295.
[45]肖斌,赵鹏大,侯景儒.现代地质统计学的新进展[J].世界地质,1999,18(9):81-91.
[46]夏邦栋.地质学概论[M].北京:高等教育出版社,1992.
[47]姚娟.云南老寨湾金矿床成矿物质来源分析及矿床成因探讨[D],成都:成都理工大学.2008.
[48]吕贻中,李保国.土壤学[M].北京:中国农业出版社,2006.
[49]夏邦栋.普通地质学[M].北京:地质出版社,1995.
[50]庞绪贵,姜相洪,战金成,朱继托,曹秀华,王增辉,崔元俊.山东黄河流域下游土壤粒度和深度试验成果[J].地质通报,2005,24(8):767-771.
[51]AllegreC.J.,LewinE..Scaling laws and geochemical distribution[J].Earth and Planetary science Letters,1995,132:1-13.
[52]宋廷山、葛金田.统计学—以excel为分析工具[J].北京:北京大学出版社,2009.
[53]邱永华,金爱武,杜亮亮,金晓春,宋艳冬.影响毛竹林冬春笋产量的因子分析——基于对遂昌县三仁乡的调查[J].竹子研究汇刊,2010,29(3),26-31.
[54]包赞.基于因子分析法评价安徽省各城市经济综合实力[J].经济论坛,2010,10,88-90.
[55]王华义,谭卫红,李春燕.钢铁企业高温作业工人血压影响因素的因子分析[J].职业与健康,2010,26(19),2171-2173.
[56]许平,雷芬丽,文晓涛,贺振.碳酸盐岩储层流体识别因子分析[J].油气地球物理,2010,8(4),1-4.
[57]谭荣波,梅晓仁.SPSS统计分析实用教程[M].北京:科学出版社,2007.