地球物理测井(简称测井)作为重要组成部分和关键技术之一，与钻探取心、岩心测试和地面地球物理技术一起构成了科学钻探的完整体系(牛一雄等，2004).近年来，我国相继开展了江苏东海中国大陆科学钻探(Pan et al.，2002)，青海湖、柴达木盆地环境科学钻探(李斌凯等，2007)，松辽盆地白垩纪科学钻探(高有峰等，2011)，汶川地震断裂带科学钻探(唐力君等，2013；李海兵等，2013)等工程，其中地球物理测井为地学研究提供了大量的原位测量数据和解释成果，从而为推断地球内部的物质组成、动态演化及驱动力，了解地球内部结构、构造和动力机理，了解与地质环境、灾害有关的信息提供了重要依据.

随着找矿空间向深部拓展，把重要成矿带和矿集区深部结构探测和成矿学研究密切结合，探索大型矿集区和巨型矿床形成的深部控制因素，开辟了成矿学研究和深部找矿的新思路(吕庆田等，2011).通过科学钻探能探测地壳精细结构，探测深部物性结构，探测深部物质成分并实验模拟，实现矿集区的“透明化”(董树文等，2014).庐枞盆地构造特征复杂、岩浆-成矿作用活跃，区内广泛发育与火山-潜火山岩及侵入岩有关的铁、铜、铅锌、明矾石矿床，是长江中下游成矿带中7个重要的矿集区之一(常印佛等，1991；翟裕等，1992；唐永成等，1998).在该地区选择重要地球物理异常或深部成矿有利区，以具有代表性矿床的成矿环境为主要探测和研究对象，部署2000 m岩心钻探孔，并开展地球物理测井工作，达到获得矿集区深部地层、物性信息，验证地球物理异常(吕庆田等，2011)的目的.

庐枞矿集区砖桥科学钻探ZK01孔地球物理测井为“深部探测技术与实验研究专项(SinoProbe)”中的研究内容，依托于项目3(SinoProbe-03)“深部矿产资源立体探测及试验研究”的“重要异常的钻探验证与金属垂向分布规律研究”课题.地球物理测井可以划分和校验地质剖面、查明目的层，研究与金属矿有成因关系的围岩蚀变以及发现金属矿，同时又是一种原位测定各种岩矿石物理性质的可靠手段.本文通过庐枞矿集区砖桥科学钻探ZK01孔的地球物理测井与井中物探研究工作，建立钻孔附近的岩石物性剖面，研究庐枞矿集区深部物性分布状况，为研究深部成矿理论提供岩矿石的地球物理信息，以助于提高对我国东部成矿背景、成矿动力学过程的认识.同时在钻孔深部发现了强放射性异常，对庐枞深部找放射性矿具有理论和实际意义.

庐枞盆地位于扬子板块北缘，郯庐断裂带东侧，是长江中下游断陷带内一东北向的不对称耳状盆地.盆地中的岩石地质单元大体可划分为两部分：一 是前寒武纪的变质结晶基底，二是以寒武系以来的沉积盖层(常印佛等，1991；任启江等，1991；胡文瑄等，1991).ZK01孔位于庐枞盆地中部，孔深2012.35 m. 钻孔2～1469.81 m主要为砖桥组火山岩，主要岩石为粗安岩、高岭石化粗安岩、黄铁矿化粗安岩、硬石 膏化粗安岩和辉石粗安岩；1488.84～1603.35 m为火山岩地层与下伏岩体的接触带；1603.35～2012.35 m 为岩体，自上而下依次为(石英)正长岩、石英二长岩、黑云母石英二长岩，相互之间为渐变过渡关系(张舒等，2014).

ZK01孔的测井数据采集工作由中国地质科学院物化探所负责，严格按照相关规范分三次进行，累计完成测井工作量1994.02 m.测井数据采集设备为英国RG公司的Micrologger II便携式测井系统、北京中地英捷物探仪器研究所的PSJ-2数字测井系统、中国地质科学院物化探所和重庆地质仪器厂研制的JGCX-1高精度井中三分量磁力仪.

