Elsevier

Geothermics

Volume 89, January 2021, 101936
Geothermics

Heat generation effects from shear friction along Xianshui river strike-slip fault in western Sichuan, China

https://doi.org/10.1016/j.geothermics.2020.101936Get rights and content

Abstract

The Xianshui river fault zone is located in the western Sichuan Plateau, on the eastern edge of the Himalayan Geothermal Belt; numerous hot springs are exposed at the surface and distributed along the fault zone, forming an area of high-temperature hydrothermal activity. The genetic mechanisms of these high-temperature geothermal systems is currently unclear, which poses an urgent scientific problem in the field of geothermal research This study uses numerical simulation methods based on geological, hydrological, geothermal, and geophysical data to study the effects of shear frictional heat on high-temperature hydrothermal activity in the research area under the conditions of different slip velocities, parameters of thermal materials, and strain rates of the Xianshui river left-slip fault zone. The results show that shear frictional heat is the main contributor of the local thermal anomalies in the high-temperature hydrothermal system of the Xianshui river fault zone. The local high temperature generated by shear friction causes the crust to partially melt at a depth of 20–25 km; this provides the heat source for the striped geothermal system, which is controlled by the Xianshui river fault zone. In addition, the surface heat flow generated by shear friction was measured at 10 Ma, along a two-dimensional section across the Xianshui river fault in an E–W direction, and it was found to be approximately 5–50 mW/m2. Furthermore, the nonlinear relationship between thermal conductivity, thermal capacity, shear stress, and temperature gives the shear frictional heat generation process "positive" and "negative" feedback effects, which maintain the dynamic equilibrium.

Introduction

The high-temperature geothermal belt in western Sichuan is one of the most active geothermal regions in China. It is located in the Tethys tectonic region, which lies in the eastern part of the Himalayan Geothermal Belt (HGB) (Fig. 1a) and is an area with intense tectonic activity. Investigation has indicated that the average value of the heat flow is higher than 80 mW/m2. Recent results of hydrogeochemistry research and drilling work (Cao et al., 2006; Guo et al., 2017; Jiao et al., 2018; Luo et al., 2017) have proven the existence of numerous high-temperature geothermal systems in the west Sichuan high-temperature geothermal belt, which can be subdivided into the Batang, Litang, and Kangding geothermal belts, according to the distribution of hot springs from the west to the east (Tang et al., 2017). All the types of high-temperature surface geothermal manifestations that occur within China can be found in this region, including boiling springs, boiling fountains, steam vents, steaming ground, and a hydrothermal exploration zone (Liao, 1990a; Liao, 1990b), which are distributed with banding structures.

High-temperature geothermal systems are closely related to volcanic and seismic activities and can act as superficial supports and constraints for corresponding geodynamic processes. Therefore, defining the origin of high-temperature geothermal systems is not only important for conducting geothermal and hydrogeological groundwork, but it is also beneficial for resolving geodynamic scientific issues. In addition, modeling and analyzing the generation process of geothermal systems is important for the evaluation and development of geothermal resources. Furthermore, the modern generation of hydrothermal mineralization has an evident relation to geothermal systems. Thus, the results of associated studies can be used as a reference in order to predict the presence of metal and mineral resources.

Although the distribution of the springs in the studied area has been thoroughly explored and detailed, the generation of the geothermal system process has not been completely apprehended. Thus, exploring the heat source is the key to analyzing the problems associated with the generation of geothermal systems. The high crustal heat beneath the HGB has been explained and discussed (Bird, 1978; England et al., 1992; Hochstein and Regenauer-Lieb, 1998). The exposed volcanic rocks, generated more than 30 Ma ago, having anomalously low 3He/4He ratios of geothermal gases that were collected from several high-temperature systems within the HGB show that the high heat flow in the HGB (including the western Sichuan geothermal systems) has an intra-crustal origin, and is not related to the intrusions of sub-crustal melts that have been occurring during the collisional period (Giggenbach et al., 1983). Research (Regenauer-Lieb, 1992) has shown that plastic deformation of the ductile crust may also have played a role in forming the HGB and the enriched radioactive elements in crust due to the crustal thickening during collision may have increased the radioactive heat flow. The uplift of the thickening crust related to the isostatic compensation has led to the elevation of the isothermal surface and increased the heat flow background value. It has also been proposed that shear heating occurs in the crustal thrust interface between the Indian and Asian plates beneath the HGB(Bird, 1978; England et al., 1992; Shi and Zhu, 1993). Furthermore, heat can be stored within magma — in the form of latent heat when partial melting is induced by fault shear frictional heat — and can affect the heat flow to some extent. Heat from the late-magmatic intrusion (<10 Ma) and released melting latent heat generated from magma solidification can also be considered. Overall, the heat source is mostly located within the crust in the western Sichuan high-temperature geothermal systems.

