山东省兰陵县草主山矿区铁矿床地质特征和找矿远景分析

2017-12-06 08:11冯玺平边丛林刘超高燕
山东国土资源 2017年12期
关键词:黑云磁铁铁矿

冯玺平,边丛林,刘超,高燕

(山东省鲁南地质工程勘察院,山东 兖州 272100)

山东省兰陵县草主山矿区铁矿床地质特征和找矿远景分析

冯玺平,边丛林,刘超,高燕

(山东省鲁南地质工程勘察院,山东 兖州 272100)

兰陵县草主山矿区属苍峄铁矿带一部分,区内赋存有全隐伏的“鞍山式”铁矿床,经勘查共圈定了南北2个主矿带,5个矿体。矿体赋存于泰山岩群山草峪组中,其形态受基底褶皱控制,具有平行延伸的特点,经过以往地质勘查工作,基本查明了-800m以浅矿体的产出位置、形态、产状、规模、分布、厚度、矿石质量等。矿体沿走向倾向尚未封闭,有继续延伸的趋势,并且矿石品位和矿体厚度未有减弱趋势,由此可得结论其-800m以深仍有很大的找矿前景,满足深部找矿的要求,总结其成矿规律、找矿标志,建立找矿模式,可为进一步勘查提供可靠的地质参考依据。

兰陵县;草主山矿区;“鞍山式”铁矿;找矿远景;成矿规律

草主山铁矿区位于兰陵县城西约25km,西距枣庄市约20km,区内公路交通发达。区内铁矿地质调查工作自20世纪30年代就有文字记载,勘查历史悠久,资源丰富,工作价值大,其成因类型被确定为沉积变质型“鞍山式”铁矿[1]。山东省鲁南地质工程勘察院先后对草主山矿区、刘家庄-幸福岭矿段、草主山南矿段开展了铁矿详查工作,并提交详查报告。结合以往勘查工作以及与鲁西地区同类型铁矿,如汶上东平铁矿带的东平县大牛铁矿进行类比分析,发现矿体深部仍有很大找矿前景[2-4]。经过深部及外围详查*山东省鲁南地质工程勘察院,山东省兰陵县草主山矿区深部及外围铁矿详查,2014年。,发现矿体沿走向倾向均未封闭,具有继续延伸的趋势,该矿区深部仍具有良好的资源潜力。

1 区域地质背景

1—第四系;2—奥陶系;3—寒武系;4—青白口系;5—南华系;6—山草峪组;7—黑云石英闪长岩;8—断层及产状;9—矿带;10—地质界线;11—角度不整合地质界线;12—平行不整合地质界线图1 兰陵县草主山铁矿区区域地质简图

2 矿区地质

2.1 地层

矿区地层发育有新太古代泰山岩群山草峪组,岩性以黑云变粒岩为主,局部夹磁铁角闪石英岩、磁铁石英角闪岩、黑云角闪片岩等。发育新元古代青白口纪土门群黑山官组,南华纪土门群二青山组、佟家庄组;古生代寒武纪长清群李官组、朱砂洞组、馒头组,地层总体走向近SN或NNW,NNE;倾向NEE或SEE,倾角8°~15°。新太古代泰山岩群山草峪组是区内沉积变质型(鞍山式)铁矿的赋存层位[9]。

2.2 构造

矿区构造特征与区域构造相似,按构造形式可分为基底构造和盖层构造2种,两者对铁矿床均有不同程度的影响。基底泰山岩群变质地层经过强烈的变质变形作用形成了一系列紧密的复式褶皱,对矿体的形态、产状等影响较大。区内主要发育轴向近东西的太白向斜,控制了区内南北2矿带的产出形态。太白向斜为一向西倾伏的不对称紧密褶皱,走向285°~295°,轴面倾向N,倾角80°,两翼倾角北陡南缓,北翼倾角65°~88°,南翼倾角33°~75°。长度大于12.0km,核部由黑云变粒岩组成,东部为盖层所覆盖。向斜枢纽自东向西平缓倾伏,倾伏角约3°~6°。向斜两翼赋存有沉积变质型铁矿带。草主山铁矿带即位于太白向斜的中东段。

区内盖层构造以断裂为主,按走向可分为近SN向、NE向2组,主要为白水牛石断裂(F6)及其派生的NE向次级断层,区内主要发育4条NE向断层,断层截穿了矿体,对矿体有不同程度的破坏。

