“地质统计学”相关问题请教

geoscience

郝老师您好，看论坛里您的帖子“闲扯地质统计学“，感觉里面信息量很大，看的有点晕了，呵呵，有几个问题想请教下
  （1）克里金插值的原则，最优无偏线性（blue），关于”最优“和”无偏“的解释，数学表达式看的有点晕，我看有些资料里面写的是”估计值与实际值的差的方差最小（最优）“，”估计值与实际值的差的数学期望为零（无偏）“，想问下，这个估计值，实际值到底是什么意思呢，是有井点的网格处的值就是实际值，没有井点的网格处算出来的值就是估计值？
（2）想问下克里金插值在petrel软件中实现的时候是按什么顺序计算的，我理解的是：输入已知井点数据—将已知数据转化成正态分布—得到不同方向下的变差函数——利用变差函数计算未知点的值（满足blue原则），这样就完成了克里金插值的过程，不知道理解对不对
  刚接触地质建模不久，很多不懂的问题，请老师不吝赐教，谢谢！

Joseph

简单点儿：
1、是的。
2、是的。克里金插值本身是通过方向性选择来补偿权重，所以不需要考虑顺序，它只按照各已知点跟未知点的距离和方向来分配权重。SGS（序贯高斯模拟）之类的过程因为后来的未知点计算借用了部分前面步骤的克里金插值结果才产生计算顺序的问题（即遍历网格路径问题）。

Geokenneth

学习了谢谢

geoscience

Joseph 发表于 2014-5-25 01:39
简单点儿：
1、是的。
2、是的。克里金插值本身是通过方向性选择来补偿权重，所以不需要考虑顺序，它只按 ...

谢谢郝老师，想继续请教几个问题
（1）您在”闲扯地质统计学“中讲到，“这样，我们就明白了，Lag Distance其实就是Lag这个过滤器的远近，换个说法大概就是Lag的个数乘以Lag Tolerance，当然为了这个过滤器的需要除了Lag Tolerance还设定了其它的尺寸参数”。lag的意思我明白就是“铅笔头儿”上每一个蓝色格子，但是lag distance和lag tolerance还是不理解，我看petrel里面lag tolerance的默认值是50%，这说明lag tolerance是基于某个值的一个变化范围吗？我把您在”闲扯地质统计学“中的与之相关的一张图和您在石油大学（华东）讲座的资料中的一张图对比了一下（如图），好像这两个参数在图上的表示不一样，郝老师能否再解释一下。
（2）”形象的说，Lag就是一个筛子一样的东西，要把不同距离的变差结果通过这样的计算简单放进不同的位置。如果我们前面总结的所谓“距离越远你猜的越不准确”的结论成立，则应该随着Lag个数的增加而看到变差变大的效应。（你知道这个过程中你在继续损失什么吗？如果你知道，你就会明白为什么我说想做精确模型的人是外行了。一笑。）“，不太理解这句话的意思。
（3）在petrel里面做数据分析的时候，是否一定需要把regression line和edit line这两根线吻合起来呢，调的时候有什么依据吗？
     谢谢老师

Joseph

geoscience 发表于 2014-5-27 17:16
谢谢郝老师，想继续请教几个问题
（1）您在”闲扯地质统计学“中讲到，“这样，我们就明白了，Lag Dist ...

1、这两张图完全一样的表达，都是用来说明术语含义的。如果一定要说不同，唯一的不同就是黑色的那张没有说明Lag tolerance的意思，它简单默认50%了。

2、如果你可以把井根据你给出的两张图来进行采样，并认为井上数据是代表着过滤网格（Lag Tolerance - 2 * Bandwidth）中心点的贡献，并以此为标准进行数据规律统计，那么你自然无法要求数据的规律性具有精确度超过这种过滤网格的精度，对吗？ 网友773377曾经批判地质统计学的这个问题，即，纯粹是万般无奈的时候的选择。虽然我不全然赞同他的观点，但是也承认如果你选择了地质统计学就要明白数据的规律性是基于这样的假设而来，那么在基于这样的规律而反推的模型是不可能具有精确性的。它甚至不表达大部分井点位置的实际规律，它表达着某种被过滤网格规范化以后的数据采样的规律性。

3、我们不强制回归曲线和适用曲线的一致，因为回归曲线毕竟是一个数学表达，尤其在数据样本不足的情况下单纯依靠数据样本本身形成的规律可能会不符合我们的整体区域地质认识。至于说参数调整的规则则恐怕会“具体问题具体分析”。

论坛同行很多，我也只是抛砖引玉希望与大家就这个平台展开讨论，所以发问并不需要特定给我。——虽然我一般看到自己知道的问题会主动作答。