测井及井中物探测量项目包括视电阻率、极化率、磁化率、纵波速度、超声成像、自然伽马、密度、井斜、井径、井温、泥浆电阻率、井中三分量磁测等10多种方法.测井数据质量满足设计和规范要求.测井资料处理采用Advanced Logic Technology公司的WellCAD4.4测井解释软件，主要包括数据预处理、图像生成和解释.

测井作为一种可在接近原始条件下获得岩层原位物理特性的物探方法，不会像实验室岩心测试，受到岩心从井下到地面的采集、搬运过程中所发生的物理、化学性质变化影响.因此，利用测井方法来获得岩石物性、识别岩性非常必要.ZK01孔钻遇岩性主要为火成岩，钻遇率较高有代表性的岩性有粗安岩、凝灰岩、粉砂岩、安山玢岩、角砾岩、正长岩和二长岩.通过常规测井曲线与钻孔岩心编录资料的对比分析，选择视电阻率、极化率、磁化率、纵波速度、自然伽马、密度测井数据进行岩性测井响应和岩性识别研究.通过直方图统计，得到各岩性测井响应值.统计结果见表 1.

各岩性测井响应的常见值能够很好地代表某岩性的特征，最大值和最小值则可以指示岩性测井响应的可能分布范围.从岩性测井响应值统计表中可以认识到：

(1)正长岩、二长岩因其正长石含量较高，其自然伽马的值也比其他几种岩性的大，常见值超过300 API.

(2)视电阻率区分岩性效果较好，并按照二长岩、正长岩、次生石英岩、黄铁矿化粗安岩、安山玢岩、粗安岩、凝灰岩的顺序依次降低.

(3)纵波速度、极化率、磁化率、密度对各岩性的测井响应区分度不高，难以从单一测井曲线上有效辨别岩性.

不同火成岩石的测井特征相互重叠，为在交会图上能直观地区分各种岩石的分界和所分布的区域，需要综合应用各种测井曲线信息.以1490 m为分界线将砖桥组和岩体分别进行分析研究，砖桥组选取粗安岩、高岭石化粗安岩、黄铁矿化粗安岩、硬石膏化粗安岩和辉石粗安岩绘制测井属性交会图(图 1)，岩体则选取了正长岩、石英正长岩、石英二长斑岩和黑云母石英二长岩(图 2).本文在岩石物性分析中应用的交会图有：自然伽马-视电阻率交会图、极化率-视电阻率交会图、密度-纵波速度交会图和极化率-磁化率交会图.

图 1

Fig. 1

图 1 ZK01孔砖桥组地层各测井参数交会图 (a)自然伽马-视电阻率交会图；(b)密度-纵波速度交会图;(c)极化率-视电阻率交会图；(d)极化率-磁化率交会图. Fig. 1 Logging data Crossplots of Zhuanqiao formation (a)GR-RT Crossplots；(b)η-RT Crossplots；(c)DEN-Vp Crossplots；(d)η-M Crossplots.

图 2

Fig. 2

图 2 岩体各测井参数交会图 (a)自然伽马-视电阻率交会图；(b)密度-纵波速度交会图；(c)磁化率-视电阻率交会图. Fig. 2 Logging data Crossplots of Bedrock (a)GR-RT Crossplots；(b)DEN-Vp Crossplots；(c)M-RT Crossplots.

从图 1可以看出，自然伽马能有效的将高岭石化粗安岩与其他岩石区分开；极化率和视电阻率配合能将大部分硬石膏化粗安岩与其他岩性区分出来；密度和纵波速度配合，能识别部分黄铁矿化粗安岩；极化率和磁化率配合能识别部分辉石粗安岩.

从图 2可以看出，自然伽马能将正长岩与其他岩石区分开；视电阻率能将石英二长斑岩与其他岩性区分开；纵波速度和密度配合，能将石英正长岩、石英二长斑岩和黑云母石英二长岩这三种岩石区分开；磁化率能将正长岩和石英正长岩与石英二长斑岩和黑云母石英二长岩区分开.

表 1 (Table 1)