(1)白水牛石断裂(F6):位于矿区东侧,走向近SN,纵贯全区,区域延伸长度达17.5km。断层倾向W,倾角60°~70°,为高角度正断层,水平断距200m左右,垂直断距达247m。断层破碎带宽15~40m,构造角砾岩发育,泥质、硅质胶结良好,局部发育有4m的断层泥。断层切穿了含矿地层及铁矿带,对矿带的连续性起明显的破坏作用。因断层陡直且距矿区较远,对矿体无影响。

(2)西石门断层(F2):断层长约3km,总体走向35°,北端与F3相交,断层倾向NW,倾角87°,断裂带宽0.5~1m,其内发育角砾岩。该断裂为西北盘上升的逆断层。位于矿区西侧,对矿体无影响。

(3)东石门断层(F3):位于矿区西端,南起枣庄断层,北至老书房延出矿区,断层走向南段为25°,北段为0°,长度>6.0km,具波状弯曲,在北段与F2交会。断层倾向NW或W,倾角45°~70°,为正断层。断层破碎带一般宽2~4m,最宽达14m左右,其内断层泥和角砾岩发育,角砾成分主要为黑云变粒岩、磁铁石英岩、页岩、灰岩等,次棱角状,少数呈压扁透镜状。该断层垂直断距为158m,切割了太白向斜和石闫背斜,使断层西侧矿体埋深增大,对矿区内矿体基本无影响。

(4)F4断层:位于矿区东南部,延伸长度约1500m,断层走向60°,倾向NW,倾角75°~82°,为高角度正断层。据钻孔揭露,断层水平断距20~50m,垂直断距15~30m,沿倾向于-300m处交会于F5断层。断层向北延伸较短,仅对南矿带东段略有错断。

(5)F5断层:位于矿区东南部,与F6呈40°斜交,为F6的次级断层。断层延伸长度约7km,走向45°~60°,倾向NW,倾角59°~63°,为高角度正断层。据钻孔揭露,断层水平断距2~60m,垂直断距40~60m。断层上部倾角较陡,为59°~63°,沿倾向向下趋于平缓尖灭,由中段图判断该断层尖灭于-600m标高之上。断层切割错断了区内矿床,对矿区东部-600m标高以浅矿体的连续性造成破坏。

2.3 岩浆岩

区内岩浆岩不发育,主要为新太古代石英闪长岩和元古代二长花岗岩。多呈脉状分布于山草峪组地层中,对区内铁矿体的分布影响较小。

3 矿体地质

3.1 矿体特征

草主山铁矿为隐伏矿床,发育南、北2条主矿带5个矿体,北(N)矿带位于矿区北部,发育N①,N②,N③3个矿体;N②,N③矿体自7线向西合为一个矿体;具有平行延伸的特点(图2)。南(S)矿带位于矿区南部,发育S①,S②两个矿体(表1)。两矿带总体走向270°~280°,顶端间距为440~550m,平行展布,相向而倾,表现为不对称的向斜构造特征(太白向斜中东段)。两矿带厚度自西向东具收缩变窄的趋势。矿带沿走向和倾向呈舒缓波状延伸,产出稳定,局部具膨胀狭缩、尖灭再现特征,沿走向向东有收敛变窄的趋势。矿体呈层状、似层状赋存于新太古代泰山岩群山草峪组变质地层中,产状与围岩片麻理产状基本一致[10],产出较稳定。矿体自上而下呈变陡趋势,深部局部具倒转现象。矿体厚度变化小,属厚度稳定型矿体,矿化连续,属品位均匀型矿体。矿体上覆盖层为青白口纪至寒武纪沉积地层,矿石成分为磁铁角闪石英岩,顶底板围岩为黑云变粒岩或黑云角闪片岩。矿体赋存于+70m~-1026m标高之间,顶部埋深0~170m。矿床平均品位TFe32.36%,mFe19.67%,矿体总厚度34.83m。

1—第四系;2—二青山组石英砂岩段;3—泰山岩群山草峪组;4—磁铁矿体及编号;5—角度不整合地质界线;6—地质界线;7—钻孔位置及编号;8—基点位置及编号图2 草主山7号勘探线剖面图

