综合物探方法在滨海县月亮湾地热资源勘查中的应用
Application of a comprehensive geophysical exploration methods in the exploration of geothermal resources in Yueliangwan, Binhai County
第一作者:
责任编辑: 王萌
收稿日期: 2022-04-22 修回日期: 2022-07-31
基金资助: |
|
Received: 2022-04-22 Revised: 2022-07-31
采用可控源音频大地电磁法和广域电磁法在江苏省滨海县月亮湾进行地热资源勘查,对可控源音频大地电磁法数据进行近场和过渡场曲线辅助校正和反演约束,获得滨海港地下地热资源电性结构信息,同时采用微动探测方法获得地下几何结构信息。综合分析3种地球物理资料进行解释结果,获得研究区地热模型,确定异常位置,随后钻井深度2 919 m,出水量2 171 m3/d,水温51 ℃,综合物探成果与地质及地热井资料高度吻合,说明运用综合物探方法进行地热勘查,提高了勘探结果的可靠性。
关键词:
This study explored the geothermal resources in Yueliangwan, Binhai County, Jiangsu Province using the controlled source audio-frequency magnetotellurics (CSAMT) method and the wide-field electromagnetic method. Through the auxiliary correction of near-field and transition-field curves, as well as the inversion based on the CSAMT data, this study obtained the electrical structure information of underground geothermal resources in the Binhai port. Meanwhile, this study acquired the information on the underground geometric structure using the microtremor exploration method. By comprehensively analyzing the interpretation results of three kinds of geophysical data, this study obtained the geothermal model of the study area and determined the locations of the anomalies. A geothermal well with a depth of 2 919 m was drilled in the study area, obtaining water yield of 2 171 m3/d with a water temperature of 51 ℃. The high consistency between the results from the comprehensive geophysical exploration and the geological and geothermal well data indicates that the comprehensive geophysical exploration method can improve the reliability of geothermal exploration results.
Keywords:
本文引用格式
王军成, 赵振国, 高士银, 罗传根, 李琳, 徐明钻, 李勇, 袁国境.
WANG Jun-Cheng, ZHAO Zhen-Guo, GAO Shi-Yin, LUO Chuan-Gen, LI Lin, XU Ming-Zuan, LI Yong, YUAN Guo-Jing.
0 引言
1 研究区地质概况
研究区位于盐城市滨海县东部,地层隶属于扬子地层区,区内未见基岩出露,地表被新生界地层所覆盖,松散层(第四系Q +新近系N) 厚度大约为300~500 m。据钻孔揭露,研究区内基岩地层由老到新依次为奥陶系(O)浅海相碳酸盐岩、志留—泥盆系(S-D)海相碎屑岩、石炭—二叠系(C-P)浅海相碳酸盐岩及碎屑岩夹含煤碎屑岩以及白垩系(K)、古近系(E)陆相碎屑岩(图1)。
图1
研究区岩石圈厚度约80~100 km,居里面(560 ℃)深度约25~30 km,大地热流值约70 mw/m2,地温梯度低于27 ℃/km。按当地2.7 ℃/100 m的地温梯度计算,排除其他活动的影响,理论上本区含水层埋深1 000 m时,地热水温可达约42 ℃;埋深2 000 m时,可达约69 ℃。区内断裂构造非常发育,其中NE向断裂F5为研究区内控热、控水构造。研究区深部为奥陶系灰岩地层,在断裂构造作用下形成的裂隙是良好的含水热储层,其断裂构造裂隙层是本研究区地热水资源勘查的最佳目标层位。地热盖层主要有第四系、新近系盐城组上段,厚度为240~380 m,以黏土与砂层间互沉积;其中黏性土隔水性强,热导率相对较低,保温性良好。基于以上分析,研究区在“源、通、储、盖”这4个方面都基本具备了形成地下热水的地热地质条件。
2 地层物性特征
收集区域电阻率、波速参数综合统计于表1,从统计结果可看出地层岩石电阻率特征:新生界为低阻,上古生界碎屑岩为低阻、碳酸盐岩为高阻;地层电阻率从上而下表现为低—中低—高—中低—高,当地层破碎时电阻率明显降低。波速特征:基岩与上覆地层波速差异明显,由浅至深波速逐渐增大;当断裂引起基岩错动或明显破碎时,将出现地层的错动或缺失,波速明显减弱。上述物性特征为本次采用电磁法、微动勘查地层结构和断裂构造,进而推断地热资源提供了良好的电性前提条件。
表1 地层岩石电阻率参数统计
Table 1
地层 | 岩性 | 电阻率/(Ω·m) | 综合特征 | 波速特征/(m·s-1) | ||
---|---|---|---|---|---|---|
标本测定 | 测井 | MT反演 | ||||
Q | 砂质黏土、粉砂、黏土、砂砾层 | 16 | 13.2 | 52.1 | 低阻层 | 900~1800 |
N | 砂岩、泥岩 | 8.7 | 8.3 | 24.1 | 低阻层 | 1900~2400 |
E2-3s | 砂岩、泥岩夹玄武岩 | 3.9 | 6.6 | 12.7 | 低阻层 | 2200~4200 |
E2d | 砂岩、泥岩 | 6.1 | 11.6 | 低阻层 | ||
E1f | 砂、泥岩夹玄武岩 | 3.9 | 5.8 | 低阻层 | ||
E1t | 砂岩、泥岩 | 5.4 | 6.9 | 低阻层 | ||
K2p | 粉砂岩 | 136 | 18.6 | 34 | 中低阻层 | 4000~4700 |
P3d | 页岩、泥灰岩 | 388 | 68 | 99.8 | 中低阻层 | |
P2l | 页岩、粉砂质泥岩、岩屑砂岩 | 282 | 37 | 96 | 中低阻层 | |
P1g | 硅质岩、页岩、泥岩、粉砂岩 | 467 | 8 | 125 | 中低阻层 | |
P1q | 灰岩 | 4059 | 684 | 226 | 高阻层 | |
C2 | 灰岩 | 2538 | 529 | 153 | 高阻层 | 5000 |
C1 | 灰岩、砂岩、泥岩 | 1950 | 517 | 165 | 高阻层 | |
D3w | 石英砂岩 | 381 | 362 | 113 | 中低阻层 | |
S2m | 石英砂岩 | 351 | 404 | 134 | 中低阻层 | |
S1f | 粉砂质泥岩 | 94 | 198 | 44 | 中低阻层 | |
O3S1g | 泥岩、粉砂质泥岩、泥质粉砂岩 | 193 | 95 | 28 | 中低阻层 | |
O3w | 泥岩 | 151 | 26 | 中低阻层 | ||
O1-2 | 灰岩、白云岩、泥灰岩 | 1134 | 330 | 115 | 高阻层 | |
![]() | 白云岩、白云质灰岩、灰岩 | 1504 | 462 | 143 | 高阻层 | |
![]() | 炭质泥页岩、灰岩、白云岩、硅质岩 | 507 | 163 | 高阻层 |
3 方法技术
由于可控源音频大地电磁法(CSAMT)施工效率高,成本低,适合大面积施工,广域电磁法(WFEM)有效探测深度大,分辨率高,而微动施工效率低,成本高,故本次研究工作以CSAMT勘查为主, WFEM及微动探测复核及评价。布置3条平行测线均垂直于断裂构造F5(图2)。根据CSAMT法成果初步确定地热井井位,并在井位附近(L2线上)分别开展WFEM法和微动探测工作,以进一步确定地热井井位。
图2
3.1 CSAMT法
3.1.1 CSAMT法技术参数
CSAMT法采用赤道偶极装置观测与场源正交的磁场水平分量Hy及与场源平行的电场水平分量Ex,计算卡尼亚视电阻率
式中:f表示频率。
图3
3.1.2 CSAMT法资料处理
本次数据处理引入了AMT资料对近场和过渡场曲线进行辅助校正和反演约束,尽可能地还原近场和过渡场数据的真实性[9]。
图4
对CSAMT数据进行近区和过渡区校正后的视电阻率曲线基本等同于AMT测量结果[11],可利用较为完善的AMT的处理方法对校正后的CSAMT数据进行反演处理。
图5a是未经近场校正的CSAMT反演断面,在勘探深度约900 m之后形成一个高阻的基底,这显然是近场数据畸变造成的数据失真,同时也掩盖了深部地层和构造信息,浅部是远区的数据,对地层和构造的反应还是真实可信的。
图5
图5c则是结合AMT视电阻率曲线,对校正后的视电阻率曲线再次进行校正后的数据反演断面。经上述方法校正后剖面整体受到近场效应的影响得到较大改善,深部地质信息更加清晰,保证了勘探深度,而且其成果与后期钻探所掌握的地层深度高度一致。
3.2 WFEM法
WFEM法是相对于传统的CSAMT及MELOS法的基础上提出的,采用更大功率的人工源发射电流,只观测电场分量,使用适合于全域的公式来计算广域视电阻率,以达到探测不同埋深地质目标体的频率域电磁测深方法[12]。
本次WFEM法采用继善高科研发的广域电磁系统,AB偶极布设长度1.2 km左右,收发距大于12 km,接收测线信号均落在AB电极中间张角60°范围内,接收点距50 m,与CSAMT工作同点位,采集频率选择为0.011 7~8 192 Hz,发射电流最大146 A,WFEM法的野外工作方式与CSAMT法类似(图3)。
研究区内的道路、工地、电线等影响了电场采集信号,致使电场有一定的畸变,采取加权正演去噪处理(图6),曲线明显变得光滑。
图6
图6
去噪处理
a—去噪前;b—去噪后
Fig.6
Denoising processing
a—before denoising;b—after denoising
图7
图7
静态校正
a—处理前;b—处理后
Fig.7
Static correction
a—before correction;b—after correction
3.3 微动探测
微动探测通过圆形台阵获取微弱振动信号中面波(瑞雷波)的频散曲线,求取台阵下方S波的速度特征,进而推测地层、结构特征[16]。
本次微动探测采用美国A-tom单站式地震仪。观测时采用的是0.1 Hz检波器进行接收,每个测点周围均布设了13个采集器,以四重圆台阵观测系统进行采集,其半径分别为375、750、1 125、1 500 m,共计布置了6个测点(图8)。
图8
图9
微动探测是在计算出频散曲线的基础上,求取台阵中心点下方的视横波波速结构,由多个台阵测点的视横波速度进行插值,获得视横波波速剖面[18]。
4 资料综合解释
4.1 CSAMT法资料解释
经对近场和过渡场曲线进行辅助校正和反演约束后,3条CSAMT成果如图10所示,由图可以看出剖面清晰地反映出了3 000 m深度范围内的地层、构造特征。
图10
结合研究区地质、物性资料,将反演电阻率断面图上反映出的3 000 m范围内的地层大致分为3个大的电性层:第一层厚度约为300 m,具低阻特征,推测是第四系和新近系(Q+N)松散层的反映;第二层厚度约为900 m,呈中低阻特征,推测是志留系(S)海相碎屑岩地层的反映;第三层厚度较大,具高阻特征,推测为奥陶系(O)浅海相碳酸盐岩夹海陆交互碎屑岩地层的反映。
断裂构造推断依据:电阻率断面图横向上出现电阻率阶跃状突变或纵向上见陡立的低阻异常一般是判断断裂构造的主要标志。反演电阻率断面图上存在这些电阻率横向突变带,被推断由断层所致[19]。
从3条CSAMT反演电阻率剖面图看,均有明显的断裂构造迹象,并且3条剖面上的断裂构造在空间位置上对应较好,推断为同一条断裂,即DF1断裂(如图10)。该断裂分别穿过L1线1 200 m,L2线的1 250 m以及L3线的1 200 m附近,均呈明显的低阻异常反映。其中,L3线经过宽大的水面及建筑工地,干扰严重且测线无法完全展布,对断裂构造反映尚有欠缺。L1、L3线大号点存在左倾断裂,为干扰引起。
4.2 WFEM法资料解释
本次广域电磁法工作在L2线1.3~2.9 km范围内开展,其解释剖面如图11所示。由图可以看出地层呈明显的层状分布:第一层深0~0.3 km为低阻层,是第四系和新近系(Q+N)松散层的反映;第二层深0.3~1.2 km为中低阻层,是志留系(S)海相碎屑岩地层的反映;第三层深1.2~3.1 km为高阻层,是奥陶系(O)浅海相碳酸盐岩地层的反映;由于该方法探测深度相对较深,可发现3.1 km以下位置存在中高阻层,推测是寒武系(Є)浅海相碳酸盐岩夹碎屑岩的反映。
图11
剖面1.65 km附近,深1 km以下电阻率呈现明显的低阻反演,推测该位置存在断裂构造,即为L2线CSAMT成果中的DF1断裂。受其影响,深部岩层破碎含水,电阻率明显降低,尤其2.5 km以下呈现大范围低阻异常,与CSAMT法对比,结合DF1构造特征,WFEM法在深2.5 km以深所反映出的地质信息更符合实际,更加可信,体现出了其大勘探深度的优势。
4.3 微动探测资料解释
图12是微动横波速度解释断面,其波速仅反映地层相对速度。结合研究区地质资料,由图可以看出微动剖面范围内速度呈现明显的层状反映,第一层速度小于470 m/s,层厚约300 m,推测是第四系和新近系(Q+N)松散层的反映;第二层速度470~1 050 m/s,层厚800 m左右,推测是志留系(S)海相碎屑岩地层的反映;第三层速度大于1 050 m/s,推测为奥陶系(O)浅海相碳酸盐岩地层的反映。在建议的DR01孔深部附近,速度呈现明显的低速,推测可能是岩层发生了破碎,富水性较强的原因所致。
图12
4.4 异常综合解释及模型特征分析
由CSAMT法推断、WFEM法以及微动探测辅助评价印证结果表明, DF1断裂构造与八滩—小街断裂(F5)对应,结合研究区地质资料分析认为,该断裂构造为该区地热通道,建议的DR01地热井储水条件较好,含水层(组)以基岩构造裂隙水和岩溶水为主;其上覆300多米的第四系及新近系(Q+N)是相对较好的地热盖层。
研究区地热水主要由大气降水提供补给,其次是弱含水层的越流补给及构造破碎带内的垂直补给。当地下水贮存于构造裂隙以及岩溶裂隙之中时能够原地静态地接收来自地球深部的热量(由自然增温形成,而且温度和深度大致呈线性关系),同时围岩中的微量元素被不断吸收。控制性断裂F5(八滩—小街断裂)为NE向的张性断裂,而且F5处于断裂的开启状态,对地下热水的运移和储存都十分有利。总结出的研究区地热模型如图13所示。
图13
图13
研究区预测地热模型概念
Fig.13
Conceptual diagram of prediction geothermal model of rearch area
地热源:总体属于地温正常场区,地热主要是通过自然增温形成,以水流体的形式表现;储水层埋藏越深,温度通常也就越高。按当地2.7 ℃/100 m的地温梯度计算[20],排除其他活动的影响,理论上本区含水层埋深1 000 m时,地热水温可达42 ℃;埋深2 000 m时,可达69 ℃。
地热通道:深部热源沿断裂或由岩浆带至浅部;水源通过断裂或裂隙运移。NE向断裂F5为研究区内控热、控水构造。伴随着长期的构造活动,不仅加强了与深部热源之间的沟通,同时也能使深部含水层的富水性有所增强[21]。
地热储层:深部为奥陶系灰岩地层,裂隙是良好的含水热储层。
盖层:本区地热盖层主要有第四系、新近系盐城组上段,厚度为240~380 m,以黏土与砂层间互沉积;其中黏性土隔水性强,热导率相对较低,保温性良好[22]。
5 钻孔验证
表2 DR01地热井遇钻情况统计
Table 2
地层 | 厚度/m | 岩性描述 |
---|---|---|
第四系+新近系(Q+N) | 0~315 | 棕色、灰黄色粉砂质黏土为主,夹细砂、中粗砂 |
志留系(S) | 315~1116 | 灰白、灰绿色砂岩、长石石英细砂岩、粉砂岩为主,夹杂色、紫色泥岩,局部含硅质条带 |
奥陶系中—上统(O2-3) | 1116~2800 | 灰白色、灰色灰岩、细晶灰岩、泥灰岩为主,局部夹泥岩、页岩,裂隙较发育 |
奥陶系下统(O1) | 2800~2919 | 浅灰—灰白色白云质灰岩 |
图14
图14
CSAMT法综合解释剖面与钻孔验证对比
Fig.14
The comparison diagram of comparative chart of CSAMT and borehole verification
本次地热流体中阳离子以 Na+离子为主,阴离子以 Cl-离子为主,地热流体的水化学类型为 Cl-Na型水。该地热水矿化度为 3 388 mg/L,按矿化度分类属咸水; pH值7.5,按酸碱度分类属碱性水;热水中所含钙、镁离子的毫摩尔当量为11.41 mmol/L,按地下水硬度分类为极硬水;该井出水水温达 51 ℃,按地热资源温度分级为低温地热资源中的温热水,可用于理疗、洗浴、采暖、温室、养殖等。
6 结语
1)CSAMT法在资料处理过程中引入了AMT资料对近场和过渡场曲线进行辅助校正和反演约束的处理手段,提高了成果对深部构造和地层的反映能力,使反演推断结果更加准确。
2)WFEM法有效探测深度明显比CSAMT法要深,虽然本次CSAMT通过数据处理,其解释成果资料达到了3 000 m,但是其深部可靠性要低于WFEM法。
3)通过以CSAMT法为主、WFEM法及微动探测为辅的综合物探方法,获取了研究区深部电性、速度特征和模型特征。分析认为:研究区在“源、通、储、盖”这4个方面都基本具备了形成地下热水的地热地质条件,NE向断裂F5(八滩—小街断裂)及其周围较发育的裂隙为本区的储水构造和热源通道,保温盖层为区内相对较厚的松散沉积层,为裂隙型带状热储类型。
4)优选了地热井1处,并成功出水。其揭示的含水破碎带位置、地层层位与综合物探成果推断高度吻合,证实了综合物探方法技术的有效性,也表明了WFEM法在有效压制干扰的同时还能保证勘探深度,优势明显,可作为地热资源勘查的首选方法。
参考文献
甘肃临泽县城区深部地热资源调查评价
[J].
Investigation evaluation and method study of deep geothermal resources in Linze County,Gansu Province
[J].
综合物探方法在嘉兴地热勘查中的应用
[J].
The application of integrated geophysical exploration to geothermal exploration in Jiaxing City
[J].
综合物探在中深层地热勘查的应用研究
[J].
Research on the application of comprehensive geophysical prospecting in middle-deep geothermal exploration
[J].
综合物探方法在安徽省岳西县溪沸地热勘探中的应用
[J].
Application of integrated geophysical exploration method to Xifei geothermal exploration in Yuexi County of Anhui Province
[J].
综合地物化方法在张家口北部地热资源调查中的应用
[J].
The application of comprehensive geophysical and geochemical method in geothermal resources investigation in northern Zhangjiakou
[J].
苏北盆地滨海隆起月亮湾地热资源特征
[J].
Geothermal resource characteristics of Yueliangwan in Binhai Uplift of North Jiangsu Basin
[J].
可控源音频大地电磁法(CSAMT)在茅山地热资源勘查中的应用
[J].
The application of CSAMT in geothermal resource exploration in Mount Maoshan,Jiangsu Province
[J].
AMT近源干扰压制方法
[J].
Method of suppressing the adjacent source interference in AMT
[J].
CSAMT 资料处理中的静校正问题
[J].
Static correction of CSAMT date
[J].
基于牛顿迭代法和遗传算法的CSAMT近场校正
[J].
Near-field correction of CSAMT data based on Newton iteration method and GA method
[J].
广域电磁法在武陵山区页岩气勘探中的探索应用——以黔北桐梓地区为例
[J].
The application of WFEM to shale gas exploration in Wuling Mountain area:A case study of Tongzi area in northern Guizhou
[J].
广域电磁法在宜春某地地热勘查中的应用
[J].
The application of wide field electromagnetic method to geothermal exploration in Yichun
[J].
基于广域电磁法的页岩气有利区预测——以渝东北巫山地区为例
[J].
The prediction of shale gas favorable area based on WFEM: A case study of Wushan area in northeast Chongqing
[J].
华北地区瑞雷面波相速度层析成像
[J].
Rayleigh wave tomography of the phase velocity in north China
[J].DOI:10.1029/JZ054i001p00067 URL [本文引用: 1]
综合物探方法在水库堤坝隐患探测中的应用
[J].
Application of comprehensive geophysical exploration technique to detecting hidden defects of reservoir dams
[J].
利用SPAC法估算地壳S波速度结构
[J].
Application of SPAC method to estimate the crustal S-wave velocity structure
[J].
二维微动剖面探测“孤石”:以深圳地铁7号线为例
[J].
Mapping spherically weathered "Boulders" using 2D microtremor profiling method:A case study along subway line 7 in Shenzhen
[J].
CSAMT在沿海围垦区地热勘查中的应用
[J].
The application of CSAMT method to geothermal exploration in coastal tideland reclamation areas
[J].
江苏地热资料类型及开发利用前景
[J].
Types of geothermal data and its development and utilization prospects in Jiangsu Province
[J].
CSAMT在淮安金湖森林公园景区地热勘查中的应用
[J].
Application of CSAMT to geothermal exploration in Jinhu forest park scenic area in Huai’an
[J].
/
〈 |
|
〉 |
