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、,():,():张云,高亮,田增霞,等 唐山马头营干热岩 井钻井工艺 中国地质调查,():,():齐晓飞,上官拴通,张国斌,等 河北省乐亭县马头营区干热岩资源孔位选址及开发前景分析 地学前缘,():,():王时林,张博明,乔海波,等 马头营凸起馆陶组低幅度构造油藏精细评价 特种油气藏,():,():,上官拴通 马头营区干热岩地热资源赋存分布特征及开发利用前景 能源与环保,():,():,许佳男 化学刺激剂对河北马头营区增强型地热储层改造作用研究 长春:吉林大学,:,赵阳升,万志军,康建荣 高温岩体地热开发导论 北京:科学出版社,:,张龙飞,张大明,潘苗苗,等 地球物理方法在干热岩不同勘查阶段中
27、的应用 中国科技信息,():,():杨冶,姜志海,岳建华,等 干热岩勘探过程中地球物理方法技术应用探讨 地球物理学进展,():,(),():雍世和,张超谟 测井数据处理与综合解释 北京:中国石油大学出版社,:,陶宏根,王宏建,傅有升 成像测井技术及其在大庆油田的应用 北京:石油工业出版社,中国地质调查 年 ,:,赵立新,陈科贵,王文文,等 声波测井新技术及应用实践 北京:石油工业出版社,:,上官拴通,孙东生,张国斌,等 唐山地区 深部地应力测量及断层稳定性分析 地质学报,():,():谢紫霄,黄中伟,熊建华,等 天然裂缝对干热岩水力压裂裂缝扩展的影响规律 天然气工业,():,():唐宜家,马天寿,陈力力,等 基于二维裂缝网络数值模拟的干热岩储层热采效率评价 天然气工业,():,():李约翰 试井 王福松,董恩环,译 北京:石油工业出版社,:,(,;,):,:;(责任编辑:常艳,王晗)