[交流] 关于学习地质统计学的一些建议 已有11人参与

最近收到很多虫友关于地质统计学方面的咨询，也都逐个回复了，但是还是觉得非常有必要单开一个帖子说点什么。其实我学习地质统计学的时间也不是很长，到目前为止才两年不到，所以说不上懂行，如果下面有说的不对的地方还请各位批评指正。 相信很多虫友接触到地质统计学都是因为要对研究对象进行空间变异分析，或者克里格插值，其中最为关键的应该就是变异函数（又称半方差函数、半变异函数）拟合了。相信很多朋友也发现了，介绍变异函数或者地质统计学的著作很多，但是几乎没有对应的商业型拟合软件。这一点主要是因为变异函数拟合存在很大的不确定性：由于原始的变异函数模型及参数未知，实际过程中我们只能通过有限的观测数据去推断，因此结果不可避免的会存在不确定性。变异函数拟合的方法有很多种，具体的参见文献《Comparison of inference methods for estimating semivariogram model parameters and their uncertainty: The case of small data sets》以及其中涉及到的一些文献。在回复一些虫友的咨询中我也提到，现在国外比较认可的拟合方法主要有加权最小二乘法、极大似然估计等，对应的程序有VARFIT、MLREML，其中VARFIT已经有虫友上传到小木虫了。这两个程序都是同一个作者编写的，关于程序的介绍都有相应的文献，对这方面感兴趣的朋友可以多阅读一下这个作者的其他文章，大多是关于变异函数拟合、不确定性分析、克里格插值方面的，相信对大家会很有帮助。另外，希望大家在进行地质统计分析之前多看看专业相关的著作，中文的推荐《地统计学概论》、《空间变异理论及应用》，英文的有《GSLIB》、《Geostatistics for environmental scientists》、《Introduction to Geostatistics: Applications in Hydrogeology》、《Geostatistics: modeling spatial uncertainty》，英文的可以重点看第一和最后一本。其中GSLIB是一套完整的体系，包括几乎地质统计的各个模块，里面的很多小程序都非常实用，学习起来也比较简单。最后，跟大家提一下GS+和ArcGIS。GS+由于傻瓜式的操作，使得用户不需要太多地了解地质统计学的相关知识就可以进行空间变异分析，尤其是融合了变异函数拟合功能，使得它广受一些学者的喜欢。但是，GS+在国外的认可度是比较低的，因为GS+在拟合变异函数时采用的是最小二乘法，其拟合出来的结果精度往往会很不理想，甚至会得出明显错误的结论，而如果没有地质统计学知识背景的话，很有可能就会将错误的结果运用到自己的论文当中（有些已经刊出的文章就发现过这个问题）。ArcGIS相比于GS+，结果的可信度就更差了（至于为什么会出现这个问题，需要查看ArcGIS中地统计学模块的相关代码，奈何自己编程水平太菜，所以就留给其他人吧），但ArcGIS胜在出图美观，因此，可以外部（比如用VARFIT、MLREML）拟合好变异函数后，将相应的模型及参数代入到ArcGIS中。

回复此楼