Na-K-Mg三角图常被用来评价水—岩平衡状态和区分不同类型的水样，是水化学分析的一个重要方法。基于Giggenbach等的Na-K-Mg三角图、Benjamin等的Na-K温标公式和Fournier等的K-Mg温标公式，利用Matlab平台研发一套可判断水—岩平衡状态和热储温度计算公式选取的独立程序，解决了水化数据联立分析的问题。以海南岛15个温泉（深井）数据为实例绘制Na-K-Mg三角图，得出海南温泉水样均为部分成熟水和未成熟水，所得结论与前人一致。利用Na-K-Mg三角图（不同Na-K温标公式）结果，算得HN14水样的热储温度与实际相符。可见基于Matlab研发的Na-K-Mg三角图绘制程序可信度高，实用性强，可用于评价水—岩平衡状态和热储温标的选取，并可便捷进行其他数据的联立分析。
Na-K-Mg triangular diagram is a very important method in hydro-chemical analysis and is often used to evaluate the water rock equilibrium state and distinguish different types of water samples. Based on the Na-K-Mg triangular diagram of Giggenbach, the Na-K temperature scale formula of Benjamin et al, and the K-Mg temperature scale formula of Fournier et al., this paper develops a set of programs usingthe Matlab platform, which can judge the water-rock equilibrium state and select calculation formulas for the heat storage temperature.The programs solved the problem of simultaneous analysis of hydration data. The Na-K-Mg triangular diagrams were drawn for the data of 15 hot springs (deep wells) in Hainan Island. The results showed that the hot spring water samples in Hainan were partially mature water and immature water, which was consistent with previous results. According to the results of Na-K-Mg triangular diagrams (different Na-K temperature scale formulas), the calculated heat storage temperature of the HN14 water sample is consistent with the actual situation. It can be concluded that the Na-K-Mg triangular diagram drawing program developed based on Matlab has high credibility and strong practicability, which can be used to evaluate the water-rock balance state and select heat reservoir temperature scale and can be used for simultaneous analysis of other data conveniently.
2022,43（5）: 111-119 收稿日期：2022-08-17
DOI：10.3969/j.issn.1003-3246.2022.05.014
基金项目：中国地震局地震科技星火计划：海南岛陆地区孕震水文地球化学背景研究（项目编号：XH202102Y）
作者简介：刘阳(1984-),高级工程师,主要从事地震地下流体监测及水化学方法研究工作。E-mail:liuyang__07@163.com
参考文献：
陈志,杜建国,周晓成,等. 2012年6月30日新源MS 6.6地震前后北天山泥火山及温泉的水化学变化[J].地震,2014,34(3):97-107.
杜建国,刘丛强.同位素地球化学在地震研究方面的作用[J].地震,2003,23(2):99-107.
何案华,王言章. Matlab在地下流体数据分析中的基础应用[J].大地测量与地球动力学,2019,39(9):966-970.
苏鹤军,张慧,史杰.祁连山断裂带中东段地下水地球化学特征研究[J].西北地震学报,2010,32(2):122-128.
孙小龙,王广才,邵志刚,等. 海原断裂带土壤气与地下水地球化学特征研究[J].地学前缘,2016,23(3):140-150.
王欣,漆继红,许模,等. Na-K-Mg三角图修正与Na-K温标选取[J]. 煤田地质与勘探,2019,47(2):121-128.
张磊,刘耀炜,任宏微,等. 水化学分析方法在地下水异常核实中的应用[J].地震,2019,39(1):29-38.
张炜斌,杜建国,周晓成,等.首都圈西部盆岭构造区地热水水文地球化学研究[J].矿物岩石地球化学通报,2013,32(4):489-496.
张颖.海南岛温泉特征及成因研究[D]. 北京:中国地质大学(北京),2019.
Benjamin T, Charles R and Vidale R.Thermodynamic parameters and experimental data for the Na-K-Ca-geothermometer[J]. J Volcanol Geoth Res, 1983, 15(1/3):167-186.
Fournier R O, SoreyM L, Mariner R H, et al. Chemical and isotopic prediction of aquifer temperatures in the geothermalsystem at Long Valley, California[J]. Journal of Volcanology and Geothermal Research, 1979, 5(1/2):17-34.
Fournier R O, PANG Zhonghe. Geothermometer:Analysis of Na-K-Mg relationship in geothermal water[J]. Foreign Geology (Beijing), 1991, (5):34-37.
Giggenbach W F. Geothermal solute equilibria.Derivation of Na-K-Mg-Ca geoindicators[J]. Geochimica et Cosmochimica Acta, 1988, 52(12):2 749-2 765.
Woith H, Wang R J, Maiwald U, et al.On the origin of geochemical anomalies in groundwaters induced by the Adana 1998 earthquake[J]. Chemical Geology, 2013, 339:177-186.
