研究区干热岩勘探主要分布在太古界, 钻井取心(图2)及录井岩屑资料显示, 本区干热岩主要岩性为角闪岩和花岗岩两大类, 还有少量变质花岗闪长岩以及混杂变质岩。其中, 花岗岩以变质角闪二长花岗岩为主, 颜色主要为浅青灰色和灰白色, 成分以斜长石、钾长石为主, 石英次之, 角闪石少量, 角闪石呈黑绿色, 夹有乳白色钙质填充物; 角闪岩成分以角闪石为主, 斜长石次之或者极少; 变质花岗闪长岩颜色为浅灰至灰绿色, 成分以斜长石为主, 角闪石次之, 含少量石英; 花岗岩和角闪岩均有变晶结构, 个别岩石夹有白色钙质填充物, 推断为裂缝内部充填物。

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 X井岩心照片(角闪二长花岗岩)Fig.2 Core photos of Well X (hornblende monzogranite)

不同岩性、矿物成分的岩石物理特征值差异比较大, 所以在进行岩性识别时, 可根据不同岩性的主要岩石组成, 各类岩石的主要矿物组分, 再依据各矿物的岩石物理特征值的差异来进行测井特征综合分析, 从而建立岩性判别模式^{[16, 17]}。例如角闪岩类、和闪长岩类的岩石由于U、Th、K含量较低, 自然伽马数值表现为低值, 但当角闪岩中含有斜长石时, 会使自然伽马数值升高。角闪石的岩石光电吸收截面指数(photoelectric absorption cross-section index, PE)数值较高, 为5.99 b/e, 同为斜长石类矿物的钠长石和钙长石PE数值差异较大, 分别为1.7 b/e和3.1 b/e, 而闪长岩和花岗岩的PE数值较低, 小于2.8 b/e, 所以各类岩石的物理特征差异与矿物含量的多少息息相关。

角闪岩的测井特征为自然伽马低值, 多低于40 API, 密度数值较高, 在2.9 g/cm³以上, 中子数值多大于10 P.U, PE大于3.5 b/e; 矿物成分相对单一的花岗岩自然伽马数值较高, 在80 API以上, 密度数值约为2.6 g/cm³, 中子数值多小于5 P.U, PE约为2.2 b/e, 当长石含量增加时, 花岗岩的各项物理特征会有所变化^[18]。研究区内各种岩石的岩石物理特征统计见表1。

表1

Tab.1

表1(Tab.1)

表1 研究区太古界非储集层岩石物理特征值 Tab.1 Petrophysical characteristics of Archean non-reservoir rocks in the study area 岩性 自然伽马/API 中子/P.U 密度/(g·cm-3) 光电吸收截面指数/(b·e-1) 声波/(μs·m-1) 电阻率 角闪岩 15~65 10~22 2.97~3.16 5.0~6.0 165~197 中等电阻率值 变质角闪二长花岗岩 50~200 -1~8 2.66~2.68 1.8~3.2 160~175 电阻率相对高值 变质花岗闪长岩 70~220 8~15 2.6~2.7 50.0~70.0 165~230 高电阻率

表1 研究区太古界非储集层岩石物理特征值 Tab.1 Petrophysical characteristics of Archean non-reservoir rocks in the study area

根据测井曲线特征、录井岩屑、钻井取心等资料综合分析识别, 结合该区地层发育的主要岩性以及不同岩性特征值的差异, 找出对岩性识别敏感的测井曲线, 如自然伽马、声波、补偿中子和岩性密度等, 应用二维或者三维交会图确定岩性(图3)。从自然伽马-补偿中子交会图(图3(a))、补偿中子-光电吸收截面指数交会图(图3(b))、补偿中子-密度交会图(图3(c))中可以看到, 不同岩性分布在不同的区域, 利用这些分布区域的差异能够很好区分岩性, 准确地识别角闪岩、花岗岩、变质角闪二长花岗岩等。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 岩性识别交会图Fig.3 Cross plot diagram of lithology identification

