基于Matlab的Durov三线图的实现

doi:10.12118/j.issn.1000-6060.2018.04.09

摘要/Abstract

摘要： 水化学分类是实现水资源可持续开发利用的重要前提之一。Durov三线图作为主要的水化学分类图因绘制复杂而限制了其在现实中的应用，且一些常见的水化学软件也不具备绘制Durov三线图的功能。Durov三线图包括原始型Durov三线图及在其基础上演变而成的扩展型Durov三线图。通过介绍两种Durov三线图的差异、组成特征和投影原理，提出了在Durov三线图中建立坐标系的方法，并推导了将水化学数据投影到Durov三线图的公式。在此基础上，阐述了基于Matlab语言编程实现原始型Durov三线图和扩展型Durov三线图的方法和流程，并以公开发表的数据绘制了两种Durov三线图。基于相同的数据比较发现，相关水化学软件绘制的原始型Durov三线图及一些公开发表的研究成果中的扩展型Durov三线图与本文所用Matlab语言绘制的Durov三线图相一致，表明该程序准确可行。本程序既可以实现原始型Durov三线图的绘制，也可以实现常见水化学软件不能完成的扩展型Durov三线图的绘制。该程序的应用，将有助于Durov三线图的高效应用，特别是扩展型Durov三线图的应用不但可以作为水化学分类的依据，同时还可以指示某些水化学过程。

关键词: 水化学, Durov三线图, Matlab, 制图

Abstract: Hydrochemical classification of water resource is one of the vital preconditions for water resource sustainable utilization.There are numerous methods with regard to hydrochemical classification for water resource,but the most effective and simplest method is the diagram method.As one of the most important diagram methods for hydrochemical classification,the Durov diagrams (including both the original type of Durov Diagram and the expanded Durov Diagram) are used widely not only because it can reveal the hydrogeochemical characteristics of water,but also because it can show the hydrochemical processes occurring within different hydrogeological systems.However,comparing with other diagrams (such as the Piper Diagram),these diagrams are very complex and difficult to plot.Few commercial software packages are designed to plot Durov diagrams,which limits researchers to use the diagrams.In this paper,we introduced both the original type of Durov Diagram and the expanded Durov Diagram,and we also compared the similarities and differences in them.To better understand the diagram,we interpreted the composition of the original type of Durov Diagram,especially the part of cation triangle,anion triangle and the square field.We also interpreted the composition of the expanded Durov Diagram,in which each of the cation and anion triangles is split into two subtriangles and a diamond,while the square field is split into nine fields.We stressed how to establish a coordinate to plot the original type of Durov Diagram and the expanded Durov Diagram respectively.And we deduced formulas to transform the hydrochemical data so that they can be projected into both the original type of Durov Diagram and the expanded Durov Diagram.Using these principles,we wrote programs for the output of both the original type of Durov Diagram and the expanded Durov Diagram based on Matlab.To test the validity of the programs,we plotted the original type of Durov Diagram with published data based on our Matlab programs.And we also plot the original type of Durov Diagram with the same data based on commercial software.The results showed that our program was as correct as the commercial software,which indicated our program is feasible.However,because none of commercial software packages we found has the function of plotting the expanded Durov Diagram,we compared the diagram plotted by our Maltlab program with diagram appearing in some classical literature that used the same data.The results also confirmed that our Matlab program is reliable.Both the original type of Durov Diagram and the expanded Durov Diagram plotted by Matlab will help researchers to interpret hydrochemical data in graphical methods conveniently.Especially,the expanded Durov Diagram provides a method to indicate hydrochemical processes.

Key words: hydrochemistry, the Durov diagrams, Matlab, plotting

中图分类号:

P342

任孝宗, 刘敏, 张迎珍, 贺志明, 朱秉启. 基于Matlab的Durov三线图的实现[J]. 干旱区地理, 2018, 41(4): 744-750.

