孤岛潜山凸起区是以古生界为基底的披覆背斜^[8], 大地构造位置^{[9, 10]}如图1所示。

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 孤岛凸起区大地构造位置示意图[9, 10]Fig.1 Tectonic location of Gudao buried hill uplift area[9, 10]

中生代, 该区整体抬升, 形成NE倾向的单面山, 潜山坡面倾角15° ~20° ; 燕山期, NE向孤南断裂和孤北断裂将孤岛潜山分割为相对独立的构造单元^[10]。潜山顶部位于Gg1孔附近, 发育面积约20 km²的“ 天窗” , 缺失上古生界与中生界。山顶两侧新生界下伏地层由古生代石炭系— 二叠系逐渐过渡为中生代侏罗系— 白垩系(图2)。除孤西、孤南及孤北断裂外, 潜山凸起区还发育规模较大的NNW向和NE向断裂, 对寒武纪— 奥陶纪碳酸盐岩裂隙岩溶的发育具有重要影响。

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 新生界下伏地层及奥陶系顶板埋深等值线图Fig.2 Depth contour of Cenozoic underlying strata and Ordovician top plate

孤西断裂: 走向NW, 倾向WS, 倾角40° ~50° (图3), 断层面较平直^{[11, 12]}。印支期, 断层开始活动, 为压性逆冲断层; 燕山中晚期, 区域应力场转变为NE-SW向拉张应力, 孤西断裂发生负反转成为正断层; 喜山期, 孤西断裂经区域扭张应力的改造, 古近纪沙河街组沉积停止活动^[13]。

图3

Fig.3

	Figure Option Fig.3 A-A' cross section of western Gudao fault'>ViewDownloadNew Window
Fig.3 A-A' cross section of western Gudao fault'>	图3 孤西断裂A-A'地质剖面图Fig.3 A-A' cross section of western Gudao fault

孤北断裂: 是孤岛潜山凸起区北缘断裂^[14], 其走向为70° ~75° , 长度约32 km, 为南盘上升、北盘下降的正断裂, 倾角约50° , 两侧奥陶系顶板落差达1.78 km(图4)。

图4

Fig.4

	Figure Option 地质剖面图 Fig.4 B-B' cross section of northern and southern Gudao fault'>ViewDownloadNew Window
地质剖面图 Fig.4 B-B' cross section of northern and southern Gudao fault'>	图4 孤北、孤南断裂B-B' 地质剖面图Fig.4 B-B' cross section of northern and southern Gudao fault

孤南断裂: 是孤岛潜山凸起区南缘断裂, 其走向为70° ~75° , 长度约42 km, 为北盘上升、南盘下降的正断裂, 倾角约45° , 两侧奥陶系顶板落差可达1.92 km(图4)。

孤岛潜山凸起区东北部地层发育较齐全; 中部潜山顶部缺失上古生界和中生界, 局部地段缺失古近系, 新近纪馆陶组呈角度不整合覆盖在奥陶系之上; 西南局部地段缺失中生界。

区内地层由老至新可划分为: 新太古代泰山岩群, 岩性为花岗片麻岩; 早古生代寒武系— 奥陶系, 岩性主要为碳酸盐岩, 厚900~1 200 m; 石炭系— 二叠系, 岩性为页岩及砂、泥岩, 夹煤层, 厚0~120 m; 中生代侏罗系— 白垩系, 岩性为砂岩和泥岩互层, 厚0~1 000 m; 新生代古近系和新近系, 岩性为砂岩与泥岩互层, 厚1 300~2 500 m, 西南薄, 东北厚; 第四系为松散碎屑岩, 厚200~250 m。

2.1.1 地层特征

2014年, 山东省国土资源厅对全省地层进行了统一划分, 将三山子组和炒米店组作为跨寒武纪— 奥陶纪的岩石地层单位, 建立了奥陶纪马家沟群, 自下而上分为东黄山组、北庵庄组、土峪组、五阳山组、阁庄组和八陡组^[15]。本文根据该地层划分方案, 以Gg2孔钻探资料为依据(图5), 由下而上, 总结孤岛凸起区古生代地层特征。

图5

Fig.5

	Figure Option ViewDownloadNew Window
	图5 Gg2孔寒武系— 奥陶系柱状图及测温曲线Fig.5 Column map and temperature curve of Cambrian-Ordovician stratigraphy at Gg2 borehole

馒头组(€ _2-3m): 包含原地层划分方案的馒头组、毛庄组和徐庄组, 岩性以紫红色、灰紫色、灰绿色页岩为主, 夹薄层灰岩、白云岩及含鲕灰岩。底部为灰色隐晶质白云岩, 含燧石结核和燧石条带, 厚约200 m。根据区域地层裂隙岩溶发育特征^[16], 推测该地层裂隙岩溶不发育, 一般不能形成热储, 但在断裂发育区段可能形成较好的热储。此外, 馒头组底部与下伏太古宇的角度不整合接触面附近可形成良好的热储。