录井岩屑、钻井取心及薄片资料显示, 研究区内干热岩储层岩性主要为角闪岩和变质花岗岩, 块状构造, 反映构造的层理或层界面相对较少。储层空间主要是裂缝和孔洞, 有效裂缝主要是高导缝、不连续的(半)充填缝和微裂缝, 局部存在孔洞。根据储集空间的发育程度, 参考孔隙发育大小和裂缝发育程度, 将干热岩储层分为I类、II类和III类(表2)。

表2

Tab.2

表2(Tab.2)

表2 研究区太古界储层分类标准 Tab.2 Classification standard of Archean reservoir in study area 储层类别 孔隙度/% 电阻率/Ω ·m 电阻率相对值 裂缝发育程度 岩性 Ⅰ 类 ≥10.0 ≤200 ≥0.3 发育 角闪岩、花岗-角闪岩 Ⅱ 类 (5.0, 10.0) ≤2 000 [0.2, 0.3) 比较发育 角闪岩、角闪花岗岩 Ⅲ 类 [2, 5.0] ≤4 000 <0.2 不发育 变质角闪二长花岗岩、花岗岩

表2 研究区太古界储层分类标准 Tab.2 Classification standard of Archean reservoir in study area

从测井解释评价结果可以看到, 其中发育I类层、II类层主要岩性为: 角闪岩、花岗-角闪岩, 少量变质角闪二长花岗岩、花岗岩; 发育Ⅲ 类层的主要岩性为: 变质角闪二长花岗岩、少量花岗岩、角闪岩。研究区太古界井的储层按厚度统计结果显示: 花岗-角闪岩类储层占27.4%、角闪-花岗岩类储层占58.9%、花岗岩占12.6%、角闪岩占1.1%。

2.2.1 常规测井曲线识别储层

不同岩性在测井曲线上反映的特征不同, 研究区测井曲线在进入到太古界时整体自然伽马数值变化较大, 角闪岩自然伽马值为低值, 花岗岩自然伽马值为中高值, 由于角闪石和长石含量的不同导致伽马数值变化较大。当干热岩储层孔隙度和裂缝不发育时, 3种孔隙度曲线特征表现为中子测井值、声波时差低值, 密度高值, 电阻率呈现明显的高值特征, 自然电位曲线比较平直。在裂缝和孔洞发育段出现电阻率低值, 孔隙度数值偏离骨架值, 向孔隙度增大方向变化, 电阻率相对于致密层会有明显的下降。

(1)Ⅰ 类层。受孔洞和裂缝发育的影响, 储层被泥浆填充, 显示为电阻率数值相对于围岩有明显下降, 下降幅度较大, 且存在明显的幅度差。3种孔隙度曲线数值偏离骨架值, 声波时差和补偿中子测井值增大, 密度数值减小, 计算的孔隙度数值在10%以上; 自然电位曲线有较明显的负异常显示。图4中的169号和171号层, 岩性为花岗闪长岩, 电阻率数值在500 Ω ·m以下, 最低到130 Ω ·m, 围岩电阻率则在20 000 Ω ·m以上, 储层电阻率明显低于围岩电阻率, 且深浅电阻率数值差异明显, 具有明显幅度差, 所以将169号和171号层评价为Ⅰ 类层。

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 X井常规测井资料显示储层特征Fig.4 Reservoir characteristics shown by conventional logging data of well X

(2)Ⅱ 类层。测井曲线特征显示为电阻率数值相对于围岩有明显的下降, 下降幅度中等, 且存在较明显幅度差。3种孔隙度曲线数值偏离骨架值, 声波时差和补偿中子测井值增大, 密度数值减小, 计算的孔隙度数值为5%~10%; 自然电位负异常幅度较 Ⅰ 类层略小。图4中167号层, 相对于 Ⅰ 类层电阻率数值略高, 幅度差也略小, 3种孔隙度曲线显示孔隙度增大。电阻率下降略小, 因此被评价为 Ⅱ 类储层。

(3)Ⅲ 类层。相校Ⅰ 类层和Ⅱ 类层, Ⅲ 类层电阻率数值较高, 但低于围岩电阻率, 相对围岩下降幅度较小, 且具有幅度差。3种孔隙度曲线略有向孔隙度增大方向变化趋势, 或者变化不明显, 显示岩性较致密, 自然电位正异常不明显, 如图4中的166号层, 电阻率相对于围岩有小幅度的下降, 3种孔隙度曲线有增大趋势, 但相对于Ⅰ 类层和Ⅱ 类层增大幅度要小。