矿带矿体编号形态规模延长(m)延深(m)均真厚度(m)厚度变化系数(%)平均品位(%)变化系数(%)TFemFeTFemFeNN①层状17759536.4448.6932.1119.8215.7035.32N②层状17769618.3447.6332.8319.0114.8036.36N③层状14639965.5947.2033.9122.7114.3730.95SS①似层状10876927.0963.0729.8117.6618.5139.50S②层状159110407.3249.2332.1919.5718.4435.42

3.1.1 N①矿体

N①矿体位于北矿带上部,其铁矿石资源量占矿床总资源量的15.73%,占草主山矿区北矿带资源量的26.78%。矿区内沿走向工程控制长度1775m,控制最大斜深953m。矿体一般厚度1.98~9.55m,最大厚度25.13m,最小厚度1.05m,平均厚度6.44m。厚度变化系数48.69%,属厚度稳定型矿体。矿体沿倾向厚度较稳定,略具中间薄、上下厚的特征。矿体单样品最高品位TFe40.38%,mFe31.78%,单工程平均品位TFe29.50%~37.41%,mFe17.17%~26.06%,矿体平均品位TFe32.11%,mFe19.82%。品位变化系数TFe15.70%,mFe35.32%,矿体矿化连续,属品位均匀型矿体。

3.1.2 N②矿体

N②矿体位于北矿带中部,N①矿体的北侧,距N①矿体10~32m。其铁矿石资源量占矿床总资源量的28.54%,占草主山矿区北矿带资源量的48.58%,为矿区的主矿体之一。矿区内沿走向工程控制长度1776m,控制最大斜深961m。矿体一般厚度1.56~16.72m,最大厚度24.37m,最小厚度1.50m,平均厚度8.34m。厚度变化系数47.63%,属厚度稳定型矿体。矿体单样品最高品位TFe40.51%,mFe32.02%,单工程平均品位TFe27.84%~35.71%,mFe12.72%~25.41%,矿体平均品位TFe32.83%,mFe19.01%。品位变化系数TFe14.80%,mFe36.36%,矿体矿化连续,属品位均匀型矿体。

3.1.3 N③矿体

N③矿体位于北矿带下部,N②矿体的北侧,距N②矿体0~9m,于3线以西与N②矿体合并为同一矿体。其铁矿石资源量占矿床总资源量的14.47%,占草主山矿区北矿带资源量的24.64%。矿区内沿走向工程控制长度1463m,控制最大斜深996m。矿体厚度一般为2.11~12.44m,最大厚度13.25m,最小厚度2.03m,平均厚度5.38m。厚度变化系数47.20%,属厚度稳定型矿体。矿体单样品最高品位TFe40.06%,mFe32.89%,单工程平均品位TFe30.85%~37.57%,mFe17.32%~27.07%,矿体平均品位TFe34.05%,mFe22.54%。品位变化系数TFe14.37%,mFe30.95%,矿体矿化连续,属品位均匀型矿体。

3.1.4 S①矿体

S①矿体位于南矿带上部,其铁矿石资源量占矿床总资源量的11.50%,占草主山矿区南矿带资源量的27.87%。矿区内沿走向工程控制长度1087m,控制最大斜深692m。矿体厚度一般为1.99~11.15m,最大厚度18.79m,最小厚度1.55m,平均厚度7.09m。厚度变化系数63.07%,属厚度较稳定型矿体。矿体单样品最高品位TFe38.39%,mFe30.87%,单工程平均品位TFe20.69%~35.07%,mFe12.29%~23.83%,矿体平均品位TFe29.81%,mFe17.66%。品位变化系数TFe18.51%,mFe39.50%,矿体矿化连续,属品位均匀型矿体。

3.1.5 S②矿体

S②矿体位于南矿带下部,S①矿体南侧,距S①矿体6~35m。其铁矿石资源量占矿床总资源量的29.76%,占草主山矿区南矿带资源量的72.13%,为矿区的主矿体之一。矿区内沿走向工程控制长度1591m,控制最大斜深1040m。矿体厚度一般为1.70~12.81m,最大厚度19.45m,最小厚度1.29m,平均厚度7.32m。厚度变化系数49.23%,属厚度稳定型矿体。矿体单样品最高品位TFe41.86%,mFe32.36%,单工程平均品位TFe18.03%~35.76%,mFe14.49%~25.93%,矿体平均品位TFe32.19%,mFe19.57%。品位变化系数TFe18.44%,mFe35.42%,矿体矿化连续,属品位均匀型矿体。

