首页 > 最新文献

Lithosphere最新文献

英文 中文
Petrogenesis and Tectonic Implication of Jurassic Granites in Central Guangdong, SE China: Constraints from Zircon U-Pb-Hf-O and Whole-Rock Geochemical and Sr-Nd Isotopic Data 中国东南部广东中部侏罗纪花岗岩的成岩过程和构造影响:锆石U-Pb-Hf-O和全岩地球化学及Sr-Nd同位素数据的制约因素
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_327
Zhiguang Lai, Yongxin Xu, Chunbo Xin, Xuewen Luo
The origin and tectonic regime responsible for the inland Jurassic granites in Southeast (SE) China remain controversial. This study presents zircon secondary ion mass spectrometry (SIMS) U-Pb ages, in situ zircon Hf-O isotopes, and whole-rock geochemical and Sr-Nd isotopic data for the Fogang and Xinxing Batholiths in central Guangdong. Mineralogical and geochemical features indicate that these granites are high-K (>4.8 wt% K2O at 72 wt% SiO2), calc-alkaline I-type granites. SIMS U-Pb analyses on magmatic zircons yield consistent ages ranging from 158 to 163 Ma, suggesting that the Fogang and Xinxing granites were emplaced in the period of 163–158 Ma. In addition, these granites have whole-rock initial Sr87/Sr86 ratios of 0.6802–0.7072 and negative εNd(t) values of −9.5 to −8.2, zircon negative εHf(t) values of −12.34 to −0.56, and high δ18O values of 7.64‰–10.08‰. The above features imply that the granites were most likely generated through the mixture of supracrustal sedimentary components with minor addition of mantle-derived magmas. Granites from the Fogang and Xinxing Batholiths in SE China should be derived from the Proterozoic crustal reworking due to asthenosphere upwelling or underplating and intrusion of mafic magmas. These Jurassic granites reflect anorogenic magmatism probably formed in an intraplate extensional setting resulted from the foundering of the flat slab beneath SE China.Granite is a primary component of continental crust, preserving abundant information about the formation, evolution, and accretion of crust, as well as interactions between the crust and mantle. Multiperiod Mesozoic granites are widely distributed in Southeast (SE) China, with a concentration in the Triassic, Jurassic, and Cretaceous, respectively Figure 1(a) [1-3]. Among them, the Nanling region is mainly characterized by Jurassic granite, while the coastal areas are dominated by Cretaceous granite (see Figures 1(a) and 1(b)) [4]. The coexistence of multiperiod rocks from different origins is of great significance for understanding the genesis of the granite, crust-mantle interaction, magma differentiation, and mixing processes [5-8]. Previous researchers have reported the geochronology, petrology, mineralogy, and geochemistry of the Nanling granites. However, there has been ongoing debate regarding their petrogenesis and tectonic mechanism.The Fogang and Xinxing Batholiths represent Late Mesozoic basements in the Nanling region, with Fogang Batholith being the largest and most representative granite basement in the region (Figures 1(c) and 1(d)) [4, 7]. Due to intense fractional crystallization, the batholiths exhibit complex geochemical characteristics, making their genetic types difficult to determine [9]. Different scholars have classified the Fogang Batholith as I-type [8, 9], A-type [5], S-type [6], or high-fractionated I-type granites [10]. Similarly, there are different views on the genetic type of the Xinxing Batholith, such as I-type [11], A-typ
中国科学院广州地球化学研究所同位素地球化学国家重点实验室选取了23CH01-1、23CH02-1和23JM03-1三个样品,采用CAMECA IMS-1280 HR SIMS(二次离子质谱法)进行锆石U-Pb定年。椭球形斑的大小约为 20 × 30 μm。使用 Plěsovice 和清湖的标准锆石校准 U-Th 含量和 Pb/U 比值。结果见在线补充表 S3。数据处理采用 Isoplot 4.15 程序[32]。锆石 O 同位素分析采用中国南京大学的 CAMECA IMS 1300-HR3 进行。电流为 2-3 nA、总冲击能量为 20 keV 的 133Cs+ 主离子束聚焦在样品表面。研究中使用了 5 × 5 μm 的光栅,并使用法向入射电子枪实现了电荷补偿。总分析时间约为 4.5 分钟。锆石 Hf 同位素分析是在南方海洋科学与工程广东实验室(珠海)对先前进行过 SIMS U-Pb 年龄和 O 同位素分析的相同锆石点上进行的。分析使用 Neptune XT/Thermo Scientific 仪器进行,该仪器配备了 NWR FemtoUC/Elemental Scientific 激光器。分析所用的光斑尺寸为 50 µm,重复频率为 5 Hz,激光功率为 ~3 J/cm2。本研究使用锆石 91500、蓬莱和普列索维采样品作为参考标准。分析程序、校准方法和参考标准由 Wu 等人[33] 和 Li 等人[34] 描述。测量的 176Hf/177Hf 比值归一化为 179Hf/177Hf = 0.7325。为了计算 εHf 值,将 176Lu 的衰变常数值设定为 1.865 × 10-11/y [35],并使用了软玉和贫化地幔现今的 176Hf/177Hf 和 176Lu/177Hf 比率(分别为 0.28277 和 0.0332,0.28325 和 0.0384)[36]。单级模型年龄是相对于贫化地幔计算的,使用的现今 176Lu/177Hf 比率为 0.28325,176Lu/177Hf 比率为 0.0384 [37]。此外,两阶段模型 Hf 年龄的计算假设大陆地壳 176Lu/177Hf 平均值为 0.015 [38]。Sr-Nd 同位素分析采用经典的两步离子交换色谱法,并在武汉样品溶液分析公司使用 Thermo Fisher Scientific Triton Plus 多收集器热电离质谱仪进行测量。有关分析程序的详细信息,请参见 Li 等人的文章[39]。用于钕同位素分析的整块岩石粉末在 HF + HNO3 + HClO4 溶解前加入混合的 149Sm-150Nd 示踪剂,然后溶解在 Savillex Teflon 旋盖胶囊中。整个过程的空白值为 100 pg Sm-Nd。同位素比值按 146Nd/144Nd = 0.7219 进行质量分馏校正。在数据收集期间,采用国际标准样品 JNdi-1 来评估仪器的稳定性。JNdi-1 Nd 标准的测量值为 143Nd/144Nd = 0.512108 ± 0.000006(n = 2,2SD)。JNdi-1 的 143Nd/144Nd 数据与之前公布的热电离质谱(TIMS)技术[39]数据吻合。从 23CH01-1 至 23CH01-10、23CH02-1 至 23CH02-10、23JM03-1 至 23JM03-10 三种花岗岩中选取了 30 个样品进行全岩主要元素和痕量元素分析。23CH01 样品的特点是 SiO2 含量高(72.89,74.46 wt%),K2O 含量高(5.42,5.99 wt%),CaO(1.04,1.41 wt%)和 Na2O(2.65,2.99 wt%)含量相对较低。它们的 K2O/Na2O 比值较高(1.86,2.16),在 SiO2 与 K2O 的关系图中属于霰石系列(图 4(a))。根据 Na2O + K2O 与 SiO2 总碱二氧化硅(TAS)的关系图,这些样品都属于亚碱性,而且都属于花岗岩(图 4(b))。此外,它们的 Al2O3 含量为 13.18%-13.75%,铝饱和度指数(ASI)(A/CNK = [Al2O3/(CaO + Na2O + K2O) mol%] )为 1.03-1.06,表明具有金属铝特征(图 4(c))。铁#数(FeOT/[FeOT + MgO])为 0.76-0.81 ,跨越了镁质花岗岩和铁质花岗岩(图 4(d))。23CH02 样品的 SiO2 含量为 73.18-75.63 wt%,K2O 含量为 4.96-5.83 wt%,CaO 和 Na2O 含量相对较低(0.62,1.64 wt%)(2.66,3.43 wt%),K2O/Na2O 比值较高(1.45,2.20),在 SiO2 与 K2O 的关系图中属于高 K(钙碱性)至霞石系列(图 4(a))。在 TAS 图(图 4(b))上,它们都表现出亚碱性花岗岩特征。它们的 ASI 指数范围很宽,为 0.92-1。
{"title":"Petrogenesis and Tectonic Implication of Jurassic Granites in Central Guangdong, SE China: Constraints from Zircon U-Pb-Hf-O and Whole-Rock Geochemical and Sr-Nd Isotopic Data","authors":"Zhiguang Lai, Yongxin Xu, Chunbo Xin, Xuewen Luo","doi":"10.2113/2024/lithosphere_2023_327","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_327","url":null,"abstract":"The origin and tectonic regime responsible for the inland Jurassic granites in Southeast (SE) China remain controversial. This study presents zircon secondary ion mass spectrometry (SIMS) U-Pb ages, in situ zircon Hf-O isotopes, and whole-rock geochemical and Sr-Nd isotopic data for the Fogang and Xinxing Batholiths in central Guangdong. Mineralogical and geochemical features indicate that these granites are high-K (>4.8 wt% K2O at 72 wt% SiO2), calc-alkaline I-type granites. SIMS U-Pb analyses on magmatic zircons yield consistent ages ranging from 158 to 163 Ma, suggesting that the Fogang and Xinxing granites were emplaced in the period of 163–158 Ma. In addition, these granites have whole-rock initial Sr87/Sr86 ratios of 0.6802–0.7072 and negative εNd(t) values of −9.5 to −8.2, zircon negative εHf(t) values of −12.34 to −0.56, and high δ18O values of 7.64‰–10.08‰. The above features imply that the granites were most likely generated through the mixture of supracrustal sedimentary components with minor addition of mantle-derived magmas. Granites from the Fogang and Xinxing Batholiths in SE China should be derived from the Proterozoic crustal reworking due to asthenosphere upwelling or underplating and intrusion of mafic magmas. These Jurassic granites reflect anorogenic magmatism probably formed in an intraplate extensional setting resulted from the foundering of the flat slab beneath SE China.Granite is a primary component of continental crust, preserving abundant information about the formation, evolution, and accretion of crust, as well as interactions between the crust and mantle. Multiperiod Mesozoic granites are widely distributed in Southeast (SE) China, with a concentration in the Triassic, Jurassic, and Cretaceous, respectively Figure 1(a) [1-3]. Among them, the Nanling region is mainly characterized by Jurassic granite, while the coastal areas are dominated by Cretaceous granite (see Figures 1(a) and 1(b)) [4]. The coexistence of multiperiod rocks from different origins is of great significance for understanding the genesis of the granite, crust-mantle interaction, magma differentiation, and mixing processes [5-8]. Previous researchers have reported the geochronology, petrology, mineralogy, and geochemistry of the Nanling granites. However, there has been ongoing debate regarding their petrogenesis and tectonic mechanism.The Fogang and Xinxing Batholiths represent Late Mesozoic basements in the Nanling region, with Fogang Batholith being the largest and most representative granite basement in the region (Figures 1(c) and 1(d)) [4, 7]. Due to intense fractional crystallization, the batholiths exhibit complex geochemical characteristics, making their genetic types difficult to determine [9]. Different scholars have classified the Fogang Batholith as I-type [8, 9], A-type [5], S-type [6], or high-fractionated I-type granites [10]. Similarly, there are different views on the genetic type of the Xinxing Batholith, such as I-type [11], A-typ","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139465141","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
Investigating the Mechanism of Strong Roof Weighting and Support Resistance Near Main Withdrawal Roadway in Large-Height Mining Face 大采高工作面主要回撤巷道附近的强顶配重和支护阻力机理研究
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_288
Junwu Du, Qingxiang Huang
Aiming at investigating the strong roof weighting when the large height mining face is nearing the main withdrawal roadway, the 52,304 working face (WF) nearly through the main withdrawal roadway mining in a colliery of Shendong coalfield was taken as the research background. The ground pressure, roof structure, and superposition effect of stress in the last mining stage were studied by field measurement, physical simulation, and numerical calculations. The obtained results demonstrated that the main roof formed the “long step voussoir beam” structure under the influence of the main withdrawal roadway. The superposition effect of the front abutment pressure of the WF and the concentrated stress of the main withdrawal roadway caused the stress asymmetrical distribution on the two sides -level hard rock straof the main withdrawal roadway, and the stability of the pillar on the mining side decreases. The initial average periodic weighting interval was 20.7 m. While the WF approaches the main withdrawal roadway, the pillar near the WF of the main withdrawal roadway collapsed, the main roof was broken ahead of the WF, and the actual roof control distance of support and the periodic weighting interval increased by 2.56 and 1.26 times the normal state, respectively. Consequently, the “static load” of the immediate roof and the “dynamic load” of the sliding unsteadiness of the long step voussoir beam increased. The structural model of the “long step voussoir beam” under the superposition of “static and dynamic load” was established concerning those results, and an expression was proposed to compute the support resistance. Meanwhile, the mechanism of strong roof weighting was revealed when the WF was nearly through the main withdrawal roadway. The research conclusion is expected to provide a guideline for the safe withdrawal of the large-height mining faces under similar conditions.To increase the withdrawal speed and yield efficacy of the working face (WF) and avoid the tense connection between face mining and entry driving, predriving double withdrawal roadway is widely used in coal mines to reinforce the withdrawal operation [1]. In this scheme, the main and auxiliary withdrawal roadways are advance driven at the stop-mining line of the WF. After the primary withdrawal roadway is connected with the WF, the reinforcements are withdrawn through the connecting entry between the primary and secondary withdrawal roadways. Consequently, the withdrawal speed of the WF increases 3–5 times compared with the traditional methods, thereby increasing the production rate and improving the mining efficiency [2, 3]. Although this method has remarkable advantages, it has some shortcomings, including low mining speed in the last mining stage, concentrated mining-induced stress field, and high roof pressure [4]. More specifically, the superposition effect of the lateral and front abutment pressure of the main withdrawal roadway and the WF near the main withdrawal roadway
在这种情况下,采空区一侧煤柱的集中应力为10.0 MPa,是正常情况下的2.5倍。当距离为 2.0 m 时,WF 的前墩压力传递到主回风巷道的煤柱壁上,峰值应力达到 11.5 MPa,WF 的应力场发生超前叠加。分析结果表明,当 WF 接近主回风巷道时,前方支护压力与主回风巷道集中应力的叠加效应显著,导致采掘侧煤柱完全垮塌。物理模拟和数值计算表明,当 WF 接近主回风巷道时,前方支护压力和主回风巷道集中应力的叠加效应显著,主回风巷道 WF 附近的煤柱垮落,WF 前方主顶板破碎。在这种情况下,WF 的风险最大。考虑到主回风巷道对 WF 的影响,WF 的顶板控制距离和周期加权间隔都有所增加。因此,主顶板呈现出 "长台阶伏梁 "结构。此外,直接顶板的 "静载荷 "和 "长台阶伏溜梁 "结构的 "动载荷 "都有所增加,由液压动力支架承担。为进一步研究末采阶段强采压的作用机理,根据物理模拟和数值计算得出的基本结论,针对大采高工作面顶板的结构特点,建立了 "静、动荷载 "叠加下的 "长台阶伏梁 "结构模型,如图 8 所示,其中 h1 和 h 分别表示直接顶和主顶板的厚度。其中,h1 和 h 分别表示直接顶和主顶板的厚度,M 和 N 分别表示主顶板的关键块体,ω θ 分别表示块体的旋转角度。此外,b 是相应的台阶高度。A、C 和 B 代表关键砌块的铰接点。T 是水平挤压力。RM 和 W 分别为支撑所承受的动荷载和静荷载。R1 为开采侧支柱的残余加固力。R0 为煤矸石对关键块 N 的加固反力,P 为液压动力支架的工作阻力。参照 "伏流梁 "结构的应力分析方法[27],由于岩块转角挤压面的高度较小,破碎关键块接触面的高度可以忽略不计。因此,WF 接近主要回撤巷道时的关键块模型可简化如下(图 9):其中,l 表示 WF 通过主要回撤巷道前的平均周期加权间隔;lz 为 WF 接近主要回撤巷道时周期加权间隔的增加长度;h 为主要顶板厚度。P1 和 P2 分别是关键区块 M 和 N 的重量及其承受的荷载。此外,QA 和 QB 分别代表铰链接头 A 和 B 处的剪力。在 C 点,关键块 M 由关键块 N 加固,而关键块 N 则由掘进巷道中的塌落矸石加固。根据关键块的平衡特性,关键块 M 在 C 点的力矩总和为 0,可用数学公式表示如下:此外,沿关键块垂直方向的结果力如公式(5)所示。将公式(1)、(3)、(4)和(5)合并可得到以下表达式:根据 "S-R "稳定性理论[28],除非满足以下不等式,否则该结构容易发生滑动失稳:其中 P1 为关键块 M 所承受的荷载,可通过以下表达式计算得出:根据 Terzaghi 的土压力理论,荷载传递系数的计算公式为。 最后,揭示了WF接近主要回撤巷道时顶板冒落和支架破碎的机理。根据物理模拟,52304 大采高工作面在无主要回撤巷道影响的情况下,主要顶板冒落形成 "阶梯伏梁 "结构,主要顶板关键块体的平均塌落角和旋转角分别为 65°和 5°。在主撤退巷道影响前,平均周期配重间隔为 20.7 m,平均支护工作阻力为 17540 kN,WF 前承压力和主撤退巷道集中应力的叠加效应导致支柱完全倒塌。同时,主顶板在 WF 前方塌陷,形成 "长台阶伏梁 "结构。因此,支护的实际顶板控制距为正常状态下的 2.56 倍,周期加权间隔为正常状态下的 1.26 倍。在此基础上,建立了 "静、动载荷 "叠加下的 "长台阶溜子梁 "结构模型,并推导出了液压动力支架在 WF 接近主撤退巷道时的合理工作阻力表达式。最后,揭示了当 WF 接近主撤退巷道时顶板冒落和支架破碎的机理。本文有望为类似条件下大采高回采工作面的安全回撤提供指导。本文收录了主要相关数据,相应作者将根据合理要求提供其他相关数据。作者声明,本文的发表不存在利益冲突。我们感谢国家自然科学基金、陕西省自然科学基础研究计划、煤炭资源精细勘查与智能开发国家重点实验室对本研究的支持。感谢学术编辑和匿名审稿人提出的善意建议和宝贵意见。
{"title":"Investigating the Mechanism of Strong Roof Weighting and Support Resistance Near Main Withdrawal Roadway in Large-Height Mining Face","authors":"Junwu Du, Qingxiang Huang","doi":"10.2113/2024/lithosphere_2023_288","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_288","url":null,"abstract":"Aiming at investigating the strong roof weighting when the large height mining face is nearing the main withdrawal roadway, the 52,304 working face (WF) nearly through the main withdrawal roadway mining in a colliery of Shendong coalfield was taken as the research background. The ground pressure, roof structure, and superposition effect of stress in the last mining stage were studied by field measurement, physical simulation, and numerical calculations. The obtained results demonstrated that the main roof formed the “long step voussoir beam” structure under the influence of the main withdrawal roadway. The superposition effect of the front abutment pressure of the WF and the concentrated stress of the main withdrawal roadway caused the stress asymmetrical distribution on the two sides -level hard rock straof the main withdrawal roadway, and the stability of the pillar on the mining side decreases. The initial average periodic weighting interval was 20.7 m. While the WF approaches the main withdrawal roadway, the pillar near the WF of the main withdrawal roadway collapsed, the main roof was broken ahead of the WF, and the actual roof control distance of support and the periodic weighting interval increased by 2.56 and 1.26 times the normal state, respectively. Consequently, the “static load” of the immediate roof and the “dynamic load” of the sliding unsteadiness of the long step voussoir beam increased. The structural model of the “long step voussoir beam” under the superposition of “static and dynamic load” was established concerning those results, and an expression was proposed to compute the support resistance. Meanwhile, the mechanism of strong roof weighting was revealed when the WF was nearly through the main withdrawal roadway. The research conclusion is expected to provide a guideline for the safe withdrawal of the large-height mining faces under similar conditions.To increase the withdrawal speed and yield efficacy of the working face (WF) and avoid the tense connection between face mining and entry driving, predriving double withdrawal roadway is widely used in coal mines to reinforce the withdrawal operation [1]. In this scheme, the main and auxiliary withdrawal roadways are advance driven at the stop-mining line of the WF. After the primary withdrawal roadway is connected with the WF, the reinforcements are withdrawn through the connecting entry between the primary and secondary withdrawal roadways. Consequently, the withdrawal speed of the WF increases 3–5 times compared with the traditional methods, thereby increasing the production rate and improving the mining efficiency [2, 3]. Although this method has remarkable advantages, it has some shortcomings, including low mining speed in the last mining stage, concentrated mining-induced stress field, and high roof pressure [4]. More specifically, the superposition effect of the lateral and front abutment pressure of the main withdrawal roadway and the WF near the main withdrawal roadway","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139758627","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
The Long-Lasting Exhumation History of the Ötztal-Stubai Complex (Eastern European Alps): New Constraints from Zircon (U–Th)/He Age-Elevation Profiles and Thermokinematic Modeling Ötztal-Stubai复合体(东欧阿尔卑斯山)的长期侵蚀历史:来自锆石(U-Th)/He 年龄-海拔剖面和热运动模型的新制约因素
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_174
Kyra Hölzer, Reinhard Wolff, Ralf Hetzel, István Dunkl
The Eastern European Alps formed during two orogenic cycles, which took place in the Cretaceous and Cenozoic, respectively. In the Ötztal-Stubai Complex—a thrust sheet of Variscan basement and Permo-Mesozoic cover rocks—the record of the first (Eoalpine) orogeny is well preserved because during the second (Alpine) orogeny, the complex remained largely undeformed. Here, new zircon (U–Th)/He (ZHe) ages are presented, and thermokinematic modeling is applied to decipher the cooling and exhumation histories of the central part of the Ötztal-Stubai Complex since the Late Cretaceous. The ZHe ages from two elevation profiles increase over a vertical distance of 1500 m from 56 ± 3 to 69 ± 3 Ma (Stubaital) and from 50 ± 2 to 71 ± 4 Ma (Kaunertal), respectively. These ZHe ages and a few published zircon and apatite fission track ages were used for inverse thermokinematic modeling. The modeling results show that the age data are well reproduced with a three-phase exhumation history. The first phase with relatively fast exhumation (~250 m/Myr) during the Late Cretaceous ended at ~70 Ma and is interpreted to reflect the erosion of the Eoalpine mountain belt. As Late Cretaceous normal faults occur at the margins of the Ötztal-Stubai Complex, normal faulting may have also contributed to the exhumation of the study area. Subsequently, a long period with slow exhumation (<10 m/Myr) prevailed until ~16 Ma. This long-lasting phase of slow exhumation suggests a rather low topography with little relief in the Ötztal-Stubai Complex until the mid-Miocene, even though the Alpine orogeny had already begun in the Eocene with the subduction of the European continental margin. Accelerated exhumation since the mid-Miocene (~230 m/Myr) is interpreted to reflect the erosion of the mountain belt due to the development of high topography in front of the Adriatic indenter and repeated glaciations during the Quaternary.Mountain belts with thick continental crust, such as the European Alps, the Himalaya, or the North American Cordillera, are formed during long-lasting plate convergence with crustal shortening by nappe stacking and folding [1-3]. Due to the isostatic uplift of the thickened crust, the internal parts of such orogens become the locus of erosion, which removes material at the Earth’s surface and leads to the cooling and exhumation of metamorphic rocks [4, 5]. Apart from erosion, another important mechanism that may cause rock exhumation and cooling is normal faulting because tectonic slip along normal faults transports rocks in their footwalls toward the Earth’s surface [6-9].To quantify the cooling history of metamorphic rocks, it is necessary to determine the temperature conditions in rocks through time, which is possible by applying geochronological methods such as Sm/Nd, Rb/Sr, or Ar/Ar dating to minerals with different closure temperatures [10-12]. The final cooling in the upper crust from temperatures of ~250°C to ~60°C can be constrained with low-temperature ther
两个模型的厚度均为 55 公里,并在模型上添加了现今的地形。由于研究区域的最高峰海拔约为 3.5 千米(图 3),因此两个模型的最大厚度均为 58.5 千米。鉴于欧洲阿尔卑斯山的局部高地形是第四纪大量冰川和河流侵蚀的结果[82, 83],我们将所有模型运行开始时的地形相对于现在的地形缩小了三分之二。在所有模型运行中,我们使用的热扩散率为 1×10-6 m2/s,辐射产热量为 3×10-6 W/m3 [84],随着深度的增加,辐射产热量减小,电子折叠长度为 20 km(表 4)[52]。模型底部的温度边界条件为 800°C,并在所有模型运行期间保持不变[52]。这些热参数和边界条件导致模型上部10千米处的初始地热梯度约为30°C/千米,并向下递减。为了根据PECUBE模拟的热历史计算冷却年龄,我们对ZFT使用了Galbraith和Laslett的扩散模型[85],对AFT使用了Ketcham的退火模型[86]。为了计算 ZHe 年龄,我们将 Reiners 等人[87] 的扩散模型调整到稍低的闭合温度(即 170°C,冷却速率为 10°C/Myr),以考虑到我们的锆石样品相当低的辐射损伤(1.28 × 1016 到 2.26 × 1017 α衰变事件/克;表 2)[88,89]。我们测量的 ZHe 年龄(以及来自 Elias [30] 的 AFT 和 ZFT 年龄)与模型预测年龄之间的不拟合度计算公式为:N 为数据点数,αi,data 为观测数据,αi,model 为预测值,σi,data 为年龄的不确定性[79]。为了找到最佳拟合模型,我们采用了两步邻域算法反演,并对参数空间进行迭代采样,以最小化观测年龄与预测年龄之间的不拟合度[80]。对于考内塔尔模型和斯图巴塔尔模型,我们总共进行了 12,880 次模型运行(即 161 次迭代,每次迭代 80 个模型)。为确保对参数空间进行充分的探索性采样,我们选择了 0.9 的高重采比率(参见 79)。为了对五个自由模型参数的不确定性进行定量约束,我们使用了边际概率密度函数[81-90].现在,我们采用反模拟的方法来再现这些冷却年龄,并以定量的方式重建三相掘出的历史。所有模型中的掘起和物质迁移都是垂直发生的。在每个模型的运行中,三个不同的掘起阶段的掘起速率都保持不变。通过大量的 12 880 次模型运行,我们探索了以下参数空间。模型从110 Ma开始,第一个掘出阶段在80-65 Ma之间结束。第一阶段的掘起速率介于 0.1 至 0.7 km/Ma 之间。第二阶段的掘起速率较低,介于0至0.05 km/Myr之间。第三掘出阶段开始于 20 至 15 Ma 之间,允许的掘出速率为 0.15 至 0.35 km/Myr。请注意,这五个自由参数是同时变化的。反演建模的结果在两个海拔剖面的三幅图中各有说明(图 6)。这些图的横轴和纵轴显示了所探索的参数空间。两个最佳拟合模型--白星表示--表明斯图拜塔尔横断面和考内塔尔横断面的掘起历史相似。在第一掘起阶段,最佳拟合模型得出的掘起速率分别为 0.27 ± 0.03 和 0.23 ± 0.03 km/Myr。这一阶段的结束(以及第二阶段的开始)分别发生在 67.9 ± 3.7 Ma(斯图巴塔尔)和 71.7 ± 3.0 Ma(考内塔尔)(图 6(a) 和(b))。第二掘起阶段的掘起速率较低,分别为~0.004 km/Myr(Stubaital)和~0.007 km/Myr(Kaunertal),而第三阶段的掘起速率分别为0.23 ± 0.14和0.23 ± 0.02 km/Myr(图6(a)和(b))。两个剖面的第三阶段开始时间几乎相同(即 16.0 ± 1.2 Ma 和 16.2 ± 1.2 Ma)。两个最佳拟合模型预测的ZHe、AFT和ZFT年龄与在Stubaital和Kaunertal海拔剖面上观测到的冷却年龄非常吻合(图7)。我们从Ötztal-Stubai复合体中部的两个海拔剖面上得到的新ZHe年龄从~70 Ma到~50 Ma不等(即在Stubaital从69 ± 3 Ma到56 ± 3 Ma,在Kaunertal从71 ± 4 Ma到50 ± 2 Ma;图4和图5)。 与来自厄兹塔尔的额外 ZHe 年龄一起,综合数据集在年龄-海拔空间中定义了一个窄带(在线补充图 S1),这表明厄兹塔尔-Stubai 复合体的中部自晚白垩世以来一直是一个连贯的地壳块。我们的两个最佳拟合热运动模型证实了这一解释,这两个热运动模型揭示了两个海拔剖面相似的掘起历史(图6),并很好地再现了几乎所有的ZHe年龄(图7)。在第一阶段,晚白垩世的隆升速度约为250 m/Myr(即Stubaital模型和Kaunertal模型的隆升速度分别为270 ± 30 m/Myr和230 ± 30 m/Myr)。这一掘起阶段一直持续到大约70Ma(即斯图巴塔尔模型为68±4Ma,考内塔尔模型为72±3Ma),随后在古新世和中新世早期出现了一个漫长的缓慢掘起阶段。在这一阶段,我们的样品通过锆石中氦的部分保留区缓慢冷却,这也解释了两个剖面中最高和最低样品之间约 15-20 Ma 的年龄差异。第二阶段的排湿速率分别为 4 ± 7 m/Myr(Stubaital)和 7 ± 6 m/Myr(Kaunertal;图 6)。尽管这些数值有很大的不确定性,但我们的模型清楚地表明,当时的掘起率很低。如果第二阶段的掘出率较高,则ZHe冷却年龄会明显变小,最高样本和最低样本之间的年龄差也会变小。因此,只有在第二掘出阶段采用小于10 m/Myr的掘出率时,我们的热运动模型才有可能再现观测到的冷却年龄。这样低的侵蚀速率可能有些令人惊讶,但我们注意到,在不同的现今山地环境中,利用宇宙成因核素测量到的类似侵蚀速率为 2-15 m/Myr[91-93]。同样值得注意的是,第二阶段的剥蚀率明显低于两个年龄-海拔图中回归线的斜率(图 5)。这一观察结果表明,将这些斜率解释为掘出率是不可取的(参见 20、52),因为很难甚至不可能证明基本假设得到满足(见第 5 节)。根据我们的模型,第三掘起阶段开始于中新世,时间为 16±1 Ma(图 6)。这一阶段的掘起速率约为 230 m/Myr,很好地再现了 P5/13 和 P4/2 样本(位于两个模型的边界内)的 AFT 年龄(图 7)。最后,我们注意到,在我们的模型中,两种蜕皮速率的变化是瞬时发生的,而在自然界中,蜕皮速率的变化更可能是渐变的,可能要经历100-200万年。通过将三个蜕皮阶段的持续时间与各自的蜕皮速率相乘,我们可以从两个最佳拟合模型中量化出蜕皮的总量。这些模型表明,斯图巴塔尔和考内塔尔地区分别有~15公里和~13公里的岩石被移走。在本节中,我们将结合东阿尔卑斯山的构造和地质演化,讨论根据热运动学模型(图8)得出的厄兹塔尔-斯图拜复合体的三期隆升史。在美利亚塔洋被东南向东倾的俯冲带吞噬之后,奥斯特阿尔卑斯基底岩层及其沉积覆盖层的堆积造成了地壳的缩短和增厚,这种现象始于早白垩世(约 135 Ma),并随着时间的推移自东向西传播[如 39、60、94]。由于大陆内部的俯冲作用,奥陶系部分岩层经历了奥陶系高P/低T变质作用[如3
{"title":"The Long-Lasting Exhumation History of the Ötztal-Stubai Complex (Eastern European Alps): New Constraints from Zircon (U–Th)/He Age-Elevation Profiles and Thermokinematic Modeling","authors":"Kyra Hölzer, Reinhard Wolff, Ralf Hetzel, István Dunkl","doi":"10.2113/2024/lithosphere_2023_174","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_174","url":null,"abstract":"The Eastern European Alps formed during two orogenic cycles, which took place in the Cretaceous and Cenozoic, respectively. In the Ötztal-Stubai Complex—a thrust sheet of Variscan basement and Permo-Mesozoic cover rocks—the record of the first (Eoalpine) orogeny is well preserved because during the second (Alpine) orogeny, the complex remained largely undeformed. Here, new zircon (U–Th)/He (ZHe) ages are presented, and thermokinematic modeling is applied to decipher the cooling and exhumation histories of the central part of the Ötztal-Stubai Complex since the Late Cretaceous. The ZHe ages from two elevation profiles increase over a vertical distance of 1500 m from 56 ± 3 to 69 ± 3 Ma (Stubaital) and from 50 ± 2 to 71 ± 4 Ma (Kaunertal), respectively. These ZHe ages and a few published zircon and apatite fission track ages were used for inverse thermokinematic modeling. The modeling results show that the age data are well reproduced with a three-phase exhumation history. The first phase with relatively fast exhumation (~250 m/Myr) during the Late Cretaceous ended at ~70 Ma and is interpreted to reflect the erosion of the Eoalpine mountain belt. As Late Cretaceous normal faults occur at the margins of the Ötztal-Stubai Complex, normal faulting may have also contributed to the exhumation of the study area. Subsequently, a long period with slow exhumation (<10 m/Myr) prevailed until ~16 Ma. This long-lasting phase of slow exhumation suggests a rather low topography with little relief in the Ötztal-Stubai Complex until the mid-Miocene, even though the Alpine orogeny had already begun in the Eocene with the subduction of the European continental margin. Accelerated exhumation since the mid-Miocene (~230 m/Myr) is interpreted to reflect the erosion of the mountain belt due to the development of high topography in front of the Adriatic indenter and repeated glaciations during the Quaternary.Mountain belts with thick continental crust, such as the European Alps, the Himalaya, or the North American Cordillera, are formed during long-lasting plate convergence with crustal shortening by nappe stacking and folding [1-3]. Due to the isostatic uplift of the thickened crust, the internal parts of such orogens become the locus of erosion, which removes material at the Earth’s surface and leads to the cooling and exhumation of metamorphic rocks [4, 5]. Apart from erosion, another important mechanism that may cause rock exhumation and cooling is normal faulting because tectonic slip along normal faults transports rocks in their footwalls toward the Earth’s surface [6-9].To quantify the cooling history of metamorphic rocks, it is necessary to determine the temperature conditions in rocks through time, which is possible by applying geochronological methods such as Sm/Nd, Rb/Sr, or Ar/Ar dating to minerals with different closure temperatures [10-12]. The final cooling in the upper crust from temperatures of ~250°C to ~60°C can be constrained with low-temperature ther","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139476961","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
Asymmetrical Microfracture Density Across an Active Thrust Fault: Evidence from the Longmen Shan Fault, Eastern Tibet 跨越活动性推断断层的不对称微裂隙密度:来自西藏东部龙门山断层的证据
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_193
Hu Wang, Peisheng Luo, Yi Liang, Dongming Li, Kaijin Li, Lin Deng, Lichun Chen
Microfracture density in fault damage zones can reflect spatial variability that decays in intensity as a function of distance from the fault, which is crucial in understanding the mechanical, seismological, and fluid-flow properties of the fault system. However, few studies explored the characteristics of fracture density between the two sides of active dip-slip faults due to rare field observations. Here, we measured and modeled microfractures across an active thrust fault associated with the 2008 Mw 7.9 Wenchuan earthquake in the Longmen Shan, eastern Tibetan Plateau. The results showed that the microfracture density at the Qingping site developed more intensely in the hanging wall than in the footwall for an exposed thrust fault, indicating an asymmetrical pattern. The hidden thrust fault at the Jushui site showed that microfractures developed more intensely in vertical planes in the hanging wall than in the footwall, whereas microfractures developed similarly in horizontal planes within the two sides, indicating a quasiasymmetrical pattern. Comparing the data at the two sites with computational modeling, we suggest that fault geometry might exert a first-order control of the asymmetrical microfracture density pattern, which is helpful for revealing different deformational behaviors of rock masses in the fault damage zones and better understanding the hanging-wall effect for evaluating seismic hazards on active thrust faults.A fault damage zone, expressed as a zone with numerous fractures surrounding a narrow fault core, has been considered to be related to coseismic loading and, therefore, has the potential to reveal the rock deformational mechanics and past earthquake rupture conditions [1-7]. Moreover, such a damage zone is expected to act as conduits, barriers, or combined conduit-barrier systems that play a fundamental role in crustal fluid flow [8-10]. Therefore, quantitative determination of characteristics of fractures in the fault damage zone is critical to understand the mechanical and seismological properties of the fault system.Geometrically, fracture density is one of the key parameters in evaluating the spatial variability that decays in intensity as a function of distance from the fault [11, 12]. Many studies have measured micro/mesofracture density on fault-perpendicular transects to show that fracture density decreases gradually away from the fault core, which can be simplified to fit either an exponential decay model [13] or a power law decay model [14, 15] in the fault damage zone. Moreover, previous studies have suggested that the characteristics of fracture density might be influenced by the amount of slip across the fault, the size of the fault, lithology, rupture processes, and movement history [8, 13, 16]. For example, Caine et al. [8] suggested that a wide damage zone may indicate the effect of more repeated seismic events with greater accumulative deformation than that of a narrow damage zone. Ostermeijer et al. [12]
此外,VTS2 的最大 MLPA、MNPA 和 DZW 似乎大于 HTS 和 VTS1,这表明微裂缝的发展在三维尺度上可能并不均匀。出现这种差异的原因可能是 VTS2 平面主要与压缩变形引起的最大主应力有关。然而,巨水遗址的数据显示,整个断层的微断裂密度非常复杂。具体来说,断层两侧 HTS 的最大 MLPA 和 MNPA 以及主要 DZW 似乎相似。然而,VTS1 的数据显示,悬壁的最大 MLPA 和 MNPA 值以及主要 DZW 值均大于底壁。虽然 VTS2 在巨水岩层脚壁的线性拟合结果在主要 DZW 方面存在一些不确定性,但悬壁的最大 MLPA 和 MNPA 均大于脚壁。因此,巨水遗址垂直面和水平面的微裂隙发育似乎有所不同。垂直平面的微裂隙密度形态不对称,突出表现为挂壁的微裂隙比底壁的更密集,而水平平面的微裂隙密度形态基本一致。此外,VTS2 的最大 MLPA、MNPA 和主要 DZW 值似乎大于 HTS 和 VTS1,这与清平矿区的情况一致。这种差异可能与断层的成熟度有关,成熟度表明断层破坏带微裂隙的发育状况[15]。例如,YBF 的规模和同震偏移要比 JGF 大得多[26-28]。以前的研究表明,断层破坏带的微裂隙密度可能受断层滑移量、断层大小、岩性、破裂过程和运动历史的影响[8, 13, 16]。这两个地点在岩性、破裂过程和运动历史方面表现出相似的特征。例如,两个地点的地层主要由石灰岩和砂岩组成(图 1(b)和 1(c));这两个断层在过去的古地震中同时发生过断裂[29],并且至少在中新世中期就已经发生过推移变形[36]。另外两个因素,如更大的滑动和断层规模,可能会导致更密集的微裂隙密度[13]。然而,上述因素都不能很好地解释两地微断裂密度的不对称模式。共震变形量可被视为反映微裂隙密度的代用指标,它可能与长期瞬态共震荷载有关[1-3, 5, 6],这有助于我们理解微裂隙密度模式的机理。两个地点的同震变形模型显示,悬壁有陡峭的收敛趋势,而底壁有平缓的收敛趋势(图 10 和 11),这与两个地点薄断面测量的非对称衰减曲线一致。模型还显示,清平遗址的同震变形量(平均约为巨水遗址的 10 倍)远大于巨水遗址。这表明,清平地块的微裂缝密度将比巨水地块更为密集,这一点也可以从两个地块的最大 MLPA 和 MNPA 以及主要 DZW 参数比较中得到证明(表 2)。换句话说,不对称的微裂隙密度在清平病区应该比在巨水病区更明显,这或许可以解释为什么巨水病区的 HTS 没有明显表现出不对称模式。然而,巨水遗址从 VTS1 到 VTS2 微断裂的不对称模式为我们进一步探索微断裂的发展提供了另一条线索。虽然巨水遗址的断层是隐伏的,但建模结果(图 11)清楚地表明,断层埋藏点(深度在 0 至 2 km 之间)上方的潜在变形界面几乎是垂直分布的,这与野外调查对取样点与断层距离的估计是一致的。
{"title":"Asymmetrical Microfracture Density Across an Active Thrust Fault: Evidence from the Longmen Shan Fault, Eastern Tibet","authors":"Hu Wang, Peisheng Luo, Yi Liang, Dongming Li, Kaijin Li, Lin Deng, Lichun Chen","doi":"10.2113/2024/lithosphere_2023_193","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_193","url":null,"abstract":"Microfracture density in fault damage zones can reflect spatial variability that decays in intensity as a function of distance from the fault, which is crucial in understanding the mechanical, seismological, and fluid-flow properties of the fault system. However, few studies explored the characteristics of fracture density between the two sides of active dip-slip faults due to rare field observations. Here, we measured and modeled microfractures across an active thrust fault associated with the 2008 Mw 7.9 Wenchuan earthquake in the Longmen Shan, eastern Tibetan Plateau. The results showed that the microfracture density at the Qingping site developed more intensely in the hanging wall than in the footwall for an exposed thrust fault, indicating an asymmetrical pattern. The hidden thrust fault at the Jushui site showed that microfractures developed more intensely in vertical planes in the hanging wall than in the footwall, whereas microfractures developed similarly in horizontal planes within the two sides, indicating a quasiasymmetrical pattern. Comparing the data at the two sites with computational modeling, we suggest that fault geometry might exert a first-order control of the asymmetrical microfracture density pattern, which is helpful for revealing different deformational behaviors of rock masses in the fault damage zones and better understanding the hanging-wall effect for evaluating seismic hazards on active thrust faults.A fault damage zone, expressed as a zone with numerous fractures surrounding a narrow fault core, has been considered to be related to coseismic loading and, therefore, has the potential to reveal the rock deformational mechanics and past earthquake rupture conditions [1-7]. Moreover, such a damage zone is expected to act as conduits, barriers, or combined conduit-barrier systems that play a fundamental role in crustal fluid flow [8-10]. Therefore, quantitative determination of characteristics of fractures in the fault damage zone is critical to understand the mechanical and seismological properties of the fault system.Geometrically, fracture density is one of the key parameters in evaluating the spatial variability that decays in intensity as a function of distance from the fault [11, 12]. Many studies have measured micro/mesofracture density on fault-perpendicular transects to show that fracture density decreases gradually away from the fault core, which can be simplified to fit either an exponential decay model [13] or a power law decay model [14, 15] in the fault damage zone. Moreover, previous studies have suggested that the characteristics of fracture density might be influenced by the amount of slip across the fault, the size of the fault, lithology, rupture processes, and movement history [8, 13, 16]. For example, Caine et al. [8] suggested that a wide damage zone may indicate the effect of more repeated seismic events with greater accumulative deformation than that of a narrow damage zone. Ostermeijer et al. [12]","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139500784","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
Data-Driven Dynamic Inversion Method for Complex Fractures in Unconventional Reservoirs 非常规储层复杂裂缝的数据驱动动态反演方法
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_347
Ruixue Jia, Xiaoming Li, Xiaoyong Ma, Liang Zhu, Yangdong Guo, Xiaoping Song, Pingde Wang, Jiantao Wang
Hydraulic fracturing is a crucial technology for enhancing the recovery of oil and gas from unconventional reservoirs. Accurately describing fracture morphology is essential for accurately predicting production dynamics. This article proposes a new fracture inversion model based on dynamic data-driven methods, which is different from the conventional linear elastic fracture mechanics model. This method eliminates the need to consider complex mechanical mechanisms, resulting in faster simulation speeds. In the model, the fracture morphology is constrained by combining microseismic data and fracturing construction data, and the fracture tip propagation domain is introduced to characterize the multi-directionality of fracture propagation. The simulated fracture exhibits a multi-branch fracture network morphology, aligning more closely with geological understanding. In addition, the influence of microseismic signal intensity on the direction of fracture propagation is considered in this study. The general stochastic approximation (GSA) algorithm is employed to optimize the direction of fracture propagation. The proposed method is applied to both the single-stage fracturing model and the whole well fracturing model. The research findings indicate that in the single-stage fracturing model, the inverted fracture morphology aligns closely with the microseismic data, with a fitting rate of the fracturing construction curve exceeding 95%, and a microseismic data fitting rate exceeding 93%. In the whole well fracturing model, a total of 18 sections were inverted. The fitting rate between the overall fracture morphology and the microseismic data reached 90%. The simulation only took 5 minutes, demonstrating high computational efficiency and meeting the needs of large-scale engineering fracture simulation. This method can effectively support geological modeling and production dynamic prediction.The world has abundant shale gas reservoir resources; however, due to the influence of reservoir rock properties, its development poses significant challenges [1-4]. Hydraulic fracturing technology can effectively enhance the physical properties of reservoirs and form complex fracture networks within the reservoir, thereby promoting oil and gas production [5-7]. In order to assess the development impact of shale gas reservoirs and devise appropriate development plans, it is necessary to establish a numerical model that is specific to the shale gas reservoir in question. Accurately describing the post-fracturing fracture morphology is crucial for model construction and subsequent flow simulations, as it is a key factor in ensuring the accuracy of model calculation results [8]. Moreover, the morphology of fractures post-fracturing is often highly complex, characterized by a network structure of fractures [9, 10]. Many existing fracture propagation models only consider a simplified quasi-three-dimensional or three-dimensional straight fracture structure. However, these mo
(为了验证本文所提模型的拟合性能,我们以实际区块模型的单级压裂模型为例。为了验证本文所提模型的拟合性能,我们以实际区块模型的单段压裂模型为例,建立了一个同时受微地震数据和水力压裂施工数据约束的裂缝反演模型,并对裂缝形态进行了动态模拟。相关模拟参数见表 1。图 3(a)显示了模型中使用的微震数据。实际水力压裂施工曲线如图 3(b)所示。从微震数据可以看出,裂缝半长约为 98 米,裂缝宽度约为 38 米,裂缝分布较为均匀,没有异常复杂区域。同时,为了获得每个微震数据点的坐标,使用 GetData 软件进行提取。考虑到每个微震数据点的半径比例,考虑了微震强度系数І。在处理构造曲线时,本文直接考虑了裂缝的传播段,忽略了井筒压力上升的过程。裂缝的实时反演结果如图 4 所示。图中描述了不同时间间隔(20 分钟、45 分钟、60 分钟、85 分钟、105 分钟和 135 分钟)的裂缝形态。从图中可以明显看出,在每个阶段,泵压的拟合率都很高,超过 95%。同时,裂缝形态与微地震数据的拟合率也很高,达到 93% 左右。模型反演的断裂长度为 95 米,断裂宽度为 35 米,与微震结果相差 3 米,这在工程应用中是可以接受的。此外,断裂形态充分解释了微震数据的分布范围,并表现出复杂的断裂形态,这是传统数值模拟方法难以实现的。这些结果有力地证明了本文提出的模型可以有效地满足工程要求。本文将该方法应用于实际区块的单井模型,并选择某油田的 W1 井作为研究对象。将该方法应用于实际区块的单井模型,选取某油田的 W1 井作为研究对象,在开发初期对该井进行压裂,并进行微震监测。微震监测结果如图 5 所示。根据 W1 井的微震数据和压裂施工数据,进行了裂缝形态反演。W1 井的模拟参数如表 2 所示,压裂形态结果如图 6 和图 7 所示。图 6 显示了带有微地震点的断裂形态,表明与微地震数据高度匹配,拟合率约为 90%。大部分误差集中在井筒附近,可能是由于井筒边界的重大影响,这在模型中很难考虑。图 7 提供了裂缝结构的清晰视图,平均反演裂缝长度为 189 米,与微震数据的偏差约为 7 米,为微震监测结果提供了合理解释。但必须注意的是,在实际压裂监测过程中,微震数据点是动态的。在本模型中,考虑的是静态微震数据点,这会对压裂方向产生重大影响。在未来的研究中,将探索纳入动态微震数据,并将其与裂缝扩展过程拟合,以实现更高精度的裂缝形态。本文建立了基于微震数据和水力压裂施工数据的多目标约束裂缝反演模型。本文建立了基于微地震数据和水力压裂施工数据的多目标约束裂缝反演模型,并利用概念和实际区块案例验证了模型计算的可靠性和准确性。主要结论如下:构建了基于微地震数据和水力压裂施工数据的多目标约束裂缝反演模型。该模型考虑了裂缝传播的多方向性以及微地震信号强度对裂缝传播方向的影响。微地震数据用于约束断裂传播方向,水力压裂施工数据用于约束断裂节点数量。通过这两种约束方法的耦合,实现了对断裂形态的动态反演。
{"title":"Data-Driven Dynamic Inversion Method for Complex Fractures in Unconventional Reservoirs","authors":"Ruixue Jia, Xiaoming Li, Xiaoyong Ma, Liang Zhu, Yangdong Guo, Xiaoping Song, Pingde Wang, Jiantao Wang","doi":"10.2113/2024/lithosphere_2023_347","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_347","url":null,"abstract":"Hydraulic fracturing is a crucial technology for enhancing the recovery of oil and gas from unconventional reservoirs. Accurately describing fracture morphology is essential for accurately predicting production dynamics. This article proposes a new fracture inversion model based on dynamic data-driven methods, which is different from the conventional linear elastic fracture mechanics model. This method eliminates the need to consider complex mechanical mechanisms, resulting in faster simulation speeds. In the model, the fracture morphology is constrained by combining microseismic data and fracturing construction data, and the fracture tip propagation domain is introduced to characterize the multi-directionality of fracture propagation. The simulated fracture exhibits a multi-branch fracture network morphology, aligning more closely with geological understanding. In addition, the influence of microseismic signal intensity on the direction of fracture propagation is considered in this study. The general stochastic approximation (GSA) algorithm is employed to optimize the direction of fracture propagation. The proposed method is applied to both the single-stage fracturing model and the whole well fracturing model. The research findings indicate that in the single-stage fracturing model, the inverted fracture morphology aligns closely with the microseismic data, with a fitting rate of the fracturing construction curve exceeding 95%, and a microseismic data fitting rate exceeding 93%. In the whole well fracturing model, a total of 18 sections were inverted. The fitting rate between the overall fracture morphology and the microseismic data reached 90%. The simulation only took 5 minutes, demonstrating high computational efficiency and meeting the needs of large-scale engineering fracture simulation. This method can effectively support geological modeling and production dynamic prediction.The world has abundant shale gas reservoir resources; however, due to the influence of reservoir rock properties, its development poses significant challenges [1-4]. Hydraulic fracturing technology can effectively enhance the physical properties of reservoirs and form complex fracture networks within the reservoir, thereby promoting oil and gas production [5-7]. In order to assess the development impact of shale gas reservoirs and devise appropriate development plans, it is necessary to establish a numerical model that is specific to the shale gas reservoir in question. Accurately describing the post-fracturing fracture morphology is crucial for model construction and subsequent flow simulations, as it is a key factor in ensuring the accuracy of model calculation results [8]. Moreover, the morphology of fractures post-fracturing is often highly complex, characterized by a network structure of fractures [9, 10]. Many existing fracture propagation models only consider a simplified quasi-three-dimensional or three-dimensional straight fracture structure. However, these mo","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140313388","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
Effects of Stress on Transport Properties in Fractured Porous Rocks 应力对断裂多孔岩输运特性的影响
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_103
Jing Ba, Jinyi Min, Lin Zhang, José M. Carcione
The nonlinear characteristics of the rock transport properties (permeability and electrical conductivity in this study) as a function of stress are closely related to the geometry of the pore space, which consists of stiff pores, microcracks, or microfractures. We consider two behaviors of the pore space, one linear and the other exponential, related to the stiff pores and microfractures, respectively, where the relation between stress and strain can be described by the Two-Part Hooke’s Model. With this model, the relations between porosity, transport properties, and effective stress (confining minus pore pressure) can be obtained and validated with the experimental data of four tight sandstones collected from the Shaximiao Formation of Sichuan Basin, southwest China. The agreement is good. At low effective stresses, the closure of cracks is the main mechanism affecting the transport properties, whose behavior is similar in terms of their parameters. Subsequently, experimental data of nine tight sandstones from the Yanchang Formation, collected from the Ordos Basin, west China, are employed to confirm the previous results, indicating that the fluid and electrical current follow the same path in the pore space.Reservoir rocks have pores, cracks, or microfractures and are generally heterogeneous [1-4]. The deformation under loading is different in stiff pores and microfractures, which affects the elastic and transport properties, especially in low-permeability rocks. Since cracks provide a permeability path for the flow of reservoir fluids [5-9], understanding of the relationships between the transport properties and effective stress is important for detecting and monitoring reservoir fluids.Previous studies revealed that the exponential function describes the behavior of permeability and conductivity as a function of effective stress [10-20]. However, an important point is to describe the behavior of the sharp decrease of these transport properties when the effective stress increases at low values, especially for low-permeability rocks [21, 22]. The power law has also been adopted to describe such variation [23-27]. For instance, Jones and Owens [28] and Walsh [29] reformulated the expression of power law. On the other hand, Kaselow and Shapiro [30] applied a four-parameter exponential equation to analyze the electrical conductivity as a function of the effective pressure.The closure of cracks with increasing effective stress leads to lower porosity, and permeability or electrical conductivity shows a similar behavior. The transport properties as a function of porosity can be studied with a power law [31, 32] or by analyzing experimental data [21]. Archie [33] established an empirical relation between the formation factor (the ratio between bulk resistivity and that of water) and porosity. Subsequently, some researchers investigated the relationships between electrical conductivity and porosity [34, 35], clay content [36-38], crack radii, aspect r
孔隙度是用氦膨胀法得到的,渗透率是用非稳态脉冲瞬态衰减技术测定的。使用阻抗-电容-电阻仪测量了盐水饱和(5% 浓度)时的电导率。孔隙度和渗透率的测量范围为 2-58 兆帕,电导率的测量范围为 5-35 兆帕。表 1 列出了这些特性。图 2 显示的是 TS4 样品的薄切片,这是一种岩性长石砂岩,大部分为孔隙接触胶结。矿物成分主要是石英,还有长石和岩石碎片,间隙物质包括粘土和方解石。图 3 和图 4 分别比较了四种致密砂岩的理论孔隙度和渗透率(公式 (8) 和 (9))以及实验孔隙度和渗透率与有效应力的关系。在低应力下,这些特性呈非线性下降,在高应力下呈线性趋势。低应力下的非线性变形可能是由于与软质部分相对应的槽状裂缝闭合造成的。TS1-TS4 样品的孔隙率(渗透率)拟合系数 R2 分别为 0.97 (0.98)、0.96 (0.96)、0.92 (0.85) 和 0.97 (0.88)。拟合参数的获取方法(见表 2)与 Liu 等人[46]和 Zheng 等人[22]的方法相同。具体来说,在高应力下,log(ke) 和 σ 之间存在线性关系,可以用 ke=ke,1exp[-βCeϕe,1Δσ] 拟合。然后,直线的斜率 -βCeje,1 可用来估算 β 值,其中 Ce 和 ϕe,1 的值是通过孔隙率-应力数据拟合得到的。ke,1的值是通过外推σ1处的直线得出的。根据 kt=α(ϕt)m 可以估算出 α 和 m 的值。表 2 显示,TS2 样品的软部孔隙率最大,而 TS3 样品的软部孔隙率最小。这表明在低应力下,TS2 的孔隙率和渗透率下降幅度最大。相反,TS3 的这些特性下降幅度最小。在高应力下,孔隙率和渗透率对数与应力呈线性关系,这表明在此范围内可以忽略软质部分的贡献。此外,高应力范围可设定为 20-58 兆帕,2-15 兆帕可定义为低应力范围,可用于软质部分的参数估算。与孔隙率和渗透率类似,电导率随应力呈非线性下降。然而,这种行为在 TS4 样品中并不明显,与软质部分相对应的常数 n 明显小于其他样品(见图 6)。TS3 和 TS4 软部件数据点较少。这是因为与 TS3 和 TS4 中的软质部分一样,微裂缝或微裂纹的长宽比可能较小,在较低的有效压力范围内,软质部分会急剧减小并逐渐接近闭合。方程 (14) 的结果与实验数据一致,四个样品的拟合系数 R2 分别为 0.99、0.96、0.99 和 0.95。在低应力时,电导率变化明显,而在高应力时,裂缝闭合,对流体流动的连续性没有贡献,呈线性下降趋势。与渗透率类似,在低应力时,TS2 样品的软部孔隙率下降幅度最大,而 TS3 样品的软部孔隙率下降幅度较小。结果表明,软质部分是流体迁移的主要路径,在低应力时不可忽视。在高应力下,log(Se)和 σ 之间存在线性关系,可用公式 (11) 拟合。然后,直线的斜率 -aCeϕe,1 可用来估算 a,其中 Ce 和 ϕe,1 的值是通过孔隙率-应力数据拟合得到的。Se,1的值通过外推σ1处的直线来确定。根据公式 (12) 可以得到软质部分的电导率,然后根据公式 (13) 可以估算出 b 和 n。表 2 显示,γt,1 只占总孔隙率的一小部分,并且随着应力的增加而急剧下降。软岩部分的体积模量 Kt(9-13 兆帕)远小于致密砂岩样本的体积模量 Kt(29-31 千兆帕)。有些参数(Ce 和 Kt)有其特定含义,可以通过实验室实验直接测量,但本文是通过拟合实验数据计算得出的。此外,Liu 等人[46] 指出了通过拟合计算得到的孔隙压缩性和软质部分体积模量等参数的合理性。
{"title":"Effects of Stress on Transport Properties in Fractured Porous Rocks","authors":"Jing Ba, Jinyi Min, Lin Zhang, José M. Carcione","doi":"10.2113/2024/lithosphere_2023_103","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_103","url":null,"abstract":"The nonlinear characteristics of the rock transport properties (permeability and electrical conductivity in this study) as a function of stress are closely related to the geometry of the pore space, which consists of stiff pores, microcracks, or microfractures. We consider two behaviors of the pore space, one linear and the other exponential, related to the stiff pores and microfractures, respectively, where the relation between stress and strain can be described by the Two-Part Hooke’s Model. With this model, the relations between porosity, transport properties, and effective stress (confining minus pore pressure) can be obtained and validated with the experimental data of four tight sandstones collected from the Shaximiao Formation of Sichuan Basin, southwest China. The agreement is good. At low effective stresses, the closure of cracks is the main mechanism affecting the transport properties, whose behavior is similar in terms of their parameters. Subsequently, experimental data of nine tight sandstones from the Yanchang Formation, collected from the Ordos Basin, west China, are employed to confirm the previous results, indicating that the fluid and electrical current follow the same path in the pore space.Reservoir rocks have pores, cracks, or microfractures and are generally heterogeneous [1-4]. The deformation under loading is different in stiff pores and microfractures, which affects the elastic and transport properties, especially in low-permeability rocks. Since cracks provide a permeability path for the flow of reservoir fluids [5-9], understanding of the relationships between the transport properties and effective stress is important for detecting and monitoring reservoir fluids.Previous studies revealed that the exponential function describes the behavior of permeability and conductivity as a function of effective stress [10-20]. However, an important point is to describe the behavior of the sharp decrease of these transport properties when the effective stress increases at low values, especially for low-permeability rocks [21, 22]. The power law has also been adopted to describe such variation [23-27]. For instance, Jones and Owens [28] and Walsh [29] reformulated the expression of power law. On the other hand, Kaselow and Shapiro [30] applied a four-parameter exponential equation to analyze the electrical conductivity as a function of the effective pressure.The closure of cracks with increasing effective stress leads to lower porosity, and permeability or electrical conductivity shows a similar behavior. The transport properties as a function of porosity can be studied with a power law [31, 32] or by analyzing experimental data [21]. Archie [33] established an empirical relation between the formation factor (the ratio between bulk resistivity and that of water) and porosity. Subsequently, some researchers investigated the relationships between electrical conductivity and porosity [34, 35], clay content [36-38], crack radii, aspect r","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139589747","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
Production Dynamic Characteristic of Fractured Wells in Multilayer Reservoirs Considering the Effect of Non-Uniform Flux Distribution 考虑非均匀流量分布影响的多层储层中裂缝井的生产动态特性
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_173
Yanzhong Liang, Bailu Teng, Wanjing Luo
Hydraulic fracturing stimulation, which improves matrix permeability and reduces production costs, has been extensively used in the exploitation of multilayer reservoirs. However, little research on the production dynamic characteristics of vertically fractured wells in stratified reservoirs has been done in the literature. The influence of flux variation along the fracture on the pressure transient behavior has been ignored in these previous works. Therefore, this paper introduces a novel semi-analytical model for fractured wells in multilayer reservoirs, in which the finite difference method is used to characterize fluid flow in the fracture and the Green’s function method is used to characterize fluid flow in the matrix. With the aid of the model, the production dynamic characteristics of fractured wells in multilayer reservoirs can be readily investigated. In addition, based on the assumption of nonuniform flux distribution along the fracture, we successfully recognize four flow regimes occurring in the pressure drop and pressure derivative curves. Following that, the influences of several parameters on the pressure dynamics and layered flux contribution are studied. The calculation results indicate that a larger storability ratio, as well as a larger permeability ratio, can increase the values of the pressure drop and the pressure derivative; the greater the fracture height, the greater the fluid flow into each layer of the fracture. During the production of this model, increasing the fracture conductivity can reduce the pressure drop and pressure derivative, which means lower flow resistance in the fracture.With the growing dependence on fossil energy, deep and ultra-deep areas have steadily evolved into the next main potentials of resource exploration and development. Recently, China has consistently discovered a huge number of deep-layer reservoirs, such as the Puguang, Tahe, Shunbei, and Anyue oilfields, showing great resource potentialities and considerable economic benefits [1]. Considering the influence of the complex sedimentary environment, most deep reservoirs are composed of several layers with different stratigraphic characteristics. Commingling production is commonly adopted to increase producing profit for stratified reservoirs in oil and gas fields [2]. Given this, extensive literature was related to the pressure dynamics of a vertical well in stratified reservoirs [3-6]. Rahman and Mattar [7] derived a new analytical solution for the commingled-layered reservoir with unequal initial pressures in the Laplace domain. Onwunyili and Onyekonwu [8] developed a coupled model that can more accurately simulate the commingled production behavior of multilayer reservoirs. Shi et al. [9] investigated the impact of the vertical inhomogeneous closed boundary radii on pressure transient behaviors of the multilayered commingled reservoir. These previous researches give us a basic understanding of the production dynamic characteristics for th
方程如下:其中,x0i,j、y0i,j 和 z0i,j 是连续声源(i, j)的位置坐标。为了计算方便,我们将方程 (4) 转换为无量纲形式。公式(5)中的无量纲参数定义见表 1。最后,将方程(5)应用于 Nf 个裂缝元素,建立 Nf 个无量纲裂缝压力与 Nf 个从基体进入裂缝的无量纲通量之间的关系。无量纲井底压力可以用 Peaceman[17] 提出的模型求解,其中井筒电导率假定为无限大。由于与井筒接触的每个裂缝元素都可能发生流体交换行为,因此可以给出一般方程:其中 pwD 为无量纲井底压力,qf-wD 为从裂缝进入井筒的无量纲通量,pnwD 为无量纲井筒压力,reqD 为无量纲等效半径。借助方程 (6),我们有 Nj 个流动方程参与线性方程组的求解。在本节中,我们可以发现已经给出了 2Nf + Nj + 1 个流动方程和 2Nf + Nj + 1 个未知数(Nf 无量纲压裂压力、Nf 从基体到压裂的无量纲通量、Nj 从压裂到井筒的无量纲通量和一个无量纲井底压力)。为了证明所提模型的正确性,我们比较了所提模型、Bennett 等人的分析模型[19]和商业软件(Eclipse)的计算结果。值得注意的是,Bennett 等人[19]的模型忽略了储层中与断裂面平行的流动,这意味着储层中发生的流动行为是二维的。图 3 展示了在 Eclipse 中建立的物理模型的俯视图和侧视图,其中红色网格通过设置较高的渗透率表示断裂。值得注意的是,属于裂缝区域的局部网格系统经过了细化,以更好地匹配实际结果。储层模型被离散为 201 × 201 × 21 个网格,每个网格的大小为 10 × 10 × 1 米。在方法验证过程中,基质、裂缝和流体的参数如下:层 1(ctm1 = 1.2 × 10-5 MPa-1,jm1 = 0.1,h1 = 10 m,km1 = 0.01 md),第 2 层(ctm2 = 1.5×10-5 MPa-1,jm2 = 0.3,h2 = 10 m,km2 = 0.03 md),kf = 1 × 104 md,Cw = 0 bbl/psi,Cf = 200 md∙m,ctf = 1.3 × 10-5 MPa-1,jf = 0.2,pi = 30 MPa,qw = 0.1 m3/d,μ = 1 mPa∙s,B = 0.985,rw = 0.05 m,xf = 20 m,zf = 20 m。如图 4 所示,建议模型的压降和压力导数与商业软件计算的结果非常吻合,这证明了建议模型的有效性。然而,Bennett 等人的模型[19]的计算结果与其他两个模型明显不同,这说明考虑非均匀流量分布的必要性。使用所提出的方法,可以很容易地构建一个两层储层,以确定内部油井生产过程中观察到的流态。随后,对生产动态特性的几个参数进行了敏感性分析,包括井筒存储系数、存储率、渗透率、裂缝高度、裂缝电导率和渗透率各向异性。本节使用的详细参数如下:无量纲压力降和压力导数曲线以对数坐标绘制,如图 5 所示,其中我们可以很容易地识别出四种不同类型的流态。(1) 井筒储存状态[10]。井筒存储状态出现在生产初期,反映了井筒内流体弹性能量的大小。压降曲线和压力导数曲线重合,斜率值为 1。 (2) 线性流动状态 [13]。线性流态是确定压裂导流系数的关键数据。在此期间,流体将以线性方式进入井筒。压力导数曲线的斜率为 1/2。 (3) 双线性流动机制[20]。顾名思义,在压裂-井筒线性流动的同时,基质中还有一股线性流进入压裂。
{"title":"Production Dynamic Characteristic of Fractured Wells in Multilayer Reservoirs Considering the Effect of Non-Uniform Flux Distribution","authors":"Yanzhong Liang, Bailu Teng, Wanjing Luo","doi":"10.2113/2024/lithosphere_2023_173","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_173","url":null,"abstract":"Hydraulic fracturing stimulation, which improves matrix permeability and reduces production costs, has been extensively used in the exploitation of multilayer reservoirs. However, little research on the production dynamic characteristics of vertically fractured wells in stratified reservoirs has been done in the literature. The influence of flux variation along the fracture on the pressure transient behavior has been ignored in these previous works. Therefore, this paper introduces a novel semi-analytical model for fractured wells in multilayer reservoirs, in which the finite difference method is used to characterize fluid flow in the fracture and the Green’s function method is used to characterize fluid flow in the matrix. With the aid of the model, the production dynamic characteristics of fractured wells in multilayer reservoirs can be readily investigated. In addition, based on the assumption of nonuniform flux distribution along the fracture, we successfully recognize four flow regimes occurring in the pressure drop and pressure derivative curves. Following that, the influences of several parameters on the pressure dynamics and layered flux contribution are studied. The calculation results indicate that a larger storability ratio, as well as a larger permeability ratio, can increase the values of the pressure drop and the pressure derivative; the greater the fracture height, the greater the fluid flow into each layer of the fracture. During the production of this model, increasing the fracture conductivity can reduce the pressure drop and pressure derivative, which means lower flow resistance in the fracture.With the growing dependence on fossil energy, deep and ultra-deep areas have steadily evolved into the next main potentials of resource exploration and development. Recently, China has consistently discovered a huge number of deep-layer reservoirs, such as the Puguang, Tahe, Shunbei, and Anyue oilfields, showing great resource potentialities and considerable economic benefits [1]. Considering the influence of the complex sedimentary environment, most deep reservoirs are composed of several layers with different stratigraphic characteristics. Commingling production is commonly adopted to increase producing profit for stratified reservoirs in oil and gas fields [2]. Given this, extensive literature was related to the pressure dynamics of a vertical well in stratified reservoirs [3-6]. Rahman and Mattar [7] derived a new analytical solution for the commingled-layered reservoir with unequal initial pressures in the Laplace domain. Onwunyili and Onyekonwu [8] developed a coupled model that can more accurately simulate the commingled production behavior of multilayer reservoirs. Shi et al. [9] investigated the impact of the vertical inhomogeneous closed boundary radii on pressure transient behaviors of the multilayered commingled reservoir. These previous researches give us a basic understanding of the production dynamic characteristics for th","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140151178","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
Lithospheric Conductivity Structure in the Middle Segment of the Jiangnan Orogenic Belt: Insights into Neoproterozoic Tectonic–Magmatic Processes 江南造山带中段岩石圈传导结构:对新元古代构造-岩浆过程的启示
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_325
Jiayong Yan, Hui Chen, Juzhi Deng, Hui Yu, Yuexin You, Yidan Wen, Min Feng
The Jiangnan Orogenic Belt (JOB) evolved from the Yangtze and Cathaysia blocks through multi-stage oceanic-continental subduction, collisional orogeny, and intracontinental deformation, which is an important region to study the formation and evolution of the South China Continent (SCC). Magnetotelluric soundings were collected along a 520-km-long northwest (NW)-trending profile across the middle segment of the JOB to explore the possible remnants of ancient tectonic–magmatic processes beneath the central SCC by combining with the satellite gravity and magnetic data. The resistivity model reveals that the crust in the middle segment of the JOB and its adjacent area is characterized by high resistivity anomalies, while the uppermost mantle is characterized as medium resistivity anomalies and separated by several subvertical, lithospheric-scale conductors. Two trans-crust anomalies of high conductivity and low density beneath the Jiujiang–Shitai Buried fault (JSBF) and Jiangshan–Shaoxing fault (JSF) extend south-eastward to the lithosphere, which are attributed to the NW and southeast boundaries of the middle segment of the JOB. The imaged NW-trending of JSF reflects the tectonic process of the JOB subducting under the Cathaysia Block. Two lower-crustal conductors also imaged beneath the Jiuling area are interpreted as the partial melting of the lower crust, which may be related to the deep southeast subduction of the Paleo-south China Ocean during 970 to 860 Ma. In addition, the trans-lithosphere high conductivity adjacent to the ancient collisional zone of the Jinning period II (ACZII) is probably related to the asthenosphere upwelling caused by the soft collision between the Yangtze and Cathaysia Blocks, which triggered the contemporaneous magmatism in the Jiuling area. This work provides a new insight into the lithospheric evolution in SCC during the Neoproterozoic.The South China Continent (SCC) is located at the junction of the Paleo-Asian Ocean, Tethys, and Pacific tectonic domains, bordered by the North China Block to the north, the Indochina Block to the south, the Qinghai-Tibet Plateau to the west, and the West-Pacific Plate to the east [1] . Its present status comes from the composite evolution of multi-stage plate tectonics in the Paleo-south China Block, making it one of the most complex geological evolution history areas since the Neoproterozoic [2, 3]. The Jiangnan Orogenic Belt (JOB) in the middle of SCC is spread in a NE-NEE direction, with the Yangtze Block on the northwest (NW) and the Cathaysia Block on the southeast (Figure 1). This area records the collisional assembly of these two ancient microplates, which is of great significance for understanding the crustal accretion, tectonic evolution in the SCC, and the breakup of supercontinent Rodinia [4, 5].Previously proposed models for the tectonic evolution of the JOB include (1) plate subduction collision model [6, 7]; (2) plume model [8]; and (3) plate-rift model [9]. The first
然而,由于 JOB 中段厚厚的沉积覆盖层,没有岩浆岩裸露出来,因此无法深入了解研究区域地下的构造-岩浆过程。因此,对 JOB 下的地下结构进行地球物理成像,是了解 SCC 构造-岩浆过程的关键。然而,对于SCC,之前报道的MT观测要么集中于整个SCC,要么主要局限于SCC的西部[12, 13]。在本研究中,我们介绍了沿穿越 JOB 中段的西北剖面进行 MT 勘测的结果。我们采用三维(3-D)反演方法确定岩石圈的电性结构,并研究了主要解析电性特征的可能起源及其构造影响。通过结合区域、地质和地球物理数据,我们打算对构造边界及其在 JOB 中段之下的深部延伸以及长江和国泰区块之间的深部过程提出新的见解。JOB分布于长江和国泰两块断块的交界处,西起广西北部,经贵州省、湖南省西部、江西省、安徽省南部、浙江省西部至浙江省北部,绵延1500多千米,宽约200千米(图1)。JOB大致可分为三部分:东北段(怀玉岩或双溪武岩系)、中段(九岭山岩系)和西南段(湖南北部和广西北部岩系),主要由中新生代低品位变质地层、新元古代花岗岩、少量黑云母岩和两套公认的蛇纹岩带组成[4,11]。其中,露头地层主要由南华火山沉积岩、未发生变质作用或低级变质作用的南华系巨厚裂谷沉积岩以及新元古代至早古生代稳定的浅海沉积序列组成。褶皱基底包括双溪乌组、溪口组、双桥山组、泗浦组、冷家溪组和板溪组,主要由一套绿斜面变质片岩和箭状岩、夹长凝灰岩、矽卡岩和基性超基性岩组成[21-23]。这里的新元古代岩浆岩的年龄主要集中在950至750Ma,在820Ma左右形成一个重要的峰值。其岩性主要为铝质花岗岩和少量硅质玄武岩、黑云母-超黑云母岩,形成了西南段的蚌洞、三房、元宝山侵入体,中段的九岭侵入体和东北段的许村、歙县、休宁侵入体。侏罗纪中、晚侏罗世以来的一系列构造运动,包括区域性深断层的发育及其向隆起和凹陷区块的剖分,主要控制了侏罗纪中、晚侏罗世以来的构造格局[22, 24, 25]。该地区的主要断层包括西北部近 NE 向分布的九江-石台埋藏断层(JSBF)和安化-罗城断层,中部近 NEE 向分布的宜丰-景德镇断层和东南部近 NE 向分布的江山-绍兴断层(JSF)[3, 23, 26]。JOB的构造变形主要包括晋宁期长江块体与国泰期块体碰撞形成的逆S形构造体系、加里东期华南断裂闭合形成的复合构造体系和印支-燕山期大陆内造山运动形成的复合构造体系[11]。目前,学术界对JOB东北段的边界位置、成因环境和区域构造演化有了全面的认识[2,4],但对中段和南段的研究程度相对薄弱,尤其是对中段的构造边界位置、深部边界延伸和岩石圈俯冲过程仍存在争议[11,13]。近年来,利用地球物理方法在SCC开展了大量深部探测研究,取得了丰硕成果。Han 等人基于横跨景德镇和温州的 MT 剖面的二维反演结果,对该地区的地壳运动进行了深入研究。
{"title":"Lithospheric Conductivity Structure in the Middle Segment of the Jiangnan Orogenic Belt: Insights into Neoproterozoic Tectonic–Magmatic Processes","authors":"Jiayong Yan, Hui Chen, Juzhi Deng, Hui Yu, Yuexin You, Yidan Wen, Min Feng","doi":"10.2113/2024/lithosphere_2023_325","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_325","url":null,"abstract":"The Jiangnan Orogenic Belt (JOB) evolved from the Yangtze and Cathaysia blocks through multi-stage oceanic-continental subduction, collisional orogeny, and intracontinental deformation, which is an important region to study the formation and evolution of the South China Continent (SCC). Magnetotelluric soundings were collected along a 520-km-long northwest (NW)-trending profile across the middle segment of the JOB to explore the possible remnants of ancient tectonic–magmatic processes beneath the central SCC by combining with the satellite gravity and magnetic data. The resistivity model reveals that the crust in the middle segment of the JOB and its adjacent area is characterized by high resistivity anomalies, while the uppermost mantle is characterized as medium resistivity anomalies and separated by several subvertical, lithospheric-scale conductors. Two trans-crust anomalies of high conductivity and low density beneath the Jiujiang–Shitai Buried fault (JSBF) and Jiangshan–Shaoxing fault (JSF) extend south-eastward to the lithosphere, which are attributed to the NW and southeast boundaries of the middle segment of the JOB. The imaged NW-trending of JSF reflects the tectonic process of the JOB subducting under the Cathaysia Block. Two lower-crustal conductors also imaged beneath the Jiuling area are interpreted as the partial melting of the lower crust, which may be related to the deep southeast subduction of the Paleo-south China Ocean during 970 to 860 Ma. In addition, the trans-lithosphere high conductivity adjacent to the ancient collisional zone of the Jinning period II (ACZII) is probably related to the asthenosphere upwelling caused by the soft collision between the Yangtze and Cathaysia Blocks, which triggered the contemporaneous magmatism in the Jiuling area. This work provides a new insight into the lithospheric evolution in SCC during the Neoproterozoic.The South China Continent (SCC) is located at the junction of the Paleo-Asian Ocean, Tethys, and Pacific tectonic domains, bordered by the North China Block to the north, the Indochina Block to the south, the Qinghai-Tibet Plateau to the west, and the West-Pacific Plate to the east [1] . Its present status comes from the composite evolution of multi-stage plate tectonics in the Paleo-south China Block, making it one of the most complex geological evolution history areas since the Neoproterozoic [2, 3]. The Jiangnan Orogenic Belt (JOB) in the middle of SCC is spread in a NE-NEE direction, with the Yangtze Block on the northwest (NW) and the Cathaysia Block on the southeast (Figure 1). This area records the collisional assembly of these two ancient microplates, which is of great significance for understanding the crustal accretion, tectonic evolution in the SCC, and the breakup of supercontinent Rodinia [4, 5].Previously proposed models for the tectonic evolution of the JOB include (1) plate subduction collision model [6, 7]; (2) plume model [8]; and (3) plate-rift model [9]. The first ","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140036800","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
Experimental Investigation on the Grouting Performance of Foam-CNT Composite Grouts in Vertical Inclined Fractures Under Flowing Condition 流动条件下泡沫-CNT 复合灌浆材料在垂直倾斜裂缝中的灌浆性能实验研究
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_341
Yaohui Liu, Fang Lv, Zebin Ouyang, Tao Wang
Grouting is an effective method to solve the problem of water inrush in tunnel and underground engineering. However, rock fractures are often simplified as horizontal and smooth fractures in most grouting studies, while studies on vertical inclined fractures are still rare. To investigate the diffusion law in vertical inclined fractures, a vertical inclined fracture grouting simulation device was developed. A new type of cement slurry with low weight and high flowing water resistance was developed by combining carbon nanotube (CNT) slurry with foamed cement. Physical simulation experiments were conducted to investigate various factors (initial flowing water, inclination angle, sand content, and grouting rate) on the sealing efficiency of grouting. Results show that the high foam content has a negative effect on the compressive strength of the slurry, and has a positive effect on the fluidity and water resistance. The optimum ratio of slurry is 30% foam content, 1.0% CNT content, 1.3 water/cement ratio, and 3% additive content. The inclination angle and inclination direction of the fracture have a great influence on the sealing efficiency of grouting. Foam-CNT composite grouts can meet the requirement of flowing water grouting in vertical inclined fractures.Water inrush is a common problem in tunneling and underground engineering. Water inrush will delay the project and result in high personal injury and property damage [1, 2]. Grouting is an effective method to solve the problem of water inrush [3]. Grouting can improve the strength and reduce the permeability of the rock by injecting slurry into the rock fractures [4-7]. Considerable progress has been made in the area of flowing water grouting in recent years [8-10]. However, the theory of grouting still cannot meet the requirements of practical engineering.Many scholars studied the sealing and diffusion law of flowing water grouting [11-14]. Sui et al. [15] investigated the effects of fracture width, initial flowing rate, grouting time, and grouting amount on the sealing efficiency of flowing water grouting through laboratory simulation experiments. Liang et al. [16] demonstrated that the inclination of fracture has a significant effect on the sealing efficiency and the diffusion law of flowing water grouting. Depending on the relation between the direction of flowing water and fracture, rock fractures can be divided into horizontal fractures, horizontally inclined fractures, and vertically inclined fractures (Figure 1). However, in most studies on flowing water grouting, rock fractures are simplified as horizontal fractures [17, 18], and studies on the vertical inclined fracture grouting are still rare.Many scholars have proved that the diffusion law of liquid in inclined fractures and horizontal fractures is completely different [19]. Graf et al. [20] proposed a numerical method to discretize inclined nonplanar two-dimensional (2D) fractures within a three-dimensional (3D) finite element grid
灌浆是解决隧道和地下工程涌水问题的有效方法。然而,在大多数灌浆研究中,岩石裂缝通常被简化为水平和光滑裂缝,而对垂直倾斜裂缝的研究仍然很少。为了研究垂直倾斜裂缝的扩散规律,开发了一种垂直倾斜裂缝注浆模拟装置。通过将碳纳米管(CNT)浆液与发泡水泥相结合,开发了一种重量轻、流动水阻力大的新型水泥浆。通过物理模拟实验研究了各种因素(初始流水、倾斜角、含沙量和注浆率)对注浆密封效率的影响。结果表明,高泡沫含量对浆液的抗压强度有负面影响,而对流动性和抗水性有正面影响。浆液的最佳配比为泡沫含量 30%、CNT 含量 1.0%、水灰比 1.3、添加剂含量 3%。裂缝的倾斜角度和倾斜方向对灌浆的密封效率有很大影响。泡沫-碳纳米管复合灌浆料可以满足垂直倾斜裂缝中流动水灌浆的要求。涌水是隧道和地下工程中常见的问题,会延误工期,造成严重的人身伤害和财产损失[1, 2]。灌浆是解决涌水问题的有效方法[3]。灌浆可通过向岩石裂缝注入浆液来提高岩石强度并降低渗透性 [4-7]。近年来,流水注浆领域取得了长足的进步 [8-10]。许多学者研究了流水注浆的封闭性和扩散规律 [11-14]。隋文等[15]通过实验室模拟实验研究了裂缝宽度、初始流速、灌浆时间和灌浆量对流动水灌浆封孔效率的影响。梁等人[16]的研究表明,裂缝的倾角对流动水灌浆的封孔效率和扩散规律有显著影响。根据流水方向与断裂的关系,岩石断裂可分为水平断裂、水平倾斜断裂和垂直倾斜断裂(图 1)。然而,在大多数流水灌浆研究中,岩石裂缝都被简化为水平裂缝[17, 18],对垂直倾斜裂缝灌浆的研究还很少。许多学者已经证明,液体在倾斜裂缝和水平裂缝中的扩散规律完全不同[19]。Graf 等人[20]提出了一种数值方法,在三维有限元网格内对倾斜非平面二维(2D)裂缝进行离散化处理,用于地下流动和运移模拟。Mustapha 等人[21]提出了一种新方法,在三维断裂介质中离散化倾斜非平面二维裂缝,用于地下流动和传输模拟。一些学者通过数值方法研究了垂直倾斜裂缝中灌浆料的封堵和扩散规律[22]。Mu 等[23]采用多方向扇形(MDS)的空间逐步法(SSM)计算单斜粗糙断裂中的浆液扩散。Li等人[24]采用数值方法研究了海水条件下具有时变粘度的C-S灌浆料在倾斜岩石裂缝中的扩散。结果表明,垂直裂缝的灌浆难度远高于水平裂缝。在实际工程中,合适的灌浆材料是灌浆成败的关键因素[25]。为满足不同工程的灌浆要求,学者们开发了不同功能的灌浆材料[26, 27],如耐久性和强度高的水泥浆[28]、流动性和泵送性强的化学浆[29-31]、灌浆性能均衡的复合浆[32]、环境污染少的环保浆[33]等。水泥基浆液是实际工程中应用最广泛的注浆材料。一些学者通过各种方法改善灌浆材料的性能,如减小粒径[34]、用粉煤灰替代部分水泥[35]、掺入添加剂[36]等。然而,对灌浆材料的研究大多还集中在水平裂缝灌浆方面,对适用于垂直倾斜裂缝的灌浆材料的研究还很少。
{"title":"Experimental Investigation on the Grouting Performance of Foam-CNT Composite Grouts in Vertical Inclined Fractures Under Flowing Condition","authors":"Yaohui Liu, Fang Lv, Zebin Ouyang, Tao Wang","doi":"10.2113/2024/lithosphere_2023_341","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_341","url":null,"abstract":"Grouting is an effective method to solve the problem of water inrush in tunnel and underground engineering. However, rock fractures are often simplified as horizontal and smooth fractures in most grouting studies, while studies on vertical inclined fractures are still rare. To investigate the diffusion law in vertical inclined fractures, a vertical inclined fracture grouting simulation device was developed. A new type of cement slurry with low weight and high flowing water resistance was developed by combining carbon nanotube (CNT) slurry with foamed cement. Physical simulation experiments were conducted to investigate various factors (initial flowing water, inclination angle, sand content, and grouting rate) on the sealing efficiency of grouting. Results show that the high foam content has a negative effect on the compressive strength of the slurry, and has a positive effect on the fluidity and water resistance. The optimum ratio of slurry is 30% foam content, 1.0% CNT content, 1.3 water/cement ratio, and 3% additive content. The inclination angle and inclination direction of the fracture have a great influence on the sealing efficiency of grouting. Foam-CNT composite grouts can meet the requirement of flowing water grouting in vertical inclined fractures.Water inrush is a common problem in tunneling and underground engineering. Water inrush will delay the project and result in high personal injury and property damage [1, 2]. Grouting is an effective method to solve the problem of water inrush [3]. Grouting can improve the strength and reduce the permeability of the rock by injecting slurry into the rock fractures [4-7]. Considerable progress has been made in the area of flowing water grouting in recent years [8-10]. However, the theory of grouting still cannot meet the requirements of practical engineering.Many scholars studied the sealing and diffusion law of flowing water grouting [11-14]. Sui et al. [15] investigated the effects of fracture width, initial flowing rate, grouting time, and grouting amount on the sealing efficiency of flowing water grouting through laboratory simulation experiments. Liang et al. [16] demonstrated that the inclination of fracture has a significant effect on the sealing efficiency and the diffusion law of flowing water grouting. Depending on the relation between the direction of flowing water and fracture, rock fractures can be divided into horizontal fractures, horizontally inclined fractures, and vertically inclined fractures (Figure 1). However, in most studies on flowing water grouting, rock fractures are simplified as horizontal fractures [17, 18], and studies on the vertical inclined fracture grouting are still rare.Many scholars have proved that the diffusion law of liquid in inclined fractures and horizontal fractures is completely different [19]. Graf et al. [20] proposed a numerical method to discretize inclined nonplanar two-dimensional (2D) fractures within a three-dimensional (3D) finite element grid","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140301026","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
Numerical Modeling of Melting Processes During Slab Break-off: Insights Into Tectonic Setting for Massif-Type Anorthosites 板岩断裂期间熔融过程的数值建模:洞察地块型正长岩的构造背景
IF 2.4 4区 地球科学 Q2 Earth and Planetary Sciences Pub Date : 2024-01-12 DOI: 10.2113/2024/lithosphere_2023_344
Qian Yuan
The concept that lithosphere detachment or break-off has long been conceived as a viable mechanism to explain prominent geological phenomena in Earth’s crust and the surface. One of the strengths of slab delamination mechanism is that it can account for the extensive magmatism in active orogenic belts due to the upwelling of the asthenosphere after the slab break-off. However, in the last 20 years, geodynamic simulations show that the inflow of the asthenosphere upon slab break-off is insufficient to cause significant melting of the overriding lithosphere adjacent to the slab window. The primary reasons include the occurrence of slab break-off at a location that is too deep to effectively heat the overriding lithospheric mantle. Another factor is the presence of a thin film of crustal material that is retained during the slab break- off, inhibiting a significant thermal perturbation within the lithosphere. In this work, we couple petrological–thermomechanical simulations with magmatic melting processes to examine the lithospheric melting and surface lithological expression associated with slab break-off. Our work shows that in the early Earth when the mantle temperature is relatively higher, shallow slab break-off can give rise to significant lithospheric melting during the development of slab break-off. Moreover, because the slab becomes weaker in the earlier hotter mantle, it may break-off prior to the stage of continental collision, thus the magmatism it induced may not give a direct constraint on the time of continental collision. Our study has implications for the interpretation of geological and tomography studies in orogenic belts. It also provides insights into reconciling conflicts between geodynamic and geological studies regarding slab break off-induced melting and magmatism.One of the most peculiar lithologies in Earth’s middle age is Proterozoic massif-type anorthosites (PMAs), a plutonic batholith-forming rock type temporally restricted to the Proterozoic [1-4]. Formally, PMAs are composed of at least 90% plagioclase feldspar accompanied by minor mafic silicates and Fe-Ti oxides [5]. PMAs are areally and volumetrically extensive, with the largest PMA being the Kunene Complex in SW Angola, which covers an area of 18,000 km2 [6]. A consensus has been largely reached on the mechanism by which anorthosites were concentrated. They formed through the accumulation of magmatic plagioclase at the top of a magma chamber due to the low density of plagioclase compared to coexisting melt [7]. However, despite their simple mineralogy and have been studied for over a century, the geodynamic setting accounting for PMAs remains hotly debated [1, 2, 8, 9].A variety of tectonic settings have been proposed for PMAs, including Andean-type continental arc, post-orogenic, anorogenic, and continental rift settings [2, 4, 10, 11]. In recent years, a new tectonic regime—slab break-off—has been adopted accounting for the origin of several PMAs in Asia based o
在这里,我们将浅层板块断裂情景作为板块断裂诱发岩浆活动的最低要求,因为如果浅层板块断裂不显示任何岩石圈熔融,我们就不需要在中层或深层板块断裂下考虑它。为此,我们开始通过改变板块速度(1-10 厘米/年)、板块屈服应力(60、200 兆帕)、长度和海洋岩石圈的年龄来探索促进浅板块断裂的条件。在板块断裂的发展过程中,我们特别关注覆岩岩石圈的熔融过程。如下文所示,在目前的地幔温度下,即使是浅层板块断裂也无法使俯冲板块附近的岩石圈地幔产生熔融。接下来,我们在地球早期地幔温度较高的情况下进行了浅层断裂模拟,结果清楚地表明,板块断裂过程中岩浆层的上涌会引起强烈的热扰动,从而熔化覆盖的地幔岩石圈。最后,由于地幔对流代码没有设计熔体分离和化学分化过程的模型,我们进行了pMELTs/pHMELTs热力学建模来模拟后续的岩浆过程。在此,我们详细描述了使用这种特殊设置的浅断裂参考模型(图 2)。由于规定的较高板块速度(10 厘米/年)以及大洋板块和大陆板块之间的薄弱区,具有 20 兆年历史的大洋岩石圈很快开始俯冲(图 2(a)和(e))。2.5 Myr 之后,整个大洋板块发生俯冲,导致两块大陆发生碰撞(图 2(b) 和 (f))。就在大陆碰撞之后和大陆地壳开始俯冲之前,俯冲板块在 50 千米深处迅速变窄,俯冲板块下形成明显的星体层流入(图 2(c) 和 (g))。随后,断裂扩大,导致俯冲板块沿俯冲面脱离(图 2(d) 和 (h))。与先前的研究[25]一致,我们发现在板块断裂后的粘性蠕变过程中,地壳薄膜持续存在。考虑与水有关的流变学并没有使结果发生重大变化。总之,尽管参考模型表现出一种浅层断裂模式,地幔通过板块窗口强劲流入,但没有观察到与板块断裂相关的岩石圈地幔熔化。在另一个模型中,我们将岩石圈地幔的屈服应力从80兆帕提高到150兆帕,更接近较高的估计值[52],我们观察到与参考模型类似的浅层断裂动力学(图3)。由于最大屈服应力较高,板块断裂发生的时间较长,但在大陆碰撞后板块断裂过程中仍表现出强烈的星体层上涌(图 3(b)和(e))。尽管如此,我们还是没有观测到板块断裂过程引起的岩石圈地幔熔融。然而,在地球早期形成 PMA 时,地幔势温度可能比现在高 150-250°C [32]。在这里,模型 3 的特征与参考案例相同,只是地幔温度比参考案例高 150°C,并且我们在开始时相应地增加了地壳地热梯度。与参考案例类似,模型开始后不久,大洋板块发生俯冲(图 4(a))。然而,由于地幔温度升高,板块变得越来越薄弱,并迅速开始变窄(图 4(b))。随着俯冲板块继续变窄,天体层地幔开始大量流入,取代了大陆板块下的变窄板块区域。同时,由于星体层温度较高,加热速度较快,星体层的上涌确实导致了凌空岩石圈地幔的部分熔化(图 4(b) 和插入部分)。随着板块颈缩的持续,当星体层地幔上涌时,会导致更显著的熔化(图4(c)和插入部分),整个熔化过程可持续近0.2 Myr(图4(d)和插入部分)。为了进一步说明岩石圈熔化是由板块断裂时的星体层上涌直接引起的,而不是由基底熔化引起的,我们将凌空岩石圈的固结温度稍微降低了20°C(1100°C)。如图 5 所示,凌覆岩石圈的显著熔化恰恰是沿着因星体层上涌而形成的弯曲区域形成的,而在凌覆岩石圈的其他位置没有观察到熔化。
{"title":"Numerical Modeling of Melting Processes During Slab Break-off: Insights Into Tectonic Setting for Massif-Type Anorthosites","authors":"Qian Yuan","doi":"10.2113/2024/lithosphere_2023_344","DOIUrl":"https://doi.org/10.2113/2024/lithosphere_2023_344","url":null,"abstract":"The concept that lithosphere detachment or break-off has long been conceived as a viable mechanism to explain prominent geological phenomena in Earth’s crust and the surface. One of the strengths of slab delamination mechanism is that it can account for the extensive magmatism in active orogenic belts due to the upwelling of the asthenosphere after the slab break-off. However, in the last 20 years, geodynamic simulations show that the inflow of the asthenosphere upon slab break-off is insufficient to cause significant melting of the overriding lithosphere adjacent to the slab window. The primary reasons include the occurrence of slab break-off at a location that is too deep to effectively heat the overriding lithospheric mantle. Another factor is the presence of a thin film of crustal material that is retained during the slab break- off, inhibiting a significant thermal perturbation within the lithosphere. In this work, we couple petrological–thermomechanical simulations with magmatic melting processes to examine the lithospheric melting and surface lithological expression associated with slab break-off. Our work shows that in the early Earth when the mantle temperature is relatively higher, shallow slab break-off can give rise to significant lithospheric melting during the development of slab break-off. Moreover, because the slab becomes weaker in the earlier hotter mantle, it may break-off prior to the stage of continental collision, thus the magmatism it induced may not give a direct constraint on the time of continental collision. Our study has implications for the interpretation of geological and tomography studies in orogenic belts. It also provides insights into reconciling conflicts between geodynamic and geological studies regarding slab break off-induced melting and magmatism.One of the most peculiar lithologies in Earth’s middle age is Proterozoic massif-type anorthosites (PMAs), a plutonic batholith-forming rock type temporally restricted to the Proterozoic [1-4]. Formally, PMAs are composed of at least 90% plagioclase feldspar accompanied by minor mafic silicates and Fe-Ti oxides [5]. PMAs are areally and volumetrically extensive, with the largest PMA being the Kunene Complex in SW Angola, which covers an area of 18,000 km2 [6]. A consensus has been largely reached on the mechanism by which anorthosites were concentrated. They formed through the accumulation of magmatic plagioclase at the top of a magma chamber due to the low density of plagioclase compared to coexisting melt [7]. However, despite their simple mineralogy and have been studied for over a century, the geodynamic setting accounting for PMAs remains hotly debated [1, 2, 8, 9].A variety of tectonic settings have been proposed for PMAs, including Andean-type continental arc, post-orogenic, anorogenic, and continental rift settings [2, 4, 10, 11]. In recent years, a new tectonic regime—slab break-off—has been adopted accounting for the origin of several PMAs in Asia based o","PeriodicalId":18147,"journal":{"name":"Lithosphere","volume":null,"pages":null},"PeriodicalIF":2.4,"publicationDate":"2024-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140313461","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