2.2.2 电成像资料识别储层

电成像资料在测量过程中, 沿井壁每0.254 cm进行一次采样, 测量阵列电极上的电流强度反映出电流流过的地层由于岩石结构或电化学性质的差异而引起的微电阻率变化。电成像资料可以直观地描述地层特征, 较其他测井具有相对明显优势^[19]。通过成像测井资料结合常规测井、地质录井等资料进行分析表明, 该区花岗岩和角闪岩, 多呈块状构造, 未见明显反映构造的层界面。高导缝、应力释放缝比较发育, 见重泥浆压裂缝, 裂缝轨迹多不完整, 被不同程度地充填, 还可见孔洞发育。

图5为该区Y井储层发育段中不同类型储层的电成像资料显示特征, 图5(a)中, 在亮色背景下的暗色团块状显示为花岗角闪岩中的孔洞发育段, 说明该段地层孔洞比较发育, 孔隙度较大, 为孔洞型储层特征; 图5(b)中可看到亮色背景下明显的暗色条纹特征, 连续的正弦线(红色)为典型的高导缝特征, 被低阻泥浆、泥质或导电矿物充填所致; 有些正弦线轨迹不完整的(黑色)是高角度的半充填缝, 还有一些微裂缝特征, 是裂缝型储层特征; 图5(c)中显示成片的暗色特征, 为溶蚀孔洞, 同时伴有连续的正弦线, 储层为孔洞-裂缝型储层。因此, 利用电成像资料可以直观的根据裂缝和孔洞的发育层段, 来进行储层划分。

图5

Fig.5

	Figure Option ViewDownloadNew Window
	图5 Y井不同类型储层电成像资料显示特征Fig.5 Electrical imaging characteristics for different types of reservoirs in Well Y

2.2.3 偶极声波测井识别储层

偶极声波测井是识别裂缝最有效的测井方法之一, 偶极声波测井主要是利用滑行波信号产生纵波、横波、斯通利波的数值大小与波形幅度变化进行地层分析。在裂缝发育带, 由于裂缝中充填的物质不同, 导致在裂缝附近地层波阻抗相差大, 会引起声波发生反射和散射, 使得地层中传播的纵波、横波、以及斯通利波的速度、幅度和波形特征有所变化^[20]。尤其是斯通利波, 在无裂缝的岩石中, 斯通利波能量变化不大, 但当地层中存在与井眼相交的裂缝时, 井内泥浆与地层中流体的相互流动会造成斯通利波能量的损失, 能量曲线负向偏移, 且变化剧烈, 并在裂缝的边界形成反射, 斯通利波波形此时会出现“人”字形状。裂缝段一般在单极变密度波具有如下特点: 条带颜色变浅, 由于波的到达时间变长, 波形向时间增大方向弯曲, 同时有波形幅度衰减、纵横波扰动错断现象。纵波、横波、斯通利波幅度衰减程度与裂缝的充填物质有关, 当流体充填裂缝之中, 充填的流体会导致地层径向波阻抗数值明显减弱, 衰减幅度也随之增大, 裂缝则为有效裂缝。

图6为X井干热岩裂缝发育段偶极声波裂缝示意图。图中3 825~3 845 m井段波形出现明显波幅度衰减和向时间增大方向弯曲, 有波形幅度衰减的响应特征, 为典型的裂缝发育段特征。3 773~3 780 m、3 863~3 875 m相应层段纵、横波扰动错断现象, 纵、横波及斯通利波有幅度衰减现象, 裂缝特征明显。

图6

Fig.6

	Figure Option ViewDownloadNew Window
	图6 研究区X井偶极声波裂缝示意图Fig.6 Fracture schematic diagram of dipole acoustic fracture in Well X in the study area

