1.   引　言

现今，世界能源危机和环境问题日益严峻，发展储藏量大且环境友好的新型替代能源逐渐受到了各国政府和科学家的重视。可以有效开发的水热型地热资源仅占已探明地热资源中的10％左右，更多的地热能储存于干热岩型地热资源中（汪集旸等, 2012; 许天福等, 2018）。作为干热岩型地热资源的有效利用手段，增强型地热系统（Enhanced Geothermal System, EGS）的开发实施需要对低孔、低渗的岩石进行改造（多吉和王贵玲, 2022）。作为中国首例EGS工程，在青海省共和县恰卜恰镇进行的EGS工程需要针对花岗岩型干热岩储层进行人工造储，而花岗岩储层内部裂缝的分布规律在场地评价中扮演着尤为重要的角色（姜晓宇等, 2020; Johnson et al., 2021）。干热岩储层中的裂隙从形成原因上可分为天然裂隙和诱发裂隙两种类型，干热岩储层改造阶段所形成的裂隙系统即为诱发裂隙，岩石物理实验表明压裂改造形成的诱发裂隙是在天然裂隙的基础上，受现今应力特征所控制，复杂缝网的形成与天然裂隙发育与否有较大的关联，施工中应在天然裂隙较丰富的层段开展水力压裂，便于实现复杂裂隙网络（周舟等, 2019; 王璜等, 2020; 刘德民等, 2021）。天然裂隙从发育规模上可以分为大尺度断裂带和小尺度裂缝，EGS工程井位部署应选择小尺度裂缝发育的区域，同时还要尽可能避开大尺度断裂带，防止导致断层活化，诱发地震（张森琦等, 2019）。花岗岩型干热岩储层内部裂缝系统分布规律复杂，裂缝的研究应当贯穿于干热岩开发的整个过程中（张雨晴等, 2016; Johnson et al., 2021; 解经宇等, 2022）。因此，如何准确刻画花岗岩储层中的裂缝系统成为EGS场地优选、压裂作业实施以及开发钻孔合理部署的关键（汪集旸等, 2012; 刘国平等, 2016; 李振宇等, 2022; 谢紫霄等, 2022）。花岗岩储层裂隙发育影响人工压裂造储裂缝延伸，是影响EGS工程规模化开发的主要因素，也是一项具有挑战性的难题（张盛生等, 2019）。

为推进深部热能（干热岩）开发利用产业化，中国地质调查局组织实施了“干热岩资源调查与勘查试采示范工程”，该工程在青海共和盆地实现了干热岩试验性发电（Zhang et al., 2022），期间应用三维地震属性融合技术成功地对花岗岩型干热岩裂缝进行了评价，本文分析其应用效果，以期为花岗岩型干热岩储层裂缝探测提供参考。

2.   工作区概况

共和盆地处于昆仑山脉与秦岭山脉之间，是古近系初形成的呈北西西方向展布的菱形断陷盆地，该盆地属于东昆仑、西秦岭造山带接合部位。按地质力学划分，属昆仑—秦岭纬向构造带与河西系构造复合的沉降带（图1，张森琦等, 2020）。盆地周边均受山前深大断裂的控制，北侧是青海南山、拉脊山断褶隆起带，南侧是河卡南山、巴吉山断褶隆起带，西为鄂拉山构造岩浆岩带，东为扎马山断褶隆起带，中间被瓦里贡山构造岩浆岩隆起带分隔出共和、贵德两盆地（李振宇等, 2022）。挽近期构造活动强烈，山体隆升，盆地大幅度下降，沉积了较厚的第四系冲洪积、河湖相堆积物和新近系、古近系湖相堆积物（图1，张森琦等, 2018）。工作区内三叠系广泛出露在盆地周边隆起带和组成盆地基底，中—下三叠统为海相沉积，上三叠统为陆相沉积。盆地内侵入岩以花岗岩、花岗闪长岩、石英闪长岩和斑状花岗岩等为主，侵入于中—下三叠统，与三叠系共同构成褶皱基底（张森琦等, 2019）。盆地内发育三级湖相阶地和多级黄河阶地，说明随青藏高原隆升，盆地一直处于间歇性抬升状态（唐显春等, 2020）。共和盆地内具有局部高大地热流及高地温特征，在盆缘边界断裂带有多处温泉出露，已在恰卜恰及扎仓寺等地多口钻井钻遇盆地基底印支期花岗岩类干热岩。盆地内主要的地热储层岩性为印支期—燕山期花岗闪长岩、二长花岗岩等（Gao et al., 2018）。