The focus of this study is the Xianshui river high-temperature geothermal systems, controlled by the Xianshui river strike-slip fault, where the Yulingong in Kangding, Yuemachaqu, and Qimeikechaqu in Daofu hot springs are located along the fault. The well was drilled within Kangding City during 2013, which has shown that temperatures at the top and bottom of the well reach 140 ℃ and 208 ℃ at 267 m, respectively; this classifies the springs as high-temperature systems and implies that the Kangding geothermal systems are possibly significant latent heat resources. Studies on the hydrogeochemical and isotopic characteristics of hot springs and gases in the Kangding geothermal area have also been conducted, and it has been found that the reservoir temperature reaches approximately 260 °C (Guo et al., 2017; Luo et al., 2017). Furthermore, the heat-storage value of the Kangding geothermal systems has also been estimated (Liao, 2018; Wang et al., 2019).

Nevertheless, it is currently unclear how a geothermal system such as Kangding is controlled by the Xianshui river strike-slip fault; therefore, this study conducts quantitative simulations of temperature evolution beneath the Xianshui river strike-slip fault using COMSOL Multiphysics software to elucidate this. The contributions of mantle heat flow, radioactive heat flow, shear frictional heat flow, block uplift heat flow, and melting latent heat flow to surface heat flow, are also respectively calculated. The simulations consider several important developments: first, fault geometry inferred from analysis of historical seismic events since 1870 is used to constrain the model structure; second, different kinematic histories of the Xianshui river strike-slip fault are tested; third, high-temperature and pressure experimental petrology physics data are used to assist in understanding the important feedback processes between temperature, shear heat generation, and the thermal structure. As a typical representation of the western Sichuan geothermal belt, the modeling results for the Kangding geothermal system may provide an insight into the deep dynamic mechanism and thermal controlling structures of the system and other analogous high-temperature geothermal systems in the Tethys tectonic region.

Section snippets

large-scale setting

The western Sichuan plateau is located at the junction of the Qiangtang, Songpan-ganzi, and Yangtze active blocks (Fig. 1a) (Copley et al., 2010; Harrowfield and Wilson, 2005; Zhang, 2008). It is part of the eastern Tibet Plateau that formed during the collision and continuous compression of the Indian and Eurasian plates (beginning approximately 55 Ma ago) (Clark and Royden, 2000; Tapponnier et al., 2001; Xu et al., 2006), which resulted in the formation of numerous deep faults extending down

Model geometric constraints

In this study, we simplified the investigation by focusing on a well-constrained, two-dimensional, vertical, cross-sectional model of the Kangding geothermal system which has a length and depth of 100 km and 30 km, respectively. The position of the cross-section is represented by the black solid line in Fig. 1. For an initial approximation, we simulated the Xianshui river fault friction heating process using a kinematic model, which is broadly consistent with geological and geophysical

Results

According to the key parameters discussed in Section 3, such as geometric constraints, kinematic parameters, and rock rheology, we conducted a series of frictional heat simulations for the Xianshui river fault zone.

Fig. 3a shows the simulated results of surface heat flow values along the section, where the horizontal coordinate of 0 m corresponds to the position of the Xianshui river fault zone, and the different colors represent the surface heat flow values calculated under different

Discussion

Unlike many previous studies, this model achieves a coupling between thermal conductivity, thermal capacity, and the temperature of rock in the crust, which share complex feedback processes between them. In this respect, the fault slip generates shear frictional heat on the contact surface of the fault, and the temperature of the surrounding rock increases. Thermal conductivity then sharply decreases while the heat capacity increases, and the thermal diffusion coefficient further decreases with

Conclusions

Based on geometric constraints, kinematic parameters, rheology, and other parameters, a 2D model of the Xianshui river fault zone is established and a series of thermal simulations are conducted based on the use of different parameter values. Combined with geological, hydrological, geothermal, and geophysical data, the heat generation effects of shear friction along the Xianshui river strike-slip fault on the high temperature geothermal activity in the studied region is also discussed. In

Funding

This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences, Grant No. XDB42020104 and the National Natural Science Foundation of China [No. 41430319, 41574074].

CRediT authorship contribution statement

