综合地球物理技术在银川盆地东缘地热研究中的应用

虎新军，陈晓晶，仵阳，安百州，倪萍

(1.宁夏回族自治区地球物理地球化学勘查院，宁夏 银川 750001； 2.中国地质大学(武汉) 资源学院，湖北 武汉 430074)

0 引言

我国地热资源类型多样，分布广泛，不同的地质条件造就了不同类型资源。根据赋存埋深和温度将地热资源划分为浅层地热资源、水热型地热资源和干热岩3个大类，再以地热资源富集赋存的关键地质要素为依据，按照源(热源和水源)—通(通道及传输)—储(储集体)—盖(盖层)等地质要素和热的传输、储集、保存、散失等地质作用的特征，将水热型地热资源进一步细分为岩浆型、隆起断裂型和沉降盆地型3亚类[1-3]。

近年来，银川盆地东缘地热勘探取得了新突破，以天山海世界地热井田最具代表性，地热类型显著区别于盆地内部“沉降盆地型”地热，具有“埋藏浅、温度高、水量大、水质好”等优点，是典型的盆地边缘“隆起断裂型”[4]，成为了本地区地热研究的热点。地热研究中的地球物理方法有多种，以地温场调查、重力勘探、磁法勘探、电法勘探、大地面波测深和人工地震勘探[5]为主流手段。早期的研究侧重于单一方法的应用分析，黄力军及胡宁等[6-7]认为花岗片麻岩、火成岩一般都具有较强磁性，通过精度高的磁法测量可确定火山岩、侵入岩及蚀变带的分布位置，间接划定地热分布范围；Pamukcu O A等[8]将自由空气重力异常与航磁异常的明显负相关关系作为了指示区域地下地热远景区的一个重要指标，李学云等[9]根据重力异常分布，发现河北汤泉地区地热热源与深部构造的酸性热液活动有关；刘长生等[10]将MT法运用到辽河凹陷西部，将埋深较大的低阻电性层确定为潜在热储层，并得到了实际钻孔数据的验证，汪琪[11]根据银川平原地区主要的电性剖面，并结合地质剖面分析研究了该地区的热储类型；付微等[12]认为微动勘探方法对深部低速层分辨率较高，对破碎带低速异常体探测效果良好，能有效圈定含水破碎带区域，有助于提高地热勘查精度，刘远[13]通过对地微动探测和大地电磁测深两种方法的研究，探索性地提出利用地微动观测结果对大地电磁数据进行约束反演的联合使用方法。随着技术的快速发展，利用多种物探方法相互结合来查明地热地质条件成为通用的做法，张宇等[14]在银川平原地热调查中利用浅层地温、视电阻率、卫星遥感图像、土壤汞质量分数、地热钻井等方法综合圈定了地热田范围；高亮等[15]在银川沙湖地区的地热研究中利用可控源音频大地电磁法结合岩性测深，确定了地热井位置，曹学刚等[16]利用MT方法辅助以钻孔资料，探讨了银川平原黄河东岸地热资源的分布特征及赋存规律；左丽琼等[17]针对地质条件复杂区的深部地热勘探，采用深孔测温测量、高精度重力剖面、微动探测、CSAMT法组成的综合物探方法，准确查明深部地热的分布特征。可以看出，针对不同地区地质条件与地热成藏模式的差异化与复杂化，物探方法的选取也需要具有很强的针对性。所以，在地热勘查与研究的实际工作中，以地质—地球物理特征为基础，合理选取一种地球物理方法的优化组合，不仅是一个重要的技术课题，而且对实际的地热勘查开发具有很强的指导意义。

本次在前期地质条件分析的基础上，利用1∶5万重力测量、可控源大地电磁测深、微动测量，反演了基底构造、确定了储盖组合、细化了盖层特征，综合分析了银川盆地东缘临河地区“传导—对流”混合型地热的成藏条件，并以此为指导，进行了地热资源远景区预测，对陶乐—横山堡冲断带的中北段地热资源勘查开发提供了先导性研究。

1 区域地质地球物理概况

1.1 区域地质概况