图  1  共和盆地及周边地质概况图（据张森琦等, 2018）

1—全新统；2—上更新统；3—中—下更新统；4—上新统；5—上三叠统；6—中—下三叠统；7—印支期花岗岩体；8—地层界线；9—断裂；10—地热钻孔；11—油井；12—热泉

Figure  1.  Geological map of the Gonghe Basin and surrounding areas (after Zhang Senqi et al., 2018)

1–Holocene; 2–Upper Pleistocene; 3–Middle–Lower Pleistocene; 4–Pliocene; 5–Upper Triassic; 6–Middle–Lower Triassic; 7–Indosinian granitic rocks; 8–Stratigraphic boundary; 9–Fault; 10–Geothermal drilling hole; 11–Oil well; 12–Hot spring

下载: 全尺寸图片 幻灯片

研究区位于共和盆地东北部，基础地质研究认为区内快速抬升减压造成的深部部分熔融体持续为干热岩体供热，盆地深大断裂沟通热储，形成局部对流型热循环通道，巨厚的放射性花岗岩体起到增温和传热导热作用，上部覆盖约2000 m厚的低热导率古近系以来的沉积盖层起到保温聚热作用（张森琦等, 2021）。

3.   研究思路

青海共和干热岩试采场地在选址阶段开展了重力、磁法、大地电磁和深反射地震剖面探测等多种地球物理勘查手段，明确了共和盆地结构特征，优选出EGS工程试采场地，该场地位于盆地中西部，远离盆地边缘大断裂。其中地震02测线（图2）显示研究区周边多条断裂断穿花岗岩基底，断层性质以西倾逆断层为主。花岗岩顶面构造整体呈西低东高的构造形态。从沉积岩层的发育情况看，西部沉积厚度较大，东部厚度较薄，自西向东地层呈楔形。为进一步明确EGS试采场地范围内裂缝发育特征，项目组进一步开展了宽频带、宽角度、高密度采样的三维地震工作。三维地震工作区面积4 km²，覆盖次数达到238次以上，面元大小10 m×10 m，道距20 m，激发点距10 m或20 m，横纵比0.77；数据处理软件使用东方地球物理公司的GEOEAST。

图  2  研究区二维地震剖面结果

Figure  2.  Result of 2D seismic profile in the study area

下载: 全尺寸图片 幻灯片

由于花岗岩型干热岩储层位于近2000 m厚的沉积盖层下，花岗岩地层内幕杂乱地震波杂乱反射，因此提取其内部裂缝带发育规律难度较大。故本文首先使用构造导向滤波和多窗口倾角扫描技术进一步提升三维地震资料的信噪比（Nguyen et al., 2014; 范存辉, 2015; 谢舟等, 2017）。通过多方位角约束的地震属性技术对断裂分布展开预测，针对裂缝簇，叠前应用最大似然体进行属性预测，叠后采取相干体、曲率体、方差体等三维地震属性预测裂缝簇的空间分布密度（Khair et al., 2015）。最后，利用多属性融合技术实现了多种分析结果的对比，最终实现对裂缝发育带的预测。

4.   地震属性裂缝预测技术

4.1   多窗口倾角扫描