3.2 矿石质量

3.2.1 矿石矿物成分

根据岩矿鉴定和野外观察,矿石中金属矿物主要有磁铁矿、磁黄铁矿、假象赤铁矿、黄铁矿及褐铁矿、黄铜矿等;非金属矿物主要有角闪石、石英、黑云母、铁闪石、透闪石、方解石,少量绿帘石、绿泥石、石榴子石、磷灰石、金红石等。

表2 矿石矿物成分

3.2.2 矿石化学成分

根据基本分析和组合分析测试结果,矿石中主要有用组分为Fe,矿床平均品位TFe32.36%,mFe19.67%。其它化学成分主要为SiO2,其次为Al2O3,CaO,MgO,As,S,P等含量较少。SiO2含量30.14%~56.91%,平均为46.16%;CaO含量1.66%~7.29%,平均为3.54%;MgO含量1.58%~7.88%,平均2.85%;Al2O3含量0.44%~15.79%,平均4.75%;S含量0.002%~3.887%,平均为0.33%;P含量0.046%~0.31%,平均0.13%;As含量(0.07~34.10)×10-6,平均为4.44×10-6。矿石中各类化学组分含量相对较为稳定。

表3 各化学成分含量

备注:As×10-6;其他×10-2

3.2.3 矿石结构构造

矿石结构为条纹状结构、片状粒状变晶结构和柱状变晶结构。矿石构造较为简单,以条带状构造为主,块状构造次之。

3.2.4 矿石类型

矿石的自然类型按主要铁矿物划分属磁铁矿石;按主要脉石矿物划分属石英型和闪石型;按矿石结构构造划分属条纹条带状铁矿石。即石英闪石型条纹条带状磁铁矿石。该矿区矿石平均品位TFe32.36%,mFe19.67%,为低品位的贫铁矿石,其工业类型属需选(贫)铁矿石[11]。

3.3 矿体围岩与夹石

3.3.1 矿体围岩

矿体赋存于泰山岩群山草峪组变质地层中,矿体顶底板围岩主要为黑云变粒岩和黑云角闪片岩,少量含磁铁黑云角闪片岩及磁铁角闪石英岩、磁铁石英角闪岩等。黑云变粒岩、黑云角闪片岩、含磁铁黑云角闪片岩的颜色、结构、构造特征与矿石差异性明显,其与矿体界限清晰;磁铁角闪石英岩、磁铁石英角闪岩的颜色、结构、构造特征与矿石基本一致,仅磁性铁含量较矿体低,其与矿体无明显界限,需根据化验结果按指标划分。南北矿带底板近矿围岩为黑云变粒岩,其中南矿带底板近矿围岩的黑云变粒岩相对其他部位同类岩性颜色较深、粒度较细。

3.3.2 夹石

该矿区矿体矿化较连续,矿体内部结构简单,夹石较少,多呈似层状、透镜状产出,工程控制夹石厚度为1.00~4.49m,岩性以含磁铁黑云角闪片岩为主、其次为磁铁角闪石英岩,少量黑云角闪片岩或黑云变粒岩。

含磁铁黑云角闪片岩品位一般为TFe14.87%~31.80%,mFe0.90%~10.93%。磁铁角闪石英岩的品位一般为TFe22.93%~30.68%,mFe6.34%~21.07%。其岩性特征、结构构造、矿物组合等与矿石基本一致,只是磁性铁含量比矿石低。

4 矿床成因及找矿标志

4.1 矿床成因

新太古代,鲁西地区沉积了巨厚的泥砂质及基性火山岩建造,局部火山喷发的硅铁质,经海水搬运分异沉淀,铁质得到了初步富集。经中温中压环境下的区域变质作用,原岩物质发生重结晶,铁质进一步富集,形成了以磁铁矿、角闪石、石英、黑云母等为代表矿物的变质岩石,即黑云变粒岩、黑云角闪片岩、磁铁角闪石英岩等,构成赋矿层位。磁铁矿体中的矿物成分与围岩大致同时沉积形成,因变质过程中未发生大规模的物质迁移及交代,原始沉积层理得以保存,而形成明显的条带状构造。

