地质统计学资料.ppt
《地质统计学资料.ppt》由会员分享,可在线阅读,更多相关《地质统计学资料.ppt(77页珍藏版)》请在课桌文档上搜索。
1、地震数据约束的地质统计学,INTRODUCTION,近几年来,石油勘探开发中研究的一个热点是在储层建模中整合(integrate)不同类型的数据。最典型的就是利用地震数据信息来约束井间孔隙度的度量。目的就是用一种或几种密集采样的二级属性(如从三维地震数据中提取出的声阻抗、振幅或旅行时等)来约束相关的一级变量(孔隙度、渗透率、深度等),从而获得相应的空间分布。在这个过程中,地质统计学中不同的方法理论得到广泛的应用。,地质统计学理论方法简介,第一部分,本世纪50年代初期,南非矿业工程师克里金(kriging)在矿山工作时观察到金属的分布在空间上并非是纯随机的,而是在空间上具有相互联系。法国巴黎矿业
2、学院马特隆教授(G.Matheron)将克里金的经验和方法上升为理论,即区域化变量理论(Regio-nalized Variable Theory)的雏形,从而为地质统计学理论体系的形成创造了条件。美国斯坦福大学应用地球科学系儒尔奈耳(A.G.Journel)教授等人在1978年出版的专著矿业地质统计学(Mining Geostatistics)对地质统计学进行了系统的叙述,并总结了地质统计学在矿业中应用的实际经验。80年以来,在石油勘探开发领域广泛的应用于储层和油藏参数的空间估计建模及非均质性分析。,地质统计学的发展史,本世纪50年代初期,南非矿业工程师克里金(kriging)在矿山工作时观
3、察到金属的分布在空间上并非是纯随机的,而是在空间上具有相互联系。法国巴黎矿业学院马特隆教授(G.Matheron)将克里金的经验和方法上升为理论,即区域化变量理论(Regio-nalized Variable Theory)的雏形,从而为地质统计学理论体系的形成创造了条件。美国斯坦福大学应用地球科学系儒尔奈耳(A.G.Journel)教授等人在1978年出版的专著矿业地质统计学(Mining Geostatistics)对地质统计学进行了系统的叙述,并总结了地质统计学在矿业中应用的实际经验。,地质统计学的定义,以区域化变量理论为基础,以变差函数为基本工具,研究那些在空间分布上既具有随机性又具有
4、结构性的自然现象的科学。,显然,按定义,凡是要研究空间分布数据的结构性和随机性,并对其进行最优无偏估计,或要模拟所研究对象的离散性、波动性或其他性质时均可应用地质统计学的理论与方法。,地质统计学的定义,地质统计学的若干基本假设及理论(一)、区域化变量,所谓区域化变量是指以空间点的三个直角坐标为自变量的随机场。当对它进行了一次观测后,就得到了它的一个实现,它是一个普通的三元实值函数或空间点函数。区域化变量的两重性表现在:观测前把它看成是随机场(依赖于坐标),观测后把它看成一个空间点函数(即在具体的坐标上有一个具体的值)。,地质统计学的若干基本假设及理论,在地质统计学研究中是用变差函数表示研究范围
5、内区域化变量的空间结构性的,要用下式 计算变差函数时,必须要有,这一对区域化变量的若干现实,而在实际中只有一对这样的现实,即在、点只能测得一对数据(因为不可能恰在同一样点上取第二个样品),也就是说,区域化变量的取值是唯一的,不能重复的。为了克服这个困难,提出了如下的平稳假设及内蕴假设。,(1)、平稳假设(stationary assumption)满足下列两个条件时,称该区域化变量满足平稳假设:1、区域化变量的数学期望是一个常数:2、在整个研究区内,区域化变量的空间协方差函数存在且平稳:,协方差平稳意味着方差及变差函数平稳,从而有关系式:,地质统计学的若干基本假设及理论(二)基本假设,(二)、
6、内蕴假设(intrinsic assumption)(只考虑变量的增量,而不考虑变量本身)满足下列两个条件时,称该区域化变量满足内蕴假设:1、在整个研究区内,随机函数的增量的数学期望为0:2、对于所有矢量的增量的方差函数存在且平稳,即:,地质统计学的若干基本假设及理论(二)基本假设,变差函数(Variogram)的定义,距离为h的变量值差的数学期望的平方的一半。,注:有时把 也称为变差函数,变差函数能够量化的描述区域化变量的空间结构变化特征,在地质统计学及储层建模中有无比重要的作用。变差函数是一个距离函数。描述不同位置变量的相似性,值越大,相关性越差。通常,值随着距离矢量的增大而增大,直到达到
7、一定值,达到其极大值,而后保持这个常数值不变。,变差函数的基本参数,变程 用来度量空间相关性的最大距离。一般来说,随样品点间距离增大,变差值趋于增大,使变差函数达到一定的平稳值时的空间距离叫做变程。当空间距离较变程大时,变差函数仍保持其平稳值。基台值 变差函数在变程处达到的平稳值 此时,其变差值应为。然而,由于诸多因素的影响,比如抽样和实验误差以及小尺度的变异,上述结论不一定正确。例如在短距离内的大变异引起间隔非常近的样品有十分不相近的值,这就导致变差函数在原点的不连续性。块金值 在原点附近非零的变差函数值。块金效应 大变异性对原点附近变差函数的影响.它通常用块金值与基台值的比表示。相对块金效
8、应常用百分比的形式,下图给出了一个实验变差函数的例子,同时对上面引进的三个参数加以说明。,变差函数的定义,实验变差函数的计算,1 基本的变差函数公式变差函数的定义式是从所有可能的样品点对计算出来的。在实际中,假设是间距为的所有点对的总数,则变差函数可以通过下式计算:式中N(h)是步长为h数据的对的数目,和 是相距为h的两点采样值。因为变差函数是一种统计方法,在一定的步长上采样点越多,变差函数的估计结果越可靠。由于取样的特性,上面的等式不能直接在实际中运用。因为取样点对很少严格满足间距为h,除非是按照规则的网格取样,为了采用此公式必须考虑如下几点:,间距h为的取样点对的数目必须根据邻域的概念来计
9、算。换而言之,在矢量h的首尾用一维线段,二维面或三维体来计算点对,意思是如果一点位于以矢量h的首端为中心的邻域,而另一点位于以矢量h的尾端为中心的邻域,这对点视为一个距离为h的取样点对。变差仅对所选择的多个空间距离h来计算,通常的作法是选择一组规则的空间间距,如h,h,h等。这就暗示着连续的变差函数用一系列的离散值来逼近。变差函数的特征,比如形状和大小,通过对照所选的空间距离的离散变差点来检查。为了获得一个非常好的变差函数,取样点对不能太少,即取样的总数应尽可能地大。少量的样品点对将会产生怪异的实验变差值,因为变差函数是根据差的平方的平均值来定义的,少量反常的样品值便会很容易地歪曲变差图。,实
10、验变差函数的计算,采用邻域的概念,可以通过如下方法来修正,用以进行变差函数的实际计算 其中xi,xj代表两个取样位置,上面等式的右边被除,这纯粹为了数学的方便,另一个变化便是邻域概念的使用。取代一个确切的距离矢量,我们仅仅要求两个样品的距离近似等于先前定的距离矢量。这种改进使它可以对所有的有效距离都能计算出变差函数。变差函数是关于h对称的,这种关系暗示着对任何特殊的方向上计算的变差函数等于其反方向的计算结果。这种对称性使它很容易在克立格方法中应用。,实验变差函数的计算,定义样品点邻域的方法不止一种。一种典型的方法就是采用图所示的截断的楔形来定义样品的邻域,对距离和方向设定一个容许范围,这种方法
11、把一个点扩展到一个面或一个体。实际上,容许范围的应用是恰当的,因为地学测量中很难设法确定取样点的确切位置。在变差函数计算中允许范围的使用实际上是合理的。邻域的概念也使变差函数值有某种程度的光滑性,这有利于建立变差函数模型。,计算实验变差函数,.,例如:从第6个延迟距离开始计算,从一个点开始,考虑距离和角度的允许范围内的所有点。,.,移到下一个结点继续计算,计算实验变差函数,计算所有的结点,然后对多个空间距离重复上面的步骤,.,计算实验变差函数,理论变差函数模型的意义,实际中,实验变差十分混乱,使得掌握区域化变量的属性和利用变差函数进行结构分析变得困难,为了根据实验变差函数获取区域化现象的主要空
12、间结构,理论变差函数是必需的。在估计网格结点时,样品点和估计点之间的空间相关性必须给定,这无法通过实验变差来产生,因为样品点和估计点间的距离可能是任何数而不是距离间隔。(这类似于回归分析中的预测问题,拟合线被用来预测与任何自变量值(不一定是观测值)相应的因变量的值。),基本变差函数模型,1、块金效应模型块金效应模型、纯块金效应模型是最简单的变差函数,它们描述了在坐标原点不连续性的现象。许多实验变差,在原点或近原点具有明显的不连续性,这意味着有一明显的跳跃现象存在在距离原点的一个小的距离范围内,变差值从零跳到一个较大的值,这种现象可以由下面的一个函数来刻划:值得注意的是块金效应通常并不作为单独的
13、基本模型来考虑,而在变差函数中当作常数来考虑,然而这个常数给出了在原点的不连续性的程度,上面的记号可理解为纯块金效应的标准化形式,块金效应的大小是通过来给出的,而也表示在原点不连续性的程度,在变差的套合模型中这个符号的方便之处可以很容易看出来。,基本变差函数模型,2、球状模型球状模型是最普遍采用的变差模型,它的标准化形式为其中a是变程,一个主要特点是在原点附近的小范围内表现出线性行为,但在大距离时变得平缓,当为变程时达到基台值。模型的另一个特点是原点的切线在变程时便达基台值,这个事实在拟合实验变差函数时非常有用,在图中用黑实线给出的便是一球状模型。,基本变差函数模型,基本变差函数模型,3、指数
14、模型另一个普遍使用的跃迁模型是指数模型:其中a是常数,这模型渐近达到它的基台值,使变差值达到基台值的95的距离a被认为是近似的变程,相似于球状模型,指数模型在原点附近是线性的,它逐步增加,当值逐渐增加到一定程度时变平缓了,原点处的切线在变程值的/附近达到基台值,在拟合实验变差时记住这一点是非常有用的。图中的点线便是指数模型。,基本变差函数模型,基本变差函数模型,4、高斯模型高斯模型是用于刻划空间连续性的另外一个跃迁模型,其定义如下:a是常数,函数渐近地达到它的基台值,定义为变程,它使变差值达到基台值的,高斯模型的特点是在原点附近表现出抛物线性质,在图中用虚线表示。这是仅有的一个含有拐点的基本变
15、差模型。,基本变差函数模型,基本变差函数模型,5、线性模型某些实验变差在整个图形上表现出线性,线性模型用于描述这种线性连续性,它没有跃迁行为,模型定义如下 这个模型没有一个确定的变程。,基本变差函数模型,1、块金效应模型、纯块金效应模型,2、球状模型,3、指数模型,4、高斯模型,5、线性模型,基本变差函数模型,空间的各向异性,一个各向异性的变差函数,随着方向的改变,它的变程或基台值也发生明显地改变。在地质统计学的文献中,定义了两种主要的各向异性:几何的和带状的。几何各向异性描述在各个不同的方向上有不同的变程,但在所有的方向上有不变的基台值。而带状各向异性,基台值随方向变化而变程不变。,图几何各
16、向异性和带状各向异性的图示,实验变差函数建模过程,(1)需要判断的是被考察的空间连续性是各向同性还是各向异性。如果是各向同性,在变差模型中仅用全向变差就可以了;如果是各向异性,模型将变得非常复杂,需要用多个步骤来完成这个过程。在各向异性变差模型中的第一步就是辨别各向异性,这通过结合定性和定量信息能够做到,例如,主要的轴,可以通过地质特征,比如地质体的延伸方向、倾向等等来确定。(2)构造一个模型,能够描绘整个变差函数特征在距离和方向上的改变。(3)变差函数的标准化 这是焦点部分,就是,以使各向异性的问题可以采用各向同性变差函数同样的方法来对待。三维的各向异性变差函数需要结合不同方向的模型,这个组
17、合模型按基台值和变程来说在所有的方向上将是一致的。这个过程称为变差函数标准化。它通过一个转换,即把各个方向变差函数简化成一个具有统一的标准化变程的公用模型来实现。其关键是转化间距以使标准化模型在相同间距条件下在所有方向上有相同的变差值。,变差函数的套合(nest)结构,所谓套合结构就是把分别出现在不同距离h上和(或)不同方向上同时起作用的变异性组合起来。每一个变差函数代表了一种特定尺度上的变异性。,其中,变差图拟合方法地质法,()根据地质、地球物理数据,确定地质体走向。()沿走向计算实验变差值;()根据所有方向的变差模型,选择具有最大变程的一个可靠变差图,那么该方向的变程就是走向的变程,该方向
18、定义为主轴方向,旋转角就是这走向的方位角;()沿与主要走向垂直的方位计算另一组变差函数。模拟所有这些变差图并选择一个具有最大变程的可靠变差模型;其方向即为第二个轴方向,旋转角为倾角,跟这两个轴垂直的方向为第三个轴的方向。()用全部的数据集合计算沿第三个轴的变差模型并确定块金效应值,计算整体的变差图和基台值。,变差图拟合方法穷举法,()计算一组实验变差图,沿每一个固定角度的方位和每一个固定的倾斜角方向,根据每一个滞后距离增量计算变差图。()把所有的实验变差图转化到一个三维坐标系上。然后,这些数据在三维块段模型中进行插值,根据实验变差图,每一个子块段都被赋予一个变差值,鉴别主轴,确定主轴方向和角度
19、。()沿着上一步所确定的主轴方向计算特殊的变差模型。()使用全部数据计算沿第三主轴变差并确定块金值,计算整体变差图并确定基台值。(5)根据三个角度具有三个变程、块金值和基台值的三维变差模型最终被构造出来。,克里金方法简介,在确定性的储层参数建模中,往往应用插值方法对空间上每个网格赋以储层参数值(孔隙度、渗透率或含油饱和度)。插值方法很多,大致可分为传统的统计学估值方法和地质统计学克里金估值方法。由于传统的数理统计学插值方法(如距离平方反比加权法)只考虑观测点与待估点之间的距离,而不考虑已知点位置之间的相互联系,即地质规律所造成的储层参数在空间上的相关性,因此插值精度相对较低。为了提高对储层参数
20、的估值精度,人们广泛应用克里金方法来进行井间插值。,克里金方法简介,克里金方法,是根据待估点周围的若干已知信息,应用变差函数所特有的性质,确定待估点周围的已知数据点的参数对待估点的贡献(即加权值),然后对待估点的未知值作出最优(即估计方差最小)、无偏(即估计误差的数学期望为0)的估计,即最佳线性无偏估计(BLUE:best linear unbiased estimator)。,克里金方法简介,克里金技术区别于其他空间估计技术(如距离平方反比加权,三次样条等)的特点,在于对变量的空间相关性进行分析,计算了局部的最优估计,并能提供出估计误差。,克里金方法简介,简单克里金普通克里金泛克里金具有外部
21、漂移的克里金协克里金指示克里金,克里金方法简介具有外部漂移的克里金,趋势模型是:,估计值:,在应用外部漂移算法时,应该满足两个条件:(1)外部变量必须在空间光滑地变化,否则可能导致线性系统不稳定;(2)在主变量的所有数据点处和要估计的位置处,外部变量都必须是已知的。,克里金方法简介协克里金(CK),普通协克里金估计的初始变量和二级变量的线性组合 估计值:,传统的普通协克里金估计的方程组如下:,克里金方法简介配置协克里金(Collocated Cokriging),配置协克里金是协克里金的一种简化形式,有的文献中称为同位协克里金,即如果二级变量密集取样时,只保留与估计点同位的二级变量。配置协克里
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地质 统计学 资料
链接地址:https://www.desk33.com/p-248030.html