AMT在新疆三屯河地区地下赋水性研究中的应用
Application of audio magnetotellurics in the study of the subsurface water-bearing properties of the Santunhe area, Xinjiang
通讯作者: 张汉雄(1990-),男,助理研究员,硕士,主要研究方向为地热地质学。Email:444219401@qq.com
责任编辑: 王萌
收稿日期: 2022-11-22 修回日期: 2023-03-2
基金资助: |
|
Received: 2022-11-22 Revised: 2023-03-2
作者简介 About authors
杨明远(1989-),男,地质矿产工程师,学士,主要从事地质矿产勘查工作。Email:
在能源、矿产等勘探领域中,因为音频大地电磁测深法(AMT)具备勘探效率高、垂向分层能力强等优点,从而得到大范围的应用。本文采用基于数据空间的三维反演算法,对新疆三屯河地区地热勘查采集的音频大地电磁数据进行了三维反演。反演结果表明,三维反演在该地区避免了测线上不均匀地质体的影响,取得了更为丰富和直观的三维地电异常。结合收集的研究区地质资料,利用AMT数据三维反演结果,对研究区地下电性特征及和热储有关的地层赋水性进行分析,并结合研究区地热成藏规律推断了几处地热成藏的有利区。
关键词:
Audio magnetotellurics (AMT) is widely used for energy and mineral explorations because of its high exploration efficiency and high vertical resolution. Using a three-dimensional (3D) inversion algorithm based on data space, this study performed 3D inversion of the AMT data collected from geothermal exploration in the Santunhe area of Xinjiang. As indicated by the inversion results, the 3D inversion avoids the influence of inhomogeneous geobodies on the survey lines in the study area and yielded very rich and intuitive 3D geoelectric anomalies. In combination with the geological data and the 3D inversion results, this study analyzed the subsurface electrical properties of the study area and the formation water-bearing properties related to geothermal reservoirs, and finally inferred several favorable areas for geothermal reservoirs from the geothermal accumulation patterns of the study area.
Keywords:
本文引用格式
杨明远, 张汉雄, 马超, 杨海磊, 朱威.
YANG Ming-Yuan, ZHANG Han-Xiong, MA Chao, YANG Hai-Lei, ZHU Wei.
0 引言
地热资源的开发利用在经济发展和环境保护相互协调的可持续发展战略中具有重要的地位[1],赋水性特征是研究地热成藏的要素之一,查明其特征对于地热勘查工作具有重要意义。音频大地电磁法使用1 Hz~10 kHz天然场源信号,具有施工效率高、勘探深度可达千米范围、垂向分层能力强等优点,特别适用于深度在千米左右的能源、资源及地下结构勘查,尤其是对浅部特别是低阻层具有较高的分辨率[2-3]。数据采集方面,AMT由于其场源是天然场,信号微弱,在矿山、城区附近极易受到环境影响,但在地广人稀的西部地区较易获取高质量的数据。资料处理方面,与常规MT处理方法并无太大区别,使用一维或者二维反演时容易受到不均匀地质体的影响[4]。目前随着反演技术的不断发展和进步,大地电磁资料的反演处理进入到三维反演阶段 [5⇓-7]。传统基于模型空间的反演算法,模型参数的数量为M,反演计算M×M数量的灵敏度矩阵,当网格参数M很大时,计算量太大,这种情况下进行三维MT反演耗时费力,不满足实用条件。采用基于数据空间的反演算法可以解决上述问题。满足观测数据的参数数量N远远小于模型参数的数量M,这样只需要计算N维矩阵,大大减少了数据计算量。Siripunvaraporn等[8]于2005年实现了该方法的三维反演。同时基于GPU并行计算的三维MT反演软件陆续出现,大幅度提高反演速度,使大尺度模型的三维MT反演成为可能[9]。本次研究是通过分析在新疆三屯河地区地热勘查获取的AMT数据的三维反演结果,对研究区地下与热储有关的地层赋水性进行了分析,结合研究区地热成藏规律推断了几处地热成藏的有利区。
1 研究区概况
1.1 区域构造
图1
图2
1.2 地层岩性
准噶尔盆地最老地层为奥陶系,最新地层为第四系未胶结黏土和砂砾堆积,南缘主要为砂泥岩互层的陆相沉积地层[11]。其中研究区剖面穿过的露头所见地层如下:
下白垩统呼图壁河组(K1h),主要岩性为灰绿色泥岩、砂质泥岩与砂岩互层,夹紫色泥岩。与上覆胜金口组(K1sh)呈整合接触。
上白垩统东沟组(K2d),主要岩性为灰棕、灰红、砖红色砾岩夹红褐色砂质泥岩、砂岩、粉砂岩,富含钙质及少量钙质结核。与上覆紫泥泉子组(E1-2z)呈平行不整合接触。
古新统—始新统紫泥泉子组(E1-2z),属河湖相紫红、褐红色为主的砂质泥岩夹灰红色砂岩,底部为砾岩或石灰质砾岩,以底砾岩出现和颜色的改变为划分标志,故有红色层之称。与上覆安集海组(E2-3a)呈整合接触[12]。
始新统—渐新统安集海组(E2-3a)为一套湖湘沉积,多为灰绿色泥岩夹泥灰岩、薄层砂岩及介壳层,含腹足类、双壳类、介形类、鱼类化石,因多为灰绿色岩石而有绿色岩系之称。与上覆沙湾组(E3N1s)呈整合接触。
1.3 电性特征
表1 研究区周边喀拉扎地区岩(矿)石物性参数统计
Table 1
岩(矿)石名称 | 个 数 | ρs/(Ω·m) | ||
---|---|---|---|---|
最大值 | 最小值 | 几何平均值 | ||
呼图壁河组紫红色泥岩 | 30 | 143.00 | 22.38 | 54.35 |
呼图壁河组灰绿色砂质泥岩 | 24 | 75.91 | 22.47 | 46.61 |
呼图壁河组灰绿色细砂岩 | 36 | 257.67 | 72.53 | 151.59 |
呼图壁河组灰绿色粉砂岩 | 24 | 504.14 | 138.26 | 294.01 |
清水河组灰绿色钙质砾岩 | 36 | 725.02 | 188.20 | 408.32 |
清水河组紫红色钙质砾岩 | 36 | 347.51 | 100.23 | 255.21 |
清水河组灰绿色砾岩 | 36 | 855.58 | 221.00 | 455.84 |
清水河组灰绿色粉砂岩 | 36 | 251.73 | 66.03 | 102.14 |
齐古组紫红色泥质粉砂岩 | 24 | 22.88 | 5.38 | 12.6 |
齐古组紫红色砂质泥岩 | 24 | 43.05 | 5.74 | 18.55 |
2 数据采集与处理
2.1 数据采集
音频大地电磁法(AMT)是基于大地电磁法(MT)发展起来的,是以地下地质体的电性差异为基础,在地面接收天然电磁信号经过地下地质体与地质结构耦合之后的电场和磁场信号,并通过对其规律变化来研究地下结构,从而探测地下目标地质体。本次研究工作(图3)野外数据采集使用加拿大凤凰地球物理公司生产研制的V8多功能电法仪,使用AMTC-30磁传感器,野外测量点距50 m,采集4分量数据(Ex、Ey、Hx、Hy),测站的布极方式主要采用“十”字型,部分测点因地形等原因,有时也采用T”字型或“L”字型(图4)。南北、东西向电极距均为40 m。结合研究区地层情况,采集时间30 min,观测的有效频率为10 000~1 Hz,根据趋肤深度公式计算后,勘探深度能够达到探测目标深度的要求。
图3
图3
研究区地质及AMT测线位置(黑框为AMT三维反演范围)
Fig.3
The study area geology and AMT survey line position map(black box is the range of three⁃dimensional inversion of AMT)
图4
由于天然电磁场信号较弱,为获取高信噪比的野外测量数据,尽量选在地形开阔、平整的地方布设测站,两个电道方向之间的地面相对高差与电极距之比小于10%,远离电磁干扰源[14]。电极埋设应避开岩体、深挖极坑、浇灌盐水增强极罐接地的耦合能力等方法,使接地电阻不大于2 kΩ。电极线选用屏蔽线,同时防止其悬空切割磁感线产生感应电流,造成干扰。探头保持水平,两水平分量的探头相互垂直,间距大于10 m,探头的埋深大于40 cm,用土埋实,正北方向为x轴,正东方向为y轴。
2.2 数据处理
首先对采集到的时间域数据进行频域变换处理后得到Zxx,Zxy,Zyx,Zyy这4个复阻抗元素以及两种极化方式(TE和TM)的视电阻率和相位数据等,在三维反演过程中实际上主要用到的数据为Zxy、Zyx复阻抗元素,TE和TM两种不同极化方式得到的视电阻率ρTE、ρTM和相位φTE、φTM数据。同时可根据反演需要视情况选用测点的其他参数数据。
图5为2个测点的实测视电阻率及相位曲线,7线50测点位于紫泥泉子组地层中,由图中可见,TE与TM两种极化方式的视电阻率曲线基本重合,说明该测点基本没受到不均匀地质体带来的静态偏移的影响。虽然由于天然场信号噪声洞的影响高频部分数据连续性较差,但中低频部分数据的连续性尚可,使用EMD、人机联作等时间域信号去噪处理后仍能得到高信噪比数据[15],从而获得更科学准确的反演结果。其中5线250测点处在背斜核部,TE与TM两种极化方式的视电阻率曲线呈平行分开态势,说明该点的视电阻率可能受浅部、受局部电性非均匀体的静态偏移或测点附近的区域性断裂、各向异性地质体等具有明显的的三维特征地下电性结构的影响。在这种情况下,一维、二维反演结果与实际地下地质结构存在着较大的出入,影响推断解释结果。采用三维反演将静态偏移效应涵盖在三维正演响应中,同时拟合TE和TM极化方式的阻抗值,可以对地表局部非均匀体进行正确反演,从而获得更为准确的三维反演结果,更好地进行地质解释。
图5
图5
实测AMT视电阻率曲线
a—5线250测点;b—7线50测点去噪处理前;c—7线50点去噪处理后
Fig.5
Measured apparent resistivity curve of AMT
a—L5-250;b—before de-noising of L7-50;c—after de-noising of L7-50
研究区测点点距较小,三维模型剖分的网格数为75×70×17=89 250。电磁响应正演过程中,为了减小边界效应带来的影响,反演模型设置过程中向外扩展一定的距离,实际建立的反演初始模型的范围x和y方向上扩展为6~7 km,z方向上按海拔高度扩展为1.5~3 km。选用59个频点中的50个频点数据,频率范围为1.72~10 400 Hz。
3 地层电性特征
图6
使用经过去噪处理后的AMT数据进行了二维和三维反演(图7、图8),其中图8a为3、5、7号测线二维反演结果的栅栏图显示,浅部存在诸多凸起状小异常的同时整个电性特征形态同已知的地质及构造形态信息差别较大,不能准确地反映地下结构及其赋水性对电性特征的影响。AMT三维反演的电性特征更符合已知地下构造形态及赋水地层对电阻率的影响(图7、图8b),AMT测线南西至北东穿过出露四套地层是呼图壁河组(K1h),东沟组(K2d)、紫泥泉子组(E1-2z)和安集海河组(E2-3a),从昌吉背斜的核部向北翼方向延伸。结合前文对区内地层电阻率值特征的分析,AMT剖面穿过的地层电阻率值特征应为,呼图壁河组为中低阻、上白垩统东沟组高阻、紫泥泉子组中高阻、安集海组低阻的电性特征。但地层的电阻率值又同地层含水饱和度有着密切的关系,实际上研究区内AMT剖面的电性特征表现为,呼图壁组(K1h)处于测线西南端,总体表现为高阻层,电阻率相对高是由于该地层多为压实的泥岩,孔隙度较小,含水饱和度低,相对于其他砂岩和含水砾岩地层,电阻率值显示为相对高值;东河组(K2d)下部的电阻率低值可能是其中孔隙度较大的砾岩、砂岩中水饱和度高引起的,为储水的有利层位;紫泥泉子组(E1-2z)电阻率总体为相对高值,部
图7
图7
AMT三维反演结果
a—AMT三维反演结果;b—-500 m深度切片、测线剖面及100 Ω·m等值面组合
Fig.7
AMT 3D inversion results
a—display diagram of AMT-3D inversion results; b—-500 m depth horizontal slice, profile slice and 100 Ω·m isosurface combined graph
图8
图8
剖面反演结果
a—二维反演结果;b—三维反演结果切片
Fig.8
Profile inversion results
a—2D inversion results;b—slice of 3D inversion results
分低阻区域可能是由孔隙度大的砾岩充填水引起的,局部的高阻可能是孔隙度较大的砾岩层钙化形成的质密的不均匀地质体引起的电阻率值变高;北东端安集海河组(E2-3a)在测线末端,电阻率阻值较大,但测线覆盖范围较小,不具有代表性。从三维电阻率反演结果上看,整体上是一个低电阻率背景,与地质情况相符,但电阻率值又呈现明显的非均匀性。存在多处低阻、高阻异常区,未有明显的水平层状结构,同图6展示的同研究区距离最相近的CJ9909剖面中粉色框区所示的昌吉背斜地层产状较陡,倾向NW向构造样式吻合,与野外实际调查情况相符。
4 赋水性特征
图9
研究区中电性特征展布形态与地层的赋水性及水流体的连通性息息相关:首先是压实的泥岩地层由于孔隙度低、阻水性好,成为挡水层;松散的砂岩、砂砾岩地层由于孔隙度大,在良好的水源补给情况下,大量的水充满在孔隙之中,这种背景下,原本高阻特征的地层反而呈现出较泥岩地层更低的电阻率。在三维反演结果中也很明显的展现这个特征;图7中的凹陷区为低于阀值的区域(低电阻率地区,主要是东沟组高阻、紫泥泉子组地层),凹陷区的南部和东部低电阻率区深度更深,西南和东南部电阻率呈现高阻特征(主要是呼图壁河组、安集海河组地层)。
三维反演的电性特征和研究区的赋水性特征有很好的对应关系。同一地层由于赋水性的强弱展现出不同的电性特征,在切片图(图8b)的地层中有明显显示,两条红色虚线中间为研究区断层分布密集的区域,属于小范围断裂带,因此其所属的紫泥泉子组(E1-2z)地层中高阻特征是由于断裂带低水饱和度高引起的。压实的泥岩地层阻挡了地表水向下渗透及地下水由南向北的连通,研究区西南和东南部因其主要为呼图壁河组、安集海河组压实泥岩地层,在三维反演结果中电阻率呈现高阻特征。研究区中南部主要为东沟组高阻、紫泥泉子组砂岩、砂砾岩地层,在三维反演结果中电阻率呈现浅部极低阻、深部低阻特征,浅部极低阻是由于其浅部地层水主要由地表径流直接补给,水量充分、矿化度高引起的。从图7可以看出研究区S向和ES向各有一个地下水的补给通道,地下水向北传递流通时,遇北部和东北部的安集海河组泥岩地层阻隔,因此地下水更多的聚集在研究区中南部。同时F3断裂的存在,高倾斜地层中的破碎带具有加速地表水和深部地下水向东沟组高阻、紫泥泉子组砂岩、砂砾岩地层的流通作用。受研究区所处的构造和地层的控制,在浅部含水的砂岩、砂砾岩地层范围又向WN方向延伸趋势,在深部呈现出向WS向延伸趋势,这同研究区所在的昌吉背斜的展布方向具有高度的一致性。在电阻率等深图上(图10)西南和东北部存在两个高低阻转换区域,基本和区内断裂在走向上重合,断裂两侧的电阻率值差异较大,不完全重合说明地下水的补给通道主要为断裂,但在具有不同孔隙度的地层中,水的流通性差别较大,因此反映出来的地下赋水性的强弱差别也较大。
图10
图10
AMT三维不同反演结果深度切片与地热有利区预测
Fig.10
Depth slice of AMT 3D inversion results and prediction of geothermal favorable areas
5 结论与建议
研究结果表明AMT三维反演结果能够准确地探测存在电性差异的地层结构,结合前人的研究成果取得了对研究区地下电性和赋水性之间相关性的认识:东沟组及紫泥泉子组等孔隙度较大的砂岩、砂砾岩地层在高含水饱和度情况下赋水性更强,其电性特征表现为低电阻率值;致密的泥岩地层由于其隔水性导致其赋水性弱,表现出相对高阻特征。在此认识基础上进而得到了区内地下赋水性特征:地下水的补给通道主要为断裂,具有不同孔隙度的地层中水的流通性差别较大,因此反映出来的地下赋水性的强弱差别也较大,研究区中南部地下赋水性强,在浅部含水的砂岩、砂砾岩地层范围有向WN向延伸趋势,在深部呈现出向WS向延伸趋势,同研究区所处的昌吉背斜的展布方向具有高度的一致性。受制于物性资料的限制,本次研究未考虑地层含水矿化度对电阻率的影响,在浅部电阻率值发生急剧变化的部位是否是由于含不同矿化度的水所引起的值得进一步的研究。
参考文献
音频大地电磁法在地热勘查中的应用以福建省宁化县黄泥桥地区为例
[J]. ,
The application of audio frequency magnetotelluric method to the geothermal exploration:A case of Huangniqiao area,Ninghua County,Fujian Province
[J]. ,
我国大地电磁测深新进展及瞻望
[J]. ,
New progress and prospect of magnetotelluric sounding in China
[J]. ,
一维大地电磁测深几种反演算法的比较研究
[J]. ,
Comparative study of several inversion algorithms for one-dimensional magnetotelluric sounding
[J]. ,
大地电磁二维反演方法对比研究
[J]. ,
Comparative study of two-dimensional magnetotelluric inversion methods
[J]. ,
大地电磁三维反演方法综述
[J]. ,
Review of 3D magnetotelluric inversion methods
[J]. ,
Three-dimensional magnetotelluric inversion:Data-spale methed
[J]. ,DOI:10.1016/j.pepi.2004.08.023 URL [本文引用: 1]
AMT三维反演在浅部岩体勘探中的应用
[J]. ,
Application of AMT 3D inversion in shallow rock exploration
[J]. ,
准噶尔盆地南缘紫泥泉子组的时空展布及成因解释
[J]. ,
Temporal and spatial distribution and genetic interpretation of Ziniquanzi Formation in the southern margin of the Junggar Basin
[J]. ,
可探源音频大地电磁测深在新疆喀拉扎地区含铀砾石层勘探中的应用
[J]. ,
Application of soundable source audio magnetotelluric sounding in the exploration of Uranium-bearing gravel strata in Kalazar,Xinjiang
[J]. .
大地电磁测深资料的噪声干扰
[J]. ,
The noise interference of magnetotelluric sounding data
[J]. ,
矿集区大地电磁噪声处理方法及其应用
[J]. ,
Noise processing methods and application study of MT in the ore concentration area
[J]. ,
昌吉背斜构造浅析
[J]. ,
Analysis of Changji anticline Structure
[J]. ,
准噶尔盆地南缘齐古背斜构造分析及三维建模
[J]. ,
Structural analysis and 3D modeling of Qigu anticline in the southern margin of the Junggar Basin
[J]. ,
准噶尔盆地南缘硫磺沟地区水文地质特征及其对煤层气富集的影响
[J]. ,
Hydrogeological characteristics and its influence on coalbed methane enrichment in the southern margin of the the Junggar Basin
[J]. ,
准噶尔盆地南缘煤系水矿化度对低煤阶煤层气的影响
[J]. ,
The influence of the salinity of groundwater in coal measures on low rank coalbed methane in the south margin of Junggar basin
[J]. ,
准噶尔盆地南缘硫磺沟煤层气富集主控地质因素及有利区优选
[J]. ,
Factors affecting coalbed methane enrichment and CBM favorable area of Liuhuanggou area in the southern Jungger basin
[J]. ,
河南平原浅层地下水总溶解固体和水化学类型的分布特征
[J]. ,
Study on distribution characteristics of TDS and hydrochemical type of shallow groundwater in Henan plain
[J]. ,
新疆地热水分布与地质构造的关系
[J]. ,
Relationship between geothermal water distribution and geological structure in Xinjiang
[J]. ,
新疆维吾尔自治区地热资源调查与评价报告
[R]. ,
Investigation and evaluation report on geothermal resources in Xinjiang Uygur autonomous region
[R] ,
/
〈 | 〉 |