研究区位于银川盆地东部，大地构造位置属柴达木—华北板块Ⅰ级构造单元、华北陆块Ⅱ级构造单元、鄂尔多斯地块Ⅲ级构造单元、鄂尔多斯西缘中元古代—早古生代裂陷带Ⅳ级构造单元、陶乐—横山堡陆缘褶断带Ⅴ级构造单元[18]。黄河断裂是银川断陷盆地与陶乐—横山堡陆缘褶断带的分界构造，为一条超壳断层[19-23]。区域巨厚新生界覆盖范围广，前新生代地层出露有限，仅在灵武东山见有白垩系地层出露(图1)。

1—第四系；2—新近系；3—古近系；4—白垩系；5—侏罗系；6—三叠系；7—二叠系；8—石炭系；9—泥盆系；10—奥陶系；11—寒武系；12—元古宇；13—英云闪长岩；14—黑云母花岗岩；15—地质界线；16—推测地质界线；17—不整合界线；18—岩相界线；19—重要断层；20—隐伏断层；21—裸露断层

1.2 研究区重磁场特征

受西侧黄河断裂的控制，重力场呈近SN向展布且具有一定宽度的带状高值异常区，反映了深部奥陶系基底隆升幅度较大，与西侧的灵武凹陷低重力异常区形成了明显的区别(图2)。航磁场呈现出与重力场截然不同的特征，灵武以北区域为平缓展布的低磁异常区，灵武以南区域分布两处形态不完整的高磁异常，幅值较低，推测其为银川地堑南部受青藏高原东北缘NE向挤压应力的作用为形成逆掩隆升条带的反映，在深部的磁性地质体整体抬升的背景下，局部隆起、聚集所引起的航磁异常(图3)。

1—重力高异常；2—重力低异常；3—断裂

1—航磁高异常；2—航磁低异常；3—断裂

2 典型地热田地热异常

据研究区内地热钻孔揭示，天山海世界地热田为银川盆地东缘典型的地热异常区，其中DRT-03钻孔井底(1 690 m)温度为64.03 ℃，DRT-05钻孔井底(1 690 m)温度为62.87 ℃，地层温度明显高于银川盆地内部ZN钻孔的地温梯度(图4)。钻孔揭露的地层岩性资料显示，纵向上，以DRT-03钻孔为例，地温梯度曲线明显分为3段，代表了3种岩性不同的地层。浅部(200～600 m)新生界古近系—新近系(E-N)以砂泥岩互层为主的地层，地温梯度由表层的1.6 ℃/100 m深逐渐增加至浅层的2.42 ℃/100 m；中部(600～1 000 m)晚古生界石炭系—二叠系(C-P)含煤的地层，地温梯度随着深度增加陡然升高至5.04 ℃/100 m，而后快速回落至3.53 ℃/100 m；深部(1 000～1 960 m)晚古生界奥陶系(O)灰岩夹泥岩隔层的地层地温梯度逐渐趋于平稳，至(1.95～2.05) ℃/100 m区间(图5)。

图4 银川盆地及东缘地区地温对比

图5 DRT-03钻孔地温梯度

由此验证本区奥陶系(O)为深部热储层，石炭系—二叠系(C-P)为中部第一盖层(隔热层)，古近系—新近系(E-N)为浅部第二盖层(保温层)。

3 综合地球物理方法的应用

3.1 重磁异常特征分析

由重磁对应分析的结果可知，天山海世界地热井DRT-03很好地落在了重力异常小波二阶细节的局部高重异常上，且位于黄河断裂附近，反映出该区域为奥陶系(O)热储层的隆升区(图6a)。对于磁异常，已见地热井位于局部高磁异常附近，推测本区深部磁性基底埋深相对较浅，为区域性深部热源的有利区(图6b)。

3.2 奥陶系构造反演

奥陶系是本区地热赋存的重要地层，反演其构造形态，是研究地热赋存的基础。

