Heat flow determination of Archean strata under the karst thermal reservoir of D01 well in Xiong'an New Area
-
摘要:研究目的
目前雄安新区范围内的热流数据多测自于新近系沉积盖层及白云岩地层，对岩溶热储下伏太古界的热流测定研究还鲜见报道，主要受限于以往的地热勘查开发层位主要在1800m以浅的蓟县系岩溶热储，揭露岩溶热储下伏太古界的钻孔稀少。2018年以来，在雄安新区牛驼镇地热田实施了D01深部地热参数井，揭穿了高于庄组岩溶热储，并揭露太古界厚度达1723.67m，给该地区岩溶热储下伏太古界的热流测定研究提供了条件。
研究方法本文基于近稳态钻孔测温及岩心热导率测试，对D01井太古界进行了热流测定及分析。
研究结果结果表明：D01井太古界呈显著的传导型地温特征，地温梯度为18.3℃/km，相比于新近系地温梯度48.6℃/km偏低。测得25块D01井太古界片麻岩岩心样品的热导率平均值为(2.41±0.40)W/(K·m)。根据钻孔测温数据及热导率测试，计算出D01井太古界2300~2700m深度的热流值为(44.1±7.0)mW/m2。同时，估算了D01井新近系400~800m深度段的热流值为84.6mW/m2，较太古界高出40.5mW/m2。
结论分析认为D01井新近系较太古界高出热流主要为高孔渗岩溶白云岩层中地下水侧向热对流作用以及牛东断裂地下水垂向热对流作用共同所致。研究结果对前人提出的地下水运移聚热模式及导热断裂聚热模式提供了数据支撑，对雄安新区深部地热资源成因机制和太古界地温场研究具有重要意义。
创新点：基于钻孔测温和岩心热导率测试，获取雄安新区岩溶热储下伏太古界结晶基底的高质量热流数据，并揭示其现今地温场特征。
Abstract:The paper is the resultof geothermal geology survey engineering.
ObjectiveThe heat flow values in Xiong'an New Area are mostly measured from the Neogene sedimentary cap and dolomite strata. The previous geothermal exploration and utilization horizons are mainly karst thermal reservoirs shallower than 1800m in the Jixian System and the boreholes exposing the Archean are scarce in Xiong'an New Area. Heat flow determinations and analyses from Archean under the karst thermal reservoir are rarely reported by previous researchers. Since 2018, D 01 deep geothermal parameter well with a depth of 3403.67m has been implemented in Niutuozhen geothermal field in Xiong'an New Area. The karst thermal reservoir of the whole Gaoyuzhaung Formation was exposed and the Archean of 1723.67m was drilled by D01 well. The D01 well provides the conditions for the heat flow study from Archean under the karst thermal reservoir in this area.
MethodsBased on borehole temperature measurement and thermal conductivity tests, the study on the heat flow from the Archean of D01 well was carried out.
ResultsThe results showed that the Archean of D01 well has significant geothermal characteristics of conduction type, and the geothermal gradient is 18.3℃/km, which is lower than the Neogene geothermal gradient of 48.6℃/km. The thermal conductivity values of 25 Archean gneiss core samples of D01 well were measured with an average value of (2.41±0.40) W/(K·m). Based on the temperature measurement curve and thermal conductivity data of 2300-2700m section of D01 well, the Archean heat flow of D01 well was calculated to be (44.1±7.0) mW/m2. The heat flow of the Neogene sedimentary cap of D01 well was estimated to be 84.6 mW/m2.
ConclusionsThe higher part 40.5 mW/m2 than Archean of D01 well was considered to be mainly contributed by the lateral convection heat transfer of groundwater in high pore-permeable karst dolomite and the vertical convection heat transfer of groundwater in the Niudong fault. The research results can provide the basic data for the previous heat accumulation mechanism of the groundwater transport and the heat conduction fault. And the results are of great significance to the genetic mechanism of deep geothermal resources and the geothermal field of Archean in Xiong'an New Area.
-
1. 引言
找矿模型，又称勘查模型，是表述目标矿床及其对一定勘查技术产生预期效应的模拟体系，利用已有地质信息及现有找矿手段，并通过找矿模型来实现成矿的预测。找矿模型的建立不仅需要总结地质直接找矿标志，还需要借助地球物理、地球化学等间接找矿标志识别成矿信息，推断矿体空间位置（陈毓川和朱裕生，1993）。陈毓川和朱裕生（1993）指出找矿模型突出的是某类矿床的基本要素和找矿过程中特殊意义的地质、物化探和遥感影像等特征及其在空间的变化情况，总结发展该类矿床的基本标志和找矿使用的方法手段。找矿模型是矿产勘查工作的实际指导，它是缩小勘查区（或靶区）甚至发现矿床的择优技术。当前，地质找矿已由浅部逐步转变为深部矿、隐伏矿，其找矿难度也逐步增大，对隐伏矿预测定位及找矿模型的建立显得尤为关键。
20世纪80年代以来，中国的地质工作者开始探索建立找矿模型的技术，来处理和分析各类地学资料进行综合分析，为普查和评价矿产资源提供依据，根据所利用的地学数据的不同来源，综合研究成果分别表述为矿田或矿床的地球物理模型、地球化学模型、地质-地球物理-地球化学模型或物理-地质模型。20世纪90年代以来，中国许多研究者对找矿模型的建模技术进行了系统的总结和推广。在寻找隐伏矿床的新方法上，多位研究者也建立了特有的找矿模型，如活动金属离子法、金属活动态法(MOMEO、(MMI) (Mann et al., 1998)、地气法(geogas) (Malmqvistet et al., 1984)。姚金炎(1990)提出综合物探、化探的隐伏矿床预测方法生物地球化学法、离子晕法、电吸附（周奇明，1996）、植物地球化学、地电化学等。
箭猪坡矿床作为五圩矿田矿集区代表性矿床之一，该矿床于20世纪80年代被发现，先后进行了普查和详查工作，经过20多年的开采，矿山资源量已不能满足矿山开采需求，急需开展地质找矿工作，增加储量。前人对该矿床的地质特征、矿床成因，流体特征等方面开展了较多的研究，普遍认为其受构造控制的岩浆期后热液脉状矿床。然而，箭猪坡矿床乃至五圩矿田的同类型矿床至今缺少找矿模式的研究。笔者所在项目组2012—2016年先后对箭猪坡矿床开展了构造控矿规律、矿体定位预测及外围探矿权找矿预测研究等工作，首次在五圩矿田箭猪坡矿床发现了石英脉型锡矿体，矿体呈脉状，倾角70°~75°，矿体平均品位：锡0.65%，铅0.56%、锑0.73%、银46 g/t，锡石主要呈棕褐色，多沿裂隙发育，呈条带状。经锡物相分析，锡以锡石为主。并对紧邻采矿权外围的花洞、板才探矿权开展了土壤地球化学测量、地电化学测量、瞬变电磁测量、可控源测量工作，大大丰富了找矿信息。本文在综合分析箭猪坡矿床地质、物探、化探特征基础上，构建了箭猪坡锑多金属矿地质-物探-化探测量综合信息找矿模型(经验证找矿效果良好)，以期为箭猪坡锑多金属矿集区热液脉型矿床及外围同类型矿床找矿勘查和评价工作提供重要参考。
2. 矿床地质特征
五圩矿田位于丹池多金属成矿带南端（图 1）。区内出露的地层为下泥盆统塘丁组(D1t)—中三叠统(T2b)；其中赋矿层位为下泥盆统塘丁组条带状绢云母泥岩夹砂岩（箭猪坡矿床），中泥盆统罗富组泥灰岩（三排洞—芙蓉厂矿床），上泥盆统榴江组、五指山组硅质岩。区域构造较为复杂，主要由五圩复式背斜及一系列以NNW向为主的张扭性断裂组成。北北西走向断裂带是五圩矿田内最发育的一组构造形迹，呈帚状构造的形态，向北北西向（箭猪坡—芙蓉厂）逐渐撒开，在花洞南东向成逐渐收敛状态，断裂带控制着区域矿化带的形成与发育。区域上未见岩浆岩出露。
![]()
图 1 五圩矿田区域地质简图T—三叠系(粉砂质泥岩、泥质粉砂岩、灰岩）；P—二叠系（灰岩、凝灰岩、硅质岩）；C—石炭系（条带灰岩、微晶灰岩、泥岩）；D—泥盆系（泥岩、灰岩、硅质岩）; 1—背斜轴；2—向斜轴；3—断裂；4—铅锌锑银矿床；5—砷汞铅锌锑银矿床；6—锑矿点；7—汞矿床（点）；8—矿田Figure 1. Regional geologic map of Wuxu mining areaT-Triassic system (silty mudstone, argillaceous siltstone, limestone); P-Permian System (Limestone, Tuff, Siliceous Rock); C-Carboniferous system (banded limestone, microcrystalline limestone, mudstone); D-Devonian System (Mudstone, Limestone, Siliceous Rock); 1-Anticlinal Axis; 2-Syncline axis; 3-Fracture; 4-Pb-Zn-Sb-Ag deposit; 5-Arsenic-mercury-lead-zinc-antimony-silver deposit; 6-Antimony ore site; 7-Mercury deposit (spot); 8-Ore field箭猪坡矿床范围内出露的地层为下泥盆统塘丁组第一段和第二段，其中以第二段为主，该段为灰黑色条带状及薄层状绢云母泥岩；矿床内发育了两组断裂，一组为北北西向张扭性断裂，第二组为北北东向扭性断裂。在矿区范围内未见到岩浆岩出露。
箭猪坡矿床矿体主要受北北西向构造破碎带和断裂的控制，呈脉状、细网脉状，条带状充填在断裂带中或节理裂隙中，沿矿体走向和倾向分支复合、尖灭再现和尖灭侧现的现象均有出现。矿石中金属矿物主要有硫化物（铁闪锌矿、闪锌矿、辉锑矿）、硫盐矿物（脆硫锑铅矿等），碳酸盐矿物(含锰菱铁矿、含铁菱锰矿)，非金属矿物主要有石英、方解石等（刘伟，2013）。箭猪坡矿床成矿阶段划分为3个阶段，即石英-硫化物阶段(Ⅰ)、石英-硫盐矿物-多硫化物阶段(Ⅱ)和锑硫化物-石英-铁锰碳酸盐阶段(Ⅲ)，第Ⅱ成矿阶段为主要成矿阶段（刘伟等，2015）。
3. 地球物理特征
本次研究工作选用的物探工作方法有瞬变电磁测量(TEM)、可控源音频大地电磁测深法（CSAMT）。物探测线布置图见图 2。
![]()
图 2 箭猪坡矿床物化探工作测线分布图1—三叠系百逢组中段（细砂岩、泥质粉砂岩）; 2—三叠系百逢组下段（含粉砂质泥岩夹凝灰质泥岩）; 3—中三叠统罗楼祖上段（粉砂质泥岩夹灰岩）; 4—中三叠统罗楼组下段（粉砂质泥岩）; 5—下二叠统茅口组（燧石条带灰岩）; 6—下二叠统栖霞组（灰岩、泥灰岩、硅质岩）; 7—上石炭统南丹组（硅质条带灰岩、微晶灰岩）8—中石炭统大埔组（白云岩、含生物硅质条带灰岩）; 9—下石炭统巴平组（泥岩、硅质泥岩）; 10—上泥盆统五指山组（灰岩）; 11—上泥盆统榴江组（硅质岩、硅质灰岩、灰岩）; 12—中泥盆统罗富组上段（泥质灰岩、生物碎屑灰岩）; 13—中泥盆统罗富组中段（泥灰岩夹泥岩）; 14—中泥盆统罗富组下段（钙质泥岩、泥灰岩）; 15—下泥盆统塘丁组四段（泥岩）; 16—下泥盆统塘丁组三段（泥岩、粉砂质泥岩）; 17—下泥盆统塘丁组二段（泥岩夹砂质泥岩）; 18—下泥盆统塘丁组一段（泥岩）; 19—实测正断层; 20—实测逆断层; 21—岩层产状; 22—地质界线; 23—大型锑多金属; 24—工作区范围; 25—可控源测量测线及编号; 26—瞬变电磁测量、地电化学测量、土壤地球化学测量测线Figure 2. Distribution of geophysical and chemical detection lines in Jianzhupo Deposit1-Middle of Triassic Baifeng Formation (fine sandstone and argillaceous siltstone); 2-Lower of Triassic Baifeng Formation (silty mudstone mixed with tuffaceous mudstone); 3-Upper of Middle Triassic Luolouzu Formation (silty mudstone mixed with limestone); 4-Lower of Middle Triassic Luolou Formation (silty mudstone); 5-Lower Permian Maokou Formation (flint banded limestone); 6-Lower Permian Qixia Formation (limestone, marl, siliceous rock); 7-Upper Carboniferous Nandan Formation (siliceous banded limestone) Microcrystalline limestone); 8-Middle Carboniferous Dapu Formation (dolomite, siliceous banded limestone); 9-Lower Carboniferous Baping Formation (mudstone, siliceous mudstone); 10-Upper Devonian Wuzhishan Formation (limestone); 11-Upper Devonian Liujiang Formation (siliceous rock, siliceous limestone, limestone); 12-Middle Devonian Luofu Formation upper member (argillaceous limestone, Bioclastic limestone); 13-Middle of Middle Devonian Luofu Formation (marl with mudstone); 14-Lower Devonian Luofu Formation (calcareous mudstone and marl); 15-Lower Devonian Tangding Formation 4 (mudstone); 16-Lower Devonian Tangding Formation 3 (mudstone, Silty mudstone); 17-Member 2 of Lower Devonian Tangding Formation (mudstone mixed with sandy mudstone); 18-Member 1 of Lower Devonian Tangding Formation (mudstone); 19-Measured normal fault; 20-Measured reverse fault; 21-Strata occurrence; 22-Geological boundary; 23-Large antimony polymetallic; 24-Scope of work area; 25-Controlled source measuring line and No; 26-Transient electromagnetic measuring line, geoelectrochemical measuring line and soil geochemical measuring line可控源音频大地电磁法(CSAMT)是在大地电磁法(MT)和音频大地电磁法(AMT)基础上发展起来的一种方法。可控源频率测深方法具有勘探深度范围大（通常可达2 km）、分辨能力强、观测效率高，兼有测深和剖面研究双重特点，是研究深部地质构造和寻找隐伏矿的有效手段，越来越多应用于受构造控制的脉状金属矿床的找矿工作。通过可控源音频大地电磁测深法（CSAMT）测量，探索箭猪坡矿床深部隐伏构造的空间形态。
磁瞬变电磁法（TEM）是基于电性差异，利用不接地回线或接地线源向地下发送一次脉冲电场，利用线圈或接地电极观测二次涡流场或电场的方法。该方法具有穿透低阻覆盖的能力，故探测深度大。在低阻围岩地区地形的起伏仅对早期道有影响，这种影响又是易于分辨的。瞬变电磁法异常响应强，多测道曲线形态简单，分层能力较强。工作效率高，同时完成剖面测量与测深工作。通过瞬变电磁测量(TEM)，查明构造和矿（化）体空间展布特征。
3.1 矿区岩（矿）石物性特征
在物探方法工作前，先对矿区内的主要岩矿石进行物性参数取样分析。矿区内出露的主要矿石划分为5类，围岩划分为7类。每类岩（矿）石进行系统的物性参数取样。所采集标本均来自新鲜露头及坑道，大小约10 cm×10 cm×10 cm，使用测试方式为强迫电流法，获得岩矿石标本的电阻率参数，使用GDD型标本测试仪，测定结果如表 1所示。根据表 1，矿区内与低阻因素有关的是断裂破碎带、矿化蚀变及炭质岩层，因此，圈定低阻异常区是本次电法工作的重点。
表 1 五圩铅锌锑矿岩石的物性特征Table 1. Physical properties of Wuxu Pb-Zn-Sb ore and rocks
（1）密度特征
区域出露的地层从下泥盆统塘丁组至中三叠统均有，是高密度层，以泥岩、碳酸盐岩为主，其岩性主要有泥岩、白云质泥灰岩、泥质粉砂岩、炭质泥岩、硅质泥岩、泥灰岩、白云岩、硅质岩、硅质灰岩、生物碎屑灰岩等，以上岩性密度值相对较高且较为稳定，一般为2.71 g/cm3左右。
根据实测及前人的资料，五圩矿区的各类岩石中，以松散的砂岩、砾岩等密度最低，一般为2.48 g/cm3，这些岩层的厚度不大；泥岩、碳酸盐岩建造为主的岩层，其厚度可达数千米，密度平均值为2.71 g/cm3，构成高密度背景；花岗岩类岩石密度平均值为2.61 g/cm3，与泥岩、碳酸盐岩为围岩的密度存在0.1 g/cm3左右的差异，当花岗岩体具有一定规模时，可以引起重力低异常。
（2）电性特征
矿床中岩石在宏观上电阻率ρ与电极化率η呈负相关，即η增大，ρ有减小的趋势。矿床中的炭质泥岩、钙质泥岩具有“低电阻率”特征，电阻率常见为5.2×101~5.1×102 Ω·m，最大为3.9×103 Ω·m，而矿床中方解石、白云岩、硅质泥岩、灰岩具有“高电阻率”特征，电阻率常见为4.1×103~7.3×105 Ω·m，最大为9.8×105 Ω·m。含矿岩石具有低电阻、高极化率、高密度特征。矿床围岩普遍具有中高电阻率、低极化率、中等密度特征。矿区内炭质泥岩普遍发育，其具有低电阻率、高极化率特征，且具有成层出现的特点，其矿石与围岩存在较明显的电性差异，为电磁法测量工作提供了有利条件。
（3）磁性特征
工作区内绝大部分岩石为无磁性或弱磁性，花岗岩呈弱磁性，一些矿石含有少量磁性矿物磁性相对较强，但仍呈弱磁性。
3.2 可控源（CSAMT）特征
结合矿床地质特征，笔者在五圩箭猪坡矿床开展了可控源剖面测量工作（位置见图 2），试验剖面有3条，测量方位线与勘探线方位一致为80°，点距40 m。投入本次工作的仪器为GDP-32II多功能电磁仪。工作选择的是GGT-30大功率（30 kAV）发射机，工作频率：DC~8 kHz；最大输出功率：30 kW；最大输出电流：30 A；最高输出电压：1000 V；关断时间：≤125 μs；稳流精度：0.1%；同时配合工作的是一个ZMG-10(10 kW)发电机。CSAMT资料反演软件是GDP-32II机器所带的带地形平面波反演软件SCS2D。
本次CSAMT工作各技术参数如下：工作频率：1~8192 Hz。发射系统：供电偶极子长度AB=1000 m；方位角117°，偶极子中心与测线中心偏差＜15°。收发距R最小为7 km，最大为12.7 km；基本电流大于6 A，1024 Hz以下频率电流大于8 A。供电线使用耐高压、大电流的电缆线，供电极则将大张铝箔置于挖好的坑中埋好，并将盐水浇于其中。接收系统：采用一个磁道带两个电道测量方式。测量电偶极子长度MN=40 m。接收电极使用不极化电极，稳固地置于土中20~30 cm，并浇入盐水。观测数据叠加次数为32~16384次，数据离差SEM < 20时，方可认定为可靠数据，并将其记录。
由图 3（317线CSAMT测深）可以看出，物探CSAMT测深曲线圆滑，其卡尼亚电阻率均方相对误差为5.783%；大于200 mrad的阻抗相位均方相对误差为8.698%，小于200 mrad的阻抗相位均方误差为40.997 mrad，测量结果为Ⅰ级精度。根据测区岩石电性特征，矿区内与低阻因素有关的是断裂破碎带、矿化蚀变及炭质岩层，由此圈定低阻异常区是本次电法工作的重点。根据电阻率剖面结果，反演的低阻异常，对应的阻抗相位高值异常，因此低阻异常可信度高。
![]()
图 3 箭猪坡矿床CSAMT测深视电阻率曲线图Figure 3. Apparent resistivity curve of CSMAT sounding in Jianzhupo Deposit以317号勘探线剖面图为例（图 4）。选取该线的3160~4520点号，作为反演剖面，长度1.36 km，317线低阻异常及推断解释见图 4，其特征如下：3840~4080号点范围存在⑤号低阻异常，电阻率低于100 Ω·m，呈带状分布，异常上窄下宽，向东倾斜，延深约1200 m，⑤号低阻异常位于下泥盆统塘丁组第二岩性段（D1t2），主要为灰—深灰色中薄层泥岩夹深灰色中薄层条带状泥岩，局部夹有薄层砂岩透镜体和砂质泥岩，岩石中普遍见有浸染状黄铁矿，该岩性段为矿区主要含矿层位。已知含矿破碎带F3及F7位于⑤号低阻异常边部或异常内，且有钻孔ZK317-1揭露，其破碎带产状与⑤号低阻异常形态一致，推测⑤号低阻异常是由含矿破碎带F3和F7引起。目前钻孔只揭露了含矿破碎带F3及F7浅部垂深200 m范围的情况，从电阻率剖面图看出，⑤号低阻异常往深部延深至标高-700 m，说明含矿破碎带F3及F7往深部仍有延伸，该异常深部仍具有较大的找矿空间。
![]()
图 4 箭猪坡矿床317号勘探线电阻率剖面与推断解释剖面图1—泥盆系下统塘丁组第三段；2—泥盆系下统塘丁组第二段；3—断裂；4—钻孔；5—电阻率异常带Figure 4. Resistivity profile and inferred interpretation profile of 317 geological profile in Jianzhupo Deposit1-Third Member of Tangding Formation of Lower Devonian system; 2-Second Member of Tangding Formation of Lower Devonian system; 3-Fault; 4-Drill holel; 5-Resisttivity anomaly band从以上结果看出，CSAMT法探测深度大，对断裂构造反应较灵敏，在五圩矿田可以取得较好的找矿效果。
3.3 瞬变电磁特征
本次TEM法采用重叠回线装置（接收线圈Rx和发送线圈Tx相重合敷设的装置），发射和接收线圈均为100 m×100 m，网度为100 m×50 m；通过试验应采用1 Hz和4 Hz作为发射基频，观测参数为归一化的感应电动势，单位为μV/A。本次使用仪器为重庆奔腾数控技术研究所研制的WTEM-1瞬变电磁仪。以箭猪坡矿床320号勘探线剖面为例（测线分布见图 2），运用一维电阻率反演（图 5），结合地质特征分析得出：通过320线多测道剖面图可见，在23~34点之间存在凸起异常，与相对低阻异常吻合较好，从形态上看，呈脉状，表明在中深部存在导电性较好的脉状地质体，反映了中深部硫化物矿（化）体的存在。从电阻率反演图中可以得出，标高200~300 m，异常带晕为一个呈现倾向东的低阻带，电阻率值一般小于100 Ω·m，多为0~50 Ω·m，倾角为50~60°，与高阻区域区别明显，且该异常带由近地表往深部有变宽的趋势。该地段构造发育，结合矿区岩矿石物性特征及矿区的地质特征推测为含金属硫化物的岩（矿）石引起的异常。通过在地表实施探槽及钻孔（ZK320-2、ZK320-3、ZK320-4、ZK320-5）验证，均揭露到含黄铁矿化的硅化破碎带或矿（化）体，已知硅化破碎带或矿（化）体，均分布在电阻率低阻区。
4. 地球化学特征
4.1 土壤地球化学特征
在箭猪坡矿区开展了土壤剖面地球化学测量工作，以箭猪坡矿床320号勘探线剖面为例，结果显示矿区土壤地球化学化学主成矿元素Pb、Zn、Sb异常与隐伏矿体（含矿破碎带）的分布范围相一致，各异常同步性较好，并且有很好的重叠性。Pb最高值为97.4×10-6，异常带Pb均值大于65×10-6；Zn最高值为255.9×10-6，异常带均值大于220×10-6；Sb最高值为154.8×10-6，异常带均值大于90×10-6；推测异常由破碎带（含矿）引起，并经钻探工作验证，在钻孔中见到真厚0.71m，Zn品位5.05%的锌矿体(钻孔ZK320-2、J144)。
4.2 地电化学特征
地电化学提取方法是将地球物理、地球化学及电化学综合交叉为一体的综合找矿方法，该方法利用电场作用，选择性提取近地表介质中的电活动态物质，通过研究电提取元素组合、含量分布及异常特征，进而提供找矿信息的一种勘查方法。而大量的找矿勘探实践证明，地电化学勘探新技术、新方法能够发现500 m以下，其中包括150~200 m浮土以下的有色金属矿床、贵金属矿床、稀有金属矿床。在箭猪坡矿床南端开展了6条地电化学剖面测量工作（测线分布见图 2），方位线与勘探线方位一致为80°，点距一般20 m，以320号勘探线剖面地电化学特征为例（图 5）：地电化学异常与隐伏矿体或含矿破碎带的分布吻合度高，即矿体的分布范围与地电化学综合异常范围一致，异常分布的规模、强度与含矿破碎带规模呈正比关系。矿体产出均较陡，倾角75°~80°，异常为多个异常点组成的“驼峰”或者“倒钟”形态，在地表矿头对应部位附近有异常峰。Pb、Zn、Sb元素异常规模大、强度高，异常同步清晰，因此Pb、Zn、Sb可以作为本区地电化学方法寻找隐伏锌矿的指示元素。当出现Pb、Zn、Sb为主异常时，则显示深部具有隐伏矿体（矿化）存在的可能性。经钻探验证，ZK320-2、ZK320-4、ZK320-5钻孔均揭露到了含锌锑的矿体及富含黄铁矿的硅化破碎带，亦说明地电化学能很好的圈定金属硫化物异常。
![]()
图 5 箭猪坡矿床320号勘探线地质-物探（TEM）-化探（地电化学）综合剖面图1—构造破碎带；2—铅锌锑矿体；3—下泥盆统塘丁组第二段；4—钻孔Figure 5. Electrogeochemical comprehensive profile of 320 exploration line in Jianzhupo Deposit1-Tectonic fracture zone; 2-Lead-zinc-antimony ore body; 3-Second member of Tangding Formation of Lower Devonian system; 4-Drill hole5. 综合找矿模型特征
根据箭猪坡矿床地质特征、地球物理特征、地球化学特征，总结了该矿床地质、地球物理，地球化学的找矿标志（表 2），并建立了矿床的地质-化探-物探找矿模型（图 6），实质是反映化探特征、物理特征与地质构造、矿体的空间关系。
表 2 箭猪坡矿床地质-地球物理-地球化学找矿标志Table 2. Geological-geophysical -geochemical prospecting criteria of JianZhuPo Deposit
![]()
图 6 箭猪坡矿床地质-地球化学-地球物理找矿模型1—条带状泥岩；2—破碎带；3—矿（化）体Figure 6. Prospecting model of geological - geochemical-geophysical in Jianzhupo Deposit1-Ribbon mudstone; 2-Broken belt; 3-Mineralized bodies从矿床成矿条件分析，箭猪坡矿床矿体赋存在断裂构造中，矿体形态呈脉状、透镜体状，空间分布具有斜列状，叠瓦状排列等。区内矿体具有明显的低阻异常，表现为低电阻率特征，且瞬变电磁反演断面图与硅化破碎带及地质内容吻合较好。
模型的地质特征：矿体受层位控制，控矿层位为下泥盆统塘丁组沉积建造，因构造变形，控矿层位位于褶皱翼部，产状一般较陡；近矿围岩蚀变发育范围有限，蚀变类型有硅化、黄铁矿化、绢云母化、绿泥石化等。矿体形态以陡产状脉状、细脉状、透镜状为主，倾角60°~80°，厚度0.1~4.5 m不等，最大10.72 m，矿石类型分致密块状和条带状；矿石以金属硫化物为主，并含硫盐矿物；矿床类型属低温热液充填型。
模型的地球物理特征：以工业矿体为中心的矿化系统，具有低电阻率，稍高密度，高极化率特征。
模型的地球化学特征：常规土壤异常较小或无异常，而地电化学异常能较清晰地反映盲矿体成矿元素的异常，且各元素间连续性较好，异常的分布范围与深部盲矿体的垂直投影基本一致，通过钻孔揭露，能较准确的反映深部埋藏盲矿体（矿化体）的分布范围。
6. 讨论
（1）通过近几年的找矿勘查，在箭猪坡矿床，矿带走向延伸控制总长度为1400 m，并根据物探可控源测量结果，经钻探工程控制，在标高为-300~-400 m）仍能见到矿体, 从而扩大了该区的找矿空间，在该矿带南、北两个延伸方向及东部矿体倾向延深方向都是寻找隐伏矿体的有利部位，勘查平面图、剖面图等资料表明矿体的倾向延深和南北走向延伸还未完全控制，尚有较大找矿空间。而且随着深部石英脉型锡矿体的发现，为该矿区增加了新的找矿方向。
（2）矿区中局部构造陡倾或直立，且构造中有炭质、硅质等围岩，具有极低或极高的电阻率，故引起了CSAMT测量反演结果出现电阻率非常大和非常小的值。
（3）受整个区围岩炭质泥岩背景影响，该区矿化不均匀，极化率值波动较大，矿区未采样激发极化法, 其还难以识别矿异常和非矿异常。
因此，对于箭猪坡矿床成因需要加强研究，勘查模型需要逐步完善，找矿方法手段的选择也需要不断探索。
7. 结论
物化探工作是目前找矿工作中的重要手段。在充分分析、深入研究物化探成果特征，选出与成矿有关的异常范围，尤其是物化探重叠异常范围，再以地质条件为基础，化探信息为先导，选取地质、化探、物探等成矿有关的综合信息标志重叠的地区或地段，是寻找成矿有利地段的保障。
（1）物探研究表明：CSAMT出现的所有低阻异常，均具有找矿意义。这类低阻异常范围、形态大致反映矿（化）脉组在空间的展布特征，为深部勘查工程布置提供了重要依据；CSAMT异常显示该区往南北两端仍具有较大的找矿空间，中深部层状低阻异常出现在较低标高，均具有较大规模，且没有完全控制，这类异常推测与炭质层或新类型矿化有关。瞬变电磁表明：矿（化）体多分布在归一化电位高值区或低电阻率区，已知构造主要分布在归一化电位较高值区域或较低电阻率区和电位值梯度带，已知矿化体与物探异常特征吻合较好。
（2）常规的化探次生晕及地电化学研究表明：地电化学异常在已知矿体上方有比较清晰的Pb、Zn、Sb元素异常，为典型的共生组合元素，它们之间也存在良好的正相关关系，说明这些元素组合可以作为找矿指示元素，在该地区利用地电化学在该地区寻找隐伏矿，方法是可行的。化探次生晕Sb、Pb、Zn等元素套和较好，异常特征相近，可以作为寻找本区矿区的主要指示元素，进而推断矿体存在的大致位置和特征。
综上所述，可控源音频大地电磁法、瞬变电磁法、土壤地球化学、地电化学在异常区和非异常区有效性明显，利用这四种集成方法对深部勘查预测具有良好的找矿效果，达到技术集成示范，值得在类似的矿床中应用。
致谢: 本次研究得到中国科学院地质与地球物理研究所姜光政博士、王朱亭博士的帮助和建议，在此表示诚挚的感谢！ -
![]()
图 1 渤海湾盆地构造单元图(a)和牛驼镇凸起基岩地质简图(b)(改自马敬业等，1990)
Figure 1. Tectonic units of Bohai Bay basin (a) and geological sketch of bedrock in Niutuozhen uplift (b)(modified from Ma Jingye et al., 1990)
![]()
图 4 新近系(a)和太古界(b)热流计算段测温数据线性拟合
Figure 4. Linear fitting of temperature measurement data of Neogene (a) and Archean (b) heat flow calculation section
![]()
图 5 DR井静井261 d后测温曲线(Jiang et al., 2016)
Figure 5. Temperature curve of DR well for 261 days after the cessation of drilling(Jiang et al., 2016)
![]()
图 6 牛驼镇凸起地热地质概念模式图(地层剖面改自王贵玲等，2018)
Figure 6. The geothermal geological conceptual model of Niutuozhen uplift (stratigraphic section modified from Wang Guiling et al., 2018)
表 1 D01井太古界片麻岩岩心热导率测试统计
Table 1 Thermal conductivity of Archean gneiss core samples in D01 well
下载: 导出CSV
-
Chen Moxiang, Huang Geshan, Zhang Wenren, Zhang Rongyan, Liu Bingyi. 1982. The temperature distribution pattern and the utilization of geothermal water at Niutuozhen basement protrusion of central Hebei Province[J]. Chinese Journal of Geology, (3): 239-252(in Chinese with English abstract).
Dai Minggang, Wang Xinwei, Liu Jinxia, Lei Haifei, Bao Zhidong. 2019. Characteristics and influence factors of geothermal resources in the starting and adjacent zone of Xiong'an New Area[J]. Chinese Journal of Geology, 51(1): 176-191(in Chinese with English abstract).
Furlong K P, Chapman D S. 2013. Heat flow, heat generation, and the thermal state of the lithosphere[J]. Annual Review of Earth and Planetary Sciences, 41(1): 385-410. doi: 10.1146/annurev.earth.031208.100051
Gong Yuling. 2003. The Thermal Structure and Thermal Evolution of Bohai Bay in East China[D]. Nanjing: Nangjing University (in Chinese with English abstract).
Guo Sasa, Zhu Chuanqing, Qiu Nansheng, Tang Boning, Cui Yue. 2020. Formation conditions and favorable areas for the deep geothermal resources in the Xiong'an New Area[J]. Acta Geologica Sinica, 94(7): 2026-2035(in Chinese with English abstract). doi: 10.3969/j.issn.0001-5717.2020.07.011
He Lijuan, Hu Shengbiao, Yang Wencai, Wang Jiyang, Yang Shuchun, Yuan Yusong, Cheng Zhenyan. 2006. Temperature measurement in the main hole of the Chinese Continental Scientific Drilling[J]. Chinese Journal of Geophysics, 49(3): 745-752(in Chinese with English abstract). doi: 10.3321/j.issn:0001-5733.2006.03.018
He Lijuan, Hu Shengbiao, Huang Shaopeng, Yang Wencai, Wang Jiyang, Yuan Yusong, Yang Shuchun. 2008. Heat flow study at the Chinese Continental Scientific Drilling site: Borehole temperature, thermal conductivity, and radiogenic heat production[J]. Journal of Geophysical Research-Solid Earth, 113: B02404.
He Lijuan, Wang Jiyang. 2021. Concept and application of some important terms in Geothermics and Geophysics such as terrestrial heat flow[J]. China Terminology, 23(3): 3-9(in Chinese with English abstract).
Jiang Guangzheng, Gao Peng, Rao Song, Zhang Linyou, Tang Xiaoyin, Huang Fang, Zhao Ping, Pang Zhonghe, He Lijuan, Hu Shengbiao, Wang Jiyang. 2016. Compilation of heat flow data in the continental area of China (4th edition)[J]. Chinese Journal of Geophysics, 59(8): 2892-2910(in Chinese with English abstract).
Jiang Guangzheng, Tang Xiaoyin, Rao Song, Gao Peng, Zhang Linyou, Zhao Ping, Hu Shengbiao. 2016. High-quality heat flow determination from the crystalline basement of the south-east margin of north china craton[J]. Journal of Asian Earth Sciences, 118, 1-10. doi: 10.1016/j.jseaes.2016.01.009
Li Weiwei, Rao Song, Tang Xiaoyin, Jiang Guangzheng, Hu Shengbiao, Kong Yanlong, Pang Jumei, Wang Jiancun. 2014. Borehole temperature logging and temperature field in the Xiongxian geothermal field, Hebei Province[J]. Chinese Journal of Geology (Scientia Geologica Sinica), 49(3): 850-863(in Chinese with English abstract).
Ma Feng, Wang Guiling, Zhang Wei, Zhu Xi, Zhang Hanxiong, Sun Zhanxue. 2021. Influence mechanism of ancient buried hill geothermal development on land subsidence[J]. Geology in China, 48(1): 40-51(in Chinese with English abstract).
Ma Jingye, Cai Hongtao, Dong Jinsheng. 1990. Investigation Report on Niutuozhen Geothermal Field, Hebei Province[R]. Hengshui: The 3rd Hydrogeology and Engineering Group of Hebei Bureau of Geology and Mineral Resources, 1-299(in Chinese).
Norden B, Förster A, Förster H, Fuchs S. 2020. Temperature and pressure corrections applied to rock thermal conductivity: Impact on subsurface temperature prognosis and heat-flow determination in geothermal exploration[J]. Geothermal Energy, 8(1): 1-19. doi: 10.1186/s40517-020-0157-0
Pang Zhonghe, Pang Jumei, Kong Yanlong, Luo Lu, Duan Zhongfeng, Yang Fengtian, Wang Shufang. 2015. Large karstic geothermal reservoirs in sedimentary basins in China: Genesis, energy potential and optimal exploitation[C]//World Geothermal Congress 2015.
Pang Zhonghe, Kong Yanlong, Pang Jumei, Hu Shengbiao, Wang Jiyang. 2017. Geothermal resources and development in Xiong'an New Area[J]. Bulletin of Chinese Academy of Sciences, 32(11): 1224-1230(in Chinese with English abstract).
Pribnow D, Williams C F, Sass J H, Keating R. 1996. Thermal conductivity of water-saturated rocks from the KTB Pilot Hole at temperatures of 25 to 300℃[J]. Geophysical Research Letters, 23(4): 391-394. doi: 10.1029/95GL00253
Rao Song, Hu Shengbiao, Zhu Chuanqing, Tang Xiaoyin, Li weiwei, Wang Jiyang. 2013. The characteristics of heat flow and lithospheric thermal structure in Junggar Basin, northwest China[J]. Chinese Journal of Geophysics, 56(8): 2760-2770(in Chinese with English abstract).
Sass J H, Lachenbruch A H, Moses T H, Morgan Paul. 1992. Heat flow from a scientific research well at Cajon Pass, California[J]. Journal of Geophysical Research, 97(B4): 5017-5030. doi: 10.1029/91JB01504
Seipold U, Huenges E. 1998. Thermal properties of gneisses and amphibolites-high pressure and high temperature investigations of KTB-rock samples[J]. Tectonophysics, 291(1/4): 173-178.
Wang Jiyang, Hu Shengbiao, Cheng Benhe, Li Tiejun, Yang Wencai, Cheng Zhenyan. 2001. Predication of the deep temperature in the target area of the china continental scientific drilling[J]. Chinese Journal of Geophysics, 44(6): 774-782(in Chinese with English abstract). doi: 10.3321/j.issn:0001-5733.2001.06.006
Wang Jiyang, Huang Shaopeng. 1988. Compilation of heat flow data in the China Continental area[J]. Chinese Journal of Geology (Scientia Geologica Sinica), (2): 196-204(in Chinese with English abstract).
Wang Guiling, Li Jun, Wu Aimin, Zhangwei, Huqiuyun. 2018. A study of the thermal storage characteristics of Gaoyuzhuang formation, A new layer system of thermal reservoir in Rongcheng uplift area, Hebei Province[J]. Acta Geoscientica Sinica, 39(5): 533-541(in Chinese with English abstract).
Wang Kai, Zhang Jie, Bai Dawei, Wu Xingang, Yue Hangyu, Zhang Baowei, Wang Xiaojiang, Zhang Kai. 2021. Geothermal-geological model of Xiong'an New Area: Evidence from geophysics[J]. Geology in China, 48(5): 1453-1468(in Chinese with English abstract).
Wang Zhuting, Zhang Chao, Jiang Guangzheng, Hu Jie, Tang Xianchun, Hu Shengbiao. 2019. Present-day geothermal field of Xiong'an New Area and its heat source mechanism[J]. Chinese Journal of Geophysics, 62(11): 4313-4322(in Chinese with English abstract). doi: 10.6038/cjg2019M0326
Wu Aimin, Ma Feng, Wang Guiling, Liu Jinxia, Hu Qiuyun, Miao Qingzhuang. 2018. A study of deep-seated karst geothermal reservoir exploration and huge capacity geothermal well parameters in Xiong'an New Area[J]. Acta Geoscientica Sinica, 39(5): 523-532(in Chinese with English abstract).
Xu Ming, Zhao Ping, Zhu Chuanqing, Shan Jingnan, Hu Shengbiao. 2010. Borehole temperature logging and terrestrial heat flow distribution in Jianghan Basin[J]. Chinese Journal of Geology (Scientia Geologica Sinica), 45(1): 317-323(in Chinese with English abstract).
Xu Ming, Zhu Chuanqing, Tian Yuntao, Rao Song, Hu Shengbiao. 2011. Borehole temperature logging and characteristics of subsurface temperature in the Sichuan Basin[J]. Chinese Journal of Geophysics, 54(4): 1052-1060(in Chinese with English abstract). doi: 10.3969/j.issn.0001-5733.2011.04.020
Yue Gaofan, Wang Guiling, Ma Feng, Zhu Xi, Zhang Hanxiong. 2021. Evaluation of fault slip probability of geothermal large-scale development: A case study of deep karst geothermal reservoir in Xiong'an New Area[J]. Geology in China, 48(5): 1382-1391(in Chinese with English abstract).
Zhao Li, Li Li. 2016. The extensional pattern and dynamics of Bohai Bay basin in Late Mesozoic-Cenozoic[J]. Geology in China, 43(2): 470-485(in Chinese with English abstract). doi: 10.3969/j.issn.1000-3657.2016.02.009
陈墨香, 黄歌山, 张文仁, 张容燕, 刘炳义. 1982. 冀中牛驼镇凸起地温场的特点及地下热水的开发利用[J]. 地质科学, (3): 239-252. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKX198203000.htm 戴明刚, 汪新伟, 刘金侠, 雷海飞, 鲍志东. 2019. 雄安新区起步区及周边地热资源特征与影响因素[J]. 地质科学, 51(1): 176-191. 龚育龄. 2003. 中国东部渤海湾盆地热结构和热演化[D]. 南京: 南京大学. 郭飒飒, 朱传庆, 邱楠生, 唐博宁, 崔悦. 2020. 雄安新区深部地热资源形成条件与有利区预测[J]. 地质学报, 94(7): 2026-2035. doi: 10.3969/j.issn.0001-5717.2020.07.011 何丽娟, 胡圣标, 杨文采, 汪集旸, 杨树春, 袁玉松, 程振炎. 2006. 中国大陆科学钻探主孔动态地温测量[J]. 地球物理学报, 49(3): 745-752. doi: 10.3321/j.issn:0001-5733.2006.03.018 何丽娟, 汪集旸. 2021. "大地热流"等地热学重要术语的概念与应用[J]. 中国科技术语, 23(3): 3-9. https://www.cnki.com.cn/Article/CJFDTOTAL-KJSY202103001.htm 姜光政, 高堋, 饶松, 张林友, 唐晓音, 黄方, 赵平, 庞忠和, 何丽娟, 胡圣标, 汪集旸. 2016. 中国大陆地区大地热流数据汇编(第四版)[J]. 地球物理学报, 59(8): 2892-2910. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201608015.htm 李卫卫, 饶松, 唐晓音, 姜光政, 胡圣标, 孔彦龙, 庞菊梅, 王建存. 2014. 河北雄县地热田钻井地温测量及地温场特征[J]. 地质科学, 49(3): 850-863. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKX201403013.htm 马峰, 王贵玲, 张薇, 朱喜, 张汉雄, 孙占学. 2021. 古潜山热储开发对地面沉降的影响机制研究[J]. 中国地质, 48(1): 40-51. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI202101004.htm 马敬业, 蔡洪涛, 董金生. 1990. 河北省牛驼镇地热田勘查报告[R]. 衡水: 河北省地矿局第三水文工程地质大队, 1-299. 庞忠和, 孔彦龙, 庞菊梅, 胡圣标, 汪集旸. 2017. 雄安新区地热资源与开发利用研究[J]. 中国科学院院刊, 32(11): 1224-1230. doi: 10.16418/j.issn.1000-3045.2017.11.007 饶松, 胡圣标, 朱传庆, 唐晓音, 李卫卫, 汪集旸. 2013. 准噶尔盆地大地热流特征与岩石圈热结构[J]. 地球物理学报, 56(8): 2760-2770. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201308025.htm 汪集旸, 胡圣标, 程本合, 李铁军, 杨文采, 程振炎. 2001. 中国大陆科学钻探靶区深部温度预测[J]. 地球物理学报, 44(6): 774-782. 汪集旸, 黄少鹏. 1988. 中国大陆地区大地热流数据汇编[J]. 地质科学, (2): 196-204. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKX198802010.htm 王贵玲, 李郡, 吴爱民, 张薇, 胡秋韵. 2018. 河北容城凸起区热储层新层系-高于庄组热储特征研究[J]. 地球学报, 39(5): 533-541. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXB201805003.htm 王凯, 张杰, 白大为, 吴新刚, 岳航羽, 张保卫, 王小江, 张凯. 2021. 雄安新区地热地质模型探究: 来自地球物理的证据[J]. 中国地质, 48(5): 1453-1468. http://geochina.cgs.gov.cn/geochina/article/abstract/20210511?st=search 王朱亭, 张超, 姜光政, 胡杰, 唐显春, 胡圣标. 2019. 雄安新区现今地温场特征及成因机制[J]. 地球物理学报, 62(11): 4313-4322. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201911026.htm 吴爱民, 马峰, 王贵玲, 刘金侠, 胡秋韵, 苗青壮. 2018. 雄安新区深部岩溶热储探测与高产能地热井参数研究[J]. 地球学报, 39(5): 523-532. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXB201805002.htm 徐明, 赵平, 朱传庆, 单竞男, 胡圣标. 2010. 江汉盆地钻井地温测量和大地热流分布[J]. 地质科学, 45(1): 317-323. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKX201001027.htm 徐明, 朱传庆, 田云涛, 饶松, 胡圣标. 2011. 四川盆地钻孔温度测量及现今地热特征[J]. 地球物理学报, 54(4): 1052-1060. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201104022.htm 岳高凡, 王贵玲, 马峰, 朱喜, 张汉雄. 2021. 地热规模化开发断层滑动概率评估——以雄安新区深部岩溶热储为例[J]. 中国地质, 48(5): 1382-1391. http://geochina.cgs.gov.cn/geochina/article/abstract/20210505?st=search 赵利, 李理. 2016. 渤海湾盆地晚中生代以来伸展模式及动力学机制[J]. 中国地质, 43(2): 470-485. http://geochina.cgs.gov.cn/geochina/article/abstract/20160209?st=search -
期刊类型引用(4)
1. 王学阳，杨言辰，刘志宏，张继武，戴台鹏. 长江中下游成矿带大金山地区综合信息找矿效果及深部找矿潜力分析. 吉林大学学报(地球科学版). 2025(01): 125-138 . 
百度学术
2. 谭杰，徐文杰，赵毅，罗达，赵俊宏，钟玉龙，陶明荣. 广西五圩矿田多金属矿成矿模式及找矿预测. 中国矿业. 2024(S1): 508-514 . 百度学术
3. 李双飞，黄鹂，陈建，朱伟，孙丽莎，王欣，彭永和，商祥鸿，邹占春，唐名鹰，高远，赵家强. 山东省五莲七宝山金铜矿床成矿模式与找矿勘查模型——来自综合物化探的证据. 地质与勘探. 2023(05): 961-973 . 百度学术
4. 李玉坤，王兴龙，廖阿托，李科. 东天山梧桐沟钨矿区综合信息找矿模型. 新疆地质. 2023(04): 530-537 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 1782
- HTML全文浏览量: 867
- PDF下载量: 1700
- 被引次数: 6