综上所述,海底基性火山岩喷发为成矿提供了丰富的成矿物质,经海水搬运分异沉淀,铁质得到了初步富集,在区域变质中温中压条件下,原岩物质发生重结晶作用,铁质进一步富集,形成以磁铁矿、角闪石、石英等为主要矿物的磁铁矿床。故矿床应属沉积变质成因铁矿[12-14]。

4.2 找矿标志

4.2.1 地层标志

泰山岩群山草峪组为区内铁矿的赋存层位,岩性以黑云变粒岩为主,局部夹磁铁角闪石英岩、磁铁石英角闪岩、黑云角闪片岩等[15]。其中磁铁角闪石英岩、磁铁石英角闪岩,为沉积变质型(鞍山式)铁矿的含矿岩石[16],黑云变粒岩与黑云角闪片岩为矿体顶底板围岩。因此,泰山岩群山草峪组为寻找“鞍山式”铁矿的有利层位;山草峪组地层中各类片岩发育层段为成矿有利部位。

4.2.2 构造控矿标志

铁矿及其赋矿层位褶皱构造比较发育,主要为紧闭顶厚褶皱,往往在向斜转折端附近矿层加厚,形成厚大矿体。由于褶皱构造发育,造成含矿带重复出现[17]。草主山铁矿带位于太白向斜的中东段。太白向斜为一向W倾伏的不对称紧密褶皱[18],走向285°~295°,轴面倾向N,倾角80°,两翼倾角北陡南缓,北翼倾角65°~88°,南翼倾角33°~75°。长度大于12.0km,核部由黑云变粒岩组成,东部为盖层所覆盖。向斜枢纽自东向西平缓倾伏,倾伏角约3°~6°。向斜两翼赋存有沉积变质型铁矿带[19-20]。

4.2.3 物探找矿标志

以往勘查期间开展了系统的高精度磁测工作。成果表明,在磁铁矿体分布地段可引起明显的高磁异常,异常规律性明显,与矿体吻合较好。ΔT强度一般为500~800nT,峰值可达1000nT以上,并伴生有负磁异常。高精度磁测对寻找隐伏铁矿体,指导探矿工程的布设,具有良好的找矿效果。因此,在该地区分布的局部高磁异常可作为寻找同类型铁矿的间接找矿标志。

5 深部及外围找矿方向

在勘查工作中,草主山铁矿钻孔工程控制标高为-800m,经工程揭露,南北两个矿带五个矿体向下均有沿走向倾向未封闭的特点(图3),具有较好找矿前景。并且沿已有磁铁矿体走向追索是寻找“鞍山式”铁矿的重要标志[21]。

1—钻孔;2—断层及编号;3—铁矿(化)体;4—不整合界线;5—泰山岩群山草峪组图3 草主山铁矿区主要勘探线联合剖面图

6 结语

草主山铁矿成矿地质条件良好,新太古代泰山岩群山草峪组的变质作用为铁矿的形成提供了丰富的物质来源,基底变质岩层中发育的褶皱构造和断裂构造控制区内“鞍山式”铁矿层的分布与形态。区内矿体厚度稳定、品位均匀,沿走向倾向均未封闭,提供了向深部找矿的可能性,并且矿体多呈带状沿走向展布,尖灭再现分支复合现象明显,故矿体深部及外围尚有较大的找矿远景。

[1] 李厚民,王登红,李立兴,等.中国铁矿成矿规律及重点矿集区资源潜力分析[J].中国地质,2012,39(3):559-580.

[2] 徐友松,罗海伟,赵金,等.山东省东平县大牛地区铁矿成因类型及找矿方向研究[J].山东国土资源,2015,31(10):17-22.

[3] 王继芳,徐然,杜显彪,等.山东省东平县大高庄铁矿床地质特征及矿床成因[J].山东国土资源,2015,31(5):31-36.

[4] 徐东来,胡艳霞.平阴县洪范池铁矿地质特征[J].山东国土资源,2008,24(1):8-9.

[5] 李洪奎,杨永波,耿科,等.山东大地构造相研究方法划分方案与基本特征[J].山东国土资源,2014,30(4):1-12.

[6] 孙茂田,安仰生,张忠涛.苍峄铁矿带地质特征及找矿远景分析[J].山东国土资源,2017,33(4):6-14.

[7] 孔庆友,张天祯,于学峰,等.山东矿床[M].济南:山东科学技术出版社,2006:310-313.

[8] 李厚民.沉积变质型铁矿的分类[J].矿床地质,2012,31(增刊):117-118.

