Characteristics of Geothermal Geology of the Gulu Geothermal Field in Tibet
-
摘要: 西藏谷露地热田于2020年成功揭露到189.2 ℃的高温地热资源,研究谷露地热田地热地质特征,对西藏地热资源勘查开发和揭示高温地热系统成因机理都具有重要的科学与指导意义.本文以构造地质调查为基础,结合最新物探、钻探、水文资料,通过解析构造‒水‒热循环系统内在联系,总结构造控热规律,并建立了谷露地热田地热系统概念模型.谷露地热田内主要发育N-S向(F1、F3)、E-W向(F2、F4)和NE向3组断裂.大气降水与冰山融水沿盆地西侧九子拉‒桑雄断裂下渗,经深循环加热以盆地西缘F1断裂为导热通道向上运移,由于受到南北两侧F2、F4断裂的阻隔作用,热水在F1和F3断裂组成的“Y”字型断裂系统内汇聚集中,形成热储,并在地表沿NE向和N-S向断裂构成的通道系统中运移排泄;盆地基底花岗岩之上覆盖的泉胶砾岩层对热储起到了良好的隔水保温作用;热田内第四系沉积很薄,热储主要赋存于基岩裂隙中,为基岩裂隙型热储.F1与F3构成的“Y”字型断裂系统是热储赋存的主要场所,谷露地热田地热资源勘查开发要以这两条断裂为重点勘查目标.Abstract: The Gulu geothermal field successfully revealed high temperature thermal reservoir of 189.2 ℃ in 2020, therefore, studying geothermal geological characteristics of the Gulu geothermal field can be significant for guiding geothermal development and helpful for revealing the high temperature geothermal system mechanism. Based on the structural geological survey, combined with the latest geophysical, drilling, and hydrological data, this paper summarizes the structural thermal control laws of the Gulu geothermal field by analyzing the internal connection of the structure-water-thermal cycle system, and establishes a conceptual model of the geothermal system. The faults of the Gulu geothermal field can be divided into 3 groups: N-S trending (F1, F3), E-W trending (F2, F4) and NE trending. Atmospheric precipitation and iceberg meltwater migrated deep along the Jiuzila-Sangxiong fault, and were heated by deep heat sources to form high-temperature geothermal fluids, and rose along the F1 fault at the edge of the basin. The hot water was blocked by the F2 and F4 faults near surface, converging into thermal reservoir, migrating and draining along the channel system which formed by the N-S and NE faults; The conglomerate cemented by sinter overlying the basement granite played a role in ensuring water and thermal insulation for thermal reservoir. The Quaternary sediments were very thin in geothermal field, and the thermal reservoir was mainly stored in the bedrock fissures. According to the results, the "Y"-shaped fault system formed by F1 and F3 faults is the main storage place for thermal reservoir, and these two faults can be regarded as the main target for deep geothermal exploration.
-
0. 引言
温盐环流(thermohaline circulation, 简称THC) 是海洋中全球尺度的循环, 温盐环流受重力作用下高纬较重水团下沉和低纬较轻水团的上升所驱动, 是一种由温度和盐度效应引起的密度梯度驱动的密度流, 其输运量约占全球大洋环流带输运的90%, 当代大洋热盐环流以大西洋输送带环流为特征, 太平洋不存在热盐环流输送带(Broecker, 1991).在大西洋中, 大西洋表层暖水向北流动, 北大西洋深层冷而咸的海水回流, 温盐环流对大西洋中纬向热量和淡水输送起着主要作用, 因此THC的强度对于大西洋热量输送和欧洲的气候有着非常显著的影响(Marotzke, 2000).
受工业化以来人类活动的影响, 大气中温室气体(以CO2为主)的浓度不断增加, 这种变化使从地球表面到太空的辐射效率降低.这是因为, 来自地球表面的大部分向外陆地辐射被大气所吸收, 并在高纬和低温下重新射出.从而导致正的辐射强迫, 使低层大气和地球表面增温, 逃逸到太空的热量减少, 形成增强温室效应.由于温室效应引起的全球变暖会改变淡水通量和地表径流的分布(Manabe and Stouffer, 1994).在全球变暖的情况下, 极地和副极地地区, 由于高纬度降水的增强, 会被输入更多的淡水(Dai et al., 2001); 副热带地区由于蒸发增强将会失去更多的淡水(Voss and Mikolajewicz, 2001).THC的演变可以由表面的浮力扰动触发.在温室气体浓度增加情景下, 全球变暖, 海表温度升高, 大气中的水汽含量增加, 淡水的经向输运加强, 即副热带地区蒸发增强, 副极地地区降雨增加, 从而降低了高低纬度的海水密度梯度, THC的强度随之减弱, 进而影响欧洲气候.
以往的很多研究工作主要注重于海表浮力扰动在海盆尺度上对THC的作用(Manabe and Stouffer, 1994; Voss and Mikolajewicz, 2001).由于THC的下沉分支——北大西洋深层水(NADW) 的形成具有很强的区域性, 主要位于北大西洋拉不拉多海, 丹麦海峡以南区域和格陵兰-冰岛-挪威海.因此分别对不同区域的海洋过程进行分析是十分必要的.本文基于德国Max-Planck气象研究所的大气海洋环流模式, 同时耦合海冰和陆面过程模式, 对地球系统的未来气候进行预测, 分析THC在温室气体浓度增加情景假定下的演变, 着重研究了北大西洋不同海区沿等密度面的海洋过程对全球变暖的响应.
1. 模式介绍、试验设计和分析方法
1.1 模式介绍
本文所采用的模式是德国马普气象所最新发展的气候模式, 海洋模式是Max-Planck Institute Ocean Model (MPI-OM, 嵌套了最新海冰模块), 大气模式是Max Planck Institute Atmosphere Model (ECHAM5).这里我们简要介绍这一气候模式的海洋和大气分量: 海洋模式采用了移动坐标曲线正交网格, 2个坐标投影点一个位于格陵兰岛, 一个位于南极.通过这样的设计, 我们达到了在北大西洋深层水生成区进行网格加密的目的, 以通过较高的分辨率来细致地刻画深层水的形成, 海洋模式水平网格精度从20 km到180 km, 平均精度1.5°×1.5°, 垂直分40层.大气模式采用T63的网格, 水平网格分辨率为1.875°×1.875°, 垂直分31层.此耦合模式没有Flux Adjustment, 这是此模式优点之一.ECHAM5较之以前应用比较广泛的ECHAM4, 主要改进是更新了可预报的气溶胶模块, 对云覆盖重新进行了参数化过程, 同时对云里面的冰和水进行了不同的过程处理, 大大提高了对降雨过程的模拟; 同时对陆面过程的参数化进行了改进, 引入了新的vegetation ratio, leaf area index和forest ratio(Roeckner, 2003).MPI-OM相对于以前版本MPI-OM的改进, 主要是海表风应力的计算改进, 考虑了局地海流对风应力的影响作用, 对上表层海洋的动力机制刻画更为细致, 同时改进了海冰模块, 对海冰的生长、消亡和输运的刻画更为细致(Marsland, 2003).
1.2 试验设计
情景场: 在数值研究中, 我们采用了IPCC (the intergovernmental panel on climate change) 排放情景特别报告中设计的CO2等温室气体情景浓度变化(图 1).在图 1中, 2000年以前的黑实线代表实测和替代资料反演所得CO2情景浓度.A2情景, 假设了一个这样的未来世界, 即经济快速增长, 全球人口快速增长, 并未采用新的能源技术, 温室气体迅速增加.A1B情景, 全球经济和人口迅速增加, 采用一种均衡的燃料使用策略(均衡是指不过分地依赖某一种能源, 其基本假设是在所有的能源供应和终端利用技术方面具有相似的改进速率), 温室气体的浓度增加相对于A2情景有所减弱.B1框架和情景系列假设全球人口在21世纪中达到高峰然后下降, 但其经济结构快速转向服务业和信息产业, 材料消耗强度减少, 清洁高效资源技术得到利用, 温室气体增加缓慢并趋于稳定, 是一个乐观的假设情景.
基于上述模式配置, 我们在不同的温室气体情景场下, 在德国气候中心的巨型计算机上对海洋大气模式进行了耦合运行.首先模式进行了500年的“spin-up”, CO2浓度固定在工业化前的水平280×10-6.然后模式开始启动, 模拟时间为从1860年(人类社会工业化前) 到21世纪末, 结果以年平均的形式输出, 下面我们将就所得的结果进行初步分析.
1.3 分析方法
把大西洋的海水以密度为标准分成三层来分析耦合模式的输出结果, 这3个密度层分别代表大西洋表层水(UW)、深层水(DW) 和底层水(BW).在大西洋里表层暖水向北流动, 深层水和底层水向南流动.随着全球变暖海洋的密度结构会发生变化, 不同密度层的密度标准随之也发生变化.同以往采用固定的密度层的密度标准的方法不同, 本文采用了一个新的方法来定义不同密度层的变化的密度标准.选取大西洋50°N作为断面, 这主要是由于THC在50°N强度最大, 便于分析.沿着50°N断面, 把经向速度场转换到密度坐标上.位势密度是由模式输出的温盐数据以2 000 m为参考面计算得到的.然后把经向速度场按面积纬向积分, 得到每个密度层的经向输运, 再把每个密度层的经向输运从表层分别累加到每一个密度层.向北累积输运最大的密度坐标点定义为σ1, 大于σ1且经向累积输运为零的密度坐标点定义为σ2.密度 < σ1是表层水, σ1 < 密度 < σ2是深层水, 密度 > σ2是底层水.由此我们得到了全球变暖下不同密度层变化的密度标准.全球变暖的情况下, 表层海水吸收的热通量相对较多, 因此增温相对于深层水和底层水比较快.表层水的位势密度在全球变暖的21世纪由1 036.72降到1 036.28.而深层水和底层水位势密度降低的幅度没有表层水那么大.3个密度层之间穿越等密度面的输运是由从海底积分到每个密度层沿等密度面输运的散度得到的.
2. CO2浓度增加情景下THC的演化
2.1 THC强度的变化
图 2a是全球平均海表温度(SST) 在不同的温室气体浓度增加情景下随时间的变化曲线(其中, 绿线表示1860年到2000年用实测和替代资料反演所得CO2情景浓度, 2000年后固定为2000年的CO2情景浓度的模拟结果).全球平均海表温度的计算, 选取全部海洋网格点, 按照网格面积对SST进行加权平均, 输出年平均数据.对2000年后的年平均数据作线性回归, 从2000年到2100年, 在B1、A1B、A2三个情景下, 全球海表面的平均温度分别上升了2.5 ℃, 3.5 ℃, 4.0 ℃, 平均年变化率分别达0.025 ℃·a-1, 0.035 ℃·a-1, 0.04 ℃·a-1, 相关系数都超过0.05的显著性水平.全球海洋环流的深层水生成区域位于北大西洋副极地地区, 在温室气体浓度升高全球变暖的情景下, 该区淡水通量增加, 温度升高, 密度减小, 进而使高低纬密度差减小, 深层水生成率降低, 致使温盐环流THC减弱.图 2b是THC强度在不同的温室气体浓度情景下的随时间的变化曲线, 从2000年到2100年的一个世纪中, 在B1、A1B、A2三个情景下, THC强度分别减弱了4 Sv、5.1 Sv、5.2 Sv, 大体相当于减弱了20%、25%、25.1%.THC强度的减弱趋势十分明显, 同时年代际振荡十分强烈.谱分析的结果显示在95%的置信度上周期为20年的振荡是显著的, 这同Delworth et al. (1993) 和Cheng (2000)的模式结果是类似的.
2.2 NADW生成率的变化
NADW的主要生成区域位于北大西洋拉不拉多海, 丹麦海峡以南区域和格陵兰-冰岛-挪威海(Broecker, 1991).为了便于有针对性的研究, 把NADW生成区域分成三部分: 拉不拉多海区域(50°~65°N, 45°W以西)、丹麦海峡以南区域(SDSR, 50°~65°N, 45°W以东) 和格陵兰-冰岛-挪威海(GIN海域).
图 3给出在A1B情景下20世纪(3a, 3c) 和21世纪(3b, 3d) 北大西洋沿等密度面输运和穿越等密度面输运.总体上来说, 随着全球变暖上层海水的副极地环流得到加强.上层海水穿越50°N的北向输运增加了大约3 Sv, 相对于20世纪大体增加了10%.上层海水通过法鲁海峡进入GIN海域的入流增加了大约1 Sv, 相当于20世纪法鲁海峡入流的10%.根据Mauritzen (1996a, 1996b)的理论, 这种入流的增加将会加强GIN海域的逆时针环流, 通过丹麦海峡和法鲁海峡的溢流水将会增加.
穿越等密度面输运的变化在不同海区情况有所不同.在拉不拉多海和SDSR海域穿越等密度面的输运随着CO2浓度的增加呈明显减弱的趋势, 在拉不拉多海穿越等密度面的输运从8.2 Sv降到5.2 Sv, 在SDSR海域穿越等密度面的输运从2.1 Sv降到1.6 Sv.在GIN海域穿越等密度面的输运有所增加, 从5.7 Sv增加到6.7 Sv.穿过丹麦海峡和法鲁海峡溢流水的总和从5.9 Sv增加到6.2 Sv.NADW是由三部分组成的: 拉不拉多海的深层对流; 丹麦海峡和法鲁海峡的溢流水; 以及SDSR海域的深层水的卷夹.在20世纪NADW的生成率大约是16.2 Sv, 其中8.2 Sv来自拉不拉多海深层对流, 5.9 Sv来自丹麦海峡和法鲁海峡溢流水, 2.1 Sv是SDSR海域的深层水卷夹的结果.在全球变暖的21世纪NADW的生成率降低到12.9 Sv, 其中拉不拉多海贡献5.2 Sv, 丹麦海峡和法鲁海峡溢流水贡献6.1 Sv, SDSR海域的深层水的卷夹贡献1.6 Sv.
2.3 NADW生成率变化的原因
正如前面2.1所讨论的, 在温室气体浓度升高全球变暖的情景下, 极地和副极地地区由于高纬度降水的增加将会被输入更多的淡水, 海表盐度降低; 同时由于全球变暖, 北大西洋海表温度(SST) 将升高.盐度和温度的共同作用导致表层密度降低, 进而拉不拉多海和SDSR海域的垂直层结更加稳定, 垂直对流和深层卷夹减弱.这就是以上两个海区NADW生成率在全球变暖情况下降低的原因.
在GIN海域NADW生成率的变化同以上两个海区是不同的, 变化趋势是升高的.在全球变暖的情况下, 副极地风场是加强的.同时法鲁海峡周围的西南风也是加强的, 因此有更多的北大西洋表层水通过法鲁海峡进入GIN海域.根据Mauritzen (1996a, 1996b)的理论, 丹麦海峡和法鲁海峡溢流水的来源是通过法鲁海峡的入流水.因此法鲁海峡北大西洋入流的增加将导致两个海峡溢流水的增加.法鲁海峡入流的增加将把更多较咸的大西洋水带到GIN海域, 因此GIN海域表层盐度增加.虽然在全球变暖的背景下GIN海域的表层温度是升高的, 但是GIN海域温度升高对密度的减小作用小于盐度升高对密度的增加作用, 最终总的结果是在GIN海域深层对流发生区域表层密度增加.表层密度的增加降低了垂直层化的稳定性, 导致GIN海域在全球变暖情况下深层对流加强.这也导致了丹麦海峡和法鲁海峡溢流水的增加.
3. 结论
本文利用德国马普气象研究所为IPCC第4次评估报告而最新发展的气候模式(ECHAM5/MPIOM), 对3种不同的温室气体排放假设(B1, A1B, A2) 进行了可靠的数值模拟.在此基础上, 本文着重就大西洋温盐环流和北大西洋深层水形成的变化, 以及北大西洋不同海区的温盐环流对温室气体浓度增加的响应, 对模拟结果进行了深入分析.研究揭示, 到21世纪末, 在3种CO2排放情景下, THC强度呈现减弱的趋势, CO2排放增加的越快, THC强度减弱的也就越快.对模式结果的分析表明穿越等密度面的输运的加强(减弱) 同上层海洋层化稳定性的减弱(加强) 密切相关.在拉不拉多海和SDSR海域, 净淡水输入(降水减掉蒸发) 的增加和表层海水的增温导致上层海洋层化稳定性加强, 抑止了这些海区的深层对流.在GIN海域, 通过法鲁海峡的北大西洋入流水的增加导致深层对流的加强和丹麦以及法鲁海峡溢流水的增加.在全球变暖的A1B情景下, NADW生成率总的变化是从16.2 Sv降低到12.9 Sv.
-
图 3 谷露地热田E-W向断裂特征
a. F2断裂地貌特征;b. F4断裂地貌特征;c. AMT剖面L1;位置见图 2a
Fig. 3. Characteristics of E-W trending faults in Gulu geothermal field
图 4 谷露地热田N-S向断裂特征
a. F1断裂错断中更新统冰水堆积层;b. F1断裂错断马里组砂岩;c. 断层面上斜向擦痕;d. 断层滑动矢量应力反演图;e. F3断裂地貌特征;f. 热水泉眼沿F3断裂线状出露;位置见图 2a
Fig. 4. Characteristics of N-S trending faults in Gulu geothermal field
图 5 谷露地热田NE向断裂特征
a. 肉红色花岗岩中共轭节理;b. 2号泉华台表面NE向裂隙,热水通道;c. 1号泉华台表面裂隙统计玫瑰花图;d. 2号泉华台表面裂隙统计玫瑰花图;e. 泉华台西侧NE向沟谷延伸通过泉华台地,位置见图 2a
Fig. 5. Characteristics of NE trending faults in Gulu geothermal field
图 6 谷露地热田热水活动与泉华分布
a.泉华台表面沸泉;b.桑曲西岸泉胶砾岩;c.桑曲河床底部红褐色泉胶砾岩;d.层状硅华;e.两级泉华台地地貌特征与地质剖面,位置见图 2a
Fig. 6. The distribution of hot water activity and sinter in Gulu geothermal field
图 8 谷露地热田测温井钻孔柱状图
位置见图 2b
Fig. 8. Histogram of boreholes of temperature measurement wells in Gulu geothermal field
-
Armijo, R., Tapponnier, P., Mercier, J. L., et al., 1986. Quaternary Extension in Southern Tibet: Field Observations and Tectonic Implications. Journal of Geophysical Research Atmospheres, 91(B14): 13803-13872. https://doi.org/10.1029/jb091ib14p13803 Bian, S., Yu, Z. Q., Gong, J. F., et al., 2021. Spatiotemporal Distribution and Geodynamic Mechanism of the nearly NS⁃Trending Rifts in the Tibetan Plateau. Journal of Geomechanics, 27(2): 178-194 (in Chinese with English abstract). Brown, L. D., Zhao, W. J., Nelson, K. D., et al., 1996. Bright Spots, Structure, and Magmatism in Southern Tibet from INDEPTH Seismic Reflection Profiling. Science, 274(5293): 1688-1690. https://doi.org/10.1126/science.274.5293.1688 Chen, Y. J., Gao, J. C., Feng, J. J., et al., 1992. A Preliminary Study on the History of Hydrothermal Activities Occurred in Southern Tibetan. Hydrogeology and Engineering Geology, 19(5): 18-21, 24 (in Chinese with English abstract). Duo, J., 2003. The Basic Characteristics of the Yangbajing Geothermal Field-A Typical High Temperature Geothermal System. Engineering Science, 5(1): 42-47 (in Chinese with English abstract). Feng, Z. J., Zhao, Y. S., Zhou, A. C., et al., 2012. Development Program of Hot Dry Rock Geothermal Resource in the Yangbajing Basin of China. Renewable Energy, 39(1): 490-495. https://doi.org/10.1016/j.renene.2011.09.005 Guo, Q. H., Wang, Y. X., Liu, W., 2007. Major Hydrogeochemical Processes in the Two Reservoirs of the Yangbajing Geothermal Field, Tibet, China. Journal of Volcanology and Geothermal Research, 166(3-4): 255-268. https://doi.org/10.1016/j.jvolgeores.2007.08.004 Guo, Q. H., Yang, C., 2021. Tungsten Anomaly of the High⁃Temperature Hot Springs in the Daggyai Hydrothermal Area, Tibet, China. Earth Science, 46(7): 2544-2554 (in Chinese with English abstract). Ha, G. H., Wu, Z. H., Liu, F., 2019. Late Quaternary Vertical Slip Rates along the Southern Yadong⁃Gulu Rift, Southern Tibetan Plateau. Tectonophysics, 755: 75-90. https://doi.org/10.1016/j.tecto.2019.02.014 He, Z. H., Yang, D. M., Wang, T. W., 2007. Geochemistry and Tectonic Setting of Miocene Granitoids near Gulou Area, Gangdise Belt. Journal of Jilin University (Earth Science Edition), 37(1): 31-37 (in Chinese with English abstract). Hu, Z. H., Gao, H. L., Wan, H. P., et al., 2022. Temporal and Spatial Evolution of Hydrothermal Alteration in the Yangbajing Geothermal Field, Xizang(Tibet). Geological Review, 68(1): 359-374 (in Chinese with English abstract). Huang, J. C., Liang, H. J., Gu, X. X., 2021. Development Situation of Geothermal Energy in China and Development Proposals in the 14th Five⁃Year Plan Period. World Petroleum Industry, 28(2): 41-46 (in Chinese with English abstract). Li, Z. Q., 2002. Present Hydrothermal Activities during Collisional Orogenics of the Tibetan Plateau (Dissertation). Chinese Academy of Geological Sciences, Beijing, (in Chinese with English abstract). Liu, J., Shen, X. H., Meng, K., et al., 2009. Preliminery Studies on the Late Quaternary Activity of the Western Margin Fault of Gulu Basin, Tibet. Earthquake, 29(3): 45-53 (in Chinese with English abstract). Liu, M. L., He, T., Wu, Q. F., et al., 2020. Hydrogeochemistry of Geothermal Waters from Xiongan New Area and Its. Earth Science, 45(6): 2221-2231 (in Chinese with English abstract). Liu, Q. S., Wu, Z. H., Ye, P. S., et al., 2005. Isotopic Dating of the Nyainqentanglha Granite and Its Significance. Acta Geologica Sinica, 79(3): 331-337 (in Chinese with English abstract). Liu, Z., 2014. The Forming Mechanism of Typical High⁃Temperature Geothermal Systems in Nimu⁃Naqu Geothermal Belt, Tibet (Dissertation). Chinese Academy of Geological Sciences, Beijing (in Chinese with English abstract). Liu, Z., Lin, W. J., Zhang, M., et al., 2014. Geothermal Fluid Genesis and Mantle Fluids Contributions in Nimu⁃Naqu, Tibet. Earth Science Frontiers, 21(6): 356-371 (in Chinese with English abstract). http://www.researchgate.net/publication/286308204_Geothermal_fluid_genesis_and_mantle_fluids_contributions_in_Nimu-Naqu_Tibet Nelson, K. D., Zhao, W. J., Brown, L. D., et al., 1996. Partially Molten Middle Crust beneath Southern Tibet: Synthesis of Project INDEPTH Results. Science, 274(5293): 1684-1688. https://doi.org/10.1126/science.274.5293.1684 Pang, Z. H., Hu, S. B., Wang, J. Y., 2012. A Roadmap to Geothermal Energy Development in China. Science & Technology Review, 30(32): 18-24 (in Chinese with English abstract). Qin, J. S., 2003. Hydrothermal Alteration and Evalution in Yangyi Geothermal Field. Journal of Taiyuan University of Technology, 34(2): 161-165 (in Chinese with English abstract). Tong, W., Zhang, M. T., Zhang, Z. F., et al., 1981. Geothermals Beneath Xizang(Tibetan) Plateau. Science Press, Beijing, 18-19 (in Chinese). Wang, J. Y., Pang, Z. H., Kong, Y. L., et al., 2020. Status and Prospects of Geothermal Clean Heating Industry in China. Science & Technology for Development, 16(S1): 294-298 (in Chinese with English abstract). Wang, S. G., Chevalier, M. L., Pan, J. W., et al., 2020. Quantification of the Late Quaternary Throw Rates along the Yadong Rift, Southern Tibet. Tectonophysics, 790: 228545. https://doi.org/10.1016/j.tecto.2020.228545 Wu, Z. H., Hu, D. G., Liu, Q. S., et al., 2005. Chronological Analyses of the Thermal Evolution of Granite and the Uplift Process of the Nyainqentanglha Range in Central Tibet. Acta Geosicientia Sinica, 26(6): 505-512 (in Chinese with English abstract). Wu, Z. H., Zhao, X. T., Wu, Z. H., et al., 2006a. Quaternary Geology and Faulting in the Damxung⁃Yangbajain Basin, Southern Tibet. Journal of Geomechanics, 12(3): 305-316 (in Chinese with English abstract). Wu, Z. H., Wu, Z. H., Hu, D. G., et al., 2006b. Holocene Seismogenic Faults along the Tanggula⁃Lhasa Section of the Qinghai⁃Tibet Railway, China. Geological Bulletin of China, 25(12): 1387-1401 (in Chinese with English abstract). Xu, S. G., Guo, Y. S., 2009. Geothermal Fundamentals. Science Press, Beijing, 19-21 (in Chinese). Xu, Z. K., Xu, S. G., Zhang, S. T., 2021. Hydro⁃Geochemistry of Anning Geothermal Field and Flow Channels Inferring of Upper Geothermal Reservoir. Earth Science, 46(11): 4175-4187 (in Chinese with English abstract). Zhang, M., Lin, W. J., Liu, Z., et al., 2014. Hydrogeochemical Characteristics and Genetic Model of Gulu High⁃Temperature Geothermal System in Tibet, China. Journal of Chengdu University of Technology (Science & Technology Edition), 41(3): 382-392 (in Chinese with English abstract). doi: 10.3969/j.issn.1671-9727.2014.03.15 Zhao, W. J., Jiang, W., Wu, Z. H., et al., 2003. Study on Deep Structure⁃Seismic⁃Geothermal Relationship and Mechanism of Yangbajing in Tibet. Chinese Academy of Geological Sciences, Beijing, 95-97 (in Chinese). Zhao, W. J., Zhao, X., Shi, D. N., et al., 2002. Progress in the Study of Deep (INDEPTH) Profiles in the Himalayas and Qinghai⁃Tibet Plateau. Regional Geology of China, 21(11): 691-700 (in Chinese with English abstract). Zhao, Y. Y., Fan, X. T., Han, J. Y., et al., 2009. Geologic and Geochemical Features and Ore Forming Process for Hot Spring Cesium Deposit of Gulu Area, Nagqu Region, Tibet, China. Geological Bulletin of China, 28(7): 933-954 (in Chinese with English abstract). doi: 10.3969/j.issn.1671-2552.2009.07.014 Zhao, Y. Y., Zhao, X. T., Ma, Z. B., et al., 2010. Chronology of the Gulu Hot Spring Cesium Deposit in Nagqu, Tibet and It Geological Significance. Acta Geologica Sinica, 84(2): 211-220 (in Chinese with English abstract). Zheng, K. Y., Zheng, F., 2020. Discussion on Prospects of Geothermal Power Generation Industry in China. Sino⁃Global Energy, 25(11): 17-23 (in Chinese with English abstract). Zheng, M. P., Wang, Q. X., Duo, J., et al., 1995. A New Type of Hydrothermal Deposit⁃Cesium⁃Bearing Geyserite in Tibet. Geological Publishing House, Beijing, 27-32 (in Chinese with English abstract). 卞爽, 于志泉, 龚俊峰, 等, 2021. 青藏高原近南北向裂谷的时空分布特征及动力学机制. 地质力学学报, 27(2): 178-194. https://www.cnki.com.cn/Article/CJFDTOTAL-DZLX202102004.htm 陈以健, 高钧成, 冯锦江, 等, 1992. 藏南的水热活动历史初探. 水文地质工程地质, 19(5): 18-21, 24. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG199205008.htm 多吉, 2003. 典型高温地热系统: 羊八井热田基本特征. 中国工程科学, 5(1): 42-47. https://www.cnki.com.cn/Article/CJFDTOTAL-GCKX200301007.htm 郭清海, 杨晨, 2021. 西藏搭格架高温热泉中钨的水文地球化学异常. 地球科学, 46(7): 2544-2554. doi: 10.3799/dqkx.2020.287 和钟铧, 杨德明, 王天武, 2007. 冈底斯带谷露区中新世花岗岩地球化学特征及构造环境. 吉林大学学报(地球科学版), 37(1): 31-37. https://www.cnki.com.cn/Article/CJFDTOTAL-CCDZ200701004.htm 胡志华, 高洪雷, 万汉平, 等, 2022. 西藏羊八井地热田水热蚀变的时空演化特征. 地质论评, 68(1): 359-374. https://www.cnki.com.cn/Article/CJFDTOTAL-DZLP202201027.htm 黄嘉超, 梁海军, 谷雪曦, 2021. 中国地热能发展形势及"十四五"发展建议. 世界石油工业, 28(2): 41-46. https://www.cnki.com.cn/Article/CJFDTOTAL-SSYY202102009.htm 李振清, 2002. 青藏高原碰撞造山过程中的现代热水活动(博士学位论文). 北京: 中国地质科学院. 刘静, 申旭辉, 孟恺, 等, 2009. 谷露盆地西缘断裂晚第四纪活动初步研究. 地震, 29(3): 45-53. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZN200903006.htm 刘明亮, 何曈, 吴启帆, 等, 2020. 雄安新区地热水化学特征及其指示意义. 地球科学, 45(6): 2221-2231. doi: 10.3799/dqkx.2019.270 刘琦胜, 吴珍汉, 叶培盛, 等, 2005. 念青唐古拉花岗岩的同位素年龄测定及其地质意义. 地质学报, 79(3): 331-337. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200503007.htm 刘昭, 2014. 西藏尼木‒那曲地热带典型高温地热系统形成机理研究(博士学位论文). 北京: 中国地质科学院. 刘昭, 蔺文静, 张萌, 等, 2014. 西藏尼木‒那曲地热流体成因及幔源流体贡献. 地学前缘, 21(6): 356-371. https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201406040.htm 庞忠和, 胡圣标, 汪集旸, 2012. 中国地热能发展路线图. 科技导报, 30(32): 18-24. https://www.cnki.com.cn/Article/CJFDTOTAL-KJDB201232015.htm 秦进生, 2003. 西藏羊易地热田的水热蚀变及地热地质意义. 太原理工大学学报, 34(2): 161-165. https://www.cnki.com.cn/Article/CJFDTOTAL-TYGY200302016.htm 佟伟, 章铭陶, 张知非, 等, 1981. 西藏地热. 北京: 科学出版社, 18-19. 汪集旸, 庞忠和, 孔彦龙, 等, 2020. 我国地热清洁取暖产业现状与展望. 科技促进发展, 16(S1): 294-298. https://www.cnki.com.cn/Article/CJFDTOTAL-KJCJ2020Z1010.htm 吴珍汉, 胡道功, 刘琦胜, 等, 2005. 念青唐古拉花岗岩热演化历史和山脉隆升过程的热年代学分析. 地球学报, 26(6): 505-512. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXB200506003.htm 吴中海, 赵希涛, 吴珍汉, 等, 2006a. 西藏当雄‒羊八井盆地的第四纪地质与断裂活动研究. 地质力学学报, 12(3): 305-316. https://www.cnki.com.cn/Article/CJFDTOTAL-DZLX200603003.htm 吴中海, 吴珍汉, 胡道功, 等, 2006b. 青藏铁路唐古拉山‒拉萨段全新世控震断裂研究. 地质通报, 25(12): 1387-1401. https://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD200612005.htm 徐世光, 郭远生, 2009. 地热学基础. 北京: 科学出版社, 19-20. 徐梓矿, 徐世光, 张世涛, 2021. 安宁地热田浅部热储水化学特征及补给通道位置. 地球科学, 46(11): 4175-4187. doi: 10.3799/dqkx.2020.401 张萌, 蔺文静, 刘昭, 等, 2014. 西藏谷露高温地热系统水文地球化学特征及成因模式. 成都理工大学学报(自然科学版), 41(3): 382-392. https://www.cnki.com.cn/Article/CJFDTOTAL-CDLG201403015.htm 赵文津, 江万, 吴珍汉, 等, 2003. 西藏羊八井深部构造‒地震‒地热关系及机理调查研究. 北京: 中国地质科学院, 95-97. 赵文津, 赵逊, 史大年, 等, 2002. 喜马拉雅和青藏高原深剖面(INDEPTH)研究进展. 地质通报, 21(11): 691-700. https://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD200211000.htm 赵元艺, 樊兴涛, 韩景仪, 等, 2009. 西藏谷露热泉型铯矿床地质地球化学特征与成矿作用. 地质通报, 28(7): 933-954. https://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD200907016.htm 赵元艺, 赵希涛, 马志邦, 等, 2010. 西藏谷露热泉型铯矿床年代学及意义. 地质学报, 84(2): 211-220. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE201002007.htm 郑克棪, 郑帆, 2020. 中国地热发电产业前景探讨. 中外能源, 25(11): 17-23. https://www.cnki.com.cn/Article/CJFDTOTAL-SYZW202011005.htm 郑绵平, 王秋霞, 多吉, 等, 1995. 水热成矿新类型: 西藏艳硅华矿床. 北京: 地质出版社, 27-32. -