Yifei Ai: Conceptualization, Formal analysis, Methodology, Software, Writing - original draft. Jian Zhang: Conceptualization, Writing - review & editing. Miao Dong: Validation. Beiyu Wang: Investigation. Gui Fang: Investigation.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References (83)

  • Y.Z. Zhang et al.

    Cooling history of the Gongga batholith: Implications for the Xianshuihe Fault and Miocene kinematics of SE Tibet

    Earth and Planetary Science Letters

    (2017)
  • G.P. Airy

    On the computations of the effect of the attraction of the mountain masses as disturbing the apparent astronomical latitude of stations in geodetic surveys

    Philos.Trans. R. Soc. London Ser. A

    (1855)
  • C.R. Allen et al.

    Field study of a highly active fault zone: The Xianshuihe fault of southwestern China

    Geological Society of America Bulletin

    (1991)
  • D. Bai et al.

    Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging

    Nature Geoscience.

    (2011)
  • D.H. Bai et al.

    Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging

    Nature Geoscience

    (2010)
  • C. Beaumont et al.

    Himalayan tectonics explained by extrusion of a low-viscosity crustal channel coupled to focused surface denudation

    Nature

    (2001)
  • P. Bird

    Initiation of intracontinental subduction in Himalaya

    Journal of Geophysical Research

    (1978)
  • M.L. Blanpied et al.

    Frictional slip of granite at hydrothermal conditions

    Journal of Geophysical Research Solid Earth

    (1995)
  • Y. Cao et al.

    Comparision of geochemical features pf warm springs between Chongqing and Kangding.

    (2006)
  • J. Chen et al.

    Surface rupture zones of the 2001 earthquake Ms8.1 west of Kunlun pass, northern Qinghai-Xizang plateau

    Quaternary Sciences. (in Chinese)

    (2003)
  • L.C. Chen et al.

    The Ms7.1 Yushu earthquake surface ruptures and historical earthquakes

    Chinese Sci Bull (in Chinese)

    (2010)
  • W. Chen et al.

    Late Cenozoic episodic uplifting in southeastern part of the Tibetan plateau-evidence from Ar-Ar thermochronology

    Acta Petrologica Sinica (in Chinese)

    (2006)
  • M.K. Clark et al.

    Topographic ooze: Building the eastern margin of Tibet by lower crustal flow: Geology, 28

    no.

    (2000)
  • A. Copley et al.

    India-Asia collision and the Cenozoic slowdown of the Indian plate: Implications for the forces driving plate motions

    Journal of Geophysical Research-Solid Earth

    (2010)
  • Q.D. Deng et al.

    Seismic activities and earthquake potential in the Tibetain Plateau

    Chinese J. Geophys. (in Chinese)

    (2014)
  • P. England et al.

    Heat-sources for Tertiary metamorphism and anatexis in the Annapurna-Manaslu region central Nepal

    Journal of Geophysical Research-Solid Earth

    (1992)
  • L.J. Gao et al.

    The study of gravity-magnetic anomaly and tectonic background in Sichuan west region

    Chinese J. Geophys. (in Chinese)

    (2015)
  • J.F. Hayford et al.
    (1912)
  • M.X. He et al.

    Deep conductivity characteristics of the southern Xianshuihe fault zone Chinese

    J. Geophys. (in Chinese)

    (2017)
  • W. Heiskanen
    (1938)
  • F. Holtz et al.

    The water-undersaturated and dry Qz-Ab-Or system revisited

    Experimental results at very low water activities and geological implications: Contributions to Mineralogy and Petrology

    (2001)
  • J.L. Huang et al.

    Pn-wave velocity structure and anisotropy in the upper mantle of sichuan, yunnan and adjacent regions

  • T. Jiao et al.

    Geochemistry of geothermal fluids with implications on the sources of water and heat recharge to the Rekeng high-temperature geothermal system in the Eastern Himalayan Syntax

    Geothermics

    (2018)
  • E. Kirby et al.

    Distribution of active rock uplift along the eastern margin of the Tibetan Plateau: inferences from bedrock channel longitudinal profiles

    Journal of Geophysical Research

    (2003)
  • E. Kirby et al.

    Late Cenozoic evolution of the eastern margin of the Tibetan Plateau: Inferences from Ar-40/Ar-39 and (U-Th)/He thermochronology

    Tectonics

    (2002)
  • H.B. Li et al.

    Co-seisimic surface rupture and dextral-slip oblique Thrusting of the Ms8.0 Wenchuan earthquake

    Acta Geologica Sinica. (in Chinese)

    (2008)
  • Z. Liao

    Yunnan-Tibet geothermal belt

    Paper read at The 6th Annual Academic Meeting of the Chinese Geophysical Society, at Wuhan, Hubei, China.

    (1990)
  • Z.J. Liao

    Yunnan-Tibet geothermal belt

    Paper read at The 6th Annual Academic Meeting of the Chinese Geophysical Society, at Wuhan, Hubei, China.

    (1990)
  • Z.J. Liao

    Thermal Springs and Geothermal Energy in the Qinghai-Tibetan Plateau and the Surroundings.

    (2018)
  • J. Luo et al.

    Geothermal potential evaluation and development prioritization based on geochemistry of geothermal waters from Kangding area, western Sichuan

    China: Environmental Earth Sciences

    (2017)
  • L.L. Luo

    Inquisition of the distribution and cause of the hot springs in western Sichuan

    Journal of Chongqing Teachers College (Natural Science Edition) (in Chinese)

    (1994)
  • Cited by (20)

    • Modeling and interpretation of geomagnetic data related to geothermal sources, Northwest of Delijan

      2022, Renewable Energy
      Citation Excerpt :

      Geothermal systems are closely related to seismic events and volcanic activity, which reflect zones of complex geodynamic processes [1].

    • Heat flux measurements and thermal potential of the Garze geothermal area in the eastern Himalayan Syntaxis

      2021, Journal of Volcanology and Geothermal Research
      Citation Excerpt :

      Moreover, the geothermal area should have received additional heat source as its heat flux to exceed that of some geothermal areas on the edge with magmatic activity. Such a replenishment heat source is probably tectonic friction heat (Wang et al., 2019; Ai et al., 2021). It is clear that the heat flux obtained by using eddy covariance tower has a large variation of about 66% on average, such as the Solfatara geothermal area (Fig. 11b).

    View all citing articles on Scopus
    View full text