研究区各地层间由浅至深存在2个密度界面，分别为新近系—古近系与石炭系—二叠系之间密度界面及石炭系—二叠系与奥陶系之间密度界面(第2密度界面(Δσ2,3))。具体地，新生界地层密度整体偏低，均值为1.91 g/cm3，与下伏的石炭—二叠系地层形成了密度差为0.68 g/cm3的第1密度界面(Δσ1,2)，深部奥陶系密度明显增大，均值为2.71 g/cm3，与上覆的石炭系—二叠系地层形成了密度差为0.12 g/cm3的第2密度界面(Δσ2,3)，上述两个密度界面直接影响了本区重力异常分布，其中第2密度界面对奥陶系基底顶面起伏形态起着决定性作用，为奥陶系基底顶面构造反演提供了密度物性前提。

基于1∶5万重力资料，利用平面帕克法密度界面反演银川断陷盆地南部奥陶系顶面构造。已知DRT-03、LS01与Ren1这3口钻孔钻遇的奥陶系顶面深度分别为791、800、788 m，与实钻深度相对比，反演深度偏差小于50 m，二者吻合，侧面印证了奥陶系储热层顶面反演结果的可靠性(表1)。

由图7可知，研究区内奥陶系隆升区域呈明显的长条状分布，南部灵武东部基本呈单一条带SN向展布，钻孔LS01位于该条带中部西侧斜坡部位，隆升幅度比较高，最高区域埋深为247 m；延伸至临河镇西南处分为两个NNE向隆升带，相较于南部隆起条带，幅度有所下降；到钻孔DRT-03西北处，两条隆起条带合二为一，并再次隆升明显，通贵乡东侧达到最高区域，埋深为485 m。上述奥陶系的隆升条带平面展布严格受黄河主断裂及其次级附属断裂的控制，是寻找天山海世界“传导—对流”混合型地热的潜力区(图7)。

图7 帕克法密度界面反演古生界基底深度

3.3 储盖层组合划分

储层为地热赋存提供了优良的空间，盖层阻断了储层中地热的逸散。

基于可控源大地电磁测量剖面WL-01剖面(CSAMT)对储、盖层叠置关系与特征进行了分析。纵向上分为4个明显的电阻率层，依次为高阻层(大于100 Ω·m)、中高阻层(25～100 Ω·m)、中低阻层(10～25 Ω·m)与低阻层(1～10 Ω·m)。高阻层无钻孔钻遇，根据区域地质特征推测为贺兰山岩群的变质岩系。经与钻孔DRT-03对比：中高阻层与储热层奥陶系相对应，为一套陆表海相沉积的含泥质灰岩，岩层裂隙发育，裂隙率为2.70%～20%，渗透率为0.01～115.4 μm2，测井温度数据显示，该段井温为38～64 ℃；中低阻层与隔热盖层石炭—二叠系对应，是泥岩、泥质粉砂岩、砂岩互层夹煤层沉积；低阻层对应保温盖层古近系—新近系的红色亚砂土、亚粘土沉积层，泥质成分含量高，此类岩石热导率较低，均小于2.3 W/(m·℃)，煤系地层热导率最低为0.61 W/(m·℃)[23](图8)。

3.4 盖层特征细化

盖层的分布特征决定了其保温效率的高低，是地热成藏的关键因素之一。

分析过DRT-03井区微动测量剖面，与钻孔分层数据对比，存在3个明显的S波速度层，深部(450～800 m)的弱低速层(1 600～2 400 m/s)为第一盖层石炭—二叠系，P波速度相对较高，层内横向波速稍有变化，推测是同层系岩性差异所致；浅部(180～450 m)的低速层(1 000～1 600 m/s)为第二盖层下段的古近系，层内速度变化不大，存在明显的局部低速异常区带，推测是砂岩与泥岩互层的表征；表层(180 m以浅)为极低速区(<1 000 m/s)，地层平坦，表明在古近系上段为胶结程度比较低的松软沉积层。综上，DRT-03井区第一、第二盖层发育完整，厚度稳定，横向变化较小，是良好的地热盖层(图9)。

图9 DRT-03井区微动反演S波速度剖面地质解释

