地学统计学地一个操作教程.docx

《地学统计学地一个操作教程.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、地学统计学一、提取数据所用的工程数据来源是老师提供的坑探取样数据和钻探取样数据，数据是已给的Access数据库表，在基于VS2010平台上设计的，并从Access中提取数据。(1)选取的工程不同，则会调用对应的不同数据。调用数据时，先要通过采样数据来选择采样数据的类型（有坑探采样和钻探采样两种类型），然后通过工程名称按钮选择相应的工程，以坑探取样下的工程”-340m1XC”为例如下图1所示；图1(2)最后通过显示工程数据按钮显示数据。其调用数据结果界面如下图2所示：图2一、程序中主要应用的算法及其关键代码在本试验中主要使用了地学统计学法中的普通

2、科里克法和杨赤中滤波法对一批沿一维分布的地质勘探（坑探和钻探两种方式）取样数据进行统计分析。(1)普通科里克法，以所选工程总长度之半为最大滞后距分别计算两个元素品位（从铁、铜、硫中任选两个）的实验变异函数值，确定各自的变程a和基台值c，变异函数的公式如下：*h=12(N-h)x=1N-h'(Zx-Zx+h)2(1)其中为该位置的变异函数值，h为滞后距，h的个数h’为滞后距数，N为样品个数，Zx为在该位置某品位的测量值。假定h为1米，则最大值为样品总数的一半取整。实现普通科里克法的关键代码如下：ArrayTOList(GradeKey);vars

3、=newdouble[GradeKey.Length/2];for(inth=1;h<=GradeKey.Length/2;h++){doublesum=0;for(intj=0;j

4、三维)逐步消除随机性变化成分，显现出基本变化成分。应用杨赤中滤波推估法对高程观测值序列进行逐遍滤波，把随机性变化和规律性变化分离出来，然后求出反映复合变量(高程观测值)变化特征的函数和数字指标作为建立高程估值协方函数的依据，使滤波与推估有机地联系起来，由此建立高程估值数学模型。其公式如下：qik=qi-1k-1+2*qik-1+qi+1k-1/4(2)r*l=i=1N(qi-q)2/N-i=1N(qik-qik)2/N(3)该算法的关键代码如下：//杨赤中滤波方法publicvoidYCZCalculation(Listdou)

5、{Listtemp=newList();temp=dou;for(inth=1;h10){YCZCalculation(temp);}}publicvoidArrayTOList(double

6、[]temp){for(inti=0;idou){doublesum=0.0;for(inti=0;i

7、nCL;}一、开始测试及分析求值分别选用坑探采样类型下的”-340m1XC”工程和钻探采样类型下的”CK4A-4”工程两个工程来进行比较分析，以两个工程中的Fe、S两种元素进行试验，确定各自的变程a和基台值c，并计算选定元素的实验随机特征函数值和随机变化系数，以及对比、分析同一工程的两个元素品位及不同工程的同一元素品位的变化特征。(1)在同一工程中的不同元素的试验坑探取样表中的”-340m1XC”工程数据的Fe品位和S品位进行分析，采用地学统计法中的普通克里格算法拟合得到的曲线和采用杨赤中滤波算法拟合得到的Fey元素曲线图如下图3所示，图中左边

8、部分为普通克里格结果下的Fe元素变程a=36.456、基台值c=421.772，右边部分为杨赤中滤波法结果下的Fe元素随机特征函数值为322.388、