张夏组(€ ₃z): 岩性为灰色厚层鲕状灰岩及微晶灰岩, 厚约190 m。据在2 523.53~2 527.18 m段采集的岩心, 灰岩内填隙物为浅色泥晶、亮晶方解石, 发育倾斜、水平及垂直裂隙, 裂隙内有方解石充填, 见缝合线构造。该段地层裂隙岩溶较发育, 可形成热储。热储的发育明显受构造控制, 一般在断裂带形成良好的热储。

崮山组(€ _3-4g): 岩性主要为薄层泥质条带灰岩, 夹竹叶状灰岩、黄绿色页岩, 厚约50 m。该段地层裂隙岩溶不发育, 一般不能形成热储, 但在断裂发育区段可能形成较好的热储。

炒米店组(€ ₄O₁c): 包括原地层划分方案的长山组和凤山组, 岩性主要为灰色泥质条带灰岩, 夹紫色竹叶状灰岩和黄绿色页岩, 厚约180 m。该段地层裂隙岩溶不发育, 一般不能形成热储, 但在断裂发育区段可能形成较好的热储。

三山子组(€ ₄O₁s): 原地层划分方案为冶里— 亮甲山组, 岩性为灰色结晶白云岩, 下段夹竹叶状灰岩, 含燧石结核, 厚约90 m。该段地层裂隙岩溶较发育, 可形成热储。热储的发育受构造控制明显, 一般在断裂带形成良好的热储。

东黄山组(O₂d): 相当于原地层划分方案的下马家沟组下段, 岩性为泥质白云岩、灰黄色角砾状灰岩, 厚约25 m。该段地层裂隙岩溶不发育, 一般不能形成热储。

北庵庄组(O₂b): 相当于原地层划分方案的下马家沟组上段, 岩性为深灰色灰岩、豹皮灰岩夹白云岩, 含少量燧石结核, 厚约90 m。该段地层裂隙岩溶较发育, 可形成良好的热储。

土峪组(O₂t): 相当于原地层划分方案的上马家沟组下段, 岩性为泥质白云岩、灰黄色角砾状灰岩, 厚约130 m。该段地层裂隙岩溶不发育, 一般不能形成热储。

五阳山组(O₂w): 相当于原地层划分方案的上马家沟组上段, 岩性为深灰色灰岩、豹皮灰岩及白云岩, 厚约70 m。该段地层裂隙岩溶较发育, 可形成良好的热储。

阁庄组(O₂g): 相当于原地层划分方案的八陡组下段, 岩性为浅灰色— 灰黄色白云岩、泥质白云岩, 局部夹角砾状灰岩, 厚约40 m。该段地层裂隙岩溶不发育, 一般不能形成热储。

八陡组( O2-3b): 相当于原地层划分方案的八陡组上段, 岩性为深灰色— 褐灰色灰岩, 夹少量泥质灰岩、白云质灰岩, 厚约40 m。该段地层裂隙岩溶较发育, 可形成良好的热储。

2.1.2 热储厚度

寒武纪— 奥陶纪, 研究区构造运动较微弱^{[17, 18, 19]}, 地层厚度稳定, 在潜山凸起区顶部Gg2孔揭露的寒武系— 奥陶系齐全, 说明在后期构造抬升中, 寒武系— 奥陶系未受强烈剥蚀。根据地层岩性及区域裂隙岩溶发育特征, 本文以灰岩、白云岩为主的地层(奥陶纪八陡组、五阳山组、北庵庄组, 寒武纪— 奥陶纪三山子组, 寒武纪张夏组)作为热储, 累计厚度约480 m。

2.2.1 裂隙岩溶发育特征

根据空隙成因不同, 裂隙岩溶可划分为不整合面风化壳孔隙、溶隙、溶孔与构造裂隙4种。不整合面风化壳发育深度约250 m。渤南洼陷渤深6孔中, 奥陶系埋深4 161~4 246.55 m区段的测试结果显示: 10件灰岩岩心样品的孔隙度为1.1%~3.3%, 平均值为1.93%; 3件白云岩岩心样品的孔隙度为12.8%~16.1%, 平均值为13.23%。埕东凸起Cg16孔(图1)中, 寒武系埋深1 584.55~1 804.96 m区段的测试结果显示: 13件灰岩岩心样品的孔隙度为1.0%~8.8%, 平均值为3.86%; 6件白云岩岩心样品的孔隙度为3.4%~23.7%, 平均值为17.2%^[20]。由此可知, 风化壳型岩溶的发育主要受地层和岩性控制, 白云岩的岩溶发育程度远大于灰岩。此外, 风化壳岩溶发育程度与地层埋深呈反比^[21]。

