首页 > 最新文献

Lithosphere最新文献

英文 中文
Development of a Large Shaking Table Test for Sand Liquefaction Analysis 开发用于砂液化分析的大型振动台试验
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-04-30 DOI: 10.2113/2024/lithosphere_2024_137
Chuanqin Yao, Huaqiao Zhong, Zhehao Zhu
This article highlights the development of a large shaking table test for sand liquefaction analysis. Two soil containers of different sizes were fabricated. The first one was small (0.87 m × 0.87 m × 1.20 m) in which the reconstitution and saturation methods could be easily tested. The dry tamping (DT) method was used to fabricate a model specimen. The subsequent field measurements suggested that the DT method provided a good distribution of sand grains in different cross sections. Before supplying the model specimen with water, carbon dioxide was flushed to replace air bubbles. This helped in obtaining a good degree of saturation, later verified by a digital moisture meter. For a given inlet water flux, the recorded pore water pressure displayed a quasi-linear trend, suggesting a good internal void system. This reconfirms the effectiveness of the DT method to yield homogeneous model specimens. The second soil container was huge (4 m × 4 m × 2 m) and used to explore liquefaction behavior in real engineering dimensions. Flexible foams were mounted on the side walls to mitigate the boundary effect. Although the boundary effect still manifested itself near the edges of the container during base shaking, half of the model specimen underwent a correct simple shear condition. For further analysis, vane shear tests were carried out before and after the liquefaction test. It was found that the intermediate layer, in general, suffered from the most severe liquefaction failure.Research activities into sand liquefaction have been conducted since the 1964 Niigata earthquake in Japan [1, 2]. In the laboratory, monotonic and cyclic triaxial tests are widely adopted to investigate liquefaction responses. As for laboratory element tests, the state parameter (considering both relative density and consolidation stress [3, 4]) and degree of saturation are two decisive indicators [5, 6] for examining the liquefaction potential. Besides, the soil fabric [7, 8] formed in different specimen reconstitution methods [9-11] has recently been proven to be another influential factor in controlling liquefaction triggering.Although triaxial tests certainly provide valuable insights into the mechanism of sand liquefaction, the understanding based on these tests is still limited by the size effect and thus only represents the liquefaction behavior of a unit soil element. This is far from representing a natural soil extent subjected to seismic loading in a semi-infinite space. Therefore, shaking table tests play an increasingly important role in the context of geotechnical earthquake engineering and contribute to improving the understanding of the liquefaction phenomenon. Many successful configurations have been presented in the literature [12-17]. Teparaksa and Koseki [18] performed a series of liquefaction tests on a shaking table to assess the effect of liquefaction history on liquefaction resistance of level ground. Ko and Chen [19] investigated the evolution of mechanical p
该方法包括在侧壁上粘贴柔性泡沫,以抑制入射波的反射。最后进行了一次完整的液化实验,以检验上述缓解策略的有效性,并检测砂土液化的敏感区域。首先给出了所用沙子的指数特性。第三部分介绍了利用小型土壤容器开发的重组和饱和方法,随后对其性能进行了详细讨论。下一节将介绍在大型土壤容器中进行的巨型振动台液化试验的实验结果。最后,我们讨论了在基础摇动过程中观察到的边界效应以及相关的现场结果。在本研究中,液化分析采用了法国贫级配参考砂 Hostun 31(HN31)。该砂的平均粒度为 D50 = 0.35 mm,均匀系数为 Cu = D60/D10 = 1.57,最小空隙率为 emin = 0.656,最大空隙率为 emax = 1.00,特殊重力为 Gs = 2.65 [20,21]。除了空隙率和固结应力[3, 22, 23]之外,沉积方法[9, 24]被认为是影响砂液化演变的另一个决定性因素。从一般意义上讲,在振动台上制作均匀的砂试样时,普遍采用气冲法这种简便的方法。然而,在大规模实验中,砂粒在重力作用下的自由下落不可避免地会在试验大厅中产生厚厚的砂尘,这对 MTS 称重传感器是有害的。此外,在本研究中,几乎不可能实现能够容纳如此大量颗粒组件的砂斗。作为替代方案,我们选择了干捣实(DT)方法,并首先使用小型容器进行检验,以检查重组试样的均匀性。Zhu 等人[25] 在不排水三轴试验中使用 DT 方法制作了元件试样,从而证实了这种方法的可重复性。为了全面了解内部均匀性,人们将叶片剪切试验作为识别 DT 方法制作的试样质量的标准[26]。在岩土工程中,叶片剪切试验是一种方便的现场方法,用于估算给定深度土壤单元的剪切强度[27]。该试验用于现场测量既快速又经济。试验采用了市售的叶片剪切设备(由英安扬公司生产),主要包括:(i) 多根可达到最大深度 3 米的延长杆;(ii) 精度大于 10%的扭矩计。因此,如图 2(a)所示,在实验室中使用 DT 方法制作了一个直径 18.7 厘米、高 40 厘米的元件试样,共分五层。控制每层的质量,使其达到相对密度为 50%的中密度状态。然后,将每一层依次放入一个分模中,并使用手持夯实机仔细夯实,以达到要求的 8 厘米厚度。叶片剪切读数在试样的中心部分进行,以尽量避免尺寸和边界效应,读数频率设定为每 5 厘米一次。从图 2(b)中可以看出,叶片剪切阻力 τf 随深度的增加呈线性增长,但在 10 厘米到 20 厘米之间略呈抛物线趋势。这一现象的合理解释是,由于深度较浅,砂等颗粒材料在较低的固结应力作用下通常具有较大的有效摩擦角。然而,这种微小的偏差并不能改变这样一个结论,即如果土柱相对均匀,则土的剪切强度应随深度的增加而呈准线性增加[26]。在小型容器中,用 DT 法重组初始高度为 60 厘米的模型砂柱,以达到 Dr = 70% 的致密状态。夯实技术包括:(i) 首先在每层表面放置一块坚固的木板,(ii) 然后落下重锤使模型试样致密,如图 3(a)所示。为了得到更好的结果,改变木板的位置重复这一过程。叶片剪切试验在五个位置(中央部分和四个角落)和四个深度(10、20、30 和 40 厘米)进行。由于在同一深度获得的叶片剪切力值的分散范围很窄(即 ±15%),因此所有实验数据都取平均值。 在理想情况下,土壤容器的侧壁应(i) 在现场条件下表现得与土体完全一样,尽可能保持应力场;(ii) 在基底震动时能够进一步改变其机械特性。由于砂土液化会导致土壤刚度和模量随着过剩孔隙水压力的增加而降低[29, 30],因此上述要求通常很难满足。尽管没有完美的技术解决方案来克服上述差异,但人们对模型容器的设计已达成普遍共识。文献中已经报道了一些例子,如 (i) 等效剪力梁容器 [31] 和 (ii) 层剪力箱 [15, 16] 等。然而,就本研究关注的规模而言,这些成功应用的成本太高,无法帮助减轻边界效应,尤其是模拟液化触发后的 "零刚度 "状态。为此,首先在土壤容器(4 m × 4 m × 2 m)中按照震动方向安装了几个 20 厘米的柔性泡沫,以抑制两侧波反射的产生,如图 10 所示。随后安装厚度为 5 毫米的定制防水土工膜,以防止孔隙水渗出模型试样。如图 11(a)所示,在沉积砂粒之前,在底部放置了一层由砾石组成的排水层,作为宏观的 "多孔石"。采用 DT 方法制作了一个 Dr 值约为 50%(4 层)的模型试样。两个接触层之间的界面被分解,以增加它们之间的内聚力。对于如此重的质量,用高频率的基础振动来对模型试样进行压实超出了现有设备的能力范围;因此,采用了动态压实机来代替人工压实。然后开始采用与小容器中相同的饱和技术(图 11(b))。如图 12 所示,本研究使用正弦波(地面加速度峰值 = 0.4 g,加载频率 f = 2 Hz)作为基础激励。输入运动在 5 至 20 秒之间进行梯形校正,以提高伺服系统的波形再现性。此外,虽然由于高压油的共振,振动台确实会出现一些轻微的运动,但在开始和结束时都分配了两个各约 5 秒钟的无运动安全时间间隔。在重组过程中,在模型试样中埋设了三个孔隙水传感器阵列(L = 左侧、ML = 左侧中部和 C = 中央),以捕捉液化演变过程。详细的仪器平面图见图 13(b)。图 14 显示了三个不同深度的过剩孔隙水压力 (EPWP) Δu 的时间历程。图中灰色不连续线表示液化的触发。在图 14(c)的底部,三条曲线几乎相互重叠。这表明底部的土壤元素经历了几乎相同的加载历史,其响应没有受到边界效应的影响。至于图 14(b)的中间部分,位于中心的曲线与位于 ML 的曲线接近。然而,它们与左角的曲线明显不同,这与底部观察到的情况相反。在图 14(a)中的顶部自由表面附近,也可以得出类似的结论。中心附近和中心处的两条曲线(C 和 ML)几乎遵循相同的趋势,与左角附近的曲线明显不同。上述现象表明,底部的边界效应非常微弱,可以忽略不计。随着深度的减小,边界效应发挥了更重要的作用,尤其是在顶部自由表面附近。这与之前文献[15, 16, 32]的结果一致。更重要的是,由于附着的柔性泡沫可以适当减少侧壁的波反射,从而保证液化分析所需的简单剪切条件,因此模型试样中心周围约一半的区域不受边界效应的影响。基底震动后,孔隙水压力立即开始消散。图 14 显示,在 t = 60 秒时,EPWP 几乎消散,ru 接近 0。 这意味着在随后的叶片剪切试验中,EPWP 消散只是次要的,尤其是在设备准备阶段持续约 1 小时之后。图 16 显示了叶片剪切读数与深度之间的关系。从图中可以看出,模型试样三分之二部分的剪切强度损失最为严重。这一现象可由以下两个独立因素解释。首先,模型试件的三分之二部分与顶部的自由排水面之间有一段距离。这无疑增强了排水状态,有利于过剩孔隙水压力 Δu 的积累。其次,该部分的覆盖层有效应力不如底层高。因此,Δu 的适度增加就足以引发砂土液化。本研究制造了两个土壤容器,以开发用于砂土液化分析的大型振动台试验。第一个容器较小,目的是寻求合适的重组和饱和方法。第二个容器很大,目的是以几乎实际的工程尺寸探索砂土液化。我们可以从获得的结果中得出以下结论:DT 方法有利于制作模型试样,因为它通常
{"title":"Development of a Large Shaking Table Test for Sand Liquefaction Analysis","authors":"Chuanqin Yao, Huaqiao Zhong, Zhehao Zhu","doi":"10.2113/2024/lithosphere_2024_137","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2024_137","url":null,"abstract":"This article highlights the development of a large shaking table test for sand liquefaction analysis. Two soil containers of different sizes were fabricated. The first one was small (0.87 m × 0.87 m × 1.20 m) in which the reconstitution and saturation methods could be easily tested. The dry tamping (DT) method was used to fabricate a model specimen. The subsequent field measurements suggested that the DT method provided a good distribution of sand grains in different cross sections. Before supplying the model specimen with water, carbon dioxide was flushed to replace air bubbles. This helped in obtaining a good degree of saturation, later verified by a digital moisture meter. For a given inlet water flux, the recorded pore water pressure displayed a quasi-linear trend, suggesting a good internal void system. This reconfirms the effectiveness of the DT method to yield homogeneous model specimens. The second soil container was huge (4 m × 4 m × 2 m) and used to explore liquefaction behavior in real engineering dimensions. Flexible foams were mounted on the side walls to mitigate the boundary effect. Although the boundary effect still manifested itself near the edges of the container during base shaking, half of the model specimen underwent a correct simple shear condition. For further analysis, vane shear tests were carried out before and after the liquefaction test. It was found that the intermediate layer, in general, suffered from the most severe liquefaction failure.Research activities into sand liquefaction have been conducted since the 1964 Niigata earthquake in Japan [1, 2]. In the laboratory, monotonic and cyclic triaxial tests are widely adopted to investigate liquefaction responses. As for laboratory element tests, the state parameter (considering both relative density and consolidation stress [3, 4]) and degree of saturation are two decisive indicators [5, 6] for examining the liquefaction potential. Besides, the soil fabric [7, 8] formed in different specimen reconstitution methods [9-11] has recently been proven to be another influential factor in controlling liquefaction triggering.Although triaxial tests certainly provide valuable insights into the mechanism of sand liquefaction, the understanding based on these tests is still limited by the size effect and thus only represents the liquefaction behavior of a unit soil element. This is far from representing a natural soil extent subjected to seismic loading in a semi-infinite space. Therefore, shaking table tests play an increasingly important role in the context of geotechnical earthquake engineering and contribute to improving the understanding of the liquefaction phenomenon. Many successful configurations have been presented in the literature [12-17]. Teparaksa and Koseki [18] performed a series of liquefaction tests on a shaking table to assess the effect of liquefaction history on liquefaction resistance of level ground. Ko and Chen [19] investigated the evolution of mechanical p","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"41 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-04-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141503461","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Quantitative Elimination of Seismic Pseudofaults and Fine Analysis of True Faults Underlying Igneous Rocks of No-Well Areas: A Case Study of Shuntuoguole Uplift in Tarim Basin 无井地区火成岩下地震伪断层的定量消除和真断层的精细分析:塔里木盆地顺图郭勒隆起案例研究
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-04-30 DOI: 10.2113/2024/lithosphere_2023_129
Ziran Jiang, Jian Jiao, Qiaomu Qi, Xingyu Deng
After multistage tectonic movement and evolution, large superimposed oil and gas basins generally developed many igneous rocks in the early rifting stages. The lithology and lithofacies of igneous rocks are complex, which is easy to lead to the distortion of the underlying migration velocity field and thus the response of seismic pseudofaults. Also, because of the obvious shielding and absorption effect of igneous rocks on seismic waves, the waveform quality of underlying strata is poor and the seismic response characteristics of faults are fuzzy. Currently, relevant studies have shown that the influence of igneous rock can be eliminated by the prestack depth migration with an accurate igneous rock velocity model. However, improving the accuracy of the velocity model needs to be corrected by well-logging data, resulting in poor applicability of the existing velocity modeling technology underlying igneous rocks without well, which is an obvious technical bottleneck. In this paper, the secondary strike-slip fault in Shuntuoguole low uplift of Tarim Basin, which has great oil and gas exploration potential but a very low degree of drilling, is selected as the research object. Aiming at difficult fault detection underlying igneous rocks caused by lack of drilling, the accuracy of fault seismic identification is improved by “interpretative fault preprocessing” and “fault sensitive attribute optimization.” In addition, through the “extreme hypothesis method” to maximize the complex migration velocity and simulate the underlying target layer distortion maximization, we realize the quantitative elimination of seismic pseudofaults. The practical application shows that this technology can determine the true and fake underlying faults quantitatively without establishing an accurate igneous rock velocity model. It is crucial not only for exploring oil and gas in the Tarim Basin’s secondary strike-slip faults but also for offering a method and technical guide for identifying faults in other basins affected by igneous rocks.Large superimposed oil and gas basins have undergone multiple periods of tectonic movement and evolution and generally experienced multiple periods of strong magmatic activity in the early stages of rifts or rifts, preserving numerous igneous rocks. As a high-velocity rock mass, igneous rock has two major impacts on the precise structural imaging of its underlying strata [1]. First, igneous rocks strongly shield and absorb seismic waves, leading to the blurring of seismic response characteristics of small structures and faults in the underlying strata. Second, the uneven distribution of thickness, the lateral and vertical variability of lithology, and the significant velocity differences between different lithologies of igneous rock bodies make it difficult to accurately describe the areal distribution, thickness, and velocity of high-velocity igneous rocks before migration imaging. This reduces the accuracy of the migration velocity field a
本文使用地震属性进行断层探测。大量实验发现,特征值一致性属性对大尺度断层敏感,曲率属性对中尺度断层敏感,自动断层提取(AFE)和似然性对小尺度断层敏感[30]。图 3 显示了不同阶段数据处理后的断层分布特征。与其他图相比,图 3(c) 对不同尺度的故障(图 3(a))进行了更清晰的描述,并显著降低了噪声(图 3(b))。从图 3(c)可以看出,大尺度断层(黄色箭头)的平面分布及其断裂边缘更加清晰,中尺度断层(粉色箭头)的平面连续性更好,小尺度断层(白色箭头)的数量明显增加。显然,所提出的断层精细预测技术能有效解决火成岩下地震断层的模糊响应特征问题。根据断层多属性融合图,对奥陶系目标层进行了平面和剖面精细断层解释,共识别出 65 条断层(图 3(d))。需要注意的是,为了提高预测精度,对原始地震数据量进行了断层增量,在突出真断层的同时,也增强了伪断层的地震响应。根据顺图郭勒隆起区二叠系火成岩的地震地质特征(图 4(a)),设计了正演模型并开展了正演模型研究。正演模型中间为高速度的英安岩,两侧为低速度的火成岩,因此模型水平方向存在岩性和速度的突变(图 4(b))。在正演建模中,采用与辉绿岩速度相近或比辉绿岩速度更低、更平滑的速度作为辉绿岩层的迁移速度,在速度突变点出现断层错动,与实际地震剖面上的断层非常相似。这证实了高速火成岩的迁移速度误差是火成岩下层产生假断层的原因(图 4(c))。此外,岩性与速度的突变点,即地层中高速火成岩厚度的突变点,也是地震伪断层的位置。众所周知,在油气勘探的早期阶段,提高火成岩的速度建模精度与降低钻井风险本质上是矛盾的。对于前者,必须首先通过钻井获得先验信息(如测井数据),以约束和校正背景速度场。对于后者,必须进行高精度速度建模,以获得精确的迁移成像数据,从而对断层和构造进行精细解释。这两项工作不可避免地要先后进行,而不能同时完成,这也是现有火成岩速度建模技术在无井地区应用的瓶颈。根据对顺图郭勒隆起主要走向-滑动断层中钻孔火成岩岩层间的岩性速度分析,认为虽然无法准确求解无井约束下高速火成岩的迁移误差真实值,但可以通过 "极端假设法 "估算迁移速度误差的范围和最大值。极限假设是一种基于数学最大化概念的技术,即在假设的极限状态下对复杂问题进行分析,使问题简单化、极端化、简明化[31]。火成岩成因复杂,岩性多样。特别是喷出面的火成岩,是岩浆沿地表流动后,在喷发及其重力的作用下逐渐凝结形成的。与爆炸面和火山沉积岩相比,喷出面的火成岩具有较高的间隔速度和较简单的岩性[32]。因此,在研究区域的不同地点,喷出面火成岩的区间速度相对稳定(表 1)。利用原始地震资料得出的高速火成岩背景速度和区域喷出面高速火成岩的区间速度,估算出迁移速度误差范围,进而确定最大误差值。
{"title":"Quantitative Elimination of Seismic Pseudofaults and Fine Analysis of True Faults Underlying Igneous Rocks of No-Well Areas: A Case Study of Shuntuoguole Uplift in Tarim Basin","authors":"Ziran Jiang, Jian Jiao, Qiaomu Qi, Xingyu Deng","doi":"10.2113/2024/lithosphere_2023_129","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_129","url":null,"abstract":"After multistage tectonic movement and evolution, large superimposed oil and gas basins generally developed many igneous rocks in the early rifting stages. The lithology and lithofacies of igneous rocks are complex, which is easy to lead to the distortion of the underlying migration velocity field and thus the response of seismic pseudofaults. Also, because of the obvious shielding and absorption effect of igneous rocks on seismic waves, the waveform quality of underlying strata is poor and the seismic response characteristics of faults are fuzzy. Currently, relevant studies have shown that the influence of igneous rock can be eliminated by the prestack depth migration with an accurate igneous rock velocity model. However, improving the accuracy of the velocity model needs to be corrected by well-logging data, resulting in poor applicability of the existing velocity modeling technology underlying igneous rocks without well, which is an obvious technical bottleneck. In this paper, the secondary strike-slip fault in Shuntuoguole low uplift of Tarim Basin, which has great oil and gas exploration potential but a very low degree of drilling, is selected as the research object. Aiming at difficult fault detection underlying igneous rocks caused by lack of drilling, the accuracy of fault seismic identification is improved by “interpretative fault preprocessing” and “fault sensitive attribute optimization.” In addition, through the “extreme hypothesis method” to maximize the complex migration velocity and simulate the underlying target layer distortion maximization, we realize the quantitative elimination of seismic pseudofaults. The practical application shows that this technology can determine the true and fake underlying faults quantitatively without establishing an accurate igneous rock velocity model. It is crucial not only for exploring oil and gas in the Tarim Basin’s secondary strike-slip faults but also for offering a method and technical guide for identifying faults in other basins affected by igneous rocks.Large superimposed oil and gas basins have undergone multiple periods of tectonic movement and evolution and generally experienced multiple periods of strong magmatic activity in the early stages of rifts or rifts, preserving numerous igneous rocks. As a high-velocity rock mass, igneous rock has two major impacts on the precise structural imaging of its underlying strata [1]. First, igneous rocks strongly shield and absorb seismic waves, leading to the blurring of seismic response characteristics of small structures and faults in the underlying strata. Second, the uneven distribution of thickness, the lateral and vertical variability of lithology, and the significant velocity differences between different lithologies of igneous rock bodies make it difficult to accurately describe the areal distribution, thickness, and velocity of high-velocity igneous rocks before migration imaging. This reduces the accuracy of the migration velocity field a","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"21 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-04-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141503520","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Method for Pressure Relief in Deep Coal Mine Roadways Using Borehole Groups and Its Application to Guqiao Coal Mine 利用钻孔群对深部煤矿巷道进行降压的方法及其在顾桥煤矿的应用
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-04-30 DOI: 10.2113/2024/lithosphere_2023_354
Yinyu Li, Li Qing, Yongshui Kang, Rui Zhang, Xiang Li, Bin Liu, Zhi Geng
With the increase in the mining depth of coal mines in China, the problem of large deformation of roadways owing to high-ground pressure has become prominent even under enhanced support systems. To reduce the high pressure on the surrounding rock, this study investigates a pressure-relief method for deep roadways using drilling borehole groups. Based on a deep roadway in the Huainan mining area of China, the influences of drilling parameters, such as borehole diameter, length, and arrangement were investigated. The results indicate that the fan-shaped arrangement of the borehole group can compensate for the dilatancy deformation of the surrounding rock. The peak stress of the surrounding rock is reduced and transferred to the inner part of the surrounding rock. Furthermore, a field experiment was conducted on an experimental roadway. The deformation of the roadway was monitored and compared with that of an adjacent roadway that did not apply the pressure-relief method. The monitoring results indicated that the deformation of the experimental roadway was significantly reduced.In China, coal is a major energy resource, which plays a dominant role in energy systems. China’s coal reserves are approximately 597 trillion tons, out of which approximately 53% are buried in the deep stratum (exceeding the depth of 1000 m) [1-3]. With the increase in the coal mining depth in recent years, an increasing number of roadways suffer from high geo-stress, which induces large deformation and failure and poses serious threats to mining safety [4-6]. Therefore, preventing large deformations of deep roadways with high geo-stress has become an important issue in coal mining. To prevent large deformation problems in deep roadways, engineers usually enhance the support system. For example, the use of high-strength and super-long bolts [6, 7], a method of bolting and shotcreting, U-steel support, grouting and floor bolting casting [3-5], and the reduction of the interval spacing of the supporting structures [4, 5]. Although some development has been achieved, large roadway deformations still frequently occur under high-ground pressure. The resistance offered by the supporting structure is extremely limited. Repeated repair is difficult and results in significant economic losses. Another approach to prevent large deformations of roadways is to release the high-ground pressure around the roadway using special measures. For example, floor grooving releases the high stress that accumulates on the floor, which is conducive to treating floor heaves [7, 8]. The rock mass within the range of pressure relief is destroyed by using high-pressure water injection softening and blasting pressure-relief methods, which reduces the elastic modulus and strength of the rock mass [5-7]. Drilling boreholes [6-9] in coal seams is beneficial for preventing coal and gas outburst accidents. Consequently, the accumulated energy on the surrounding rock surface decreases, leading to the release of
[20]通过建立模型实验,研究了深部巷道中钻孔直径对围岩力学性能的影响,并获得了合理孔距值的计算公式。Mishra 等人[21]提出了一种现场监测的钻孔泄压方法,用于评估单孔开挖巷道。Liu 等人[22]利用数值模拟分析了不同钻孔布置对围岩结构的破坏和高应力传递。Lan 等人[23]和 Guo 等人[24]模拟分析了钻孔直径、钻孔间距和煤应力对卸压的影响。Ma等人[25]和Song等人[26]都研究了巷道降压孔与螺栓的联合支护机理。然而,很少有研究关注深部煤矿巷道利用钻孔群降压方法的现场实验[27-29]。本研究通过数值模拟和现场实验研究了钻孔群对深部巷道稳定性的影响。提出了一种扇形钻孔组布置方式。此外,还在实验巷道上进行了应用实验。淮南矿区顾桥煤矿位于中国安徽省(见图 1)。淮南矿区顾桥煤矿位于中国安徽(见图 1),设计年生产能力为 1000 万吨。东翼巷道埋深超过 1000 米,原位应力测试最大值约为 32 兆帕。岩体主要由泥岩、砂质泥岩、粉砂岩和煤层等软岩组成。因此,尽管采用了加强支护方法,巷道仍经常出现巨大变形。严重的底板隆起和顶板下陷经常发生,需要反复维修。本研究选择了一条东翼带状巷道作为利用钻孔组进行降压方法的实验巷道。实验段长约 60 米(见图 2)。试验巷道上覆煤层 11-2,下伏煤层 8,11-2 煤层和 8 煤层无开采活动。该段上覆煤层 11-2 和下伏煤层 8 倾角较缓,煤层倾角约为∠2-4°。水文地质条件简单,主要水害为砂岩裂隙水,正常来水量为1-3 m³/h。监测点1#、2#布置在实验巷道内,监测点A1#、A2#布置在铁路巷道内,与皮带巷道相邻,不实施降压方案(见图2),地质条件和支护系统相似。开挖方法采用人工空气锤钻孔分层爆破法。巷道断面宽 5.6 米,高 4.6 米。岩层主要由砂质泥岩、粉砂岩和细粉砂岩组成。岩层倾角约为 3-8°。之前开挖的地质条件类似的巷道,虽然采用了加强支护方法,但巷道表面变形严重,开挖后 3 个月内巷道表面最大位移超过 500 毫米(见图 3)。假定泄压后,开挖引起的峰值应力会降低,促进高应力向围岩的深层传递,降低围岩发生大变形和破坏的风险。利用弹塑性力学理论建立了钻孔力学模型(见图 5)。为便于分析,假定这是一个平面应变问题。使用极坐标系是为了简化钻孔内壁抵抗均匀分布的径向应力 P1 的变形能力。应力平衡方程可表示为 [30, 31]:其中,σr 和 σθ 分别代表径向应力和圆周应力,r 为极坐标半径。 在降压应力和时间效应的影响下,钻孔壁岩体会发生坍塌,堵塞钻孔[23]。对实际钻孔坍塌条件的研究表明,钻孔形成后,原岩层中会形成孔洞效应,导致钻孔周围应力重新分布。岩层的原始应力有平衡和不平衡之分。沿钻孔中心轴线的周围岩体可分为四个区:破碎区、塑性区、弹塑性区和弹性区。在钻进过程中,钻具会扰动井壁周围的岩石,从而形成断裂带,岩石大多呈块状分布。由于钻孔的降压作用,周围的断裂发展成最不稳定的状态。由于巷道围岩的应力和时间的影响,断裂区域的岩体可能会发生崩落或坍塌,导致钻孔堵塞和堵塞。总之,在围岩中钻孔会形成一个降压区,破坏原岩层的应力平衡,重新分配钻孔周围的应力。钻孔过程会破坏岩体,降低岩体强度,导致孔周围出现裂缝。这会导致破碎区煤体和岩体稳定性变差,孔壁岩体坍塌,在卸压应力作用下堵塞孔口,最终达到消耗应力的目的。指定角度钻孔的目的是降低纵向应力集中区的应力,方便钻机施工。通过数值模拟和现场实验,研究了钻孔组对深部巷道稳定性的影响。通过数值模拟和现场实验,研究了钻孔组对深部巷道稳定性的影响,提出了一种扇形布置的钻孔组,并将其应用于顾桥煤矿的实验巷道。主要结论如下:高应力与低强度之间的矛盾是深部巷道产生较大变形的关键原因。数值模拟和现场监测结果表明,合理的钻孔降压措施能有效降低巷道围岩的高应力。数值模拟和现场试验结果表明,本研究提出的降压钻孔设计方案能有效减缓巷道大变形速率,降低巷道围岩变形。现场对比监测结果表明,两侧围岩变形可减少约 35.53%,顶板下垂可减少约 33.19%,底板隆起可减少约 35.85%。合理的钻孔卸压可有效降低集中应力,可作为防治深部高应力软弱围岩大变形灾害的辅助手段。数值模拟和现场监测结果表明,合理的钻孔降压措施可有效降低巷道围岩的高应力。现场对比监测结果表明,两侧围岩变形可减少约 35.53%,顶板下垂可减少约 33.19%,底板隆起可减少约 35.85%。合理的钻孔卸压可有效降低集中应力,可作为防治深部高应力软弱围岩大变形灾害的辅助手段。我们感谢淮南矿业(集团)有限责任公司在数据收集和岩石观测方面提供的支持。钻孔塑性区半径方程的推导见补充材料。文件中给出了钻孔轴向应力和轴向应变的表达式。确定了钻孔周围的应力边界条件。然后推导出 σθ、σr 和 σL 的表达式。钻孔临界屈服应力公式是通过获得钻孔周围应力公式并结合相关的 Tresca 屈服准则得出的。
{"title":"Method for Pressure Relief in Deep Coal Mine Roadways Using Borehole Groups and Its Application to Guqiao Coal Mine","authors":"Yinyu Li, Li Qing, Yongshui Kang, Rui Zhang, Xiang Li, Bin Liu, Zhi Geng","doi":"10.2113/2024/lithosphere_2023_354","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_354","url":null,"abstract":"With the increase in the mining depth of coal mines in China, the problem of large deformation of roadways owing to high-ground pressure has become prominent even under enhanced support systems. To reduce the high pressure on the surrounding rock, this study investigates a pressure-relief method for deep roadways using drilling borehole groups. Based on a deep roadway in the Huainan mining area of China, the influences of drilling parameters, such as borehole diameter, length, and arrangement were investigated. The results indicate that the fan-shaped arrangement of the borehole group can compensate for the dilatancy deformation of the surrounding rock. The peak stress of the surrounding rock is reduced and transferred to the inner part of the surrounding rock. Furthermore, a field experiment was conducted on an experimental roadway. The deformation of the roadway was monitored and compared with that of an adjacent roadway that did not apply the pressure-relief method. The monitoring results indicated that the deformation of the experimental roadway was significantly reduced.In China, coal is a major energy resource, which plays a dominant role in energy systems. China’s coal reserves are approximately 597 trillion tons, out of which approximately 53% are buried in the deep stratum (exceeding the depth of 1000 m) [1-3]. With the increase in the coal mining depth in recent years, an increasing number of roadways suffer from high geo-stress, which induces large deformation and failure and poses serious threats to mining safety [4-6]. Therefore, preventing large deformations of deep roadways with high geo-stress has become an important issue in coal mining. To prevent large deformation problems in deep roadways, engineers usually enhance the support system. For example, the use of high-strength and super-long bolts [6, 7], a method of bolting and shotcreting, U-steel support, grouting and floor bolting casting [3-5], and the reduction of the interval spacing of the supporting structures [4, 5]. Although some development has been achieved, large roadway deformations still frequently occur under high-ground pressure. The resistance offered by the supporting structure is extremely limited. Repeated repair is difficult and results in significant economic losses. Another approach to prevent large deformations of roadways is to release the high-ground pressure around the roadway using special measures. For example, floor grooving releases the high stress that accumulates on the floor, which is conducive to treating floor heaves [7, 8]. The rock mass within the range of pressure relief is destroyed by using high-pressure water injection softening and blasting pressure-relief methods, which reduces the elastic modulus and strength of the rock mass [5-7]. Drilling boreholes [6-9] in coal seams is beneficial for preventing coal and gas outburst accidents. Consequently, the accumulated energy on the surrounding rock surface decreases, leading to the release of","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"51 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-04-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140829516","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Influence of Rock Fabric on Physical Properties of Shale Oil Reservoir Under Effective Pressure Conditions 有效压力条件下岩石结构对页岩油藏物理性质的影响
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-04-30 DOI: 10.2113/2024/lithosphere_2023_338
Yu Ma, Suping Yao, Ning Zhu, Huimin Liu, Junliang Li, Weiqing Wang
The physical properties of shale oil reservoirs under overburden pressure are of great significance for reservoir prediction and evaluation during exploration and development. Based on core, thin section, and SEM observations, as well as test data such as XRD, TOC, and porosity and permeability under pressure conditions, this study systematically analyzes the variation of physical properties of different lithofacies shales in the Jiyang depression and the influence of rock fabric on the physical variation under pressure. The porosity and permeability of shale samples significantly decrease under pressure. According to the phased reduction in porosity and permeability, the pressurization process is divided into three pressure stages: low pressure (<8 MPa), medium pressure (8–15 MPa), and high pressure (>15 MPa). The reduction of porosity is fastest in the low-pressure stage and slowest in the medium-pressure stage. The reduction of permeability is fastest in the low-pressure stage and the slowest in the high-pressure stage. The rock fabric has a significant impact on porosity and permeability under pressure conditions. The permeability of laminated shale and bedded shale is higher than that of massive shale under pressure, and the permeability loss rate is lower than that of massive shales. Especially under lower pressure, the difference can be 10–20 times. In addition, the reduction rate of porosity and permeability under pressure is negatively correlated with felsic minerals content, which is positively correlated with carbonate minerals content and clay minerals content. The contribution of clay minerals to the porosity reduction rate is dominant, followed by carbonate minerals. The contribution of carbonate minerals to the permeability reduction rate is dominant, followed by clay minerals. The TOC content has no significant impact on the porosity and permeability of shales under pressure in the study due to the low maturity.With the change in global energy structure, shale oil and gas has become the core growth point of China’s oil and gas resources [1-4]. In the past decade, a series of important progresses have been made in the exploration and development of shale oil and gas in China, including breakthroughs in the exploration of shale oil in Junggar Basin, Ordos Basin, Jianghan Basin, Songliao Basin, and Bohai Bay Basin [5-8]. However, due to the heterogeneity of shales and the complexity of geological conditions in China, the prediction and evaluation of shale oil reservoirs still face many challenges [3, 9, 10].Many studies have shown that rock fabric, such as laminated structure and mineral composition, has a significant influence on the pore development and physical properties of shale oil reservoirs [10-14]. However, most of these studies were conducted under unpressurized conditions, and there are some errors with the formation conditions, which affect the prediction and evaluation of shale oil desserts. To recover the real physical
页岩油藏在覆盖层压力下的物性对于勘探和开发过程中的储层预测和评价具有重要意义。本研究基于岩心、薄片和扫描电镜观察,以及XRD、TOC、压力条件下孔隙度和渗透率等测试数据,系统分析了济阳凹陷不同岩性页岩的物性变化,以及岩石结构对压力作用下物性变化的影响。在压力作用下,页岩样品的孔隙度和渗透率明显降低。根据孔隙度和渗透率的阶段性降低,加压过程分为三个压力阶段:低压(15 兆帕)。孔隙度在低压阶段降低最快,在中压阶段降低最慢。渗透率的降低在低压阶段最快,在高压阶段最慢。在压力条件下,岩石结构对孔隙度和渗透率有很大影响。在压力作用下,层状页岩和层状页岩的渗透率高于块状页岩,而渗透率损失率低于块状页岩。特别是在低压条件下,两者相差可达 10-20 倍。此外,压力下孔隙度和渗透率的降低率与长石矿物含量呈负相关,与碳酸盐矿物含量和粘土矿物含量呈正相关。粘土矿物对孔隙度降低率的贡献占主导地位,碳酸盐矿物次之。碳酸盐矿物对渗透率降低率的贡献最大,其次是粘土矿物。随着全球能源结构的变化,页岩油气已成为我国油气资源的核心增长点[1-4]。近十年来,我国页岩油气勘探开发取得了一系列重要进展,包括准噶尔盆地、鄂尔多斯盆地、江汉盆地、松辽盆地、渤海湾盆地等页岩油勘探取得突破[5-8]。然而,由于页岩的异质性和中国地质条件的复杂性,页岩油藏的预测和评价仍面临许多挑战[3, 9, 10]。许多研究表明,层状构造和矿物组成等岩石结构对页岩油藏的孔隙发育和物性有重要影响[10-14]。然而,这些研究大多是在无压条件下进行的,与地层条件存在一定误差,影响了页岩油甜点的预测和评价。为了恢复地层条件下的真实物理参数,一些学者进行了加压实验,发现页岩的孔隙度和渗透率都随着有效压力的增加而逐渐降低,并提出了指数函数、二项式等函数模型来描述加压过程中孔隙度和渗透率随压力的变化[15-22]。然而,对于单个页岩样本而言,在相同压力条件下,不同样本之间的孔隙度和渗透率差异较大[17-19],造成这些差异的原因尚不清楚,有待进一步探讨。本研究选取了济阳凹陷的 Es3x-Es4s 页岩。通过岩石学、矿物学、有效压力下的孔隙度和渗透率实验以及 FE-SEM 观察,探讨了有效压力下页岩储层的物理特征以及岩石结构对孔隙度和渗透率的影响,为预测和评价济阳凹陷古新统页岩油藏提供依据。济阳凹陷位于中国渤海湾盆地东南部,由南向北依次为东营、惠民、沾化和车镇四个子凹陷(图 1(a))[23]。在近十年的勘探中,东营、沾化等子凹陷多层次、多类型页岩油取得重大突破,显示出济阳凹陷页岩油良好的勘探前景[24-26]。目前,Es4s、Es3x 和 Es1 是济阳凹陷页岩油的主要勘探区间(图 1(b)),其特点是碳酸盐岩含量高,平均含量超过 50%[23、27、28]。作为济阳凹陷页岩油的主要勘探区带,Es4-Es3x有许多工业油流井和采样井。 本研究的样品采集自 BYP5 井、F201 井和 FYP1 井(图 1(a)),基本信息见表 1。首先,从岩心样品中钻出一个与基底平行的塞状样品(直径 25 毫米,长 60 毫米)。然后,从塞子的一侧切割出一个 2 厘米高的微型塞子,用于薄片和扫描电子显微镜观察。剩余的塞子样品用于测量有效压力下的孔隙度和渗透率。在完成孔隙度和渗透率测量后,将塞子样品粉碎,用于总有机碳和矿物成分分析。使用蔡司显微镜观察厚度为 30 μm 的岩石薄片,以明确层理的发育特征和矿物分布。使用高分辨率场发射扫描电子显微镜(蔡司 Sigma 500)观察孔隙的形态、大小和分布。使用配备 J&M MSP200 显微光度计的蔡司显微镜对样品表面(10 × 10 μm)进行手工研磨和氩离子抛光(Leica EM 3X),以获得更好的图像质量。总有机碳 (TOC) 含量是用 Vario MICRO 立方体元素分析仪测定的。测试前,将页岩样本研磨至小于 200 目,并与稀 HCL 反应 72 小时以去除碳酸盐矿物。X 射线衍射(XRD)分析是使用 Dmax IIIa 衍射仪和铜 X 射线管(40 千伏,30 毫安)进行的。粘土矿物的分离和 XRD 测试按照中国石油行业标准(SY/T 5163-2018)进行。上述实验在南京大学表层地球化学教育部重点实验室进行,有效压力条件下的孔隙度和渗透率采用 PoroPerm-200 全自动孔隙度-渗透率仪进行测量。超压条件的实现主要依靠岩心保持系统,该系统以气体为加压介质,模拟水力压裂加压,从而实现超压条件下孔隙渗透率的在线测试。实验过程中,加压气源采用高压空气,测试采用高纯度氦气。为避免滑移效应对页岩样本渗透率的影响,在加压过程中保持驱动压力恒定。实验中设定的有效压力点在加压过程中保持不变。实验中设定的有效压力点分别为 2.41、4.48、6.55、8.62、10.69、12.07、20、30 和 50 兆帕。在每个压力点保持测量 30 分钟,稳定后测量相应压力下的孔隙度和渗透率。过压实验前,页岩样品需要用二氯甲烷和甲醇(93:7)萃取(72 小时),以去除样品中的残余油。根据层理厚度,结合岩芯和薄片观察,将页岩样品划分为层理页岩(2 厘米)。除块状页岩外,层状页岩和层状页岩的层理均较发育,表现为浅色层理和深色层理交替出现(图 3)。浅色层理的主要矿物为方解石,泥晶方解石和疏松方解石之间漂浮着少量石英、长石和白云石;深色层理的矿物成分较为复杂,以粘土矿物为主,还含有一些方解石、白云石、石英、长石和有机质(OM)。层状矿物成分的不同导致层状页岩和层状页岩中矿物分布的强烈异质性。表 2 列出了页岩样本的平均玻璃光泽反射率(Ro%)、总有机碳含量和矿物成分。本研究中的页岩样本热成熟度较低(0.78% Ro 至 1.12% Ro),总有机碳含量适中(2.85% 至 6.26 wt%,平均 3.51 wt%)。矿物包括方解石(7-69 wt%,平均 37 wt%)、粘土矿物(4-30 wt%,平均 22%)、石英(10-34 wt%,平均 19 wt%)、长石(4-22 wt%,平均 11 wt%)、白云石(3-22 wt%,平均 10 wt%)和少量黄铁矿(1-3 wt%,平均 15 MPa)。
{"title":"Influence of Rock Fabric on Physical Properties of Shale Oil Reservoir Under Effective Pressure Conditions","authors":"Yu Ma, Suping Yao, Ning Zhu, Huimin Liu, Junliang Li, Weiqing Wang","doi":"10.2113/2024/lithosphere_2023_338","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_338","url":null,"abstract":"The physical properties of shale oil reservoirs under overburden pressure are of great significance for reservoir prediction and evaluation during exploration and development. Based on core, thin section, and SEM observations, as well as test data such as XRD, TOC, and porosity and permeability under pressure conditions, this study systematically analyzes the variation of physical properties of different lithofacies shales in the Jiyang depression and the influence of rock fabric on the physical variation under pressure. The porosity and permeability of shale samples significantly decrease under pressure. According to the phased reduction in porosity and permeability, the pressurization process is divided into three pressure stages: low pressure (&lt;8 MPa), medium pressure (8–15 MPa), and high pressure (&gt;15 MPa). The reduction of porosity is fastest in the low-pressure stage and slowest in the medium-pressure stage. The reduction of permeability is fastest in the low-pressure stage and the slowest in the high-pressure stage. The rock fabric has a significant impact on porosity and permeability under pressure conditions. The permeability of laminated shale and bedded shale is higher than that of massive shale under pressure, and the permeability loss rate is lower than that of massive shales. Especially under lower pressure, the difference can be 10–20 times. In addition, the reduction rate of porosity and permeability under pressure is negatively correlated with felsic minerals content, which is positively correlated with carbonate minerals content and clay minerals content. The contribution of clay minerals to the porosity reduction rate is dominant, followed by carbonate minerals. The contribution of carbonate minerals to the permeability reduction rate is dominant, followed by clay minerals. The TOC content has no significant impact on the porosity and permeability of shales under pressure in the study due to the low maturity.With the change in global energy structure, shale oil and gas has become the core growth point of China’s oil and gas resources [1-4]. In the past decade, a series of important progresses have been made in the exploration and development of shale oil and gas in China, including breakthroughs in the exploration of shale oil in Junggar Basin, Ordos Basin, Jianghan Basin, Songliao Basin, and Bohai Bay Basin [5-8]. However, due to the heterogeneity of shales and the complexity of geological conditions in China, the prediction and evaluation of shale oil reservoirs still face many challenges [3, 9, 10].Many studies have shown that rock fabric, such as laminated structure and mineral composition, has a significant influence on the pore development and physical properties of shale oil reservoirs [10-14]. However, most of these studies were conducted under unpressurized conditions, and there are some errors with the formation conditions, which affect the prediction and evaluation of shale oil desserts. To recover the real physical ","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"23 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-04-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140932826","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Optimal Inversion Method Based on Joint Waveform Inversion and Least Squares Reverse Time Migration 基于联合波形反演和最小二乘反向时间迁移的最佳反演方法
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-04-30 DOI: 10.2113/2024/lithosphere_2023_361
Kai Zhang, Yipeng Xu, Zhenchun Li, Zilin He, Yiming Pan
Joint Waveform inversion (JWI) uses the results of reflection waveform inversion (RWI) as the initial model for full waveform inversion (FWI). Compared with the FWI, the JWI method can obtain more information about the structure of the subsurface medium. The reason is that the reflected waveform inversion can invert the long wavelength component in the middle and deep areas. In JWI, reflected waveform inversion is used to obtain the reflected wave information in the simulation record by demigration, which is computationally more expensive than FWI; the least squares reverse time migration (LSRTM) also obtains the reflected wave information in the simulated record by demigration. In order to effectively use the reflected wave information brought by the high computational amount of reverse migration in JWI, this paper proposes a simultaneous inversion method of JWI and LSRTM (JWI-LSRTM). This method can simultaneously perform an iterative update of the subsurface medium velocity of JWI and the migration imaging of LSRTM, which improves the calculation data utilization rate of each forward and inversion process. In the model test, the effectiveness of the method is proved.Reflected waveform inversion (RWI) is a technique that obtains reflected wave information from simulated records through reverse migration. It can invert long wavelength components in the middle and deep layers [1-3], but its computational cost is higher than full waveform inversion (FWI). FWI is a high-precision inversion method [4-6] with the potential to provide accurate models of subsurface media parameters. Joint waveform inversion (JWI) [7] leverages the results of reflection waveform inversion (RWI) as an initial model for FWI. Compared to FWI, JWI can retrieve more information about the subsurface media structure [8]. This is because RWI can obtain long wavelength components in the middle and deep layers [9, 10]. Ren (2019) proposed an adaptive JWI method that automatically switches between RWI and FWI by adjusting the weight factor with the number of iterations and allowable errors, without manually pausing the switch [11]. This approach addresses the limitations of traditional waveform inversion methods and improves the efficiency and accuracy of subsurface media modeling.The LSRTM [12, 13] is based on the Born approximation, and the reflection coefficient is solved by many iterations with the known background velocity. LSRTM is also used to obtain the reflected wave information from simulated records.In the above three methods, the simulated data and observed data are inverted using the generalized least squares method to obtain the corresponding gradient. However, the computational cost of forward modeling the wave equation and reverse migration is substantial, accounting for at least 90% of the total computation time in these three methods. As a result, the inversion cycle is often prolonged [14, 15]. Now there are a variety of programming techniques (MPI, Openmp, Open
联合波形反演(JWI)将反射波形反演(RWI)的结果作为全波形反演(FWI)的初始模型。与全波形反演相比,联合波形反演方法可以获得更多有关地下介质结构的信息。这是因为反射波形反演可以反演中深部的长波长分量。在 JWI 中,反射波形反演是通过反演来获取模拟记录中的反射波信息,其计算成本比 FWI 高;最小二乘反向时间迁移(LSRTM)也是通过反演来获取模拟记录中的反射波信息。为了有效利用 JWI 反演计算量大带来的反射波信息,本文提出了一种 JWI 和 LSRTM(JWI-LSRTM)同步反演方法。该方法可同时对 JWI 的地下介质速度和 LSRTM 的迁移成像进行迭代更新,提高了各正演和反演过程的计算数据利用率。反射波反演(RWI)是一种通过反向迁移从模拟记录中获取反射波信息的技术。它可以反演中深层的长波长成分[1-3],但其计算成本高于全波形反演(FWI)。全波形反演是一种高精度反演方法[4-6],有可能提供地下介质参数的精确模型。联合波形反演(JWI)[7] 利用反射波形反演(RWI)的结果作为 FWI 的初始模型。与 FWI 相比,JWI 可以获取更多有关地下介质结构的信息[8]。这是因为 RWI 可以获得中深层的长波长成分[9, 10]。Ren(2019)提出了一种自适应 JWI 方法,该方法通过调整权重系数,根据迭代次数和允许误差在 RWI 和 FWI 之间自动切换,无需手动暂停切换[11]。这种方法解决了传统波形反演方法的局限性,提高了地下介质建模的效率和精度。LSRTM[12, 13]基于 Born 近似,在已知背景速度的情况下,通过多次迭代求解反射系数。上述三种方法都是利用广义最小二乘法对模拟数据和观测数据进行反演,得到相应的梯度。然而,波方程正向建模和反向迁移的计算成本很高,至少占这三种方法总计算时间的 90%。因此,反演周期往往会延长[14, 15]。现在有多种编程技术(MPI、Openmp、Openacc、GPU)或优化方法(L-BFGS 法、伪牛顿方向校正混合共轭梯度法等)可以缩短这一周期,还有一些方法可以提高计算数据的利用率(FWI 成像、RWI 和 LSRTM 联合反演)[16, 17]。为了更有效地利用 JWI 高计算解译带来的反射波信息和波方程正演模拟数据,本文提出了一种 JWI 和 LSRTM 同步反演方法。该方法可以同时更新 JWI 中地下介质的速度和 LSRTM 的迁移剖面,提高了计算数据在正演和反演过程中的利用率。本文以较容易实现的声学为例,通过模型试验和噪声数据证明了所提方法的有效性。以下是声波的波方程和反演方程:其中 x 和 z 为空间坐标,t 为时间,p0 为背景场,p 为干扰场,vp 为 p 波速度。式(2)中,dcalxr,t;xs 为模拟数据,dcalrefxr,t;xs 为散射数据,dobsxr,t;xs 为观测数据,m 为反射系数,其表达式为:mimage(x)=Δv(x)vp(x),散射数据可通过式(2)计算。其邻接方程和梯度公式可通过邻接法获得。其中,T 为最大记录时间 dcal,m 为次表层介质参数,分别为模拟数据(由方程 dobs(1) 求得)和观测数据。
{"title":"Optimal Inversion Method Based on Joint Waveform Inversion and Least Squares Reverse Time Migration","authors":"Kai Zhang, Yipeng Xu, Zhenchun Li, Zilin He, Yiming Pan","doi":"10.2113/2024/lithosphere_2023_361","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_361","url":null,"abstract":"Joint Waveform inversion (JWI) uses the results of reflection waveform inversion (RWI) as the initial model for full waveform inversion (FWI). Compared with the FWI, the JWI method can obtain more information about the structure of the subsurface medium. The reason is that the reflected waveform inversion can invert the long wavelength component in the middle and deep areas. In JWI, reflected waveform inversion is used to obtain the reflected wave information in the simulation record by demigration, which is computationally more expensive than FWI; the least squares reverse time migration (LSRTM) also obtains the reflected wave information in the simulated record by demigration. In order to effectively use the reflected wave information brought by the high computational amount of reverse migration in JWI, this paper proposes a simultaneous inversion method of JWI and LSRTM (JWI-LSRTM). This method can simultaneously perform an iterative update of the subsurface medium velocity of JWI and the migration imaging of LSRTM, which improves the calculation data utilization rate of each forward and inversion process. In the model test, the effectiveness of the method is proved.Reflected waveform inversion (RWI) is a technique that obtains reflected wave information from simulated records through reverse migration. It can invert long wavelength components in the middle and deep layers [1-3], but its computational cost is higher than full waveform inversion (FWI). FWI is a high-precision inversion method [4-6] with the potential to provide accurate models of subsurface media parameters. Joint waveform inversion (JWI) [7] leverages the results of reflection waveform inversion (RWI) as an initial model for FWI. Compared to FWI, JWI can retrieve more information about the subsurface media structure [8]. This is because RWI can obtain long wavelength components in the middle and deep layers [9, 10]. Ren (2019) proposed an adaptive JWI method that automatically switches between RWI and FWI by adjusting the weight factor with the number of iterations and allowable errors, without manually pausing the switch [11]. This approach addresses the limitations of traditional waveform inversion methods and improves the efficiency and accuracy of subsurface media modeling.The LSRTM [12, 13] is based on the Born approximation, and the reflection coefficient is solved by many iterations with the known background velocity. LSRTM is also used to obtain the reflected wave information from simulated records.In the above three methods, the simulated data and observed data are inverted using the generalized least squares method to obtain the corresponding gradient. However, the computational cost of forward modeling the wave equation and reverse migration is substantial, accounting for at least 90% of the total computation time in these three methods. As a result, the inversion cycle is often prolonged [14, 15]. Now there are a variety of programming techniques (MPI, Openmp, Open","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"1 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-04-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141253434","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Oil Production Rate Forecasting by SA-LSTM Model in Tight Reservoirs 利用 SA-LSTM 模型预测致密油藏的石油生产率
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_197
Denghui He, Yaguang Qu, Guanglong Sheng, Bin Wang, Xu Yan, Zhen Tao, Meng Lei
The accurate forecasting of oil field production rate is a crucial indicator for each oil field’s successful development, but due to the complicated reservoir conditions and unknown underground environment, the high accuracy of production rate forecasting is a popular challenge. To find a low time consumption and high accuracy method for forecasting production rate, the current paper proposes a hybrid model, Simulated Annealing Long Short-Term Memory network (SA-LSTM), based on the daily oil production rate of tight reservoirs with the in situ data of injection and production rates in fractures. Furthermore, forecasting results are compared with the numerical simulation model output. The LSTM can effectively learn time-sequence problems, while SA can optimize the hyperparameters (learning rate, batch size, and decay rate) in LSTM to achieve higher accuracy. By conducting the optimized hyperparameters into the LSTM model, the daily oil production rate can be forecasted well. After training and predicting on existing production data, three different methods were used to forecast daily oil production for the next 300 days. The results were then validated using numerical simulations to compare the forecasting of LSTM and SA-LSTM. The results show that SA-LSTM can more efficiently and accurately predict daily oil production. The fitting accuracies of the three methods are as follows: numerical reservoir simulation (96.2%), LSTM (98.1%), and SA-LSTM (98.7%). The effectiveness of SA-LSTM in production rate is particularly outstanding. Using the same SA-LSTM model, we input the daily oil production data of twenty oil wells in the same block and make production prediction, and the effect is remarkable.The forecasting of the oil and gas production rate is one of the most important and effective evaluation indicators for measuring the success of reservoir development, and it plays a crucial role in dynamically predicting the oil and gas production rate during the development process. However, due to the geological factors of the reservoir and the construction factors during the development process, oil and gas production rate forecasting has become more complex, and the dynamic characteristics cannot be well described, resulting in the subsequent production rate forecasting being affected [1-4]. There are various methods for forecasting oil and gas production rate, including the Arps decline method explored by the production rate decline law, the analytical model method based on the permeability law and the material balance equation, and the numerical simulation method based on the geological model constructed using geological data [5, 6]. Conventional oil and gas production rate dynamic forecasting generally uses numerical simulation methods, which can comprehensively consider various geological factors, wellbore interference, and the impact of multiphase flow on oil and gas well production rate. However, for unconventional reservoirs such as tight oil res
油田产量的准确预测是每个油田能否成功开发的关键指标,但由于油藏条件复杂、地下环境未知,高精度的产量预测是一个普遍面临的挑战。为了找到一种耗时少、精度高的生产率预测方法,本文提出了一种混合模型--模拟退火长短期记忆网络(Simulated Annealing Long Short-Term Memory network,SA-LSTM)。此外,还将预测结果与数值模拟模型输出结果进行了比较。LSTM 可以有效地学习时间序列问题,而 SA 可以优化 LSTM 中的超参数(学习率、批量大小和衰减率),以达到更高的精度。通过将优化后的超参数导入 LSTM 模型,可以很好地预测石油日产量。在对现有生产数据进行训练和预测后,使用三种不同的方法预测了未来 300 天的石油日产量。然后使用数值模拟对结果进行验证,比较 LSTM 和 SA-LSTM 的预测效果。结果表明,SA-LSTM 可以更高效、更准确地预测石油日产量。三种方法的拟合精度如下:油藏数值模拟(96.2%)、LSTM(98.1%)和 SA-LSTM(98.7%)。SA-LSTM 在生产率方面的效果尤为突出。油气产量预测是衡量油藏开发成功与否最重要、最有效的评价指标之一,对开发过程中油气产量的动态预测起着至关重要的作用。然而,由于油气藏的地质因素和开发过程中的施工因素,油气产量预测变得更加复杂,动态特征不能得到很好的描述,导致后续的产量预测受到影响[1-4]。油气产量预测的方法多种多样,包括由产量递减规律探索出的 Arps 递减法、基于渗透率规律和物料平衡方程的解析模型法、基于利用地质资料构建的地质模型的数值模拟法等[5, 6]。传统的油气产量动态预测一般采用数值模拟方法,可以综合考虑各种地质因素、井筒干扰、多相流对油气井产量的影响等。然而,对于致密油藏等非常规储层,由于储层的岩性、横向连通性、垂直分布、微观结构快速变化、强异质性等问题,在准确拟合生产率历史和预测未来生产率方面存在挑战[7, 8]。随着人工智能技术的不断发展,人工智能逐渐在油气藏开发过程中发挥着非常重要的作用,使得利用机器学习技术解决油气藏开发中的一些复杂问题成为可能[9-12]。递归神经网络(RNN)被提出用于油气产量预测,同时其局限性也暴露出来。RNN 只能将上一序列的生产率数据与当前数据连接起来。当生产率数据与预测数据之间的间隔很小时,RNN 具有很好的学习能力和匹配性。但是,当引入复杂情况时,如在生产率过程中进行注水和能量补充阶段,导致生产率量发生变化时,由于这些相关信息,神经单元之间的间隔变得非常大,RNN 就会失去学习能力,无法进行准确有效的预测。RNN 的本质容易受到短序列生产率数据的影响。一旦生产率数据的序列增加,RNN 就很难将信息从早期的生产率数据转移到后期的生产率数据。在确定了 RNN 在长序列生产率预测中的不准确性后,有人提出了 RNN 的一种变体,称为长短时记忆(LSTM)。LSTM 最初由 Hochreiter 和 Schmidhuber 提出[13, 14],后经 Graves [15]改进和推广,使 LSTM 在各种应用中得到广泛应用。 从图 6 可以看出,随着迭代次数的增加,模型收敛速度在早期阶段较慢。然而,当迭代次数超过 50 次时,模型收敛速度就会增加,当迭代次数达到 150 次时,就会接近完全收敛。通过研究学习率值,可以推断出在优化范围内,当学习率在 0.0006 和 0.0008 之间时,模型收敛速度较快。超出这一范围的学习率被排除在外,因为它们在优化过程中的损失值一直较高。随着 SA 温度的降低,学习率值保持在 0.0002-0.0008 的范围内。随着迭代的进行,计算了不同学习率下的损失值。图 7 显示,损失值的变化随着图 6 中学习率的增加而减小,表明两者之间存在反比关系。因此,经过全局优化后发现,当学习率为 0.0007 时,相应的迭代次数为 232 次,最小损失值为 0.00167。为了找到衰减超参数的最佳值,通过计算模型的损失测试了不同的值。随着迭代的继续,模型在第 100 次迭代时几乎完全收敛(图 8)。根据 SA 选择有效数据的原则,衰变超参数的优化主要集中在 0.9 至 0.94 之间。在 0.9-0.94 的范围内,损失值随着衰减值的增大而减小。当学习率为 0.0007 时,能使模型损失最小化的最佳衰减值为 0.924,相应的损失值为 0.00161。在优化了学习率和衰减的超参数后,在保持其他超参数不变的情况下,通过在可控范围内随机扰动来测试不同的批量大小。表现最好的批次规模被保留下来。从图 10 中可以看出,模型在第 150 次迭代时完全收敛,最佳批量大小是均匀分布的。起初,在 SA 的高温阶段,由于收敛速度较慢,较大的批次规模作为新的解决方案是可以接受的。在确定了上述学习率和衰减的最佳值后(图 11),算法根据最小损失值接受了 15 的全局最佳批量大小。在将三个优化的超参数应用于 LSTM 模型后,我们将优化后的模型与未优化的模型进行了比较。为了评估两个模型的性能,我们比较了它们的训练和验证损失以及预测准确率。优化前后的训练和实际损失值如图 12 所示。图中左侧显示的是超参数优化前 LSTM 模型的训练损失结果。可以看出,随着迭代次数的增加,损失值逐渐趋近于 0。图中右侧显示的是超参数优化后的 LSTM 训练损失结果。一般来说,验证集损失(真值损失)低于训练集损失。表 3(a)和表 3(b)显示了模型优化前后对训练数据和预测数据的损失水平。表 3(a) 和表 3(b) 显示了优化前后模型在训练和预测数据方面的损失水平。具体来说,优化后的 LSTM 模型的训练损失比之前小得多,这导致后续预测过程中的损失更小。根据损失值的差异,可以得出结论:超参数优化后的 LSTM 模型能够更有效地预测日采油率。为了预测未来的采油率,我们根据现有的采油率数据对油藏数值模拟模型(图 13)进行了修改,以确保多口虚拟垂直井的日采油率等同于单口压裂水平井的日采油率,且注采同步。 经过对比和验证,可以确定 SA-LSTM 模型在拟合和预测致密油藏的日产油量方面具有很高的精度。经过大量的调试和属性修改,常规油藏的数值模拟精度达到 96.2%。但这一过程需要大量的时间。SA 搜索可以有效地剔除不良参数,为 LSTM 模型找到最合适的参数。该方法操作简单,搜索时间短,可获得全局最优超参数。经过比较和验证,可以确定 SA-LSTM 模型在拟合和预测致密油藏日产油量方面表现出较高的精度。本研究得到了中海油(中国)有限公司 "海洋油田显著提高石油采收率关键技术 "重大项目(KJGG2021-0501)的支持。
{"title":"Oil Production Rate Forecasting by SA-LSTM Model in Tight Reservoirs","authors":"Denghui He, Yaguang Qu, Guanglong Sheng, Bin Wang, Xu Yan, Zhen Tao, Meng Lei","doi":"10.2113/2024/lithosphere_2023_197","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_197","url":null,"abstract":"The accurate forecasting of oil field production rate is a crucial indicator for each oil field’s successful development, but due to the complicated reservoir conditions and unknown underground environment, the high accuracy of production rate forecasting is a popular challenge. To find a low time consumption and high accuracy method for forecasting production rate, the current paper proposes a hybrid model, Simulated Annealing Long Short-Term Memory network (SA-LSTM), based on the daily oil production rate of tight reservoirs with the in situ data of injection and production rates in fractures. Furthermore, forecasting results are compared with the numerical simulation model output. The LSTM can effectively learn time-sequence problems, while SA can optimize the hyperparameters (learning rate, batch size, and decay rate) in LSTM to achieve higher accuracy. By conducting the optimized hyperparameters into the LSTM model, the daily oil production rate can be forecasted well. After training and predicting on existing production data, three different methods were used to forecast daily oil production for the next 300 days. The results were then validated using numerical simulations to compare the forecasting of LSTM and SA-LSTM. The results show that SA-LSTM can more efficiently and accurately predict daily oil production. The fitting accuracies of the three methods are as follows: numerical reservoir simulation (96.2%), LSTM (98.1%), and SA-LSTM (98.7%). The effectiveness of SA-LSTM in production rate is particularly outstanding. Using the same SA-LSTM model, we input the daily oil production data of twenty oil wells in the same block and make production prediction, and the effect is remarkable.The forecasting of the oil and gas production rate is one of the most important and effective evaluation indicators for measuring the success of reservoir development, and it plays a crucial role in dynamically predicting the oil and gas production rate during the development process. However, due to the geological factors of the reservoir and the construction factors during the development process, oil and gas production rate forecasting has become more complex, and the dynamic characteristics cannot be well described, resulting in the subsequent production rate forecasting being affected [1-4]. There are various methods for forecasting oil and gas production rate, including the Arps decline method explored by the production rate decline law, the analytical model method based on the permeability law and the material balance equation, and the numerical simulation method based on the geological model constructed using geological data [5, 6]. Conventional oil and gas production rate dynamic forecasting generally uses numerical simulation methods, which can comprehensively consider various geological factors, wellbore interference, and the impact of multiphase flow on oil and gas well production rate. However, for unconventional reservoirs such as tight oil res","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"21 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139462812","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Chronological, Petrogenetic, and Tectonic Significance of Paleoproterozoic Continental Crust within the North China Craton during the Global Tectono-Magmatic Lull 全球构造-岩浆低潮期华北克拉通内古新生代大陆地壳的年代学、岩石学和构造意义
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_234
Heng Liu, Lei Liu, Dexian Zhang, Inkyeong Moon, M. Santosh, Yanyan Zhou, Tianyang Hu, Shisheng Kang
The 2.45–2.20 Ga period during the early Paleoproterozoic era is considered to have witnessed a global “Tectono-Magmatic Lull (TML)” and thus marks a relatively quiescent period. Our study unveils a 2.45–2.20 Ga magmatic suite from the Xiong’ershan area in the southern North China Craton, offering some key constraints on localized active tectonics during the TML. Zircon U-Pb dating shows Paleoproterozoic ages for the meta-basalt (2.31, 2.28 Ga), Na-rich meta-andesite (~2.33 Ga), tonalite-trondhjemite-granodiorite (TTG) gneisses (2.36, 2.30 Ga), K-rich granodiorite (~2.29 Ga), and monzogranite (2.33, 2.27 Ga). The meta-basalts geochemically and petrographically belong to calc-alkaline basalts and show distinctive Nb, Ta, and Ti contents and primitive mantle normalized patterns from different places in the Xiong’ershan area. Combined with their enriched εHf(t) values, the magmas were derived from subduction-related enriched mantle sources within a convergent plate boundary. The meta-andesites display high MgO content (average 4.5 wt%) and Mg# (44–57), strongly fractionated rare-earth pattern, calc-alkaline affinity, and negative Nb, Ta, and Ti anomalies. The TTG gneisses are of high SiO2 type (>62 wt%), high (La/Yb)N (17.5, 39.2), and Sr/Y (50.2, 104.3) and mostly display positive Eu anomalies and high-pressure type. Zircons from these rocks show a relatively narrow range of δ18O isotope values (5.35‰, 6.79‰) with εHf(t) isotope characteristics (−9.3, −3.3), suggesting derivation from partial melting of a thickened mafic lower crust. The youngest K-rich granodiorite and monzogranite show high K2O/Na2O ratios (0.65, 2.45). Variable molar ratio Al2O3/(CaO+Na2O+K2O) (A/CNK) and low zircon εHf(t) values suggest that the K-rich granitoids formed from the partial melting of different levels of crust. The presence of meta-basalt to andesite assemblages and diverse intermediate to felsic magmatic rocks implies magmatic activity within a convergent plate boundary tectonic environment with potential influence from plume-triggered extensional processes, supported by evidence of slab rollback and upwelling of mantle material.After 2.5 Ga, the globe has witnessed a relatively quiescent period for over 200 million years in terms of active plate tectonics, referred to as the “Tectono-Magmatic Lull (TML, 2.45–2.20 Ga),” with no significant continental crust growth or major orogenesis [1-6]. In this regard of the geological processes of TML, Silver and Behn [7] suggested stagnation of the global subduction system leading to a decrease in volcanic activity and continental growth, Condie et al. [1] referred to unusual period as a crustal age gap, while Spencer et al. [8] referred to it as a TML. At the Archean/Proterozoic boundary (2.50 Ga), the Earth underwent significant episodic evolution and transformation in the early Paleoproterozoic period [1, 9, 10].There are controversial opinions about the tectonic evolution of the Precambrian era. Cawood et al. [11], Palin
古近代早期的2.45-2.20 Ga时期被认为是全球 "构造-岩浆低潮期(Tectono-Magmatic Lull,TML)",因而是一个相对静止的时期。我们的研究揭示了华北克拉通南部熊耳山地区的2.45-2.20 Ga岩浆岩组,为TML时期局部活跃构造提供了一些关键的制约因素。锆石U-Pb定年显示了元玄武岩(2.31, 2.28 Ga)、富Na元安山岩(〜2.33 Ga)、黑云母-特长闪长岩-花岗闪长岩(TTG)片麻岩(2.36, 2.30 Ga)、富K花岗闪长岩(〜2.29 Ga)和独居岩(2.33, 2.27 Ga)的古新生代年龄。元玄武岩在地球化学和岩石学上属于钙碱性玄武岩,在熊耳山地区的不同地方表现出不同的Nb、Ta、Ti含量和原始地幔正常化模式。结合其富集的εHf(t)值,这些岩浆来源于汇聚板块边界内与俯冲有关的富集地幔源。元安山岩显示出较高的氧化镁含量(平均 4.5 wt%)和镁#(44-57),强烈的稀土分馏模式,钙碱性亲和性,以及负的 Nb、Ta 和 Ti 异常。TTG片麻岩具有高SiO2类型(>62 wt%)、高(La/Yb)N(17.5,39.2)和Sr/Y(50.2,104.3),并大多显示出正Eu异常和高压类型。这些岩石中的锆石显示出相对较窄的δ18O同位素值范围(5.35‰,6.79‰)和εHf(t)同位素特征(-9.3,-3.3),表明这些锆石来源于增厚的岩浆岩下地壳的部分熔融。最年轻的富钾花岗闪长岩和单斜花岗岩显示出较高的 K2O/Na2O 比值(0.65,2.45)。Al2O3/(CaO+Na2O+K2O)摩尔比(A/CNK)的变化和锆石εHf(t)值的降低表明,富钾花岗岩是由不同层次的地壳部分熔融形成的。元玄武岩到安山岩组合以及多种多样的中间岩浆岩到长英岩的存在,意味着岩浆活动是在板块交汇的构造环境中进行的,并可能受到羽状触发的延伸过程的影响,板块后退和地幔物质上涌的证据也支持了这一点。在 2.5 Ga 之后,全球板块构造活动经历了 2 亿多年的相对平静期,被称为 "构造-岩浆平静期(Tectono-Magmatic Lull,TML,2.45-2.20 Ga)",没有明显的大陆地壳增长或重大造山运动[1-6]。关于 TML 的地质过程,Silver 和 Behn[7]认为全球俯冲系统的停滞导致火山活动和大陆生长的减少,Condie 等[1]将不寻常时期称为地壳年龄间隙,而 Spencer 等[8]则将其称为 TML。在Archean/Proterozoic边界(2.50 Ga),地球在Paleoproterozoic早期经历了显著的偶发演化和转变[1, 9, 10]。Cawood等[11]、Palin和Santosh[12]以及Tang等[4]认为,全球板块构造体系开始运行于中新世至新元古代。此外,Campbell和Griffiths[13]提出,地幔势温在TML期间显著下降。来自花岗岩和碎屑沉积物的锆石U-Pb年龄分布表明,在TML期间(2.365, 2.235 Ma)有一个异常强烈的最小值,尽管没有发现公认的年龄差距[2]。然而,一些研究报告称火成岩和碎屑锆石位于TML窗口内,例如在巴西米纳斯造山带、加拿大阿罗史密斯造山带和华北克拉通(NCC)南部发现的火成岩和碎屑锆石[14]。值得注意的是,约 2.45-2.20 Ga 这一引人入胜的时期似乎在促进从 Archean 到古近代超大陆的过渡中发挥了关键作用。一些研究发现,在 2.45 至 2.20 Ga 期间,至少有二十四个陨石坑或陨石块发生了构造岩浆活动,包括造山活动、被动边缘以及岩浆活动,如绿岩带和黑云母-透闪石-花岗闪长岩(TTG)。这些岩浆活动在从阿基坦构造体系的分散陨石坑向现代板块构造体系的转变过程中起到了至关重要的作用[8]。以往的研究表明,板块构造类似于现代的高角度深俯冲,在大约 3.0 Ga 时就已经开始,并在新元古代到古近代早期在全球范围内盛行[18, 19],这可以从全球金属生成事件中得到证明[6]。Condie等人[1]汇编了来自碎屑沉积岩和岩浆岩组的全球锆石U-Pb年龄,提供了2.45Ga时期全球岩浆活动持续约200-250Ma的证据。俯冲带和断裂带的出现是板块构造的一个显著特征[16]。
{"title":"Chronological, Petrogenetic, and Tectonic Significance of Paleoproterozoic Continental Crust within the North China Craton during the Global Tectono-Magmatic Lull","authors":"Heng Liu, Lei Liu, Dexian Zhang, Inkyeong Moon, M. Santosh, Yanyan Zhou, Tianyang Hu, Shisheng Kang","doi":"10.2113/2024/lithosphere_2023_234","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_234","url":null,"abstract":"The 2.45–2.20 Ga period during the early Paleoproterozoic era is considered to have witnessed a global “Tectono-Magmatic Lull (TML)” and thus marks a relatively quiescent period. Our study unveils a 2.45–2.20 Ga magmatic suite from the Xiong’ershan area in the southern North China Craton, offering some key constraints on localized active tectonics during the TML. Zircon U-Pb dating shows Paleoproterozoic ages for the meta-basalt (2.31, 2.28 Ga), Na-rich meta-andesite (~2.33 Ga), tonalite-trondhjemite-granodiorite (TTG) gneisses (2.36, 2.30 Ga), K-rich granodiorite (~2.29 Ga), and monzogranite (2.33, 2.27 Ga). The meta-basalts geochemically and petrographically belong to calc-alkaline basalts and show distinctive Nb, Ta, and Ti contents and primitive mantle normalized patterns from different places in the Xiong’ershan area. Combined with their enriched εHf(t) values, the magmas were derived from subduction-related enriched mantle sources within a convergent plate boundary. The meta-andesites display high MgO content (average 4.5 wt%) and Mg# (44–57), strongly fractionated rare-earth pattern, calc-alkaline affinity, and negative Nb, Ta, and Ti anomalies. The TTG gneisses are of high SiO2 type (&gt;62 wt%), high (La/Yb)N (17.5, 39.2), and Sr/Y (50.2, 104.3) and mostly display positive Eu anomalies and high-pressure type. Zircons from these rocks show a relatively narrow range of δ18O isotope values (5.35‰, 6.79‰) with εHf(t) isotope characteristics (−9.3, −3.3), suggesting derivation from partial melting of a thickened mafic lower crust. The youngest K-rich granodiorite and monzogranite show high K2O/Na2O ratios (0.65, 2.45). Variable molar ratio Al2O3/(CaO+Na2O+K2O) (A/CNK) and low zircon εHf(t) values suggest that the K-rich granitoids formed from the partial melting of different levels of crust. The presence of meta-basalt to andesite assemblages and diverse intermediate to felsic magmatic rocks implies magmatic activity within a convergent plate boundary tectonic environment with potential influence from plume-triggered extensional processes, supported by evidence of slab rollback and upwelling of mantle material.After 2.5 Ga, the globe has witnessed a relatively quiescent period for over 200 million years in terms of active plate tectonics, referred to as the “Tectono-Magmatic Lull (TML, 2.45–2.20 Ga),” with no significant continental crust growth or major orogenesis [1-6]. In this regard of the geological processes of TML, Silver and Behn [7] suggested stagnation of the global subduction system leading to a decrease in volcanic activity and continental growth, Condie et al. [1] referred to unusual period as a crustal age gap, while Spencer et al. [8] referred to it as a TML. At the Archean/Proterozoic boundary (2.50 Ga), the Earth underwent significant episodic evolution and transformation in the early Paleoproterozoic period [1, 9, 10].There are controversial opinions about the tectonic evolution of the Precambrian era. Cawood et al. [11], Palin","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"41 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140072420","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Influence of Hydraulic Conditions on Seepage Characteristics of Loose Sandstone 水力条件对松散砂岩渗流特性的影响
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_275
Mingtao Jia, Quan Jiang, Qi Xu, Xuebin Su
To investigate the impact of hydraulic conditions on the seepage characteristics of loose sandstone, this study employed optimized methods to prepare loose sandstone samples. Subsequently, seepage experiments were conducted under different injection pressures, flow rates, and flow volumes. The permeability, porosity, particle size distribution, and other parameters of the rock samples were obtained. By analyzing the response of seepage characteristics to pore and particle size characteristics, the influence of different hydraulic conditions on the seepage characteristics of loose sandstone was explored. The results indicated that improvements in the parameters of hydraulic conditions had different effects on various rock samples. For rock samples with developed seepage channels, increasing the value of each hydraulic condition parameter could expand the channels and discharge particles, and improve permeability. For rock samples with a larger number of small pores, increasing each hydraulic condition parameter caused particles to crack under pressure, drove particles to block holes, and thus reduced permeability. In this experiment, the permeability parameter had a significant positive response to the proportion of pores larger than 0.1 µm and a significant negative response to the proportion of particles smaller than 150 µm.In the fields of oil extraction and solution mining, loose sandstone is a common resource-bearing rock mass, and the seepage characteristics of this type of rock are directly related to the process design and efficiency of resource extraction [1-5]. During the processes of extraction and injection in rock formations, problems often occur, such as increased rock permeability leading to imbalanced extraction and injection or decreased permeability leading to low mining efficiency [6-8]. In the above situations, a common countermeasure in industrial practice is to use agents or equipment to adjust the permeability of rock strata to improve production. For example, in oil exploitation, dispersed gel particles with a certain particle size are used to temporarily plug the high permeability area to better drive the production of the reservoir [9, 10], and unblocking agents are used in situ to unblock clogged channels to improve the efficiency of uranium leaching [11-13]. Implementation of the above methods provides immediate improvement in the seepage effect, and the cost can be controlled, effectively solving problems related to abnormal seepage in rock formations. However, it is important to note that while the solutions for abnormal seepage are relatively effective, they are still confined to a reactive, postevent stage. The time and costs associated with these solutions continue to impede the enhancement of production efficiency. The objective of this study is to address these issues by starting from the mechanisms of permeability changes. It delves into the influence of hydraulic conditions on the seepage characteristics of loo
孔隙分布特征参数可以有效揭示岩石的孔隙发育与渗流情况,并据此推导出渗流特征的变化特征与机理[26-29],而颗粒特征的分析则侧重于粒度分布对孔隙结构的影响以及游离颗粒对渗流通道的影响[29, 30]。研究表明,孔隙特征和颗粒特征基本上决定了岩石的渗透性,而且这两种特征是相互关联的[31]。因此,将对岩石孔隙分布和粒径分布特征的研究与从微观角度进行的渗流实验相结合,可以充分揭示松散砂岩的水力条件与渗流特征之间的关系。然而,在实验室中,由于松散砂岩独特的岩性特征,对其渗流特征的微观研究明显不足。松散砂岩的内聚力主要来自非饱和岩粒以及细粒土之间的毛细作用 [32,33]。这种作用力稳定性差,在饱和状态下就会失效,因此松散砂岩很难切割成型,岩石颗粒在水中容易崩解扩散,破坏渗水管道。此外,天然松散砂岩由于颗粒分布随机性大,容易产生较大的实验误差。因此,在实验室中对松散砂岩进行的渗流实验很少。一些学者采用数值模拟或大型岩样整体渗流的方法来获取渗流特性的变化[34-38]。但这类方法只能用于评价渗透率变化等参数,无法直观地获得实验前后的孔隙度、颗粒分布等关键参数,难以解释松散砂岩渗流特性变化的原理。近年来,一些研究采用侧向包裹或特殊固定装置保护松散砂岩样品[39-41],成功实现了小规模渗流实验,其岩样保护方法值得借鉴。但从相关研究的数据来看,这些实验的成功得益于砂岩样品的良好完整性。但这些实验后,岩样的渗出端出现了明显的颗粒剥落现象,导致实验误差较大,难以进行更多的实验。为了研究不同水力条件影响松散砂岩孔隙度和颗粒特征,进而改变其渗流特性的规律和机理,设计了一种特殊的优化松散砂岩制备方法,用于制备岩石样品。在不同注入压力、流速和流量条件下,对岩石样本进行了渗流实验。实验前后通过渗流仪、核磁共振波谱仪、粒度分布仪等仪器获得了渗透率、孔隙度、粒度分布等关键参数,并综合分析了水力条件、渗流特征、孔隙度分布特征、粒度分布特征的响应和原理。最后,总结了不同水力条件对松散砂岩渗流特征的影响关系和影响原理,为工程实践提供了理论参考,有助于提高含松散砂岩资源渗流相关作业的效率。本实验中,松散砂岩样品 A 和 B 分别取自新疆 A 型铀矿地下 400 米地层和内蒙古 B 型铀矿地下 500 米地层。两种岩石样品的矿物相显微鉴定结果见表 1。样品 A 的干密度为 1.69 g/cm³,样品 B 的干密度为 1.62 g/cm³。渗流实验采用海安县石油科学研究仪器有限公司生产的 HKY-1 长芯渗流监测仪进行,如图 2(a)所示。
{"title":"Influence of Hydraulic Conditions on Seepage Characteristics of Loose Sandstone","authors":"Mingtao Jia, Quan Jiang, Qi Xu, Xuebin Su","doi":"10.2113/2024/lithosphere_2023_275","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_275","url":null,"abstract":"To investigate the impact of hydraulic conditions on the seepage characteristics of loose sandstone, this study employed optimized methods to prepare loose sandstone samples. Subsequently, seepage experiments were conducted under different injection pressures, flow rates, and flow volumes. The permeability, porosity, particle size distribution, and other parameters of the rock samples were obtained. By analyzing the response of seepage characteristics to pore and particle size characteristics, the influence of different hydraulic conditions on the seepage characteristics of loose sandstone was explored. The results indicated that improvements in the parameters of hydraulic conditions had different effects on various rock samples. For rock samples with developed seepage channels, increasing the value of each hydraulic condition parameter could expand the channels and discharge particles, and improve permeability. For rock samples with a larger number of small pores, increasing each hydraulic condition parameter caused particles to crack under pressure, drove particles to block holes, and thus reduced permeability. In this experiment, the permeability parameter had a significant positive response to the proportion of pores larger than 0.1 µm and a significant negative response to the proportion of particles smaller than 150 µm.In the fields of oil extraction and solution mining, loose sandstone is a common resource-bearing rock mass, and the seepage characteristics of this type of rock are directly related to the process design and efficiency of resource extraction [1-5]. During the processes of extraction and injection in rock formations, problems often occur, such as increased rock permeability leading to imbalanced extraction and injection or decreased permeability leading to low mining efficiency [6-8]. In the above situations, a common countermeasure in industrial practice is to use agents or equipment to adjust the permeability of rock strata to improve production. For example, in oil exploitation, dispersed gel particles with a certain particle size are used to temporarily plug the high permeability area to better drive the production of the reservoir [9, 10], and unblocking agents are used in situ to unblock clogged channels to improve the efficiency of uranium leaching [11-13]. Implementation of the above methods provides immediate improvement in the seepage effect, and the cost can be controlled, effectively solving problems related to abnormal seepage in rock formations. However, it is important to note that while the solutions for abnormal seepage are relatively effective, they are still confined to a reactive, postevent stage. The time and costs associated with these solutions continue to impede the enhancement of production efficiency. The objective of this study is to address these issues by starting from the mechanisms of permeability changes. It delves into the influence of hydraulic conditions on the seepage characteristics of loo","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"174 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139589630","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Petrogenesis of the Early Paleozoic Dioritic–Granitic Magmatism in the Eastern North Qilian Orogen, NW China: Implications for Tethyan Tectonic Evolution 中国西北部祁连山东麓早古生代闪长岩-花岗岩岩浆岩的成岩作用:对哲罗纪构造演化的影响
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_297
Fei Xue, Fan Yang, Weidong Ren, M. Santosh, Zesheng Qian, Yin Huang, Zijian Tan
The North Qilian Orogen witnessed the opening, subduction, and closure of the Proto-Tethys Qilian Ocean and the post-subduction of multiple exhumation events from Late Neoproterozoic to Early Paleozoic. The Early Paleozoic dioritic–granitic magmatic suites, prominently exposed in the eastern North Qilian Orogen, offer valuable insights into the evolution of the Proto-Tethys Ocean. However, their petrogenesis, magma source, and tectonic evolution remain controversial. Here, we investigate the Leigongshan, Zhigou, and Dalongcun intrusions and present geochronological, geochemical, and isotopic data, aiming to refine the comprehension of their timing and petrogenesis, which will contribute to understanding the tectonic evolution of the Proto-Tethys Ocean. Zircon U-Pb dating reveals mean ages of 471–427 Ma for these intrusions, consistent with compiled formation ages of dioritic–granitic intrusions in the eastern North Qilian Orogen, indicating close temporal links with the tectonic evolution of the Proto-Tethys Ocean during the Early Paleozoic. The studied magmatic rocks could be categorized into two major types: granitoids and diorites. The granitoids are majorly I-type granitoids that are generated through partial melting of the mafic lower crust and fractional crystallization at the middle-upper crust, with the involvement of mantle-derived materials. The diorites underwent limited crustal contamination and fractionation of hornblende, plagioclase, and some accessory minerals. They were derived mainly from the mixture of fertile mantle and reworked crustal components, with minor contributions from subduction-related slab fluids and sediment melts. In addition, all the studied Early Paleozoic dioritic–granitic intrusions (ca. 471–427 Ma) formed within subduction-related arc settings. Combined with the tectonic evolution of the Early Paleozoic Qilian orogenic system, we interpret these Cambrian to Silurian dioritic–granitic intrusions as tectonic responses to the subduction (ca. 520–460 Ma) and closure (~440 Ma) of the Proto-Tethys Ocean, whereas the Devonian Huangyanghe intrusion witnessed the final stage of extensional collapse of the Qilian orogenic system at ca. 400–360 Ma.The Tethyan orogenic belt, a significant continent–continent collisional belt in the world, preserves records of oceanic subduction, continental collision, and extensional collapse [1-3]. This belt is divided into the Proto-Tethys (Early Paleozoic), the Paleo-Tethys (Late Paleozoic-Early Mesozoic), and the Neo-Tethys (Late Mesozoic-Cenozoic) stages [4-6]. Originating from the breakup of the Rodinia supercontinent, the Proto-Tethys continued to expand in the Cambrian [3]. Subsequently, the Proto-Tethys started to shrink and closed during the assembly of North China and Siberia–Kazakhstan Cratons during the Late Silurian [7]. The Qilian orogenic belt is the pivotal segment of the Central China Orogen and witnessed the subduction and collision processes during the closure of the
87″;海拔高度:2808米)。这些岩石主要呈灰绿色至灰色,具有中粗粒纹理和块状结构(图 3)。此外,这些岩石有时会被第四纪沉积物覆盖,并发生微弱的蚀变(图 3(g))。研究样本可分为两类:闪长岩(BL-1/2、ZG-2/1 和 ZG-2/2)和石英闪长岩(BL-1/1、YHW-3/1 和 YHW-3/2)。显微镜下,闪长岩显示出多合成孪晶斜长石(50%-55%)、强多色亚方角闪石(25%-35%)、正方斜长石(5%-10%)和正方石英(5%-10%)的组合,以及锆石、磷灰石和磁铁矿等辅助矿物(图 3(b) 和 (f))。相比之下,石英闪长岩显示出近方斜长斜长石(40%-45%)、近方角闪石(30%-35%)、近方辉石(5%-10%)、方形石英(5%-10%),以及锆石、磷灰石和钛铁氧化物等辅助矿物(图 3(d) 和 (h))。锆石晶粒在中国河北省廊坊市豫能地质矿产分离勘查中心采用标准重力和磁力分离技术进行分离,然后在一台双目显微镜下手工挑选,并固定在环氧树脂盘上进行抛光,露出内部。在北京地质分析有限公司,使用一台显微镜和配备加坦化学发光探头的扫描电子显微镜(JSM510)拍摄了透射、反射和阴极射线(CL)图像,并将其用于锆石 U-Pb、Lu-Hf 同位素和微量元素分析。锆石 U-Pb 定年和微量元素分析在兰州大学中国西部(甘肃省)矿产资源重点实验室同步进行,使用的仪器是 Agilent 7900× 电感耦合等离子体质谱仪(ICP-MS)和 Analyte 193 nm 激光烧蚀系统。分析光斑的尺寸设定为 ~30 μm,并使用锆石 Plešovice 91500 和 NIST610 作为标准,以监测精度并计算 207Pb/206Pb 比值的校正因子。中国地质调查局天津中心使用激光烧蚀多收集器ICP-MS(LA-MC-ICP-MS)结合NEW WAVE 193 nm FX激光对锆石Lu-Hf同位素进行了分析。分析光斑尺寸为 50 µm,激光重复频率为 10 Hz,激光能量为 100 mJ。根据 176Lu 衰变常数(1.865 × 10-11 y-1)、贫化地幔与现今 Hf 同位素比值(176Hf/177Hf = 0.28325 和 176Lu/177Hf = 0.0384)、平均大陆地壳(176Lu/177Hf = 0.015),计算了 Hf 模型(TDM)和两级 Hf 模型(TDM2)年龄,并进行了详细分析。全岩地球化学分析在兰州大学中国西部(甘肃省)矿产资源重点实验室进行。主要元素采用 Agilent720Axial ICP 光学发射光谱仪测定,分析误差小于 1.0%。痕量元素和稀土元素采用 Agilent 7700× ICP-MS 分析,分析误差小于 10%。GSP-2 标准用于监控整个分析过程,以确保所分析元素数据的准确性。北祁连山东麓闪长岩-花岗岩侵入体(如毛藏寺、黄羊河、雷公山等)的锆石 U-Pb 平均年龄采用等值线法绘制。采用 ArcGIS 软件 Geostatistical Analyst 模型中的反距离加权插值法绘制等值线图。样品BL-1/1、ZG-2/1和YHW-3/1中的锆石颗粒无色透明,呈八面体至近正方体,多为长棱镜体或短棱镜体。它们的尺寸范围为长 80-250 μm,宽 60-150 μm,长宽比为 4:1 至 1:1。在 CL 图像下,除少数锆石晶粒具有扇形分带纹理外,大多数晶粒都呈现振荡分带纹理,并具有异晶核(图 4(a)-(c))。来自样品 BL-1/1 的 21 个锆石显示 Th 和 U 含量以及 Th/U 比率分别为 5.0-42.8、31.5-255.5 和 0.16-0.33 ppm。 1至-0.3),两阶段钕模型年龄为1.26-1.20 Ga[35],进一步证实了岩浆源中有地壳再加工成分的参与。所研究的闪长岩和已发表的志沟、老虎山侵入体的Mg#值从32.87到75.93不等(平均值=51),大部分超过了40(在线补充表S5)。这意味着岩浆来源主要来自地幔。此外,闪长岩的铌/镱比值在 0.63-1.51 之间,明显高于富集的大洋中脊玄武岩(0.02)[67],这表明来源于肥沃的地幔[68]。此外,这些闪长岩的 Zr/Nb 比值(范围:4.23-43.02)远高于大洋岛屿玄武岩(OIB)(5.8)[67],Nb/Ta 比值(5.31-21.87)低于亚寒带玄武岩。87)低于亚软玉Nb/Ta比(3 wt.%,Rb/Sr < 1,Eu/Eu*为0.8-1.1,而S型花岗岩的Na2O < 3 wt.%,Rb/Sr > 1,Eu/Eu* < 0.6 [82]。我们研究的石英闪长岩与典型的 A 型花岗岩(图 8(d)-(e))明显不同,因为它们不含钠镁铁矿物,Zr 含量[83, 84]和 FeOt/MgO 比率[85]较低,Sr 含量较高[85, 86]。与大洋长花岗岩截然不同的矿物和地球化学特征支持了排除 M 型花岗岩的观点[76]。有几个证据可证实这些花岗岩具有 I 型地球化学特征:(1) 不含富铝矿物,而含有闪石和辉石(图 3);(2) A/CNK值大多低于 1.1(图 8(c));(3) 地球化学亲和性与 I 型花岗岩基本一致(Na2O > 3 wt.%,Rb/Sr < 1,Eu/Eu* 为 0.8-1.1),(4)I 型花岗岩区域内的地球化学图(图 8(f)),以及(5)I 型花岗岩在元素判别图中的变化趋势(图 11(e)和(f))。结合地球化学特征,北祁连造山带东部的花岗岩主要进一步划分为I型钙碱性特征(图8(b)-(f))。钙碱性I型花岗岩通常是通过(a)源于地幔的玄武质岩浆的碎裂结晶[87-90],(b)经类地幔岩浆改造的沉积物质的再加工[91],或(c)有或没有加入源于地幔的黑云母岩浆的黑云母下地壳的部分熔融[92, 93]形成的。花岗岩(石英闪长岩)并非源自闪长岩岩浆,这一点可从以下几个方面得到证明:(1) 从 La/Sm 对 La 图和 Zr/Nb 对 Zr 图(图 13(c)-(d))中可以看出,所研究的石英闪长岩和混合花岗岩并没有经历明显的分块结晶,但部分熔融可能起了关键作用、(2) 根据建模结果(图 12(a)-(b)),BL-1/1 样品与 BL-1/2 样品的成分相似,但两者的主要元素成分差异很大,表明石英闪长岩不可能来自闪长岩岩浆的碎裂结晶。所研究的石英闪长岩(BL-1/1、YHW-3/1)的εHf(t)值为4.2-15.8,TDM2在1171-442 Ma之间(图7(a)-(b)),表明其来源来自地幔或幼壳成分。这与祁连造山带前寒武纪基底岩石和早古生代侵入体的再加工物质有关[94-96]。此外,Mg#值还可以判别岩浆演化过程中是否掺入了地幔衍生物质[97]。由下地壳岩浆岩部分熔融形成的熔体,其初始 Mg# 值通常低于 40,但加入地幔衍生成分后,其 Mg# 值会增加[98]。这些研究岩石的 Mg# 值介于 37.30 至 42.80 之间(在线补充表 S4),表明来源中含有地幔成分。此外,这些样品也被绘制在偏闪长岩到偏黑长岩的部分熔融区域内,以及闪长岩区域内(图 12(c)-(d)),这进
{"title":"Petrogenesis of the Early Paleozoic Dioritic–Granitic Magmatism in the Eastern North Qilian Orogen, NW China: Implications for Tethyan Tectonic Evolution","authors":"Fei Xue, Fan Yang, Weidong Ren, M. Santosh, Zesheng Qian, Yin Huang, Zijian Tan","doi":"10.2113/2024/lithosphere_2023_297","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_297","url":null,"abstract":"The North Qilian Orogen witnessed the opening, subduction, and closure of the Proto-Tethys Qilian Ocean and the post-subduction of multiple exhumation events from Late Neoproterozoic to Early Paleozoic. The Early Paleozoic dioritic–granitic magmatic suites, prominently exposed in the eastern North Qilian Orogen, offer valuable insights into the evolution of the Proto-Tethys Ocean. However, their petrogenesis, magma source, and tectonic evolution remain controversial. Here, we investigate the Leigongshan, Zhigou, and Dalongcun intrusions and present geochronological, geochemical, and isotopic data, aiming to refine the comprehension of their timing and petrogenesis, which will contribute to understanding the tectonic evolution of the Proto-Tethys Ocean. Zircon U-Pb dating reveals mean ages of 471–427 Ma for these intrusions, consistent with compiled formation ages of dioritic–granitic intrusions in the eastern North Qilian Orogen, indicating close temporal links with the tectonic evolution of the Proto-Tethys Ocean during the Early Paleozoic. The studied magmatic rocks could be categorized into two major types: granitoids and diorites. The granitoids are majorly I-type granitoids that are generated through partial melting of the mafic lower crust and fractional crystallization at the middle-upper crust, with the involvement of mantle-derived materials. The diorites underwent limited crustal contamination and fractionation of hornblende, plagioclase, and some accessory minerals. They were derived mainly from the mixture of fertile mantle and reworked crustal components, with minor contributions from subduction-related slab fluids and sediment melts. In addition, all the studied Early Paleozoic dioritic–granitic intrusions (ca. 471–427 Ma) formed within subduction-related arc settings. Combined with the tectonic evolution of the Early Paleozoic Qilian orogenic system, we interpret these Cambrian to Silurian dioritic–granitic intrusions as tectonic responses to the subduction (ca. 520–460 Ma) and closure (~440 Ma) of the Proto-Tethys Ocean, whereas the Devonian Huangyanghe intrusion witnessed the final stage of extensional collapse of the Qilian orogenic system at ca. 400–360 Ma.The Tethyan orogenic belt, a significant continent–continent collisional belt in the world, preserves records of oceanic subduction, continental collision, and extensional collapse [1-3]. This belt is divided into the Proto-Tethys (Early Paleozoic), the Paleo-Tethys (Late Paleozoic-Early Mesozoic), and the Neo-Tethys (Late Mesozoic-Cenozoic) stages [4-6]. Originating from the breakup of the Rodinia supercontinent, the Proto-Tethys continued to expand in the Cambrian [3]. Subsequently, the Proto-Tethys started to shrink and closed during the assembly of North China and Siberia–Kazakhstan Cratons during the Late Silurian [7]. The Qilian orogenic belt is the pivotal segment of the Central China Orogen and witnessed the subduction and collision processes during the closure of the","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"13 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139951432","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Late Cenozoic Cooling History of the Xigaze Fore-Arc Basin along the Yarlung–Zangbo Suture Zone (Southern Tibet): New Insights from Low-Temperature Thermochronology 沿雅鲁藏布断裂带(藏南)的日喀则前弧盆地的晚新生代冷却史:低温热年代学的新见解
IF 2.4 4区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_211
Shida Song, Zhiyuan He, Wenbo Su, Linglin Zhong, Kanghui Zhong, Stijn Glorie, Yifan Song, Johan De Grave
The Tibetan Plateau is currently the widest and highest elevation orogenic plateau on Earth. It formed as a response to the Cenozoic and is still ongoing collision between the Indian and Eurasian plates. The Xigaze fore-arc basin distributed along the Indus–Yarlung suture zone in southern Tibet preserves important information related to the late Cenozoic tectonic and topographic evolution of the plateau. In this study, apatite fission track (AFT) thermochronology was carried out on twelve sandstone samples from the middle segment of the Xigaze basin and additionally on four sedimentary rocks from the neighboring Dazhuka (Kailas) and Liuqu Formations. Inverse thermal history modeling results reveal that the fore-arc basin rocks experienced episodic late Oligocene to Miocene enhanced cooling. Taking into account regional geological data, it is suggested that the late Oligocene-early Miocene (~27–18 Ma) cooling recognized in the northern part of the basin was promoted by fault activity along the Great Counter thrust, while mid-to-late Miocene-accelerated exhumation was facilitated by strong incision of the Yarlung and Buqu rivers, which probably resulted from enhanced East Asian summer monsoon precipitation. Sandstone and conglomerate samples from the Dazhuka and Liuqu Formations yielded comparable Miocene AFT apparent ages to those of the Xigaze basin sediments, indicative of (mid-to-late Miocene) exhumation soon after their early Miocene burial (> ~3–4 km). Additionally, our new and published low-temperature thermochronological data indicate that enhanced basement cooling during the Miocene prevailed in vast areas of central southern Tibet when regional exhumation was triggered by both tectonic and climatic contributing factors. This recent and widespread regional exhumation also led to the formation of the high-relief topography of the external drainage area in southern Tibet, including the Xigaze fore-arc basin.Orogenic belts are dominant topographic features on Earth and are characterized by high tectonic activity and high elevations. They provide the best natural laboratory to study the coupling between tectonics, erosion, and climate [1-4]. In these regions, negative feedback between fast denudation and high elevation causes enhanced erosion that, in turn, tends to reduce the topography. The collision between India and Asia led to the formation of the Tibetan Plateau, which stands ~4–5 km high over a region of ~3 million km2 (Figures 1(a) and (b)). The southern Tibetan Plateau (i.e., the southern Lhasa terrane—Tethyan Himalaya) is characterized by a high-elevation and low-relief landscape (i.e., “flat” highland) [5]. Structurally, several ~E–W trending large-scale thrust faults and a series of ~N–S striking normal faults are well developed in the southern Lhasa terrane [6-8]. Furthermore, the west-to-east flowing Yarlung river runs through the southern Tibetan Plateau, its source is high in western Tibet, and it cuts through the Namche Barwa
日喀则前弧盆地所有碎屑岩样品的AFT中心年龄都明显小于其白垩纪沉积年龄(表1和表2),这与之前在其他采样剖面的观察结果一致[21,24]。考虑到这些样品的秩方概率 P(χ²)大多大于 10%,且离散度较低,它们一定达到了比 AFT 总退火温度更高的温度(通过深埋或由于邻近的冈底斯岩浆活动产生的加热),从而完全重置了 AFT 系统。根据同样的标准,邻近的大竹卡地层和柳曲地层的样品在出露前也经历了完全的AFT热退火。在对另外一部分样品进行252Cf辐照和/或重离子轰击(针对有足够数量晶粒的样品)后,有六个样品(X28、33、35、36、38和39)产生了超过30条可测量的封闭轨道。由于自发轨迹密度非常低,在 X30、31 和 37 样品中只发现了约 20 条封闭轨迹。在轨迹长度方面,样品 X27 的平均轨迹长度值非常长,约为 15.4 μm,但这仅仅是基于四个封闭轨迹的测量结果。所有建模样品的预期模型[71]合并在图 4(b)中(并标明 95% 的概率范围区间),每个样品的详细单个热模型见在线补充图 S3(HeFTy 建模结果也在此提供,以供验证)。如上所述,我们意识到所有模型都是基于少于 50 个封闭的裂变轨迹,这与通常的做法不符;因此,模型的约束性较弱,在解释时应采取必要的谨慎态度。为了更定量地描述冷却历史,我们根据经验将冷却速率(在我们研究区域的大陆碰撞带背景下)~5-10 °C/Ma和>10 °C/Ma分别定义为中度冷却和快速冷却[13, 27, 77]。其中样品X28在渐新世晚期(~27 Ma)进入APAZ,在中新世早期(~18 Ma)冷却到~50 °C,冷却速率为~7.8 °C/Ma(图4(b))。X33、35 和 38 号样品的曲线一般显示在 ~20-9 Ma 期间通过 APAZ 冷却,冷却速度为中等到快速(图 4(b))。样品 X36 和 39 的 AFT 年龄较小,但它们的冷却历史却相当。它们在中新世晚期(约10-5 Ma)停留在APAZ的时间相对较短,冷却速度较快,约为12 °C/Ma(图4(b))。值得注意的是,与 QTQt 模型相比,HeFTy 模型为上述 6 个样本提供了相似的冷却路径(在线补充图 S3)。另外三个样品(X30、31 和 37)在中新世(~21-8 Ma)也显示了类似的中度冷却。尽管这三个模型是基于更低的封闭轨道测量值,即 ~3-4 km)。这一观察结果与之前在冈底斯浴成岩南缘狭长的开拉斯地层其他地段进行的热年代学研究结果基本一致[11, 100],这些研究结果表明开拉斯砾岩在中新世中期(~17-15 Ma)迅速隆升。我们还强调,雅鲁藏布江及其支流同时发生的侵蚀可能促进了砾岩的掘起。在主取样剖面的西侧,来自同向聚合柳曲地层的两块砾岩的AFT中心年龄分别为~15和~10Ma(表2)。它们比该地层大致的沉积年龄(~20-19 Ma [55, 101])要年轻。根据非晶AFT、磷灰石和锆石(U-Th)/He测年,Li等人[101]认为,由于印度河-雅鲁藏布江缝合带中古雅鲁藏布江的侵蚀,柳曲砾岩在~12-10Ma开始隆升。在Li等人的研究中[101],非晶AFT年龄峰值(大部分或完全退火的晶粒)在~11Ma,所有AHe单晶粒年龄都在~8-6Ma范围内。本研究分析了日喀则前弧盆地中部及邻近的大竹卡(开拉斯)地层和柳曲地层的沉积岩,显示了不同的晚渐新世至中新世加速冷却事件(图4(b))。新生代晚期(如渐新世晚期至中新世晚期)的基底冷却在北面的冈底斯浴成岩,尤其是其南缘也有广泛的记录[102]。在这一时期,拉萨南部地层经历了一段活跃的构造和岩浆活动。一方面,区域性断层,如冈底斯推力和大逆冲推力,分别在 ~27-17 Ma 和 ~23-17 Ma 期间开始活跃[78, 103, 104]。
{"title":"Late Cenozoic Cooling History of the Xigaze Fore-Arc Basin along the Yarlung–Zangbo Suture Zone (Southern Tibet): New Insights from Low-Temperature Thermochronology","authors":"Shida Song, Zhiyuan He, Wenbo Su, Linglin Zhong, Kanghui Zhong, Stijn Glorie, Yifan Song, Johan De Grave","doi":"10.2113/2024/lithosphere_2023_211","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_211","url":null,"abstract":"The Tibetan Plateau is currently the widest and highest elevation orogenic plateau on Earth. It formed as a response to the Cenozoic and is still ongoing collision between the Indian and Eurasian plates. The Xigaze fore-arc basin distributed along the Indus–Yarlung suture zone in southern Tibet preserves important information related to the late Cenozoic tectonic and topographic evolution of the plateau. In this study, apatite fission track (AFT) thermochronology was carried out on twelve sandstone samples from the middle segment of the Xigaze basin and additionally on four sedimentary rocks from the neighboring Dazhuka (Kailas) and Liuqu Formations. Inverse thermal history modeling results reveal that the fore-arc basin rocks experienced episodic late Oligocene to Miocene enhanced cooling. Taking into account regional geological data, it is suggested that the late Oligocene-early Miocene (~27–18 Ma) cooling recognized in the northern part of the basin was promoted by fault activity along the Great Counter thrust, while mid-to-late Miocene-accelerated exhumation was facilitated by strong incision of the Yarlung and Buqu rivers, which probably resulted from enhanced East Asian summer monsoon precipitation. Sandstone and conglomerate samples from the Dazhuka and Liuqu Formations yielded comparable Miocene AFT apparent ages to those of the Xigaze basin sediments, indicative of (mid-to-late Miocene) exhumation soon after their early Miocene burial (&gt; ~3–4 km). Additionally, our new and published low-temperature thermochronological data indicate that enhanced basement cooling during the Miocene prevailed in vast areas of central southern Tibet when regional exhumation was triggered by both tectonic and climatic contributing factors. This recent and widespread regional exhumation also led to the formation of the high-relief topography of the external drainage area in southern Tibet, including the Xigaze fore-arc basin.Orogenic belts are dominant topographic features on Earth and are characterized by high tectonic activity and high elevations. They provide the best natural laboratory to study the coupling between tectonics, erosion, and climate [1-4]. In these regions, negative feedback between fast denudation and high elevation causes enhanced erosion that, in turn, tends to reduce the topography. The collision between India and Asia led to the formation of the Tibetan Plateau, which stands ~4–5 km high over a region of ~3 million km2 (Figures 1(a) and (b)). The southern Tibetan Plateau (i.e., the southern Lhasa terrane—Tethyan Himalaya) is characterized by a high-elevation and low-relief landscape (i.e., “flat” highland) [5]. Structurally, several ~E–W trending large-scale thrust faults and a series of ~N–S striking normal faults are well developed in the southern Lhasa terrane [6-8]. Furthermore, the west-to-east flowing Yarlung river runs through the southern Tibetan Plateau, its source is high in western Tibet, and it cuts through the Namche Barwa ","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":"9 1","pages":""},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139463005","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
Lithosphere
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1