多窗口倾角扫描技术通过搜索各点不同的倾角窗口，得到地震数据的倾角数据体，局部反射界面的成像效果优化，提高了倾角和方位角的计算精度（Marfurt et al., 1998; 印兴耀等, 2014）。多窗口倾角扫描技术还可利用特定方位角进行倾角扫描，更好地识别特定方位走向的断裂，用于估算局部反射层的不连续性，有助于微小断层的解释（应玉双, 2016）。研究区内三维地震资料的有效扫描倾角为0~180°，以15°为间隔将三维地震数据体计算生成倾角数据体，由于研究区多发育北西向和北东向断裂，垂直于断裂走向进行倾角扫描能够更清楚地识别该断裂，最终通过对比发现，30°方位角数据体对断裂的反映更为清晰（图3）。

图  3  基于多倾角窗口扫描前后计算的相干时间切片（据应玉双, 2016）

Figure  3.  Coherent time slices computed before and after multi−window inclination scanning (after Ying Yushuang, 2016)

下载: 全尺寸图片 幻灯片

4.2   多方位角相干技术

多方位角相干技术利用地震波在不同方向传播时的差异性，通过方向滤波和相干性分析，提取具有特定方向和频率特征的地震信号，可以有效地提高地震数据的分辨率和信噪比（Bahorich and Farmer, 1995; 刘传虎, 2001）。将多窗口倾角扫描的系列结果进行相干处理，形成多方位角相干数据体。在地震解释基础上，对方位角相干数据体提取沿花岗岩顶界面相干切片，勘查结果显示预测的裂隙带展布与构造解释的类型发育带方向基本一致，研究区内大断裂并不发育（图4）。

图  4  多方位角相干体属性切片图

Figure  4.  Sliced image of multi−azimuth coherence volume attribute

下载: 全尺寸图片 幻灯片

4.3   构造导向滤波

由于花岗岩界面与上覆沉积岩有较大的波阻抗差异，使花岗岩顶面形成了一个地震波反射信号遮挡面，对地震反射信号有屏蔽作用，使得花岗岩内幕反射信号变弱，反射地震剖面受噪声影响较大，需要通过解释性处理，降低噪音影响，突出裂缝信息，以达到提高裂缝预测精度的目的。构造导向滤波是应用效果较好的预处理方法之一，它在分析同相轴的延伸方向及变化点（断点）的基础上，沿地震同相轴进行平滑，平滑操作不跨越断点，能够使断点间的同相轴变得更加连续，同时使断点显得更加突出，十分有利于裂缝预测。原始地震资料虽然能够反映一些小断裂或者裂缝（图5a)，但由于地震有效信号弱，但断点并不清晰。根据该区构造特征选择合理的参数进行构造导向滤波，得到如图5b所示的地震剖面，滤波后地震剖面断点变得更加清晰、干脆，有利于裂缝预测。

图  5  构造导向滤波前后对比图（a，原始地震数据体，b，构造导向滤波后数据体）

Figure  5.  Comparison of before and after structural-oriented filtering (a, original seismic data; b, seismic data after structural-oriented filtering)

下载: 全尺寸图片 幻灯片

4.4   锐化曲率体属性

地震曲率体属性通过计算层面的弯曲程度定量裂缝产生的可能性，作为一种常规手段在裂缝发育区预测应用较为广泛（Roberts, 2001; 王雷等, 2010）。本研究在计算曲率体的基础上，利用曲率锐化技术，将构造导向滤波后的地震数据体先计算传统曲率体，在传统曲率体的基础上利用锐化手段生成新的高精度锐化曲率属性体。通过曲率锐化不但能保持与传统曲率体预测的裂缝展布规律一致，还能更清晰地识别小断裂。从锐化前后的对比结果中可见（图6）：锐化后的属性切片对于大断裂及破碎带的整体平面展布是相似的，但是对于小尺度裂缝的预测更为丰富和清晰。裂缝发育带曲率值较大（红色），研究区发育有高曲率值的裂缝，裂缝与断裂带走向基本一致，沿北东和北西走向呈条带状分布。曲率体对花岗岩顶界面处裂缝发育带准确度较高，花岗岩体内部裂缝发育带需要结合其他属性才能够获得。

图  6  曲率体锐化前后对比图

Figure  6.  Comparison diagram of curvature body of before and after sharpening contrast