REN Xiao-zong, LIU Min, ZHANG Ying-zhen, HE Zhi-ming, ZHU Bing-qi. Plotting Durov Diagram based on Matlab[J]. 干旱区地理, 2018, 41(4): 744-750.

参考文献

[1] 石虹.浅谈全球水资源态势和中国水资源环境问题[J].水土保持研究,2002,9(1):145-150.[SHI Hong.The state of global water resource and the problems of water environment in China[J].Research of Soil and Water Conservation,2002,9(1):145-150.]
[2] 任国玉,姜彤,李维京,等.气候变化对中国水资源情势影响综合分析[J].水科学进展,2008,19(6):772-779.[REN Guoyu,JIANG Tong,LI Weijing,et al.An integrated assessment of climate change impacts on China's water resources[J].Advances in Water Science,2008,19(6):772-779.]
[3] 龙爱华,徐中民,王新华,等.人口、富裕及技术对2000年中国水足迹的影响[J].生态学报,2006,26(10):3358-3365.[LONG Aihua,XU Zhongmin,WANG Xinhua,et al.Impacts of population,affluence and technology on water footprint in China[J]. Acta Ecologica Sinica,2006,26(10):3358-3365.]
[4] 刘佳骏,董锁成,李泽红.中国水资源承载力综合评价研究[J].自然资源学报,2011,26(2):258-269.[LIU Jiajun,DONG Soucheng,LI Zehong.Comprehensive evaluation of China's water resources carrying capacity[J]. Journal of Natural Resources,2011,26(2):258-269.]
[5] 朱海勇,陈永金,刘加珍,等.塔里木河中下游地下水化学及其演变特征分析[J].干旱区地理,2013,36(1):8-18.[ZHU Haiyong,CHEN Yongjin,LIU Jiazhen,et al.Variation and evolution of groundwater chemistry in the middle and lower reaches of the Tarim River[J].Arid Land Geography,2013,36(1):8-18.]
[6] 潘峰,张清寰,何建华.甘肃董志塬地区第四系地下水补给环境与水化学特征演化[J].干旱区地理,2014,37(1):9-18.[PAN Feng,ZHANG Qinghuan,HE Jianhua.Groundwater recharge environment and geochemistry evolution of the Quaternary aquifer in the Dunzhiyuan region,Gansu Province[J].Arid Land Geography,2014,37(1):9-18.]
[7] 张人权,梁杏,靳孟贵,等.水文地质学基础[M].北京:地质出版社,2011:67.[ZHANG Renquan,LIANG Xing,JIN Menggui,et al.General Hydrogeology[M].Beijing:Geological Publishing House,2011:67.]
[8] LLOYD J W,HEATHCOTE J A.Natural inorganic hydrochemistry in relation to groundwater:An introduction[M].Oxford:Clarendon Press,1985:130-133.
[9] CHADHA D K.A proposed new diagram for geochemical classification of natural waters and interpretation of chemical data[J].Hydrogeol J,1999,7(5):431-439.
[10] 孙亚乔,钱会,张黎,等.基于矩形图的天然水化学分类和水化学规律研究[J].地球科学与环境学报,2007,29(1):75-79.[SUN Yaqiao,QIAN Hui,ZHANG Li,et al.Natural-water hydrochemistry classification and hydrochemistry rule research based on rectangle hydrochemistry diagram[J].Journal of Earth Sciences and Environment,2007,29(1):75-79.]
[11] PIPER A M.A graphic procedure in the geochemical interpretation of water analyses[J].American Geophysical Union Transactions,1944,25:914-923.
[12] 任孝宗.浑善达克沙地天然水体的水化学及水-沙相互作用研究[D].北京:中国科学院大学,2014.[REN Xiaozong.The ion chemistry of natural water and water-sand interaction in the Hunshandake Sandy Land[D].Beijing:The University of Chinese Academy of Sciences,2014.]
[13] 石培泽,马金珠,赵华.民勤盆地地下水地球化学演化模拟[J].干旱区地理,2004,27(3):305-309.[SHI Peize,MA Jinzhu,ZHAO Hua.Simulation of geochemical evolution of groundwater beneath the Minqin Basin,Gansu Province[J].Arid Land Geography,2004,27(3):305-309.]
[14] 王礼恒,董艳辉,宋凡,等.甘肃石油河流域地下水补给来源与演化特征分析.干旱区地理,2017,40(1):54-61.[WANG Liheng,DONG Yanhui,SONG Fan,et al.Recharge sources and hydrogeochemical properties of groundwater in the Shiyou River,Gansu Province[J].Arid Land Geography,2017,40(1):54-61.]
[15] 邵天杰,赵景波,董治宝.巴丹吉林沙漠湖泊及地下水化学特征[J].地理学报,2011,(5):662-672.[SHAO Tianjie,ZHAO Jingbo,DONG Zhibo.Water chemistry of the lakes and groundwater in the Badain Jaran Desert[J].Acta Geographica Sinica,2011,(5):662-672.]
[16] 朱秉启,杨小平.塔克拉玛干沙漠天然水体的化学特征及其成因[J].科学通报,2007,52(13):1561-1566.[ZHU Bingqi, YANG Xiaoping.The ion chemistry of natural water and its origins in the Taklamakan Desert,Xinjiang,China[J].Chinese Science Bulletin,2007,52(13):1561-1566.]
[17] 宋献方,李发东,于静洁,等.基于氢氧同位素与水化学的潮白河流域地下水水循环特征[J].地理研究,2007,(1):11-21.[SONG Xianfang,LI Fadong,YU Jingjie,et al.Characteristics of groundwater cycle using deuterium,oxygen-18 and hydrochemistry in Chaobai River Basin[J].Geographical Research,2007,(1):11-21.]
[18] 宋献方,刘相超,夏军,等.基于环境同位素技术的怀沙河流域地表水和地下水转化关系研究[J].中国科学(D辑:地球科学),2007,(1):102-110.[SONG Xianfang,LIU Xiangchao,XIA Jun,et al.Transformation relationship between surface water and groundwater using environmental isotopes in Huaisha River Basin[J].Science China(D:Earth Sciences),2007,(1):102-110.]
[19] AL-BASSAM A M,AWAD H S,AL-ALAWI J A.Durov Plot:A computer program for processing and plotting hydrochemical data[J].Ground Water,1997,35(2):362-367.
[20] AL-BASSAM A M,KHALIL A R.DurovPwin:A new version to plot the expanded Durov diagram for hydro-chemical data analysis[J].Computers & Geosciences,2012,42:1-6.
[21] 苏金明,傅荣华,胡卸文.基于MATLAB的地学可视化实现[J].地质灾害与环境保护,2009,20(1):83-86.[SU Jinming,FU Ronghua,HU Xiewen. The implementation of geo-visualization based on MATLAB[J].Journal of Geological Hazards and Environment Preservation, 2009, 20(1):83-86.]
[22] DUROV S A.Natural waters and graphic representation of their compositions[J].Dokl Akad Nauk SSSR,1948,59:87-90.
[23] ZAPOROZEC A.Graphical Interpretation of water-quality data[J].Ground Water,1972,10(2):32-43.
[24] BURDON D J,MAZLOUM S.Some chemical types of ground-water from Syria[M].Tehran:UNESCO,Paris,1958.
[25] LLOYD J W.The hydrochemistry of the aquifers of North-Eastern Jordan[J].Journal of Hydrology,1965,3:319-330.
[26] AL-BASSAM A M,KHALIL A R,KASSEM O M K.Automatic analysis reporting for hydro-chemical data using updated durovpwin program[J].Journal of Agricultural Science and Technology B,2013,3:757-768.
[27] PIPER A M.A graphic procedure in the geochemical interpretation of water-analyses[J].Transactions American Geophysical Union,1944,25(6):914-928.