Hydraulic Fracturing Characteristics and Fracture Mechanism of Granite under High Temperature and Pressure
-
摘要:
地热能作为一种分布广泛的可再生资源,已成为人类未来能源安全的重要保障。干热岩型地热能因其储量巨大和清洁、高效、可再生的特点具有重大的开采价值。目前干热岩开采通常需要通过水力压裂来制造渗流换热通道,因此,储层水力压裂特征和破裂机制研究对于干热岩资源的高效开采具有重要意义。笔者以常见的干热岩储层——花岗岩为研究对象,研究了不同围压、温度下水力压裂特征;分析了压裂过程中声发射特征演化规律;阐明了裂隙扩展规律及水力压裂影响因素;揭示了水力压裂破裂机制和裂隙扩展准则。实验结果表明:①常温和高温下花岗岩的水力压裂曲线均可以分为管孔充水、管孔增压、水力压裂、裂隙扩展4个阶段。②压裂后的裂隙均沿着管孔的深度方向扩展,随着围压的增加,裂隙长度增大,裂隙路径更为复杂,裂隙扩展类型包括沿晶断裂和穿晶断裂。③水力压裂特征受围压和温度的影响明显,温度一定时,围压越大,破裂压力越大,二者的关系近似呈线性;围压一定时,温度越高破裂压力越大,但这种相关性并非线性。④花岗岩水力压裂主导破裂机制为张拉破裂,但裂隙仍存在剪切破裂特征;破裂压力受围压、温度和岩石材料性质的多重影响。研究结果可为花岗岩热储层的压裂设计提供支撑。
Abstract:Geothermal energy, as a widely distributed renewable resource, has become a crucial safeguard for humanity's future energy security. Enhanced geothermal systems (EGS), particularly those involving hot dry rock (HDR), hold significant exploitation value due to their vast reserves and characteristics of cleanliness, high efficiency, and renewability. Currently, the extraction of HDR typically requires hydraulic fracturing to create permeable flow channels for heat exchange. Therefore, research on the hydraulic fracturing characteristics and fracture mechanisms of reservoirs is of great importance for the efficient exploitation of HDR resources. This study focuses on granite, a common HDR reservoir, to investigate the hydraulic fracturing characteristics under varying confining pressures and temperatures. The evolution of acoustic emission (AE) characteristics during the fracturing process was analyzed, the fracture propagation patterns were elucidated, the influencing factors of hydraulic fracturing were clarified, and the fracture mechanisms and crack propagation criteria of hydraulic fracturing were revealed. The experimental results demonstrate that: (1) The hydraulic fracturing curves of granite under both room temperature and high temperature can be divided into four stages: borehole water filling, borehole pressure increase, hydraulic fracturing, and fracture propagation. (2) Post-fracturing cracks predominantly propagate along the borehole's depth direction. As the confining pressure increases, the crack length extends, and the fracture path becomes more complex. The fracture propagation modes include intergranular and transgranular fractures. (3) Hydraulic fracturing characteristics are significantly influenced by confining pressure and temperature. At a constant temperature, the breakdown pressure increases with higher confining pressure, exhibiting an approximately linear relationship. Under constant confining pressure, the breakdown pressure rises with increasing temperature, though this correlation is nonlinear. (4) The dominant fracture mechanism in granite hydraulic fracturing is tensile failure, though shear failure characteristics are also observed. The breakdown pressure is influenced by multiple factors, including confining pressure, temperature, and the material properties of the rock. The findings of this study can provide theoretical support for the fracturing design of granite-based geothermal reservoirs.
-
Keywords:
- hot dry rock (HDR) /
- hydraulic fracturing /
- fracture propagation /
- fracture mechanism
-
新疆东天山铜镍成矿带是中国重要的铜镍矿集区(Feng et al., 2018; 师震等,2019),分布有香山、黄山、图拉尔根、葫芦等多个铜镍硫化物矿床,其空间分布可划分为两个带,分别受康古尔大断裂和黄山大断裂控制。位于黄山东至图拉尔根成矿带之间的浅覆盖区,是黄山断裂成矿带之间的空白区,成矿潜力巨大。新疆地矿局物化探大队、中国国土资源航空物探遥感中心、中国地质调查局西安地质调查中心等单位在该区域开展过1∶5万综合物化探调查和1∶5万、1∶20万重力、航磁调查以及1∶5万矿产地质调查等。新疆地质矿产勘查开发局第六地质大队根据区域航空磁测成果,在大黄山浅覆盖区发现多个航磁异常,利用综合物探方法推测出类似黄山东的基性-超基性岩体,并在钻孔中发现了含铜镍矿超基性岩体,初步估算矿体规模达到中型。尽管取得了一定的进展,但是该区域依然存在很大的找矿潜力,因此依据大黄山的找矿成果,对大黄山一带浅覆盖区内的物探异常进行筛选,发现红石岗南区域存在1∶5万航磁异常与重力异常套合较好,同时前人研究表明,该浅覆盖区域北侧的红石岗岩体存在形成铜镍硫化物矿床的潜力(王志福等,2012; 王亚磊等,2017; Feng et al., 2018),目前对红石岗南区域调查研究程度不高,该区存在大规模第四系覆盖,无岩体出露,1∶5万航磁数据表明在该区的负异常背景场上叠加少量局部正磁异常,异常形态多呈椭圆状,磁异常幅值较小,因此很难评估该区域存在的隐伏异常体,难以对其成矿潜力进行评价。因此本文在附近区域矿床及岩体研究的基础上,开展高精度地面磁测,对该调查区磁异常数据采用磁异常化极、垂向一阶导数、归一化总梯度法、二维剖面模拟和三维磁异常反演方法进行研究,实现对区内隐伏异常体和断裂构造的初步勘探研究,确定了异常体的数量、磁化强度、埋深、产状要素、规模及断裂构造的位置,为探测隐伏异常体和后续在该区域所开展的勘探工作,提供了可信的地球物理依据和勘查实例。
1. 区域地质和地球物理背景
东天山红石岗南调查区位于新疆维吾尔自治区哈密市东南约160 km,整体上大地构造位于准噶尔洋壳板块与塔里木陆壳板块的聚合处,中天山隆起和吐哈盆地之间(王亚磊等,2017),区内地层出露主要以泥盆系和石炭系为主,少量奥陶系、志留系和侏罗系(Xiao & Wang, 2017)。红石岗南调查区在康古尔-黄山韧性剪切带内(宫辰,2020),处于东天山铜镍成矿带东段,附近存在多处铜镍矿床,如图拉尔根、葫芦、黄山、黄山东、香山铜镍矿等,红石岗镁铁-超镁铁质岩体位于土墩-镜儿泉断裂和黄山-镜儿泉断裂之间,而红石岗南调查区位于黄山-镜儿泉断裂以南,整体处于镜儿泉矿集区(图1a)。在黄山-康古尔断裂的北部,主要岩层为梧桐窝子组深灰色-灰绿色的海相喷发岩,而南侧则出露干墩组深灰-灰黑色浅变质的硅质岩、泥质岩等。区内铜镍矿和镁铁-超镁铁质岩体主要分布于这两套地层中(王志福等,2012)。调查区内侵入岩以花岗岩类为主,镁铁-超镁铁质岩为辅,镁铁-超镁铁质岩主要为晚石炭世-早二叠世(王亚磊等,2017),铜镍矿床也主要赋存于这些岩体中,整体呈南西-北东向分布,年代学研究表明矿集区内铜镍硫化物矿床均形成于早二叠世(韩宝福等,2004; Mao et al., 2016; 师震等,2019)。
图 1 (a) 东天山地区主要铜镍矿床及镁铁-超镁铁岩体分布图(据Mao et al., 2014修改);(b) 东天山地区航磁异常图;(c) 红石岗南地区航磁异常图;(d) 红石岗南地区布格重力异常图Figure 1. (a) The distribution map of important copper-nickel deposits and mafic-ultramafic intrusion in the Eastern Tianshan Mountains ; (b) Aeromagnetic anomaly map of the Eastern Tianshan Mountains; (c) Aeromagnetic anomaly map of the southern Hongshigang area; (d) Bouguer gravity anomaly map of the southern Hongshigang area区域航磁调查成果显示(图1b)黄山-镜儿泉断裂以北存在诸多正磁异常,异常幅值大于200 nT,其中位于黄山和香山附近,夹于土墩-镜儿泉断裂和黄山-镜儿泉断裂的区域存在一片高磁区域,异常最强处幅值高于800 nT。该区域存在的矿床普遍位于土墩-镜儿泉断裂和黄山-镜儿泉断裂附近,矿床附近区域磁异常整体小于200 nT。雅满苏断裂形成连续的串珠状磁异常(图1b),而雅满苏断裂以北,大范围新生代沉积物覆盖区,磁异常很低。南侧中北天山北缘断裂东北部存在条带状高磁异常带,异常幅值高于200 nT。红石岗南调查区附近磁异常相对香山和黄山区域整体较弱(图1b~图1c),磁异常呈扁椭圆状,长轴方向约40°,长约8 km,宽约4 km,异常值140~250 nT,极大值282 nT,北侧伴生似圆形正磁异常,长约3.5 km,异常极大值约252 nT(图1c)。整个调查区对应在高重力和低重力异常之间的梯度带上,等值线呈密集束状,梯度变化率4.4 mGal/km(图1d)。
整个红石岗南调查区表层覆盖新生代沉积物(图1a),沉积层较厚,地形平缓,受控于区域性主干断裂带所派生的次级断裂或裂隙系统,可能为镁铁-超镁铁质岩浆的就位与成矿提供了良好的导矿和储矿空间(王亚磊等,2017)。由于岩浆活动是矿产形成的关键性因素之一(杨大欢等,2022; 肖丹等,2022),因此红石岗南调查区可能会存在铜镍硫化物矿床赋存于镁铁-超镁铁质岩体中。
2. 岩矿石磁性特征
由于调查区内覆盖层较厚,且针对该区域还未开展任何钻探工作,因此需要使用调查区附近岩石样品的磁性数据参考,根据前人的研究,区内岩石出露的主要岩石类型为闪长岩、辉长岩和苏长岩等(王亚磊等,2017; 师震等,2019),因此统计了附近铜镍矿床岩石标本的体积磁化率、剩磁和Q比(图2,表1),其中Q比是柯尼希斯贝格比(Koenigsberger ratio, Q比)(Koenigsberger, 1938),是剩余磁化和感应磁化的比值,$ Q={M}_{R}/{M}_{i}=\mathrm{N}\mathrm{R}\mathrm{M}[A/m]/\kappa \left[SI\right]\times H[A/m] $,其中$ {M}_{R} $表示剩余磁化,$ {M}_{i} $表示感应磁化,NRM是自然剩余磁化强度(Natural remanent magnetization),$ \kappa $是岩石样品的磁化率,H是当地地磁场强度,在文中计算时使用的地磁场强度为45.68 A/m (
57406.2 nT)。Q比大于1表明磁异常场主要是由剩余磁化引起的,Q比小于1主要是由感应磁化引起的,Q比等于1表明感应磁化和剩余磁化起同样的作用。一般而言大陆喷出岩具有较高的Q比(Parkinson & Barnes, 1985),但具有商业价值的铁矿床几乎都具有低Q比(Jahren, 1965)。虽然这些统计数据并不能代表所有岩性单元的真实值,但根据区内标本的实测值可以为磁异常的正反演提供有效的约束。表 1 区域岩矿石磁性范围统计表Table 1. Statistical table of magnetic range of rock and ore岩石类型 磁化率(10−5 SI) 剩磁(10−3 A/m) Q比 备注 区域变质岩 14-72 13-25 0.4-3.91 包括变粒岩、片岩等 沉积岩 50-101 50-100 1.10-4.44 包括砂岩、粉砂岩 花岗岩 19-50 10-50 1.15-2.22 图拉尔根、黄山东样品 闪长岩 12-458 11-150 0.17-2.22 图拉尔根、图拉尔根钻孔、黄山东样品 凝灰岩 40-290 11-50 0.15-2.22 图拉尔根、图拉尔根钻孔、黄山东样品 辉长岩 91-165 14-150 0.29-2.02 图拉尔根钻孔、黄山东样品 辉石岩 37-806 11-226 0.29-0.65 图拉尔根钻孔样品 橄榄岩 2160 -3610 196-792 0.20-0.45 黄山东、图拉尔根岩心、土墩井中、黄山井中样品 矿化橄榄岩 7470 -17800 1290 -12200 0.13-1.91 黄山东、土墩井中、黄山井中样品 铜镍矿石 3560 -18800 417- 7910 0.23-1.88 图拉尔根岩心样品 统计结果表明(图2),该区域的岩矿石磁化率的常见值处于101× 10−5~105 × 10−5 SI之间,剩磁常见值则为101× 10−3~105 × 10−3 A/m。橄榄岩、矿化岩石和铜镍矿石普遍具有高磁化率和高剩磁,与围岩存在显著的磁性差异(图2),和矿化有关的这些岩石磁化率高于103 × 10−5 SI,绝大多数样品的Q比小于1,这表明在矿化区域,岩矿石的感应磁化可能主导着区域磁异常。而调查区附近的区域变质岩、沉积岩和花岗岩的磁性偏弱(图2和表1),闪长岩、凝灰岩、辉长岩和辉石岩稍强,但磁化率整体小于103 × 10−5 SI,剩磁强度小于300 × 10−3 A/m。该区域所有岩矿石的Q比处于0.1~10之间,大部分样品的Q比小于1(图2),这表明该区域的围岩磁性可能也是由感应磁化主导的,当然Q比也与岩矿石中的磁性矿物的浓度和粒度有关(刘隆等,2021),需要进一步研究。因此,根据图2和表1的结果,含矿超基性岩体的磁化率普遍大于
6000 × 10−5 SI,因此在后续二维剖面模拟和三维反演中,假设异常体平均磁化率为6000 × 10−5 SI进行物性约束。3. 数据采集与处理
2022年由中国地质调查局西安地质调查中心在红石岗南调查区(图1)进行了100 m × 40 m网格的地面高精度磁测,测区约80 km2。数据采集和处理流程见图3,在获得地磁场总强度T后,通过日变改正和基于国际地磁参考场IGRF13th模型(Wardinski et al., 2020)的正常场改正(IGRF改正)之后得到地磁场异常$ \mathrm{\Delta }T $,整个调查区的研究区总磁场强度、磁倾角和磁偏角的平均值分别为57 335.79 nT、63.56°和−0.35°,磁测总精度为2.44 nT。
对磁异常数据进行化极处理,以消除斜磁化影响(Baranov, 1957; 骆遥,2013),并进行垂向一阶导数计算,更好地确定磁性体边界。而为了确定异常体的场源信息,对化极后的磁异常数据剖面进行归一化总梯度法计算。归一化总梯度法通过在跨越磁场源的截面上构造一个特殊变换的场(总归一化梯度),使其可以用来探测源位置(Elysseieva, 2019)。另一方面,为了估算到异常体的深度并圈定异常体形态,使用了Geosoft公司的Oasis Montaj软件中的2D GM-SYS建模模块,创建一个假设的地质模型,并计算磁响应(Xiangjin et al., 2017; Ekwok et al., 2019, 2022)。用于建模的两条剖面A’-A和B’-B数据来源于原始磁异常剖面数据,剖面走向分别为315°和334°,模拟中使用真实地形数据,观测高度为实际仪器高度,约2 m。
最后进行三维反演和可视化,在本研究中为了解决正则化平滑反演难以揭示真实的地质情况,模糊的边界不能很好地反映地质结构(Vanzon, 2006; Sun , 2014; Utsugi, 2019)和平滑度反演通常会低估恢复的磁化率的值(Sun, 2015)的情况,所以采取稀疏范式的三维反演方法(Meng, 2018)。在正则化中使用混合$ {\mathit{L}}_{p} $范数解决反演问题。利用迭代重加权最小二乘法(Iterative Reweight Least Square, IRLS)进行离散和评估,改写最小的模型组分,通过有限差分算子代替梯度项等步骤,得到最终的正则化函数。$ {\mathit{L}}_{p} $范数正则化是高度非线性的,随着迭代过程的进行和阈值趋近于0,正则化函数的重点关注模型值的范围,直到迭代持续到算法达到预定义的收敛标准。开源框架SimPEG(Cockett et al., 2015)已实现混合$ {\mathit{L}}_{p} $范数反演,本文中的三维模型反演工作均在SimPEG中进行,反演中使用地形数据为真实地形数据。
通过对磁异常数据的处理、二维剖面模拟、三维反演及可视化,可以确定异常体的数量、磁化强度、埋深、产状要素等。进一步结合区域地质资料和物性资料,判断异常体的岩性、规模和形态等(图3)。
4. 结果
4.1 平面磁异常解释
原始磁异常等值线图表现磁异常带总体呈现近SW—NE走向,异常中心存在高值,可达220 nT,异常高值呈条带状(图4a),主要异常位于调查区南侧,呈现近SW—NE走向。化极后磁异常形态突出,圈定了5个主要的磁异常(图4b)。C-1异常位于红石岗南调查区西部,整体异常幅值较小,幅值约210 nT,呈现近SW—NE走向,为条带状正磁异常(图4b)。因此该处异常体为近SW—NE走向,北侧等值线下降较缓慢,因此异常体埋藏可能较深并北倾。C-2异常是整个调查区的主要异常,位于调查区中心,呈现近SW—NE走向,椭圆状正磁异常,异常幅值约为160 nT(图4b),因此异常体为近SW—NE走向。通过对等值线异常分析,其正值范围很大,西北侧等值线下降缓慢,而东南侧下降很快,表明异常体向西北侧倾斜且埋藏很深。C-3异常位于调查区北部,为圆状正磁异常,幅值约160 nT(图4b)。C-4异常则处于调查区东北部,异常整体走向为近SW—NE走向,异常幅值约110 nT,东北侧等值线下降缓慢,而西南侧下降很快并出现较低值,表明异常体向东北侧倾斜且埋藏很深(图4b)。C-5异常幅值与C-4异常接近,但异常等值线西侧等值线下降缓慢,而东北侧下降很快并出现较低值,所以异常体可能向西侧倾斜(图4b)。同时化极磁异常显示在调查区北侧和南侧存在串珠状磁异常,与存在的断裂构造相关,其中北侧异常走向为南西-北东(图4b-c),垂向一阶导数更清晰地展示了该处异常(图4c蓝色箭头所指推测断层),可能是黄山-镜儿泉断裂。南部异常(图4c紫色箭头所指推测断层)走向一致,可能为区内次级断裂。垂向一阶导数结果也显示出其他5个异常体的顶面局部隆起处的地质边界(图4c)。
图 4 (a)调查区原始磁异常平面等值线图;(b)调查区磁异常化极平面等值线图;(c)调查区化极磁异常垂向一阶导数;(d)剖面A’-A和B’-B的磁异常归一化总梯度计算结果Figure 4. (a) Contour map of the RTP magnetic anomaly in the survey area; (b) Contour map of the RTP magnetic anomaly in the survey area; (c) The first vertical derivative of the RTP magnetic anomaly in the survey area; (d) Normalized total gradient of magnetic anomalies for profiles A’-A and B’-B选取的红石岗南调查区磁异常数据两条剖面A’-A和B’-B(图4a显示剖面位置)应用归一化总梯度法进行计算,极值最大值对应的深度和埋藏异常体的中心深度相关。从图4d的结果来看穿越C-1异常体的剖面A’-A下方异常中心深度约200 m,C-1异常体中心对应的埋藏深度可达200 m,穿越C-2异常体的剖面B’-B下方异常深度与剖面A’-A类似,因此C-2异常体埋藏深度可达200-300 m,但规模大于C-1异常体。
4.2 二维剖面模拟和三维反演模型的地球物理解译
进行二维剖面模拟的2条剖面A’-A和B’-B沿着东南-西北方向(图4a),贯穿调查区西侧C-1异常体和中心地区C-2异常体,剖面走向分别为315°和334°(图4f),模拟误差均小于2.62 nT(图5)。在跨越C-1异常体的剖面A’-A观察到一个规模相对较小的磁性异常体,由于上层新生代沉积物覆盖层的磁性接近0,背景磁化率设置为1× 10−5 SI,剖面A’-A的磁异常主要来源于C-1异常体。C-1异常幅值较小,假设异常体磁化率为
6000 × 10−5 SI,得到的二维总磁强度模型显示在约240 m的深度存在约120 m厚度的异常体(图5a)。而跨越C-2异常体的剖面B’-B观察到一个规模较大的磁性异常体,假设异常体磁化率为6000 × 10−5 SI时,获得的二维总磁强度模型显示在约280 m的深度存在约160 m厚度的异常体(图5b)。二维剖面模拟的结果与归一化总梯度法确定的场源深度基本是一致的。三维反演模型和可视化处理更好地揭示了隐伏异常体的形态特征和分布规律(图6),恢复的观测数据(图6a)和实际数据拟合的很好(图4a)。在红石岗南调查区,三维反演异常体的最大磁化率为
10000 × 10−5 SI,假设剩余磁化与当前地磁场方向一致时,将剩余磁化强度转换为磁化率,根据图2中磁性强度最强的矿化样品计算,该区域样品磁化率最大值为17000 × 10−5 SI,反演恢复的磁化率结果在这个范围内,因此是可信的。通过反距离插值法对反演结果进行处理,生~成三维可视化模型(图6b~图d)。可观察到异常体走向为近SW—NE向(图6b~图d),存在5个主要的异常体,异常体磁化率接近6000 × 10−5 SI,主要分布在红石岗南调查区中部和北部,形态接近长轴圆柱体,中部为调查区最大异常体(即C-2异常体),北部存在多个较小的异常体,异常体整体向北倾,总体上异常体深度可能位于新生代沉积物覆盖层之下,约200~300 m之间。图 6 (a)三维反演恢复的磁异常;(b)三维反演磁化率模型200 m深度(海拔高度1180 m)切片;(c)三维反演磁化率模型264 m深度(海拔高度1116 m)切片;(d)红石岗南调查区隐伏异常体三维可视化图像,显示阈值0.06 SIFigure 6. (a) Prediction of magnetic anomalies; (b) 200 m depth (Altitude1180 m) slice of 3D inversion susceptibility model; (c) 264 m depth (Altitude1116 m) slice of 3D inversion susceptibility model; (d) Three-dimensional visualization image showed a threshold value of 0.06 SI of the underground anomalous body in the southern Hongshigang area需要注意的是,反演时模型区域进行适当的扩大以减小边界效应,对反演数据进行补空,边界处未出现畸变(图6b~图6d),但反演结果垂向分辨率有待提高,存在一定的细节缺失。整体上,多种方法分析磁异常的结果对调查区隐伏异常体的水平和垂直位置限定相一致。
5. 找矿潜力分析
通过在红石岗南调查区开展的地面高精度磁测,并利用磁异常化极、归一化总梯度法、二维剖面模拟、三维磁异常反演确定了红石岗南调查区异常体的数量、磁化强度、埋深、产状要素及规模。化极数据显示存在5个正磁异常,走向普遍为近SW—NE向,异常体向东北侧倾斜且埋藏很深(图4a、图4b),北部和南部的串珠状异常指示调查区内存在2条近SW—NE向的断裂(图4b-c),其中北部为前人推断的黄山-镜儿泉断裂,南部为调查区内次级断裂(图4c)。利用归一化总梯度法确定了异常体的埋深可能达到200~300 m(图4d),二维模拟和三维可视化模型进一步对异常体的形态进行了精细刻画,异常体平均磁化率值约
6000 × 10−5SI,走向主要为近南西-北东向并向北侧倾斜。在对红石岗南异常体的岩性进行判断时,需要考虑该区域存在大规模新生代沉积物的覆盖,根据对北山裂谷带梧桐窝子组、甘泉组沉积物(He et al., 2021)和黄山~镜儿泉侏罗系煤窑沟群底砾岩(邵行来,2012)的磁性测试,其磁化率平均值小于100 × 10−5 SI,几乎没有磁性,因此红石岗南调查区的异常体与新生代沉积物和石炭统沉积岩无关。尽管该区域存在泥盆统火山岩,但除泥盆系下统大南湖组第四亚组凝灰岩、中基性火山凝灰岩显著高值(
13000 × 10−5 SI)外,泥盆系下统大南湖组的凝灰岩、火山角砾、粉砂岩大部分磁化率为100 × 10−5 SI和1000 × 10−5 SI量级,且此区域范围内泥盆统火山岩主要出露在土墩~镜儿泉断裂以北(图1b),所以红石岗北和红石岗的磁性异常体不是泥盆统火山岩。因此在红石岗区域异常体可能与闪长岩、花岗岩类及镁铁~超镁铁质岩体相关(图1b)。根据表1的统计结果,图拉尔根和黄山东的闪长岩的磁化率为12× 10−5~458 × 10−5 SI,剩磁量级为11× 10−3~150 × 10−3 A/m,远小于二维模拟和三维反演得到的异常体平均磁化率6000 × 10−5 SI。而花岗岩类的磁化率为19× 10−5~50 × 10−5 SI,剩磁量级为10× 10−3~50 × 10−3 A/m(表1),也远小于异常体平均磁化率6000 × 10−5 SI,东天山成矿带其他区域的闪长岩和花岗岩类也存在类似的磁性(邵行来等,2010; 刘璎等,2011; 惠卫东等,2011; 2012;Xiao, 2017; 吴功成,2018; 侯朝勇等,2021)。因此,红石岗南调查区的异常体不是闪长岩和花岗岩类。根据表1的结果,橄榄岩和矿化样品具有显著的高磁化率和剩磁值,结合区域地质情况,推测红石岗南调查区的隐伏异常体可能是镁铁-超镁铁质岩体,其平均磁化率约为
6000 × 10−5 SI,该区域新生代沉积物覆盖层厚度约50-140 m,下石炭统沉积岩厚度较大(李彤泰,2011),岩体埋深约200~300 m,走向主要为近SW—NE向,岩层向北侧倾斜。同时前人推断黄山-镜儿泉断裂通过红石岗南调查区(图1b),化极磁异常和垂向一阶导数的结果(图4b~图4c)指示了黄山-镜儿泉断裂断层的存在,同时发现了调查区内与黄山-镜儿泉断裂走向一致的次级断裂(图4c)。区域断裂构造意味着在基底内的薄弱地带,其为侵入体创造了侵入条件,超镁铁质岩等岩石的侵入对于矿床形成具有重要的控制作用,当超镁铁质岩侵入到地下时,其高温和高压作用会促进地下矿物的熔融和矿床的形成,同时改变地下流体的成分和运移路径,影响矿床的分布(Yurichev , 2017; 石煜等,2022),已在黄山-镜儿泉断裂带上发现多个铜镍矿床(图1a)。因此,在红石岗南地区可能存在镁铁-超镁铁岩和岩浆对中泥盆统火山岩和石炭统沉积岩(图1b)地广泛入侵,源区由岩石圈地幔和软流圈地幔组成,岩浆入侵期间遭受地壳物质的混染(王亚磊等,2017; 师震等,2019)。而侵位于干墩组的镁铁-超镁铁质岩体普遍规模较大,且含矿性较好(尹希文,2015; 宋谢炎等,2022),因此红石岗南调查区存在的超镁铁质岩体具有较好的找矿潜力。6. 结论
通过对红石岗南的高精度地面磁测数据的处理,利用磁异常化极、垂向一阶导数、归一化总梯度法、二维剖面模拟、三维磁异常反演等方法,发现红石岗南调查区呈现出5个磁性隐伏异常体,其平均磁化率约为
6000 × 10−5 SI,埋深大约在200~300 m范围内。异常体主要走向为近SW—NE向,并且显示出岩体向北侧倾斜的特征。磁异常反演结果揭示了该区域存在规模较大的异常体和延伸较长的断层,其中调查区中心存在C-2主要磁异常体,处于浅覆盖层之下,厚约160 m,与南部推测断层接壤。异常体可能由磁化率较高的镁铁-超镁铁质岩体构成。而推断出的2处断层具有明显与区域大构造走向一致的特征,与地质背景相符合,为解译黄山东至图拉尔根成矿带之间的浅覆盖区地质构造提供了重要线索。磁异常特征和反演结果显示,红石岗南区域可能存在潜在的镁铁-超镁铁质岩体,这些岩体具有较高的磁化率。此外,推测的断层也为铜镍矿床形成提供了有利条件。因此,红石岗南浅覆盖区内具有进一步寻找基性-超基性岩型铜镍矿床的潜力。
致谢:感谢各位审稿专家提出的宝贵意见!
-
-
甘一雄, 吴顺川, 任义, 等. 基于声发射上升时间/振幅与平均频率值的花岗岩劈裂破坏评价指标研究[J]. 岩土力学, 2020(7): 41. Gan Yixiong, Wu Shunchuan, RenYi et al. Evaluation indexes of granite splitting failure based on RA and AF of AE parameters[J]. Rock and Soil Mechanics,2020(7):41.
郭亮亮. 增强型地热系统水力压裂和储层损伤演化的试验及模型研究[D]. 吉林大学, 2016. Guo Liangliang. Test and model research of hydraulic fracturing and reservoir damage evolution in Enhanced Geothermal System[D]. Jilin University, 2016.
胡剑, 苏正, 吴能友, 等. 增强型地热系统热流耦合水岩温度场分析[J]. 地球物理学进展, 2014, 29(3): 1391−1398. Hu Jian, Su Zheng, WuNengyou et al. Analysis on temperature fields of thermal-hydraulic coupled fluid and rock in Enhanced Geothermal System[J]. Progress in Geophysics,2014,29(3):1391−1398.
李全贵, 邓羿泽, 胡千庭, 等. 煤岩水力压裂物理试验研究综述及展望[J]. 煤炭科学技术, 2022, 50(12): 62−72. Li Quangui, Deng Yize, Hu Qianting et al. Review and prospect of coal rock hydraulie fracturing physical experimental research[J]. Coal Science and Technology,2022,50(12):62−72.
李馨馨, 李典庆, 徐轶. 地热对井系统裂隙岩体三维渗流传热耦合的等效模拟方法[J]. 工程力学, 2019, 36(7): 10. LI Xinxin, LI Dianqing, XU Yi. Equivalent simulation method of three-dimensional seepage and heat transfer coupling in fractured rock mass of geothermal-borehole system[J]. Engineering Mechanics,2019,36(7):10.
雷治红. 青海共和盆地干热岩储层特征及压裂试验模型研究[D]. 吉林大学, 2020. Lei Zhihong. Study on the characteristics of hot dry rock reservoir and fracturing test model in the Gonghe Basin, Qinghai Province[D]. Jilin University, 2020.
马兵. 热固耦合下储层岩石的地质力学特性及损伤演化规律研究[D]. 重庆大学, 2019. Ma Bing. Study on Geomechanical Behavior and Damage Evolution of Reservoir Rock under Thermo-mechanical Coupling[D]. Chongqing University, 2019.
谢紫霄, 黄中伟, 熊建华, 等. 天然裂缝对干热岩水力压裂裂缝扩展的影响规律[J]. 天然气工业, 2022(042-004). Xie Zixiao,Huang Zhongwei,Xiong Jianhua et al. Influence of natural fractures on the propagation of hydraulie fractures in hot dry rock[J]. Natural Gas Industry,2022(042-004).
熊峰. 裂隙岩体非线性渗流特性及水热耦合模拟研究[D]. 武汉大学, 2020. XiongFeng. Study on nonlinear flow behaviors and coupled hydro-thermal simulation of fractured rock masses[D]. Wuhan University, 2020.
Chen Y, Nagaya Y, Ishida T. Observations of Fractures Induced by Hydraulic Fracturing in Anisotropic Granite[J]. Rock Mechanics Rock Engineering,2015(4):48.
Deng J Q, Lin C, Yang Q, et al. Investigation of directional hydraulic fracturing based on true tri-axial experiment and finite element modeling[J]. Computers and Geotechnics,2016,75:28−47. doi: 10.1016/j.compgeo.2016.01.018
Fan T-G, Zhang G-Q. Laboratory investigation of hydraulic fracture networks in formations with continuous orthogonal fractures[J]. Energy,2014,74:164−173. doi: 10.1016/j.energy.2014.05.037
He J, Lin C, Li X, et al. Initiation, propagation, closure and morphology of hydraulic fractures in sandstone cores[J]. Fuel,2017,208:65−70. doi: 10.1016/j.fuel.2017.06.080
Kumari W G P, Ranjith P G, Perera M S A, et al. Hydraulic fracturing under high temperature and pressure conditions with micro CT applications: Geothermal energy from hot dry rocks[J]. Fuel,2018,230:138−154. doi: 10.1016/j.fuel.2018.05.040
Kumari W G P, Ranjith P G, Perera M S A, et al. Experimental investigation of quenching effect on mechanical, microstructural and flow characteristics of reservoir rocks: Thermal stimulation method for geothermal energy extraction[J]. Journal of Petroleum Science and Engineering,2018,162:419−433. doi: 10.1016/j.petrol.2017.12.033
Lin C, He J, Li X, et al. An Experimental Investigation into the Effects of the Anisotropy of Shale on Hydraulic Fracture Propagation[J]. 2017.
Lin C, He J, Li X, et al. An Experimental Investigation into the Effects of the Anisotropy of Shale on Hydraulic Fracture Propagation[J]. Rock Mechanics and Rock Engineering,2017,50(3):543−554. doi: 10.1007/s00603-016-1136-4
Mao R B, Feng Z J, Liu Z H, et al. Laboratory hydraulic fracturing test on large-scale pre-cracked granite specimens[J]. Journal of Natural Gas Science and Engineering, 2017: 278−286.
Sebastian B, Michael M, Ferdinand S C, et al. Hydraulic and Sleeve Fracturing Laboratory Experiments on 6 Rock Types, Andrew P B, John M, Rob J, editor, Effective and Sustainable Hydraulic Fracturing, 2013: Ch. 20.
Shao H, Kabilan S, Childers M I, et al. Environmentally friendly, rheoreversible, hydraulic fracturing fluids for enhanced geothermal systems[J]. Geothermics,2018,58(Mar.):72.
Stephen, Rassenfoss. Digging Up New Information On What Fractures Really Look Like[J]. Journal of Petroleum Technology,2018,70(3):38−39.100.
Tan P, Jin Y, Han K, et al. Analysis of hydraulic fracture initiation and vertical propagation behavior in laminated shale formation[J]. Fuel,2017,206:482−493. doi: 10.1016/j.fuel.2017.05.033
Zhuang L, Kim K Y, Jung S G, et al. Effect of Water Infiltration, Injection Rate and Anisotropy on Hydraulic Fracturing Behavior of Granite[J]. 2018(3): 1-15.