下载: 全尺寸图片 幻灯片

4.5   裂隙玫瑰图分析

裂缝玫瑰图分析方法常用来描述特定线状特征的方位分布，便于直观分析裂隙分布规律，通过在层切片上选取主测线和联络测线组成的网格化矩形分析时窗来产生玫瑰分布图。在每个时窗内，将每个样点根据其方位最小曲率映射到玫瑰图上的花瓣上，然后相加并按比例形成玫瑰图，最终成为了玫瑰图数据体。每个玫瑰分布值均位于原始数据分析时窗的中心位置，这样就形成了表达裂缝分布的属性体，显示裂缝分布的密度和方位。分析结果显示（图7），研究区内裂缝方向以北东向和北东向为主，与GH-01井钻井揭示的裂缝方向一致。攻瑰图属性预测的裂缝发育区范围与曲率体、蚂蚁追踪以及相干体属性具有相似的特征，在裂缝发育区玫瑰图具有明显的裂缝发育方向；在裂缝不发育区，攻瑰图没有明确的裂缝发育方向。

图  7  玫瑰图1500 ms和2300 ms时间切片图

Figure  7.  Rose diagrams of time slice images at 1500 ms and 2300 ms

下载: 全尺寸图片 幻灯片

4.6   属性融合技术

属性融合技术就是利用数学运算把多种属性结合在一起，实现在同一色板模式下，同时显示多种属性的功能。以多方位角相干体和曲率锐化体成果按等比例关系属性融合（图8），把单属性对于断裂带和裂缝型储层的多解性降低，提高研究区的断裂识别精度。

图  8  花岗岩顶面沿层曲率与相干体属性融合成果图

Figure  8.  Fusion result map of the top surface of granite along the layer curvature and coherence attribute

下载: 全尺寸图片 幻灯片

5.   讨　论

共和盆地花岗岩干热岩储层内通过利用曲率体、相干体属性等融合而成的三维地震融合属性，很好的预测了裂缝的发育特征。横向上从花岗岩顶面沿层曲率与相干体属性融合解译裂缝发育带结果来看（图9），研究区内花岗岩发育区内发育裂隙带，裂隙带分布方向以北东和北西向为主，与该区历史主应力方向基本对应。纵向上对比成像测井裂缝发育图（图10），成像测井结果揭示井周发育中高角度和网状缝构造缝，该井位于裂缝较为发育的区域。通过地震裂缝预测与钻井结果对比（图10），达到了干热岩裂缝探测的评价要求。三维地震属性体综合分析结果显示研究区内裂缝较为发育，裂缝发育带主要沿北东向展布，裂缝发育带以北东和北西向展布，深度越深裂缝越发育。预测裂缝发育带分布情况与钻井资料吻合。

图  9  综合解释裂缝发育带分布图

Figure  9.  Distribution map of fracture zone interpretation

下载: 全尺寸图片 幻灯片

图  10  深度2370~2375 m范围内三维地震解释裂缝发育程度与钻孔成像测井对比图

Figure  10.  Comparison between 3D seismic interpretation fracture and borehole imaging logging in the depth range of 2370–2375 m

下载: 全尺寸图片 幻灯片

6.   结　论

（1）共和盆地干热岩靶区内裂缝发育带分布规律通过属性融合后精度大幅度提高，降低了裂缝预测的多解性。

（2）不同属性分析技术的应用效果，曲率属性对花岗岩顶界面裂缝带判别效果较好，但对花岗岩内裂缝带发育情况在相干体属性显示结果优于曲率体，利用属性融合技术可以将两种属性分析结果优势互补，对花岗岩型干热岩储层裂缝带探测具有较好的效果。

（3）通过上述几种裂缝预测技术在花岗岩型干热岩裂缝探测中的应用可以发现，研究区内不发育大型断裂，小尺度裂缝带较为发育，发育方向以北东向和北西向为主。

（4）在研究区预测出的裂缝发育带对比常规测井和成像测井结果，裂缝带预测结果基本一致。