干热岩地热储层分类标准是以常规测井解释的岩性、电性特征和物性参数为基础, 同时依据偶极声波测井和电成像资料解释的裂缝发育程度做出综合评价。孔隙度的增大和相对于围岩电阻率的下降可以表征储层的发育程度, 鉴于该区前期未有明确的储层分类标准, 参照类似储集空间类型的储层划分依据, 建立研究区的储层分类标准。分类标准采用孔隙度-电阻率交会和自然伽马—电阻率相对值交会的方法进行划分(图7), 本分类标准在后期地热投产井中应用效果好, 证明适合本区域的储层评价。

图7

Fig.7

	Figure Option ViewDownloadNew Window
	图7 研究区干热岩储层分类标准交会图Fig.7 Cross plot diagrams for classification standard of hot dry rock reservoir in study area

参考自然伽马曲线, 依据孔隙度、电阻率、电阻率相对值等测井信息确定该区的储层分类标准: I类层孔隙度φ≥10.0%, 电阻率相对值≥0.3; II类层5.0%<φ<10.0%, 0.2≤电阻率相对值﹤0.3; III类层 2%≤φ≤5.0%, 电阻率相对值﹤0.2(表2)。

地应力的大小和方向是干热岩开发注采中注采井网部署、水力压裂设计和诱发地震评估等方面的重要基础数据^[21]。通常在对干热岩进行水力压裂时, 水力裂缝于井底裸眼处起裂, 并沿平行于最大水平主应力方向扩展^{[22, 23]}, 地应力方向的确定对于确定地热井网分布很关键, 可采用电成像资料的泥浆压裂缝或者诱导缝来判断地应力方向, 也可采用横波在传播过程中的变化进行判断。快慢横波的分离程度能够反映地层各向异性的大小, 快横波方位为地层各向异性的方向, 因此可以利用偶极声波的横波测井资料分析地层的各向异性, 判断地应力的大小。该区的成像资料和偶极声波资料显示: 整体上地应力以SWW—NEE向(方位约80° ~260° )为主, 由此判断最大水平主应力方向为SWW—NEE向或近EW向。

干热岩地热储层构成的首要条件是地层具有高温特性, 地层温度参数尤为重要。要从井筒内获得准确的地层温度需要钻井液静置7 d或者更长时间, 使井筒温度与地层温度平衡, 此时利用下测方式测量井温即可认为是地层温度, 但这种方法时效低, 很少采用。目前通常采用的方法是当地热井完钻时泥浆停止循环后, 采用时间推移的方式获取同一深度的多次井温测量数据来外推求取地层的温度, 此方法和常规电测井同时进行, 不需要长时间关井恢复, 可以节约时间成本。

假定油藏为一无限大、均质和各向同性的油藏, 可以近似用地层压力恢复的原理计算井底温度, 恢复满足下述关系^[16]

T=T0-K·logΔtt+Δt 。(1)

式中: T₀为地层原始温度, ℃; T为测井在井底测量时间t的井底最高温度, ℃; t为泥浆循环时间, 可以从钻井记录上查出, h; Δt为泥浆停止循环后开始起钻时间计时开始, 直到某次测井下到井底的时间为止, h; K为系数, 无量纲。

根据多次测量的井底温度数据, 利用外推法做T-Δt/(t+Δt)交会, 求得Δt趋近无穷大(即Δt/(t+Δt)=1)时的井底温度即为该井井底的地层温度^[24]。如图8所示, 可以推算X井底地层温度为155 ℃。

图8

Fig.8

	Figure Option ViewDownloadNew Window
	图8 X井外推法地层温度分析Fig.8 Formation temperature analysis by extrapolation method in Well X

沿着最大水平主应力方向设计X井与Y井一组生产井, 电成像资料显示: X井诱导缝特征较明显, 其走向反映测量段内主要为近SEE—NWW(方位约110° ~290° )向; Y井诱导缝走向SEE—NWW向(方位约100° ~280° )(图9)。通过常规测井资料和成像测井解释可以看出, 两井在同一储层裂缝走向基本一致。经过两次压裂后投产, 产出井水温高达100 ℃, 证明压裂后裂缝连通性好。

图9

Fig.9

	Figure Option ViewDownloadNew Window
	图9 X井与Y井解释评价Fig.9 Interpretation and evaluation result of Well X and Well Y