表 1 ZK01孔主要岩性测井响应值统计表 Table 1 Log response values of different lithologies in borehole ZK01 岩性 自然伽马 (API) 视电阻率 (Ωm) 极化率 (%) 纵波速度 (m·s-1) 磁化率 (×10-4SI) 密度 (g·cm-3) 粗安岩 极大值 476 3785 44.21 6425 2158 3.12 极小值 71 4 0.19 1600 837 1.27 常见值 180 170 1.14 5800 900 2.79 黄铁矿化 粗安岩 极大值 472 5899 10.94 6764 1669 3.11 极小值 51 65 0.09 4458 853 2.52 常见值 190 2050 0.73 5800 1030 2.82 次生石英岩 极大值 857 4692 2.71 6668 1165 3.08 极小值 94 558 0.77 4608 943 2.65 常见值 230 4500 1.23 5750 950 2.76 安山玢岩 极大值 233 363 1.19 6015 1050 2.87 极小值 173 224 0.67 5471 949 2.78 常见值 205 350 0.79 5780 950 2.81 凝灰岩 极大值 329 760 6.21 6254 1240 3.11 极小值 123 39 0.35 2834 947 2.29 常见值 225 80 1.74 5400 960 2.82 正长岩 极大值 2789 11187 9.10 6934 1249 3.14 极小值 69 509 0.58 4739 937 1.63 常见值 731 5000 1.47 5800 945 2.66 二长岩 极大值 4032 12265 — 6536 2075 3.07 极小值 69 1183 — 5435 941 2.32 常见值 405 6200 — 5880 1219 2.75

表 1 ZK01孔主要岩性测井响应值统计表 Table 1 Log response values of different lithologies in borehole ZK01

人工识别岩性准确性高，但耗时相对较长，将支持向量机方法运用于岩性识别，可提高工作效率.支持向量机是建立在统计学习理论VC维理论和结构风险最小化原理基础上的机器学习方法，在解决小样本、非线性和高维模式识别问题中表现出许多特有的优势，并在很大程度上克服了“维数灾难”和“过学习”等问题(丁世飞等，2011).

本文优选了220个具有代表性的岩心作为训练样本，其中粗安岩40个，黄铁矿化粗安岩40个，高岭石化粗安岩35个，正长岩35个，安山玢岩30个，凝灰岩40个，提取它们的声波时差、自然伽马、视电阻率、密度测井数据建立了4维6类的样本空间.为了消除各特征量纲不同所带来的影响，对学习样本的测井参数进行归一化处理，将每一组测井值归一到[-1，1](邓呈祥等，2015).

本次岩性识别中，支持向量机的优化核函数参数γ为151.9852，惩罚因子C为9.1105.如图 3所示，为250~440 m井段人工识别综合分层和支持向量机判别分层的结果对比，总体上判别结果接近(邓呈祥等，2015).

图 3

Fig. 3

图 3 支持向量机的判别结果 Fig. 3 The profile of discrimination results of SVM

庐枞矿集区砖桥科学钻探ZK01孔岩性复杂，通过交会图等技术对测井响应特征分析，提取了不同岩性不同测井参数的特征.在此基础上对全孔进行了测井岩性识别，建立了钻孔钻遇主要岩性的测井解释岩性剖面，如图 4所示.测井推断本孔的地层分布124层，包含40种岩性.其中，粗安岩有15层，共计526.6 m，主要分布在1000 m以上的部位；黄铁矿化粗安岩有10层，共计319.2 m，分布在ZK01孔610~1440 m的部位；高岭石化粗安岩有16层，共计301.4 m，主要分布在1200 m以上的部位；凝灰岩有6层，共计101.6 m，主要分布在ZK01孔的中间部位；安山玢岩有4层，共计28.5 m，主要分布 在1600 m以下的部位；正长岩有5层，共计248.5 m，主要分布在1480 m以下的部位；二长岩全分布在1790 m以下.

图 4

Fig. 4

图 4 ZK01孔的测井识别岩性剖面 Fig. 4 The lithologic profile discriminated with borehole logging data in ZK01