[9] 王世进,万渝生,宋志勇,等.鲁西泰山岩群地层划分及形成时代----锆石SHRIMP U-Ph测年的证据[J].山东国土资源,2015,31(10):17-22.

[10] 曾广湘,吕昶,徐金芳,等.山东铁矿地质[M].济南:山东科学技术出版社,1998:30-35.

[11] 邵厥年,陶维屏.矿产资源工业要求手册[M].北京:地质出版社,2010:84.

[12] 周军,葛均,段红梅.鲁西台上铁矿床地质特征及矿床成因[J].山东国土资源,2010,26(10):6-10.

[13] 安仰生,韩廷宝,刘邦君,等.山东苍山王埝沟铁矿床地质特征及找矿方法探讨[J].山东国土资源,2008,24(7-8):44-47.

[14] 郝兴中,杨毅恒,李英平,等.山东苍峄铁矿带预测模型[J].吉林大学学报,2013,43(4):1136-1142.

[15] 宋明春,王沛成.山东省区域地质[M].济南:山东省地图出版社,2003:30-31.

[16] 郭宝奎,熊玉新,孙斌,等.山东省苍山县凤凰山铁矿铁品位统计分析[J].山东国土资源,2012,28(11):13-16.

[17] 于学峰,张天祯,王虹,等.山东矿床成矿系列[M].北京:地质出版社,2015:88-89.

[18] 张增奇,刘明渭.山东省岩石地层[M].济南:山东科学技术出版社,1998:30-35.

[19] 李银真,黄文院,刘建稳,等.山东省苍山县鲁城地区铁矿成矿规律及深部找矿研究[A]//中国地质学会2013年学术年会论文摘要汇编[C].2013.

[20] 刘玉强,李洪喜,黄太岭,等.山东省金、铁、煤矿床成矿系列及成矿预测[M].北京:地质出版社,2004:99-118.

[21] 郝兴中,杨毅恒,李英平,等.综合找矿方法在覆盖区的应用:以山东省单县大刘庄铁矿勘查为例[J].吉林大学学报:地球科学版,2013,43(2):641-648.

AnalysisonGeologicalCharacteristicsandProspectingPotentialofIronDepositsinCaozhushanMiningAreainLanlingCountyofShandongProvince

FENG Xiping, BIAN Conglin, LIU Chao, GAO Yan

(Lunan Geo-engineering Exploration Institute, Shandong Yanzhou 272100, China)

Caozhushan mining area is a part of Cangyi iron ore belt. Hidden Anshan type iron deposit occurs in this area. Through explorations, two main ore belts in north and south and five ore bodies have been circled. Ore bodies occur in Shancaoyu formation of Taishan group. Its forms are controlled by folds in the base with the characteristics of extending parallelly. After previous geological exploration, the position, shape, scale, occurrence, distribution, thickness, ore quality of ore bodies in the depth of -800m have been identified. Ore bodies have not closed along the trend and tendency. It has the trend of extending continuously, and ore grades and thickness do not have the trend of cutting. It is concluded that there is a big prospecting prospect in the depth of -800m. Through summarizing metallogenic regularity and prospecting criteria, prospecting model has been established. It will provide reliable geological references for further exploaration in this area.

Lanling county; Caozhushan mining area; Anshan type iron deposit; ore prospecting; metallogenic regularity

2017-04-06;

2017-06-14;

王敏

冯玺平(1995—),女,山东菏泽人,助理工程师,主要从事金矿、铁矿地质矿产勘查工作;E-mail:1561718028@qq.com

P618.31

A

冯玺平,边丛林,刘超,等.山东省兰陵县草主山矿区铁矿床地质特征和找矿远景分析[J].山东国土资源,2017,33(12):17-23.FENG Xiping,BIAN Conglin,LIU Chao,etc.Analysis on Geological Characteristics and Prospecting Potential of Iron Deposits in Caozhushan Mining Area in Lanling County of Shandong Province[J].Shandong Land and Resources,2017,33(12):17-23.

猜你喜欢
黑云磁铁铁矿
漂流计划
漂流计划
大红山铁矿找矿前景分析
黑云白云
漫画与幽默
冀东南部铁矿区重力异常特征
男妖女妖知多少
无轨斜坡道在大红山铁矿中的应用
玩磁铁
用磁铁治病,由不信到信等