地热资源地震勘探方法综述

  全球能源问题日益突出。一是供需矛盾，截止２０１０年，２０％的世界人口仍然无法用电。二是能源 需求增长与温室效应的矛盾，２０１０年初级能源需求增加５％，与之相关的二氧化碳排放量为３０．４ｔ，相比２００９年增长了５．３％，可能引起严重的温室效应。三是资源浪费，由化石燃料补贴而引起严重的浪费性消费，２０１０年增加约４０９０亿美元［１］。四是不稳定因素，能源分布不均，中东和北非等国家的动荡易引起全球性能源危机。 而地热资源作为可再生无污染资源，将对发展低碳经济起到极其重要的作用。《新能源报告》指出，地热资源的能量是全球常规能源总能量的４５万倍，是全球煤资源的１．７亿倍。此外，我国地热资源较为丰富，中国主要盆地地热资源量为２．４９６×１０２２Ｊ，折合标准煤８．５３２×１０１１ｔ，并且中国仅重点城市的浅层地热能就达到２．７８×１０２０Ｊ，每年可利用量为２．８９×１０１２ｋＷｈ［２］，因此我国地热资源的开发和利用 前景广阔。 地热资源勘探的地震方法分为主动地震和被动地震两类方法。被动地震勘探是利用密集地震台阵来监测岩层中微地震活动所产生的地震波，然后运 用地震学方法来反演微地震的活动特征，并反演出研究区三维纵横波速度、纵横波速度比、泊松比等弹性参数分布的一种勘探方法。被动地震勘探方法在国外地热资源勘探中已得到普遍应用，而国内对该方法研究尚不广泛。主动地震是指高分辨率的反射地震勘探方法。主动地震成本相对高，因此在地热资源普查阶段仍然依靠被动地震，而对于评价区的地热储层，主动地震才是最可靠的资料。 地震方法作为一种超深且高精度的勘探方法，在一定程度上弥补了重磁电方法勘探的不足。它已在我国大多数地热田展开了应用，例如西藏、天津、河北、云南等省市地热田都采用了三维地震勘探技术［３－４］，并且山东省地热田还采用高分辨率地震勘探技术［５］。但是，国内地热田仅仅利用主动地震方法来研究地质构造，解决断裂识别的问题，而如何有效利用地震信息解决热储层预测、裂隙通道精细刻画、地层水识别等关键问题尚是空白。另一方面，被动地震方法在国内地热田中的应用严重不足。也就是说，与美国、日本相比，我国地震技术在地热资源勘探中利用得远远不够，发展速度也相对滞后。在美国和冰岛等国家，地震方法早已成为地热资源勘探开发的主流手段，不仅利用叠后地震资料进行导水的断裂的识别，还综合利用叠前、叠后以及多波信息精细刻画和描述热储层特征，明显提高了地热勘探的精准度。

 

  １ 地热储层地震响应特征

 

  １．１ 地震波速度衰减 随着地层温度的升高，岩石的物性特征会发生 很大变化，其中最明显的变化是：随着温度的升高，岩石的纵波速度降低，横波速度也逐渐减小且趋于零［６］。Ｍａｋｋｙ等［７］通过岩石物理实验进一步印证了 这一结论：岩石的地震波速度随温度升高而降低，水的地震波速度也随温度升高而降低。此外，由于地热储层中蒸汽和裂缝的存在，将导致弹性波速度进一步衰减。探开发的主流手段，不仅利用叠后地震资料进行导水的断裂的识别，还综合利用叠前、叠后以及多波信息精细刻画和描述热储层特征，明显提高了地热勘探的精准度。

 

  １ 地热储层地震响应特征

 

  １．１ 地震波速度衰减 随着地层温度的升高，岩石的物性特征会发生 很大变化，其中最明显的变化是：随着温度的升高，岩石的纵波速度降低，横波速度也逐渐减小且趋于零［６］。Ｍａｋｋｙ等［７］通过岩石物理实验进一步印证了 这一结论：岩石的地震波速度随温度升高而降低（图１ａ），水的地震波速度也随温度升高而降低。此外，由于地热储层中蒸汽和裂缝的存在，将导致弹性波速度进一步衰减。

 

  １．２ 地层吸收衰减 当地层含流体后，在地震资料上振幅和频率都有所反映：纵横波振幅能量严重衰减，地震波高频成分也大幅度衰减，且含蒸汽与含水时也存在一定的差别，含蒸汽时振幅、频率衰减更严重，有时在地震剖面上形成“空白带”反射，因此地层吸收系数和地震衰减属性是进行地热储层描述和流体预测的重要参数。

 

  １．３ 弹性参数的影响 非地热储层、含水地层以及含蒸汽地层的纵横波速度比、泊松比等弹性参数有较大差异。在国外很多地热田中，都对纵横波速度比做过详细研究。总结这些研究可知，水占主导的地热储层（例如美 国ＥａｓｔＭｅｓａ和墨西哥ＣｅｒｒｏＰｒｉｅｔｏ地热田）纵横波速度比高，达到１．５５～１．６８。而气态水占主导的地热田（例如美国Ｇｅｙｓｅｒｓ和Ｃｏｓｏ地热田）纵横波速比低。因此可以根据纵横波速度比识别地热储层，并进一步确定储层中的流体性质。

 

  １．４ 裂缝发育特征 断裂构造对地热资源分布有直接的控制作用。上地幔岩浆热液沿深大断裂上升至地壳上部，将浅层的地下水加热，这些被加热的地下水沿着断裂通道上升。断裂和裂缝为地下高温水、蒸汽和热量提供通道，因此地热储层中断裂和裂缝通道十分关键。据统计，具有价值的热储一般主要位于褶皱部位以及断裂构造发育部位，因此受断裂和裂缝影响，热储类型主要分为：孔隙裂缝型热储、岩溶裂缝型热储、构造裂缝型热储等。可以说，地热勘探的另一个重点就是寻找裂缝发育带。罗马尼亚ＭｏｅｓｉａｎＰｌａｔ?ｆｏｒｍ地热区的地热勘探中，在Ｃｒａｉｏｖａ?Ｂａｌｓ?Ｓｌａｔｉｎａ、Ｖｉｄｅｌｅ和Ｌｅｏｒｄｅｎｉ?Ｒｅｃｅａ三个地区的断层附近都找到了８０℃的地层水，证实了断层、裂缝与地热储层的紧密关系［８］。

 

  ２ 被动地震勘探技术

 

  ２．１ 速度成像 由于地热异常区内蒸汽和热水的上涌，时常导致地面震动以及局部的蒸汽爆炸。因此利用高精度微震仪可探测这些微震，从而圈定热储位置、上涌通道以及分布范围等［９］。被动地震方法成本低，在地热田的普查阶段发挥着重要作用。 基于被动地震数据的旅行时反演在地热勘探中 应用较多［１０］。该方法是利用地震波旅行时数据，反演地下地质体的速度模型，找出速度异常体。Ｌａｒ?ｄｅｒｅｌｌｏ－Ｔｒａｖａｌｅ地热田就是采用旅行时反演对速度进行成像从而来圈定热储位置。该地热田以蒸汽为 主，是第一个用于发电的地热田。Ｍａｔｔｅｉｓ等［１１］在该项目中，采用ｖｐ、ｖｓ、ｖｐ／ｖｓ、ｖｐ×ｖｓ四种属性来研究储层内流体分布特征，不仅识别出地热储层，还指出了热储的流体类型。由于低ｖｐ／ｖｓ异常一般对应含蒸汽地热储层，所以图３剖面中Ｌａｇｏ井和Ｍｉｎｉｅｒａ井底部均为含蒸汽地层。而高ｖｐ／ｖｓ异常对应（并且对应低ｖｐ×ｖｓ）含水地层，因此剖面中浅层为含水层。

 

  ２．２ 随机噪声成像 随机噪声成像方法是利用从随机地震噪声中提取短周期面波，由这些数据可以得到面波的频散曲线，并利用频散曲线建立横波速度剖面，从而寻找速度异常的地热储层［１２］。该方法在探测下地壳和上地幔之间的构造中应用较为广泛［１３－１４］。目前基于随机噪声成像的研究在新西兰、南非、西班牙、韩国、 日本、美国等多个国家和地区都有所应用［１４］。随机噪声主要指来源于天气变化、大气压力、海 浪、潮汐、人类活动和机器产生的噪声。这些在地表探测到的微振动分为体波和面波，而７０％的是面波［１５］。尽管这些随机噪声的振幅和波形随时间和地点不断变化，但是统计结果表明：一般情况下在某一特定时间和地点，随机噪声的振幅和波形相对比较稳定。

 

  图４给出了该方法的简要原理及步骤：①提取随机噪声垂直分量：从多组台站记录到的持续时间较长的随机噪声中提取出垂直分量数据，该数据中主要包含瑞利面波信息；②噪声压制：根据面波频带范围，滤除非面波信号，并采用谱白化、振幅归一化等处理进一步压制无效信息；③计算互相关：在面波的频带范围内，对处理后不同台站的随机噪声计算互相关［１６－１７］。④计算频散曲线：采用时频分析方法 （ＦＴＡＮ）由互相关结果计算得到瑞利面波的相速度频散曲线。⑤面波成像：利用计算得到的频散曲线进行面波成像，得到瑞利面波相速度模型。很多学者将面波成像方法单独地进行过详细介绍［１４］。⑥横波反演：通过使得实际计算的频散曲线和频散模型曲线差异最小，得到横波速度模型的最优解［１７］。 最后，根据面波速度和横波速度结果，圈定出热储位置和范围。Ｙａｎｇ等［１８］采用地震随机噪声成像方法研究了 美国加州Ｃｏｓｏ地热田的地质构造。图５是Ｃｏｓｏ地热田不同周期的相速度的扰动平面。不同周期的相速度扰动能反映不同深度的热储特征。在这项研究中，周期为３ｓ的瑞利面波能较好地反映２ｋｍ深的热储特征，而１０ｓ的瑞利面波对１０ｋｍ深的热储反映尤为明显，对２０ｋｍ的也能较好反映。小的相速度扰动对应于热储区域。

 

  国内部分地热田也采用了该方法，并取得了预期效果。Ｘｕ［１２］通过随机噪声成像研究了江苏吴江的地热勘探项目。通过随机噪声成像最终得到横波剖面（图６）。从横波剖面上看，浅层速度横向上连续，说明该地层沉积环境相对稳定；深部存在低速度异常，说明存在断层和裂缝发育带以及含水层或石灰岩洞穴等地质体。此外，也得到了地热储层深度、厚度等重要信息。

 

  此外，国内外通常采用多种地球物理方法综合预测，从而提高预测的准确度和精度。美国犹他州西南的ＣｏｖｅＦｏｒｔ地热田项目同时采用了面波数据和局部天然地震数据［２０］。宽频带的随机噪声数据由美国国家台阵记录得到，利用该数据得到相速度频散曲线，从而用于面波成像；局部地下活动的微地震数据则由当地１０个台站记录得到，利用双精度的层析成像方法得到纵波和横波速度模型。加利福利亚北部的Ｇｅｙｓｅｒｓ浅层地热储层探测也采用了微地震数据和面波数据的联合反演方法［２１］。此外，综合面波、局部微地震数据以及重磁电数据综合分析得到横波速度体的方法也是非常常用且有效的方法［２２］，多种方法能相互印证，相互补充，使预测结果更为合理。

 

  ２．３横波分裂技术

 

  横波分裂技术一致被认为是探测地热储层裂缝密度和方向行之有效的方法［２３］。由于从微地震事件中同样也能接收到快、慢横波信息，因此完全可以将横波分裂技术应用到被动地震之中。

 

  横波传播经过裂缝时分裂为快横波和慢横波，快横波平行于裂缝发育方向，慢横波垂直于裂缝方向，当裂缝密度很大时，导致慢横波存在较大延迟。裂缝密度越大，快慢横波时间差异越大。因此，快慢横波的时间延迟可以用于探测裂缝发育密度，而快横波可用于探测裂缝发育方向［２４］。利用横波分裂预测裂缝的理论和方法国外很多学者已经做过研究。

 

  Ｓｉｌａｓ等［２９］在ＣａｓａＤｉａｂｌｏ地热田勘探中采用的就是基于微地震数据的横波分裂预测裂缝方法。首先，通过研究快横波和慢横波的振动方向，得到裂缝发育方向；然后用快慢横波的时间差计算得到裂缝发育密度。图７是ＣａｓａＤｉａｂｌｏ区横波分裂预测的裂缝发育密度深度切片。红色表示高裂缝密度发育。从图中可以得出，该区有两个有利的地热储层。

 

  一个是ＣａｓａＤｉａｂｌｏ区１～３ｋｍ深度处的北边和东部，一个是２．５ｋｍ深度的远东部（Ｍａｍｍｏｔｈ机场附近）的。 被动地震成本低，能有效预测热储位置和分布范围，但是由于被动地震数据接收的站点少，数据覆盖范围小等不足，导致预测结果的纵横向分辨率都较低，所以基于被动地震的勘探方法一般局限于地热田的普查阶段，而选区评价和井位优选时必须借助主动地震３ 主动地震勘探 利用地表上有限个台站数据进行分析，难以实现对大面积地热区进行高精度评价研究。因此，针对较大区域的地热区，需要借助三维主动地震勘探方法。主动地震勘探覆盖面积广，信息量大，能够提供更丰富的地下信息。

 

  ３．１ 地震衰减属性 在地热储层中，由于热液、蒸汽等高温流体存在，导致地震波振幅、频率吸收衰减严重。因此很多学者尝试建立地热田的三维吸收衰减系数模型，通过吸收衰减系数模型找到热储［３０－３１］。Ｗｕ等［３２］通过计算纵波的三维吸收衰减系数Ｑ模型来研究加州Ｃｏｓｏ地热储层，而Ｓａｎｄｅｒｓ等［３３］则采用横波衰减属性研究该区热储。此外，频率属性也是预测地热资源的有效手段。 Ｗｅｉ等［３４］在加利福利亚南部Ｗｉｓｔｅｒ地热田采 用了有效频带宽度属性预测高温流体分布。由于地热田高温流体的存在，导致频率衰减，地震有效带宽变窄，故有效带宽的低值区指示高温流体分布范围。因此，由图８预测出高温流体主要分布在工区西南部。 图８ Ｗｉｓｔｅｒ地热田纵波有效频带宽度平面示意（红色为 高值，紫色为低值） ［３４］３．２ 构造精细成像 地热资源主要受地质构造因素控制，与火山结构或断层活动高度相关联，这些因素导致地热资源分布极其不均匀［３５］。另外，地质构造位置和大小共同决定了地热储集体的大小［３６］。因此，地震资料的构造成像效果对地热储层预测尤为关键。 在意大利Ｔｕｓｃａｎｙ地热田勘探中，尽管常规处理结果（图９ａ）大部分资料能够满足勘探要求，但是在石灰岩地层中由于岩石的吸收衰减降低了地热储层的成像质量，故信噪比急剧降低，所以Ｃａｓｉｎｉ［３６］在前期处理和成像过程中，改进了资料品质，提高了成像精度。如图９ｂ为改进后的地震资料，Ｍ＿１井目标层段的信噪比和层位精度得到明显提高。

 

  ３．３ 三维多波多分量地震勘探 基于被动地震数据的横波分裂技术对预测地热储集体，特别是裂缝体系具有很好的效果。为了更好地服务于地热能源的开发与建设，提高勘探精度，近年来一些国家开展了地热区多波三维地震勘探。受美国能源部委托，德州大学奥斯汀分校勘探地球 物理实验室２００９年开始开展多波地热储集体的表征研究。现在该实验室已采集到了Ｓｏｄａ和Ｗｉｓｔｅｒ两块地热区多波地震资料，开展了一些研究工作。三维多波多分量的优势主要体现在以下两个方面。一是，由于ＰＰ波穿过热流体衰减严重，降低资料信噪比，而ＰＳ波受热流体影响小，因此ＰＳ波的应用会显著提高地热区的地震解释精准度。二是，利用ＰＰ波和ＰＳ波计算的弹性参数（例如ｖｐ／ｖｓ）更加可靠，对热流体分布指示更准确；此外，利用横波分裂技术能有效识别储层裂缝。总之，多波多分量三维地震勘探比基于被动地震数据的多分量勘探精度和分辨率更高，它将是地热地震勘探的重要发展方向。

 

  图１０和图１１是Ｗｉｓｔｅｒ地热田三维多波多分量地震资料应用实例。图１０是Ｗｉｓｔｅｒ地热田Ｃａｎｅ? ｂｒａｋｅ／Ｏｌｌａ／Ｄｉａｂｌｏ层段与Ｄｅｇｕｙｎｏｓ层段地震解释剖面，图１０ａ为ＰＰ波，图１０ｂ为ＰＳ波。对比两者虚线椭圆框内不难发现，ＰＰ波受浅层热流体影响，反射同相轴被削弱，资料信噪比急剧降低，而ＰＳ波受热流体影响小，虚线框内资料信噪比仍然较高，更有利于地震解释工作的开展。另外，由ＰＰ波和ＰＳ波计算得到的ｖｐ／ｖｓ数据体比常规反演得到的ｖｐ／ｖｓ数据体可靠性更高。图１１是利用ＰＰ波和ＰＳ波计算得到的Ｃａｎｅｂｒａｋｅ／Ｏｌｌａ／Ｄｉａｂｌｏ层与Ｄｅｇｕｙｎｏｓ层之间的ｖｐ／ｖｓ平面分布，由图可知Ｗｉｓｔｅｒ地热田西南部ｖｐ／ｖｓ低于２．６（红色区域），因此该区域为最有利的地热储层。

 

  ４ 结论及认识 （１）在未来的地热资源勘探开发应用中，地震方法优势巨大。 （２）地热勘探最主要的目标是寻找裂缝和断裂 发育区，而裂缝和断裂能够造成地震波弹性参数明显的变化，地震方法是捕捉这些变化最有效的手段。 因此，使用地震进行地热勘探将会明显提高勘探精度。 （３）被动地震方法理论成熟，成本较低，但是分 辨率低，适合地热田的初期开发；主动地震方法密度大，高信噪比和高分辨率，适合地热田的选区评价和井位优选；将二者有效结合，能够提高地热资源的探测精准度。 （４）地震多波技术，对描述裂缝和断裂特征有 显著效果，是地热勘探未来发展的重要方向。