从测井识别岩性的情况来看，ZK01孔附近在1480 m以上主要为喷出岩，1480 m以下为深层侵入岩，这显示了该地区同时经历了火山喷发活动与深部岩浆活动(周涛发等，2010；贾丽琼等，2014)；岩石与矿化在喷出岩的矿化区域中立体分布，则说明该地区在火山喷发结束之后，深部岩浆活动也逐渐停止，伴随的热液成矿作用(熊欣等，2014a)导致了粗安岩的硬石膏化、黄铁矿化和高岭石化(张荣华等，2010).

ZK01孔井深20~1480 m处的黄铁矿化现象比较常见，均发生在粗安岩当中.粗安岩的黄铁矿化和高岭石化、硬石膏化均有明显的测井响应，如图 5所示.硬石膏化粗安岩和高岭石化粗安岩的磁化率小于1000×10^-4SI，远小于粗安岩和黄铁矿化粗安 岩.硬石膏化粗安岩的电阻率主要分布在100~700 Ωm 之间，高岭石化粗安岩的电阻率主要分布在60~200 Ωm之间，硬石膏化粗安岩的电阻率减小，而高岭石化电阻率变化不大；根据黄铁矿化磁化率特征，790~840 m处的黄铁矿化粗安岩较900~960 m处的黄铁矿化粗安岩的矿化程度高.

图 5

Fig. 5

图 5 ZK01孔700~1000 m测井曲线图 Fig. 5 The curves of borehole logging in ZK01(700~1000 m)

黄铁矿化粗安岩集中在井深610~1440 m的部位，测井曲线中声波时差减小、密度增大、视电阻率增大；高岭石化粗安岩分布在1200 m以上部位，测井曲线中声波时差增大、密度减小、视电阻率减小、磁化率减小.

超声成像测井中超声波脉冲的井壁反射信号振幅反映井壁的声阻抗分布：地层与井液声阻抗差小，反射波信号振幅小，在成像测井振幅图像上颜色较深(邹长春等，2014).如图 6所示，当井穿过裂缝时，裂缝表现为一条深色的正弦曲线：如图 6a是250 m处的多个平行低角度缝；图 6b是ZK01孔中常见的垂直缝(倾角大于75°).

图 6

Fig. 6

图 6 ZK01孔超声波成像测井 (a)平行低角度缝；(b)垂直缝；(c)井壁扩孔. Fig. 6 The intensity image of BHTV logging in ZK01 (a)Parallel low-angle fractures;(b)Vertical fracture;(c)Borehole breakouts.

ZK01孔裂缝较为发育，共拾取了裂缝3200余条，平均值为1.6条/m，其中300～400 m和800～ 900 m深度范围内裂缝尤为发育(图 7a)，而在1100 m 以下相对不发育.

图 7

Fig. 7

图 7 ZK01孔裂缝参数随深度的变化规律 (a)密度;(b)倾角;(c)倾向. Fig. 7 Statistics of density(a)，dip(b) and inclination(c)of fractures in ZK01

全井段裂缝倾角大小整体上呈现出随深度增加而增大的趋势，且主要分布在50°～80°之间(图 7b)；裂缝倾向分布不集中，在0°～360°之间均有一定量的分布(图 7c)，裂缝倾向18~550 m主要分布于0°～140°、240°～360°，550～1000 m主要分布在150°～300°，1000～1400 m分布较为均匀，1400～1964 m主要分布于0°～150°与270°左右.综上，裂缝倾角以高角度为主，倾向分布不集中，随深度变化明显.

钻孔中，泥浆压力与地应力的不同会导致井壁上的应力释放(如井壁崩落)，从而产生双对称分布的井壁崩落或扩孔.通过井壁崩落方向可以判断地应力方向，根据井眼崩落形成的椭圆井眼特征，椭圆的长轴方向为井眼崩落方向，短轴方向为最大水平 主应力(Nie et al，2013)，如图 6c所示，1038～1042 m 的井段上，崩落方向大致为90°和270°，从图 6c的崩落方向可以推断出，最大主应力方向大致沿南北走向.

放射性测井是确定井中铀矿含量、厚度及深度的主要手段.通过测井资料获得地层中铀含量分布并找出非矿夹层，为铀储量计算提供基础数据并指导铀矿开采(骆淼等，2008).