碳酸盐岩中的白云石对方解石或文石进行交代, 发生白云岩化作用, 在晶间或晶内形成溶隙、溶孔, 该类孔洞的连通性差, 一般不易形成裂隙岩溶热储有效空隙。

断裂对裂隙岩溶热储的形成具有关键作用, 孤岛潜山地区经历了印支期、燕山期— 喜山期多期构造运动, 断裂发育^[22]。在主要断裂带附近可形成良好的裂隙岩溶热储, 如位于断裂带附近的Gg1孔、Gg4孔的单井涌水量> 1 000 m³/d^[23]。

2.2.2 热储裂隙岩溶率

潜山顶部上古生界与中生界缺失区, 热储埋深较浅, 风化壳孔隙较发育, 且断裂发育, Cg16孔测试结果表明, 热储的平均裂隙岩溶率为10.53%; 潜山外围上古生界与中生界发育区, 热储埋藏较深, 风化壳孔隙不发育, Bs6孔测试结果表明, 热储的平均裂隙岩溶率为7.58%。潜山斜坡带中生界缺失区, 热储的平均裂隙岩溶率为9.05%。

目前, 研究区尚无以寒武纪— 奥陶纪碳酸盐岩裂隙岩溶热储为目标层的地热勘探孔, 石油勘探过程中在该热储的Gg1、Gg2和Gg4孔中发现了可供开采的地热资源, 其中Gg1和Gg4孔单井涌水量> 1 000 m³/d, Gg2孔单井涌水量48 m³/d^[24]。Gg1孔曾建有地热疗养院, 目前已废弃。

分析Gg2孔测温曲线形态(图6), 垂向上, 曲线明显呈2个直线段分布: 热储顶板以上的盖层, 温度曲线斜率大; 热储层, 温度曲线斜率小, 地层升温率为1.91 ℃/100 m。

图6

Fig.6

	Figure Option ViewDownloadNew Window
	图6 Gg2孔测温曲线图[25]Fig.6 Temperature logging curve of well Gg2[25]

确定热储温度时, 有勘探孔实测温度的地区, 以实测温度值作为热储平均温度; 无勘探孔实测温度的地区, 以热储中点处估算的温度作为热储平均温度。热储温度采用两段法估算, 推导的热储温度估算公式为

t_r=0.035 5h+20.39+0.019 1H/2 , (1)

式中: t_r为热储平均温度, ℃; h为盖层累计的厚度, m; H为寒武系— 奥陶系平均厚度, m。

3.1.1 计算方法

采用热储法^[26]计算寒武系— 奥陶系裂隙岩溶热储中蕴藏的热资源量, 计算式为

Q=A· d· [ρ _s· C_s· (1-φ )+ρ _w· C_w· φ ]· (t_r-t₀) , (2)

式中: Q为热资源量, J; A为热储面积, m²; d为热储厚, m; ρ _s为热储岩石密度, kg/m³; C_s为岩石比热容, J/(kg· ℃); φ 为热储孔隙度, 无量纲; ρ _w为地热水密度, kg/m³; C_w 为地热水比热容, J/(kg· ℃); t_r为热储平均温度, ℃; t₀为下限利用温度, ℃。

3.1.2 利用方向与下限温度

孤岛潜山凸起区是孤岛油田的主产区, 地热资源丰富, 浅埋的馆陶组下段热储地热资源可满足住宅供暖需求^[27]。由于寒武系— 奥陶系裂隙岩溶热储埋藏深、温度高, 本文将地热发电确定为开发利用方向。地热发电已有100 a以上的历史^[28], 自1967年俄罗斯建立世界第一座双循环地热电站以来^[29], 中、低温地热发电技术取得了进步; 美国阿拉斯加契那发电厂利用57 ℃的地热流体进行发电^[30]。20世纪70年代初, 我国建成了7个中、低温地热发电厂, 目前仍在运行的地热电站有: 广东丰顺县邓屋地热电站, 利用地热水温度92 ℃, 发电功率300 kW; 湖南宁乡县灰汤地热电站, 利用地热水温度98 ℃, 发电功率300 kW; 江西宜春县温汤地热电站, 利用地热水温度67 ℃, 发电功率100 kW^[31]。综合考虑, 孤岛潜山凸起区地热资源的下限利用温度t₀确定为70 ℃。地热发电潜力计算时, 其热效率系数取7%^[32]。

3.1.3 计算单元划分

根据区内寒武系— 奥陶系埋藏条件及上古生界、中生界缺失情况, 将研究区划分为11个计算分区(图7)。

图7

Fig.7

	Figure Option ViewDownloadNew Window
	图7 地热资源计算分区Fig.7 Zoning for geothermal resource calculation

