地质统计学与随机建模原理4-1实例:序贯高斯模拟
- 格式:ppt
- 大小:2.57 MB
- 文档页数:17
多点地质统计学原理、方法及应用概述及解释说明1. 引言1.1 概述本文旨在探讨多点地质统计学的原理、方法及应用,为读者提供一个全面了解该领域的概述。
多点地质统计学是一门研究如何有效地利用多变量数值以及空间数据进行地质分析和预测的学科。
它通过综合多种数据,包括物理测量数据、遥感图像数据和野外调查数据等,来实现对不同地质现象和过程的建模与研究。
1.2 文章结构本文按照以下结构组织内容:首先介绍多点地质统计学的基本原理,包括其定义与概念、基本假设以及原理解释。
随后,针对多点地质统计学的方法进行详细阐述,探讨数据收集与预处理、变量选择和缺失值处理以及统计模型拟合与优化算法应用等关键步骤。
接下来,我们将通过具体案例研究来展示多点地质统计学在矿产资源评估与勘探、地下水资源管理与保护以及石油勘探与开发中的应用实践。
最后,在结论部分对全文进行概括总结,并展望未来多点地质统计学研究的发展方向。
1.3 目的本文旨在全面介绍多点地质统计学的原理、方法及应用,以帮助读者对该领域有一个清晰的认识。
通过阐述基本原理和方法,读者可以了解多点地质统计学在地质分析和预测中的重要性。
此外,通过具体案例的引入,读者将能够更好地理解多点地质统计学在实际问题中的应用价值和潜力。
最后,通过对未来研究方向的展望,读者可以获得一些启示,并为自己在该领域开展研究提供参考。
2. 多点地质统计学原理2.1 定义与概念多点地质统计学是一种广泛应用于地质科学领域的统计学方法。
它通过对多个地点上的地质数据进行收集、分析和解释,旨在揭示地下资源的分布规律和空间变异性。
多点地质统计学基于一系列假设和方法,能够提供可靠的预测结果和决策依据。
2.2 基本假设在多点地质统计学中,存在几个基本假设:- 空间自相关假设:相邻位置上的地质现象存在关联性,即一个位置的观测值可能受到相邻位置观测值的影响。
- 空间平稳假设:在整个研究区域内,不同位置上的地质变量具有类似的变异性。
基于序贯高斯条件模拟的GEDI数据联合Landsat8反演森林地上生物量罗绍龙;舒清态;余金格;胥丽;杨正道【期刊名称】《林业科学研究》【年(卷),期】2024(37)3【摘要】[目的]单一遥感技术估测森林生物量存在较大局限性,本研究旨在利用多源遥感协同技术互补激光雷达和光学遥感的优势,提高生物量估测精度。
[方法]以星载激光雷达GEDI和光学遥感Landsat8数据为主要信息源,采用序贯高斯条件模拟方法实现GEDI光斑数据由“点”到“面”的空间扩展,结合地面138块生物量调查样地,利用随机森林回归方法估测云南省香格里拉云冷杉林的地上生物量。
[结果](1)采用序贯高斯条件模拟方法对GEDI光斑点进行空间插值,模拟的12个生物物理指标在空间分布上呈现出随机性、破碎化的特征,这与森林的空间分布聚集性非常相似,参与建模的9个指标OEC均大于0.90;(2)利用单一Landsat8光学遥感数据和地形因子构建的随机森林模型精度为:R2=0.82,RMSE=35.51 t·hm^(-2),P=0.77;Landsat8数据协同星载激光雷达GEDI数据构建的随机森林模型精度为:R2=0.86,RMSE=32.11 t·hm^(-2),P=0.80,模型精度明显提升;(3)利用多源遥感技术估测的香格里拉2019年云冷杉林地上的生物量总量为37 042 605.68 t,平均生物量为123.28 t·h m^(-2)。
[结论]基于地统计学的序贯高斯条件模拟方法考虑到研究对象的空间异质性、能克服一定的平滑效应,用于实现激光点由“点”到“面”的空间扩展是可行的。
星载激光雷达GEDI与光学遥感Landsat8协同的多源遥感数据可有效填补单一遥感数据源的缺陷,提高森林生物量的估测精度,能为激光雷达联合光学遥感估测大范围、全覆盖的森林生物量提供参考。
【总页数】12页(P49-60)【作者】罗绍龙;舒清态;余金格;胥丽;杨正道【作者单位】西南林业大学林学院【正文语种】中文【中图分类】X87;S757.2【相关文献】1.联合多点地质统计学与序贯高斯模拟的随机反演方法2.基于Landsat8影像的汝城县森林地上生物量遥感估算研究3.GEDI不同算法组数据反演森林最大冠层高度和生物量精度比较4.联合光学和合成孔径雷达数据的太平湖森林地上生物量反演研究5.深度学习方法下GEDI数据的天然云杉林地上生物量反演因版权原因,仅展示原文概要,查看原文内容请购买。
储层建模Reservoir Modeling 吴胜和金振奎黄沧钿陈崇河编著石油工业出版社内容提要本书系统地阐述了储层建模的基本理论、方法和技术。
从提高建模精度的角度出发,强调应用高分辨率层序地层学及地质概念模式约束建模过程,倡导等时建模、成因控制相建模以及相控建模的基本原则,并系统地介绍了确定性建模与随机建模的基本原理、各种方法技术及其地质适用性,从实用出发系统介绍了目前国际上储层建模的主要软件。
本书可供从事石油地质、开发地质及油藏工程等工作的科研人员参考,也可作为油气地质专业研究生的教学参考书。
前言随着油气田勘探开发的不断深入,储层研究以建立定量的三维储层地质模型为目标,这是储层研究向更高阶段发展的体现。
科学的油藏评价、油藏工程管理以及三维油藏模拟均要求三维储层地质模型,即储层属性的三维定量分布模型。
储层三维建模是国外80年代中后期开始发展起来的储层表征新领域。
十多年来,这一学科领域的发展十分迅速,数学地质及计算机工作者致力于发展各种建模方法特别是各种随机模拟方法,而储层地质工作者则研究各种建模方法的地质适用性,并力求改进方法以建立符合地质实际的三维储层地质模型。
本书作者近年来一直致力于储层建模的教学和科研工作,曾对多种油气储层类型,包括冲积扇、河流、扇三角洲、辫状河三角洲、正常三角洲、滨岸相、潮坪相等储层进行过三维建模研究,并与美国、挪威、法国、英国、澳大利亚等国家的相关专家进行广泛的交流。
为了使更多的油气地质工作者了解并掌握储层建模方法并切实应用于油气田勘探开发实际,有必要系统介绍储层建模的基本理论、方法和技术。
全书共分五章,第一章储层建模概论,介绍储层建模的目的意义、储层模型类型、储层建模流程及建模策略,论述了确定性建模与随机建模、等时建模、成因控制储层相建模以及相控建模的基本原则;第二章储层地质模式,介绍三维储层建模所必须了解的各种沉积环境的储层结构模式;第三章确定性建模技术,介绍三种主要的确定性建模方法,即储层地震地层学方法、储层沉积学方法及地质统计学克里格方法;第四章随机建模技术,介绍了随机模拟的基本原理、算法和随机模型,以及主要的随机模拟方法及其地质适用性,包括标点过程、序贯高斯模拟、截断高斯模拟、序贯指示模拟、分形模拟等;第五章储层建模软件,介绍目前国际上主要的储层建模软件。
多点地质统计学随机建模方法原理详细教程多点地质统计学(Multiple-Point Geostatistics,简称MPGS)是一种用于地质建模的统计学方法,旨在综合考虑多个地质属性之间的空间关系,可以用于模拟地质体结构和属性的空间分布。
下面是一个详细的MPGS建模方法的教程。
1.数据收集和准备首先,需要收集和准备地质数据。
这些数据可以包括钻孔数据、采矿数据、地球物理数据等。
数据应该包括多个不同属性的测量结果。
2.数据预处理对收集的数据进行预处理是为了消除异常值、填充缺失值和准备数据用于建模。
这些步骤可以包括数据清洗、插值等。
3.定义模型网格创建一个用于建模的三维网格,通常由正交的网格单元组成。
网格的尺寸和边界应根据实际问题的要求进行选择。
4.模式提取在做MPGS建模之前,需要从数据中提取出具有空间一致性和相关性的模式。
这可以通过模式提取算法实现,如基于模拟退火算法的直方图匹配。
5.模式匹配在模型建模过程中,需要通过模式匹配找到与已知数据最相似的地质模式。
这可以通过计算模式之间的相似性指标,如多点统计函数(MPS)实现。
6.模式合成一旦找到与已知数据相似的地质模式,可以根据模式之间的空间关系来生成新的地质模式。
这可以通过使用概率或变异性模型来实现。
7.模型重建利用已生成的地质模式,可以在模型网格单元上对地质属性进行插值,以重建地质体的结构和属性分布。
这可以使用插值方法,如克里金插值、逼近法等。
8.模型评估和修正完成模型重建后,需要评估模型的性能并根据需求对模型进行修正。
可以利用模型与实际数据之间的比较以及其他准则来评估模型的准确性和合理性。
9.模型应用完成最终的地质建模后,可以将模型应用于相关的地质问题,如矿产资源评估、地质风险评估等。
以上是MPGS建模方法的详细教程。
这种方法在地质建模中广泛应用,可以提供更准确和全面的地质属性分布信息,对于地质资源开发和管理具有重要意义。
序贯高斯模拟基本原理
嘿,朋友们!今天咱来唠唠序贯高斯模拟基本原理。
这东西啊,就像是搭积木一样有趣!比如说吧,想象一下你面前有一堆五颜六色的积木,你想怎么搭就怎么搭,最后搭出你心目中的城堡,序贯高斯模拟差不多也就是这个道理啦!
它首先呢,会去收集一堆数据,就像我们收集各种颜色的积木一样。
这些数据可能来自于地质探测啦、环境监测啦等等好多地方。
然后呢,它可神奇了,它会利用这些数据来构建一个模型,就如同我们用积木搭建出城堡的大致轮廓。
接着呀,它会一点一点地填充细节,就好像给城堡加上窗户、门这些小部件一样。
在这个过程中,它不是死板地进行的哦,而是充满了灵活性,可会“随机应变”啦!比如在一个区域数据比较密集的时候,它就能很准确地模拟;要是数据比较稀疏的地方呢,它也能发挥想象,给出一个合理的推测。
这多有意思啊!
比如说在石油勘探领域,序贯高斯模拟就能大显身手呀!勘探人员可以用它来模拟地下油藏的分布情况,说不定就能找到大油田呢,那可就发财啦,
哈哈!在环境保护中,也可以用它来模拟污染物的分布,这样就能更好地制定治理方案啦。
序贯高斯模拟就是这样一个神奇又实用的工具!它能让我们从一堆看似杂乱无章的数据中,找出规律,构建出一个有意义的模型。
是不是超级厉害?所以呀,千万可别小瞧了它,它真的能给我们的生活和工作带来很多惊喜和帮助呢!。
几种常用的反演方法综述一、稀疏脉冲反演(C onstrained Sparse Spike Inversion)1、原理:①首先假设地下地层的波阻抗模型所对应的反射系数序列模型是稀疏的,即由起主导作用的强反射系数序列和具有高斯背景的弱反射系数序列叠加而成。
②将地震记录与子波进行稀疏脉冲反褶积得到地层反射系数,一般是使用最大似然反褶积求得一具有稀疏特性的反射系数序列Ri。
根据①的假设可以导出最小目标函数:R(K)为第一个采样点的反射系数,M 为反射层数, N为噪音变量的平方根,L 为采样总数,ƛ根据目标函数,对每一道,从上到下推测反射系数的位置点,判断反射系数的幅值大小。
如此反复迭代修改每个反射系数的位置和幅度,使最后的修改误差最小符合似然比值的判别标准即可,这样就完成了一道的反褶积,得到该道的反射系数的分布。
③通过最大似然反演导出波阻抗Zi 反演公式为Zi=Zi-1*[(Ri+1)/Ri].具体的计算方法是稀疏脉冲序列每次建立的反射系数为一个脉冲,然后在地震资料中提取子波与初始反射系数进行褶积,得到一个初始合成地震记录,并用此合成地震记录与实际地震纪录作对比得到他们之间的残差,利用这个残差的大小来修改反射序列中脉冲的个数再次进行褶积运算,得到新的合成地震记录,再与实际地震资料对比,就这样循环迭代,直到残差达到最小,最后得到一个与实际地震资料达到最佳逼近的合成地震记录,获得宽频带的反射系数。
图1 稀疏脉冲反演每次建立反射序列为一个脉冲,增加脉冲进行循环迭代约束稀疏脉冲反演采用的是一个快速约束趋势的反演算法,约束条件主要是波阻抗趋势和地质控制,而波阻抗趋势又是由解释层位和断层来控制的,从而可以把地质模式融入进去得到一个宽带的结果,恢复地质信息中缺少的低频和高频成分。
约束稀疏脉冲反演的最小误差函数是:第二项为原始地震道与合成地震道的均方差的总和;第三项为趋势协调的补偿i 是地震道样点号;di是原始地震道;Si是合成地震记录;ri 为地震道采样点的反射系数;ti是波阻抗趋势;Zi是地震道采样点的波阻抗值,介于井约束的最大和最小波阻抗之间;ɑ是趋势最小匹配加权因子,一般情况下ɑ=1;p、q是L 模因子,一般情况下p =1,q=2是调节或平衡因子,与信噪比大小有关。
原文《SIMPAT: Stochastic Simulation with Patterns》作者:G. Burc ArpatSIMPAT算法的细节该算法很先进,与petrel中已有的SNESIM算法完全是两个档次,相信不久的将来它会被加入到petrel中。
SIMPAT redefines the geostatistical modeling as an image construction problem。
1 前言在地球科学领域,序贯模拟是应用最广泛的随机模拟技术之一。
序贯模拟最早用于传统的基于变差函数的地质统计学算法,如序贯高斯模拟,序贯指示模拟(Deutsch and Journel, 1998)。
90年代早期,多点地质统计学在扩展的序贯模拟中引入了训练图像的概念代替变差函数Guardiano and Srivastava, 1993)。
后来,Sterbell于2000年提出了具有里程碑意义的SENSIM 算法用于实际油藏的多点地质统计学模拟(Caers et al., 2003; Strebelle et al., 2002)。
文本,在扩展的序贯模拟框架下,探索出了一种不同的利用训练图像的方法。
通常大家对训练图像进行多点统计关系的扫描以获取条件概率,但是,本文将从训练图像寻找地质模式的过程视为图像处理。
找到了这些地质模式之后就可以进行序贯模拟,在模拟过程中,利用局部的模式相似性判别而不是局部的条件概率作为赋值依据。
该方法建立在多级网格系统之上,考虑地质模式的多尺度迭代运算,其中的尺度关系紧密耦合。
该算法在离散和连续变量中都能良好运行,同时可以吻合许多种局部的条件概率数据,如测井和地震数据。
这种算法叫SIMPAT。
2004年Casers 和Arpat的文章讨论了该方法的基本原理,本文意在研究SIMPAT算法的技术细节,包括SIMPAT算法涉及的基本概念、算法的细节,以及具体实例。
2 概念为了简介,该部分讨论SIMPAT 算法的相关概念时采用含有大量的例子,在第三部分讨论了多种类别的离散和连续变量的模拟。
基于地质统计先验信息的随机地震反演叶端南;印兴耀;孙瑞莹;王保丽【摘要】基于地质统计先验信息的随机地震反演方法是一种基于蒙特卡洛的非线性反演方法.在贝叶斯理论框架下,通过序贯高斯模拟方法(sequential Gaussian simulation,SGS)和逐渐变形算法(Gradual Deformation Method,GDM)得到基于地质统计学的先验信息,然后构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到反演问题的解.与确定性反演结果相比,该方法能够有效地融合测井资料中的高频信息,提高反演结果的分辨率.数值模拟试验表明:本方法的反演结果与理论模型吻合较好,具有较高的分辨率;序贯高斯模拟采用一种新的逐点模拟方式,并结合GDM,有效提高了随机反演的计算效率.【期刊名称】《物探化探计算技术》【年(卷),期】2015(037)003【总页数】7页(P341-347)【关键词】地质统计先验信息;贝叶斯理论;高分辨率;计算效率【作者】叶端南;印兴耀;孙瑞莹;王保丽【作者单位】中国石油大学(华东)地球科学与技术学院,青岛 266580;中国石油大学(华东)地球科学与技术学院,青岛 266580;中国石油大学(华东)地球科学与技术学院,青岛 266580;中国石油大学(华东)地球科学与技术学院,青岛 266580【正文语种】中文【中图分类】P631.4石油勘探和开发已经从简单的构造识别转向复杂构造、薄储层以及老油区剩余油的研究。
目前工业目标地质体的厚度一般在10m甚至10m以下,而实际的反演方法较难区分10m以下的薄层顶底[1],所以提高分辨率以反演薄层信息在储层识别和油藏描述中起到了重要的作用。
随着地球物理研究的深入,非线性反演方法更适用于复杂隐蔽性储层的研究。
地球物理反演中普遍存在非线性和不适定性的问题,若采用线性方法来反演非线性反演问题通常较难得到其真值[2],而随机反演能够在空间相关性和井的约束下,模拟出地震频带以外的信息,分辨率高于常规的确定性反演方法。
多点地质统计学随机建模方法原理详细教程多点地质统计学随机建模是一种应用于地质领域的统计学建模方法,它主要用于处理地质参数在空间上的变化规律。
该方法的原理基于地质参数的随机性和空间相关性,通过构建具有地质属性的随机模型,可以模拟地质现象的空间分布。
具体而言,多点地质统计学随机建模方法主要包括以下几个步骤:1.数据准备:收集与地质参数相关的数据,例如岩性、厚度、含矿物质等。
要求数据具有一定的地质意义和空间分布规律。
2.变量描述:对收集到的数据进行统计分析,包括计算均值、方差、协方差等统计指标,以描述地质参数的分布特征。
3.变量变换:根据地质参数的实际特征,对数据进行变换,使其符合正态分布、对数正态分布或其他分布类型。
4.空间相关性建模:通过计算地质参数之间的空间相关性,可利用协方差函数、变差函数或半方差函数等,建立地质参数之间的空间相关模型。
5.随机模拟:根据变量的统计特征和空间相关模型,结合随机数生成算法,生成符合实际情况的具有随机性和空间相关性的地质参数数据。
6.模型验证:对生成的地质模型进行验证,比较随机模拟结果与实际数据的吻合程度。
可以使用统计指标如均值、方差、协方差等进行对比分析。
7.地质模型应用:根据随机模拟结果,可以进一步进行岩层插值、矿产资源评估和地质灾害风险评估等相关研究及应用。
总的来说,多点地质统计学随机建模方法是将统计学原理应用于地质参数的空间分布建模,通过对地质参数的统计特征和空间相关性的建模,生成具有随机性和空间相关性的地质模型。
这种方法可以提供地质领域研究的基础数据和分析手段,为地质灾害风险评估、资源勘探和环境评价等问题提供科学依据。
序贯高斯模拟在矿石品位估计中的应用研究刘占宁;宋宇辰;孟海东;于晓燕【摘要】随机模拟是地质统计方法的重要内容.在矿石品位估计方法中克里格方法作为一种无偏估计方法,常被用于矿石品位的估计.但克里格法估值存在平滑效应.作者在分析序贯高斯模拟和普通克里格法基本原理的基础上,运用序贯高斯模拟方法和普通克里格方法对某铁矿体内全铁(TFe)品位进行估计,给出了品位估计结果模型.研究从勘探线方向、垂直勘探线方向和竖直方向分别计算变差函数,对球状模型、指数模型、高斯模型的变差函数拟合效果进行了优选,结果表明球型模型拟合效果最好.针对序贯模拟和克里格品位估值效果进行了分析,结果显示:序贯高斯模拟结果在品位分布形态上更接近样品品位分布形态,其平滑效应更小;克里格方法估计与序贯高斯模拟方法相比仅在品位均值方面更接近样本品位均值.因此,认为序贯高斯模拟方法可以更好地刻画矿体内品位分布状态.%Stochastic simulation is an important method of geological statistics.The Kriging method is a unbiased estimation method,ften used for ore grade estimation with smoothing effect.Based on analysis of principle of Sequential Gaussian Simulations and Ordinary Kriging method they are used to estimate the TFe grade of a iron ore body and the grade model is given.The variation function is calculated from the direction of the exploration line,the vertical exploration line direction and the vertical direction.The variation function is fitted with the spherical model,the exponential model and the Gaussian model.The fitting effect of he spherical model is the best.The result show that distribution of the Sequential Gaussian Simulation is closer to the distribution of sample grade,with smaller smoothing pared withSequential Gaussian method only mean grade value of the Kriging method is closer to to that of the sample.Therefore,the authors believe that the Sequential Gaussian Simulation can better characterize the distribution of ore in the ore body.【期刊名称】《地质找矿论丛》【年(卷),期】2018(033)001【总页数】7页(P149-155)【关键词】序贯模拟;变差函数;矿石品位;克里格法【作者】刘占宁;宋宇辰;孟海东;于晓燕【作者单位】内蒙古科技大学矿业研究院,内蒙古包头014010;内蒙古科技大学矿业研究院,内蒙古包头014010;内蒙古科技大学矿业研究院,内蒙古包头014010;内蒙古科技大学矿业研究院,内蒙古包头014010【正文语种】中文【中图分类】P628.20 引言条件模拟技术是地质统计学中继克里格估计技术之后,迅速发展的一个新方法[1]。
地质统计学反演方法在煤厚及其夹矸预测上的应用作者:薛海军高阳来源:《地球》2013年第12期[摘要]在岩性解释方面,目前比较广泛使用的是稀疏脉冲反演,但其纵向分辨能力有限,在煤田勘探中,煤层厚度仅有几米的厚度,夹矸也只有几十公分的厚度,使用稀疏脉冲反演方法无法满足煤厚及其夹矸预测和描述的需要。
本文以有着高分辨率的地质统计学反演方法为基础,来论述地质统计学反演方法在煤厚及其夹矸预测上的应用。
[关键词]稀疏脉冲反演地质统计学随机反演煤厚预测岩性反演煤层夹矸[中图分类号] P628+.2 [文献码] B [文章编号] 1000-405X(2013)-12-77-20引言在纵向分辨率上,如何提高煤层的高精度预测和描述能力,是当前煤田勘探一个重要的研究问题。
三维地震资料在横向上分辨能力较高,在垂向上分辨能力较低。
为此在1992年Bortoliand就提出了地质统计反演方法,近年来随着计算机技术的快速发展,该方法得到越来越广泛的应用。
地质统计反演是以地质信息(包括地震、钻井、测井等)为基础,应用随机函数理论和地质统计学方法(变差函数分析、直方图分析、相关分析等),结合传统的地震反演技术,在每一个地震道(或多个地震道)产生多个可选的等概率反演结果的一种地震反演方法。
地质统计反演的应用主要有3个方面:随机地震反演+序贯高斯模拟系列;随机地震反演+序贯指示模拟系列;递推反演+随机模拟系列。
1方法原理简述在反演过程中,首先统计随机变量的直方图分布,统计随机变量的变差函数,确定随机变量在不同方向上的变程;然后,针对每一种实现,在地震道位置根据计算出来的反射系数人为合成地震道,比较合成道与原始地震道之间的误差,达到要求的精度后输出反演结果,主要步骤包括地质统计学分析、随机模拟和随机地震反演等步骤,具体流程如图1所示。
2实际资料应用以陕西某煤田工区为例,该工区的区域构造为平缓的南倾单斜,目标煤层厚度为1m-8m。
由于本研究区地层发育较为平稳,使用1.5km×1.5km的变程参数。
如何让地质模型更准确关于如何做好建模,每个人都有自己的看法,有人从理论上一套一套的,有人从操作上一套一套的。
那么什么样的模型是一个可靠的地质模型呢,众所周知有一个不是答案的答案:逼近地下真实的模型。
那么问题就来了,怎样来逼近地下真实呢?地下的东西看不见摸不着,谁知道到底如何。
没有调查就没有发言权,看看大伙都是怎么致力于提高模型准确度的吧。
上cnki,上spe,上aapg,国内国外期刊,sci,ei,核心,哪怕非核心都看一看。
结论就出来了,关于如何提高模型的可靠程度,几十年来大概集中在以下四个方面:1、建模算法2、地质认识3、约束建模4、模型优选一、建模算法算法有很多,成熟的(商业化的软件中),不成熟的(没有商业化的,高校与油田发表的文章中多见)。
不管是哪种算法,从数学角度来看,既然是随机模拟,每个模拟的概率是相同的,因此每个模型都是合理的。
但从油田生产角度来看,这些模型无一不经过一而再、再而三的调整,因此这些模型又都是不准确的。
因此,这些模拟实现仅在某种程度上符合储集层的数学特征,逼近其“数学真实”而非“地质真实”。
关于算法也是很多技术人员很了解的内容,这里不谈令人眼花缭乱的公式,只谈应用,认清本质,比会推导公式更重要。
先讨论下大家最为熟悉,用得也最多的序贯指示到底有什么能耐。
首先看个实例效果吧,以基于变差函数的序贯指示模拟为例。
图a是研究区的岩相分布图(蓝色为砂,白色为泥),图b,c是序贯指示两次随机模拟实现。
统计、对比发现,模拟结果砂岩发育长轴与井点统计参数基本一致(图e),而砂体形态及部分细节刻画上与原始相图差别较大,甚至相图上部分泥岩分布区域模拟出连片砂岩。
可见,地质统计学中,变差函数仅仅描述了数据的空间变异程度,而不同形态的地质体可能拥有相同的变异程度,图a,b,c就是这种情况,即井点统计变差函数模拟出来的地质体形态在服从输入数据统计的前提下,可以有无限多个,不确定性极大,但却鲜有与实际吻合,通常需反复多次调整。