ZK01孔在1500～1900 m井段发现高自然伽马异常，且异常层厚度较大.在该铀矿化段上下围岩 的自然伽马大多处于180API左右水平，而在1500～1900 m 多见高达2000API左右的尖峰状自然伽马异常.1500～1900 m井段自然伽马测井曲线及其测井岩性识别见图 8.

图 8

Fig. 8

图 8 ZK01孔1500～1900 m测井识别岩性图 Fig. 8 The lithologic profile discriminated with borehole logging curves of ZK01(1500～1900 m)

数据采集采用的自然伽马仪器标定铀含量灵敏度为877.88API相当于万分之一铀含量(0.01%eus)^-1，针对放射性异常，对铀当量进行了计算.本文仅对铀当量大于万分之一的异常进行了解释，通过对ZK01孔壁铀当量曲线统计，共发现21处异常，见表 2.其中，铀当量最高井段在1897.39～1897.93 m，厚度0.54 m，铀当量0.039%；单层最 厚异常井段在1603.96～1643.95 m，厚度39.99 m，铀当量0.011%，见图 9.解释异常累计厚度93.02 m.

表 2 (Table 2)

表 2 地层铀当量计算结果 Table 2 The uranium equivalent value calculated in formations 异常 序号 起始深度 (m) 终止深度 (m) 厚度 (m) 铀当量 (%) 岩性 异常 序号 起始深度 (m) 终止深度 (m) 厚度 (m) 铀当量 (%) 岩性 1 1499.76 1506.06 6.30 0.018 正长岩 12 1709.51 1710.19 0.68 0.011 石英正长岩 2 1506.06 1508.37 2.31 0.013 正长岩 13 1713.83 1714.60 0.77 0.018 石英正长岩 3 1514.12 1520.67 6.55 0.018 正长岩 14 1728.03 1728.86 0.83 0.016 石英正长岩 4 1523.87 1525.32 1.45 0.015 正长岩 15 1731.22 1733.60 2.38 0.013 石英正长岩 5 1527.42 1532.99 5.57 0.015 正长岩 16 1786.25 1787.63 1.38 0.011 二长岩 6 1571.16 1573.37 2.21 0.015 正长岩 17 1795.67 1796.48 0.81 0.020 二长岩 7 1575.19 1579.17 3.98 0.015 石英正长岩 18 1840.79 1842.21 1.42 0.012 二长岩 8 1603.96 1643.95 39.99 0.011 粗安岩(捕虏体) 19 1870.71 1871.64 0.93 0.015 二长岩 9 1659.89 1672.32 12.43 0.010 石英正长岩 20 1881.99 1882.88 0.89 0.020 二长岩 10 1701.74 1702.32 0.58 0.020 正长岩 21 1897.39 1897.93 0.54 0.039 二长岩 11 1705.56 1706.58 1.02 0.011 石英正长岩 合计 1499.76 1897.93 93.02

表 2 地层铀当量计算结果 Table 2 The uranium equivalent value calculated in formations

图 9

Fig. 9

图 9 1603.96～1643.95 m铀当量解释曲线 Fig. 9 The curve of uranium equivalent value(1603.96~1643.95 m)

为研究放射性异常的原因，选取异常井段14件样品进行质谱和光谱分析，得到岩石样品53种元素含量，表 3为选择的相关主要元素含量统计表.在放射性异常地段采集的岩石样品化学分析结果显示，与正长岩类岩石平均化学组成相比，大多数元素并没有明显的异常显示，但是与放射性有关的几个元素U、Th和K₂O却出现较明显且连续的异常，U、Th两元素异常吻合较好，多数样品中两元素含量的富集程度及其变化趋势一致.此外出现异常的还有Cu元素.

表 3 (Table 3)