综上所述，综合运用1∶5万重力、可控源大地电磁测深、微动测量资料，分析银川盆地东缘天山海地热田“传导—对流”混合型地热的成藏地质条件，效果良好。

4 地热资源远景区预测

以综合分析的地热成藏地质条件为指导，在1∶5万重力数据处理分析的基础上，圈定本区地热资源远景区(图10)。圈定原则有3条:一是为奥陶系顶面隆升区带，二是上覆一定厚度的石炭—二叠系与新近—古近系地层，三是处于自然保护区红线范围之外。

(a)重力小波2阶逼近场 (b)前第四系地层分布 (c)地热资源远景预测区

整体上，由重力小波2阶逼近场可知(图10a)，奥陶系储层的分布严格受黄河断裂及其次级附属断裂控制，分为灵武市—黄草坡一线的南部及横山堡—通贵乡一线的北部。南部地区基本全部位于白芨滩自然保护区，并且大面积出露白垩系宜君组与古近系清水营组，缺少优质的保温盖层(图10b)，因此不具备良好的储盖组合；北部地区覆盖有第四系全新统地层，其下伏的保温盖层(新近系—古近系)、隔热盖层(二叠系—石炭系)与储热层(奥陶系)地层完整，邻近黄河断裂，具备导热、导水通道，是寻找天山海世界“传导—对流”复合型地热的远景区。

以重力异常反演的奥陶系基底隆起区为底图，叠合白芨滩自然保护区与前第四系地层出露区，划定了本区地热远景的预测(图10c)，具体分为3个远景区(Ⅰ区、Ⅱ区与Ⅲ区)。Ⅰ区呈NNE向展布的不规则片状，面积为68.51 km2，地表被第四系完全覆盖，奥陶系储热层顶界面埋深最浅处位于钻孔DRT-03西北处4.6 km附近，埋深约430 m；Ⅱ区呈NNE展布的条带状分布于临河镇南北，面积为34.6 km2，地表覆盖第四系，奥陶系储热层埋深最浅处位于临河镇正N向2.9 km附近，埋深为1 120 m；Ⅲ区呈SN向展布的片状分布于横山堡北部，面积为44.23 km2，地表覆盖第四系，奥陶系储热层埋深最浅约为860 m，与Ⅰ区、Ⅱ区相比较，Ⅲ区远离黄河断裂，是否有良好的导热通道为该远景区提供深部热源需要进一步研究，建议布设MT电法剖面，分析低阻带、低速带，分析该区导热导水通道，明确热储模式，并布设地热钻孔验证物探方法的有效性。

5 结论

1)针对银川盆地东缘天山海世界“传导—对流”混合型地热资源的研究，1∶5万重力、可控源大地电磁测深及微动测量剖面资料是有效的资料基础，DRT-03井的成功钻遇验证了物探资料应用效果。

2)研究区地层纵向存在2个差异明显的密度界面，深部奥陶系与中深部的石炭—二叠系形成了密度差为0.12 g/cm3的第2密度界面(Δσ2,3)，为奥陶系基底顶面构造反演提供了密度物性前提。利用平面帕克法反演的奥陶系储热层顶面深度与DRT-03、LS01与Ren1这3口钻孔钻遇深度相吻合。反演结果显示出奥陶系隆升区域呈明显的长条状分布，平面展布严格受黄河主断裂及其次级附属断裂的控制。

3)本区域电性与波速的差异性为储盖组合划分提供了依据。结合WL-01剖面可控源大地电磁测深和微动测量的结果与地热钻孔DRT-03相对比，可知深部奥陶系储热层呈现中高阻层(25～100 Ω·m)的特征，中部石炭—二叠系作为隔热盖层表现出中低阻层(10～25 Ω·m)、弱低速层(1 600～2 400 m/s)的响应，浅部的低阻层(1～10 Ω·m)、低速层(1 000～1 600 m/s)则是保温盖层古近系—新近系的反映。

4)在银川盆地东缘地热储盖层组合的地球物理模型总结的基础上，预测了3处“传导—对流”混合型地热资源的远景区，为该区地热资源规模性开发利用提供了科学依据。