3.1.4 地热资源量

采用式(1)估算热储平均温度, 采用式(2)估算孤岛潜山凸起区寒武系— 奥陶系裂隙岩溶热储中可供发电的热量为804× 10¹⁶J(表1)。按30 a服务年限、热效率系数7%, 估算的孤岛潜山凸起区寒武系— 奥陶系裂隙岩溶热储的发电潜力为595 MW, 为大型地热田。

表1

Tab.1

表1(Tab.1)

表1 研究区寒武系— 奥陶系裂隙岩溶热储可供发电的热量 Tab.1 Geothermal resource calculation results of Cambrian-Ordovician fractured heat reservoir in the study area 计算分区 面积/km2 顶板埋 深/m 温度/℃ 裂隙岩溶 率/% 热量/ (1016J) Ⅰ 1 14.43 1 900 86 10.53 29.50 Ⅰ 2 2.66 2 000 102 10.80 Ⅰ 3 2.89 1 700 80 3.69 Ⅱ 1 1.39 2 100 105 9.05 6.20 101.00 Ⅱ 2 18.78 2 300 112 42.90 Ⅱ 3 10.69 2 000 102 8.95 Ⅱ 4 1.43 2 500 120 15.90 Ⅱ 5 1.97 2 900 134 5.26 Ⅱ 6 1.48 1 900 98 Ⅲ 1 11.80 2 300 112 7.58 62.90 Ⅲ 2 50.59 3 400 152 517.00 合计 804.00 注: Ⅰ 1区采用Gg1孔井口温度; Ⅰ 3区采用Gg4孔井口温度。

表1 研究区寒武系— 奥陶系裂隙岩溶热储可供发电的热量 Tab.1 Geothermal resource calculation results of Cambrian-Ordovician fractured heat reservoir in the study area

3.2.1 地热水化学成分

以Gg1孔为例, 对寒武系— 奥陶系裂隙岩溶热储中地热水的水质进行评价。Gg1孔地热水的矿化度为17.38 g/L, 阳离子以Na⁺为主, 阴离子以Cl^-为主, 水化学类型为Cl-Na型(表2)。

表2

Tab.2

表2(Tab.2)

表2 Gg1孔地热水化学成分 Tab.2 Chemical components of geothermal water at well Gg1 阳离子浓度/(mg· L-1) 阴离子浓度/(mg· L-1) 矿化度/ (g· L-1) pH值 水化学类型 K+ Na+ Ca2+ Mg2+ Cl- S O42- HC O3- 161.05 5 620 676.4 85.40 10 048 < 3.0 515.70 17.38 6.8 Cl-Na

表2 Gg1孔地热水化学成分 Tab.2 Chemical components of geothermal water at well Gg1

由表2可知, Gg1孔地热水中的S O42-含量极低, 说明地热水处于还原环境, 且可能含有油气资源, 导致其经历了较强的脱硫酸作用。

3.2.2 质量评价

地热水的腐蚀性与结垢性是影响地热发电的主要因素。地热水结垢性评价采用PHREEQC软件计算主要矿物的饱和指数。结果表明, 在不同温度下, Gg1孔地热水文石、方解石和白云石的饱和指数均> 0(表3), 处于过饱和状态, Gg1孔地热水具有结垢趋势。

表3

Tab.3

表3(Tab.3)

表3 Gg1孔主要矿物不同温度下饱和指数计算结果 Tab.3 Saturation index of major minerals at different temperatures of Gg1 温度/℃ 饱和指数 文石 方解石 白云石 25 0.50 0.64 0.77 60 0.81 0.94 1.53 70 0.92 1.03 1.68 80 1.02 1.13 1.81 90 1.13 1.23 1.90

表3 Gg1孔主要矿物不同温度下饱和指数计算结果 Tab.3 Saturation index of major minerals at different temperatures of Gg1

地热水的Cl^-含量占阴离子摩尔当量的百分比为97.1%, 远大于25%, 可采用拉申指数(I_L)对地热水的腐蚀性进行评价^[33], 评价公式为

I_L=(C_Cl+ CSO4)/Q_ALK , (3)

式中: I_L为拉申指数(I_L< 0.5, 不腐蚀; 0.5≤ I_L< 3, 轻微腐蚀; 3.0≤ I_L≤ 10, 中等腐蚀; I_L> 10, 强腐蚀); C_Cl为氯化物浓度, mg/L; CSO4为硫酸盐浓度, mg/L; Q_ALK为总碱度, mg/L; Cl、SO₄与ALK均以等当量的CaCO₃表示。

计算表明, Gg1孔地热水的拉申指数I_L为33.53, 具有强腐蚀性。此外, 地热水中含少量油气资源, 开发过程中应注意对油气资源的回收利用。