表 3 ZK01孔岩心样本中主要元素含量统计表 Table 3 The main element contents of core samples in ZK01 样本编号 深度 /m Cu /(μg·g-1) Zn /(μg·g-1) Ag /(μg·g-1) Th /(μg·g-1) K2O /(%) U /(μg·g-1) Al2O3 /(%) CaO /(%) Fe2O3 /(%) H-120 1503 5.5 4.6 34.7 191.6 0.11 34.3 13.36 1.90 0.70 H-121 1506 11.2 5.8 29.6 146.0 0.11 19.7 13.77 1.18 1.77 H-122 1518 5.5 6.8 41.3 196.1 0.22 50.9 14.25 1.27 1.27 H-123 1515 4.2 7.3 35.9 169.4 0.27 47.4 15.19 0.91 0.68 H-124 1528 131.5 30.6 65.0 40.6 3.47 12.9 15.03 4.93 4.63 H-125 1532 317.8 34.0 85.5 35.3 5.11 17.2 15.89 4.25 4.72 H-127 1622 90.9 31.4 97.0 217.0 5.40 95.6 15.34 1.31 2.18 H-126 1624 332.0 29.3 105.3 136.4 5.50 27.6 15.19 2.95 3.86 H-129 1640 92.2 32.8 76.0 109.3 6.17 34.7 14.85 1.41 3.71 H-128 1644 7.1 83.2 104.0 13.6td> 8.41 2.1 16.22 4.57 8.55 H-132 1662 236.2 37.0 60.1 153.3 3.90 12.3 12.89 0.72 5.41 H-131 1666 929.7 51.3 122.0 280.0 2.79 19.7 14.78 0.75 13.15 H-130 1668 148.6 48.3 66.7 57.5 8.37 10.3 16.29 1.43 3.39 H-133 1896 50.7 71.8 96.7 179.8 5.18 70.4 12.87 7.84 2.60

表 3 ZK01孔岩心样本中主要元素含量统计表 Table 3 The main element contents of core samples in ZK01

正长岩和二长岩的自然伽马平均值明显高于其他岩性，分别为731 API、405 API，按给定强度法计算相当于铀当量0.83×10^-2%和0.46×10^-2%.说明该孔中正长岩和二长岩的铀当量是地壳中相应岩石克拉克值得10倍以上.正长岩和二长岩中存在铀元素富集区的可能性较大.

钻孔内铀钍矿化属于高温阶段矿化，钻孔中铀钍矿物主要为铀钍石和钛油矿及少量的晶质油矿(熊欣等，2013，2014b).正长岩和二长岩富含铀钍资源，有可能成为区内铀矿床的母岩.除应加强在其内外接触带中铀资源勘查外，尤其要注意在正长岩和二长岩成岩后的热液温水活动区的勘查，热液温水活动有可能将母岩中的铀和其他元素溶解和淋滤出来，进而汇集到裂隙中沉淀成矿.

庐枞地区目前已发现品级较富的铀矿，均位于正长岩体外接触带(曹达旺和向铭，2012).ZK01孔深部放射性异常的发现，推动了庐枞地区与盆地内部碱性正长岩有关的Fe-Cu-U矿化研究.这类与正长岩有关的矿化具有其自身鲜明的特征，是长江中下游成矿带中一个新的矿化类型(张舒等，2014).

(1)ZK01孔的主要岩性为粗安岩、凝灰岩、安山玢岩、正长岩和二长岩，从上到下依次为喷出岩和侵入岩，该地区的火山喷发活动与深部岩浆活动大致同时进行，伴随着正长岩中的钍铀矿化和粗安岩中的黄铁矿化、高岭石化和硬石膏化.通过分析喷出岩和侵入岩地层的测井响应特征，完成主要岩性的人工识别与支持向量机判别，建立了钻孔岩性剖面；并根据超声成像测井资料推断南北走向为该地区的最大水平主应力方向.

(2)黄铁矿化粗安岩集中在井深610～1440 m的部位，测井曲线中声波时差减小、密度增大、视电阻率增大；高岭石化粗安岩分布在1200 m以上部位，测井曲线中声波时差增大、密度减小、视电阻率减小、磁化率减小.

(3)正长岩中的铀钍矿化产生的放射性异常在 1500～1900 m分布，其自然伽马值最高达到2000 API，铀矿化累积厚度93.02 m.本类铀钍矿化品位低，但可能会有一定规模，甚至可能有很大的储量，作为未来的利用对象仍值得注意.