多点地质统计学
- 格式:docx
- 大小:581.50 KB
- 文档页数:3
多点地质统计学python(原创版)目录1.多点地质统计学概述2.Python 在多点地质统计学中的应用3.多点地质统计学的具体方法4.Python 实现多点地质统计学的优势5.总结正文【1.多点地质统计学概述】多点地质统计学是一种研究地下岩石、矿产等地质信息的学科,它通过对大量地质数据进行统计分析,从而揭示地下结构的规律和特征。
在地质勘探、矿产资源开发等领域具有重要的应用价值。
【2.Python 在多点地质统计学中的应用】Python 作为一种广泛应用于数据处理和分析的编程语言,在地质统计学领域也有着广泛的应用。
Python 提供了丰富的库和工具,如 NumPy、SciPy、Pandas 等,可以方便地进行地质数据的处理、分析和可视化。
【3.多点地质统计学的具体方法】多点地质统计学主要包括以下几种方法:(1)克里金插值:是一种基于空间相关性的插值方法,可以预测地下岩石、矿产等地质信息的空间分布。
(2)序贯模拟:是一种基于随机过程的地质模拟方法,可以模拟地下结构的演化过程。
(3)多元统计分析:包括主成分分析、聚类分析等方法,可以对地质数据进行降维、分类等处理。
【4.Python 实现多点地质统计学的优势】(1)强大的数据处理能力:Python 可以方便地处理大规模、复杂的地质数据,满足多点地质统计学的需求。
(2)丰富的库和工具:Python 提供了丰富的库和工具,如 NumPy、SciPy、Pandas 等,可以方便地进行地质数据的分析和可视化。
(3)灵活的编程环境:Python 作为一种高级编程语言,具有语法简洁、易学易用等优点,可以提高地质统计学的工作效率。
【5.总结】多点地质统计学是一门研究地下岩石、矿产等地质信息的学科,Python 作为一门广泛应用于数据处理和分析的编程语言,在地质统计学领域有着广泛的应用。
多点地质统计学原理、方法及应用概述及解释说明1. 引言1.1 概述本文旨在探讨多点地质统计学的原理、方法及应用,为读者提供一个全面了解该领域的概述。
多点地质统计学是一门研究如何有效地利用多变量数值以及空间数据进行地质分析和预测的学科。
它通过综合多种数据,包括物理测量数据、遥感图像数据和野外调查数据等,来实现对不同地质现象和过程的建模与研究。
1.2 文章结构本文按照以下结构组织内容:首先介绍多点地质统计学的基本原理,包括其定义与概念、基本假设以及原理解释。
随后,针对多点地质统计学的方法进行详细阐述,探讨数据收集与预处理、变量选择和缺失值处理以及统计模型拟合与优化算法应用等关键步骤。
接下来,我们将通过具体案例研究来展示多点地质统计学在矿产资源评估与勘探、地下水资源管理与保护以及石油勘探与开发中的应用实践。
最后,在结论部分对全文进行概括总结,并展望未来多点地质统计学研究的发展方向。
1.3 目的本文旨在全面介绍多点地质统计学的原理、方法及应用,以帮助读者对该领域有一个清晰的认识。
通过阐述基本原理和方法,读者可以了解多点地质统计学在地质分析和预测中的重要性。
此外,通过具体案例的引入,读者将能够更好地理解多点地质统计学在实际问题中的应用价值和潜力。
最后,通过对未来研究方向的展望,读者可以获得一些启示,并为自己在该领域开展研究提供参考。
2. 多点地质统计学原理2.1 定义与概念多点地质统计学是一种广泛应用于地质科学领域的统计学方法。
它通过对多个地点上的地质数据进行收集、分析和解释,旨在揭示地下资源的分布规律和空间变异性。
多点地质统计学基于一系列假设和方法,能够提供可靠的预测结果和决策依据。
2.2 基本假设在多点地质统计学中,存在几个基本假设:- 空间自相关假设:相邻位置上的地质现象存在关联性,即一个位置的观测值可能受到相邻位置观测值的影响。
- 空间平稳假设:在整个研究区域内,不同位置上的地质变量具有类似的变异性。
储层多点地质统计学随机建模方法摘要:多点地质统计学使用训练图像代替变差函数,将更多的地质资料整合到储层建模过程中,使得最终模型更加符合地质认识。
随着研究的不断深入,越来越多的地质工作人员开始熟悉这一方法,凭借自身的独特优势,多点地质统计学将在储层建模领域占得重要的一席。
关键词:多点地质统计学训练图像储层建模一、多点地质统计学与训练图像基于变差函数的传统地质统计学随机模拟是目前储层非均质性模拟的常用方法。
然而,变差函数只能建立空间两点之间的相关性,难于描述具有复杂空间结构和几何形态的地质体的连续性和变异性。
针对这一问题,多点地质统计学方法应运而生。
该方法着重表达空间中多点之间的相关性,能够有效克服传统地质统计学在描述空间形态较复杂的地质体方面的不足。
多点地质统计学的基本工具是训练图像,其地位相当于传统地质统计学中的变差函数。
对于沉积相建模而言,训练图像相当于定量的相模式,实质上就是一个包含有相接触关系的数字化先验地质模型,其中包含的相接触关系是建模者认为一定存在于实际储层中的。
二、地质概念模型转换成图像训练地质工作人员擅于根据自己的先验认识、专业知识或现有的类比数据库来建立储层的概念模型。
当地质工作人员认为某些特定的概念模型可以反映实际储层的沉积微相接触关系时,这些概念模型就可以转换或直接作为训练图像来使用。
利用训练图像整合先验地质认识,并在储层建模过程中引导井间相的预测,是多点地质统计学模拟的一个突破性贡献。
可以将训练图像看作是一个显示空间中相分布模式的定量且直观的先验模型。
地质解释成果图、遥感数据或手绘草图都可以作为训练图像或建立训练图像的要素来使用。
理想状态下,应当建立一个训练图像库,这样一来建模人员就可以直接选取和使用那些包含目标储层典型沉积模式的训练图像,而不需要每次都重新制作训练图像。
三、多点模拟原理进行多点模拟,需要使用地质统计学中的序贯模拟。
但是,多点模拟与传统的基于变差函数的两点模拟是不同的。
多点地质统计学随机建模方法原理详细教程多点地质统计学(Multiple-Point Geostatistics,简称MPGS)是一种用于地质建模的统计学方法,旨在综合考虑多个地质属性之间的空间关系,可以用于模拟地质体结构和属性的空间分布。
下面是一个详细的MPGS建模方法的教程。
1.数据收集和准备首先,需要收集和准备地质数据。
这些数据可以包括钻孔数据、采矿数据、地球物理数据等。
数据应该包括多个不同属性的测量结果。
2.数据预处理对收集的数据进行预处理是为了消除异常值、填充缺失值和准备数据用于建模。
这些步骤可以包括数据清洗、插值等。
3.定义模型网格创建一个用于建模的三维网格,通常由正交的网格单元组成。
网格的尺寸和边界应根据实际问题的要求进行选择。
4.模式提取在做MPGS建模之前,需要从数据中提取出具有空间一致性和相关性的模式。
这可以通过模式提取算法实现,如基于模拟退火算法的直方图匹配。
5.模式匹配在模型建模过程中,需要通过模式匹配找到与已知数据最相似的地质模式。
这可以通过计算模式之间的相似性指标,如多点统计函数(MPS)实现。
6.模式合成一旦找到与已知数据相似的地质模式,可以根据模式之间的空间关系来生成新的地质模式。
这可以通过使用概率或变异性模型来实现。
7.模型重建利用已生成的地质模式,可以在模型网格单元上对地质属性进行插值,以重建地质体的结构和属性分布。
这可以使用插值方法,如克里金插值、逼近法等。
8.模型评估和修正完成模型重建后,需要评估模型的性能并根据需求对模型进行修正。
可以利用模型与实际数据之间的比较以及其他准则来评估模型的准确性和合理性。
9.模型应用完成最终的地质建模后,可以将模型应用于相关的地质问题,如矿产资源评估、地质风险评估等。
以上是MPGS建模方法的详细教程。
这种方法在地质建模中广泛应用,可以提供更准确和全面的地质属性分布信息,对于地质资源开发和管理具有重要意义。
多点地质统计学python -回复"多点地质统计学python"是一种基于Python编程语言的地质统计学方法,它通过使用多点统计技术来分析地质数据和建立地质模型。
多点地质统计学python可以用于地质建模、矿产资源评估、油气勘探等地质领域。
本文将逐步介绍多点地质统计学python的基本原理、数据处理步骤以及应用案例。
第一步:理解多点地质统计学的基本原理多点地质统计学是一种从多个样本点中提取信息进行推断的统计学方法,它考虑了地质数据在空间上的相关性和变异性。
这种方法的基本原理是通过对多个样本点进行统计分析,揭示地质变量的分布和空间结构。
多点地质统计学方法在地质学领域得到广泛应用,可以实现对地质数据的建模和预测。
第二步:了解多点统计技术在python中的实现Python是一种功能强大且易于使用的编程语言,具有丰富的数据处理和分析库,例如NumPy、Pandas和Matplotlib等。
这些库提供了处理大型地质数据集和执行统计分析的工具。
在python中,可以使用这些库来实现多点地质统计学方法,例如多点统计算法、样本点选择和空间变异性分析等。
第三步:处理地质数据集在使用多点地质统计学python之前,首先需要准备一个地质数据集。
这个数据集包含了地质变量的观测值,例如地层厚度、属性值或矿石品位等。
地质数据集通常以表格形式存储,可以使用Pandas库将数据导入到python环境中进行处理。
第四步:分析地质数据集在有了地质数据集之后,可以使用多点地质统计学python来分析数据。
首先,可以使用NumPy库计算地质数据的统计指标,例如均值、方差和协方差等。
这些统计指标可以用来揭示地质变量的中心趋势和空间关系。
接下来,可以使用多点统计算法来分析地质数据的变异性和空间结构。
其中,常见的多点地质统计学方法包括变差函数、半变异函数和变异权重等。
这些方法可以用来评估地质数据的变异性,并揭示地质变量之间的空间关联性。
多点地质统计学在点坝内部构型三维建模中的应用刘可可;侯加根;刘钰铭;史燕青;柳琳;唐力;高兴军;周新茂【摘要】首次将多点地质统计学应用于点坝内部三维建模。
基于60 m超小井距资料,统计点坝内部夹层发育特征,利用统计结果人机交互绘制训练图像,定量表征了夹层的厚度、倾角、频率、密度、水平间距等信息。
选取典型点坝,将单井解释夹层沿着夹层面垂直投影至点坝顶面,结合点坝沉积样式,获取建模过程中旋转数据体,以此来表征夹层走向信息。
利用多点地质统计学Snesim算法进行点坝内部侧积夹层三维建模,与序贯指示模拟方法在同一点坝和旋转数据体基础上建立的点坝内部三维构型模型进行对比分析。
研究结果表明,序贯指示模拟方法建立的点坝内部三维构型模型,虽然能够在一定程度上表征点坝内部夹层特征,但在表征夹层的连续性方面效果不佳,并且由于缺少训练图像的约束,只能定性展示夹层发育情况,无法达到定量刻画的程度。
而利用多点地质统计学方法建立的模型,点坝内部夹层受控于训练图像及旋转数据体的双重约束,能够定量再现夹层的发育规模和产状,精确表征点坝内部夹层的几何形态与空间结构。
结合水淹信息表明,此多点方法建立的点坝内部三维构型模型对于剩余油分布研究有指导意义。
%For the first time,an application of Multiple-point geostatistics in 3D internal architecture modeling of point bar is proposed in this paper.Data from dense well patterns ( well spacing of 60 m) offer us the statistical distribution charac-teristics of interbeds within the point bar.Training images are produced by interactive processing of the statistical results for quantitative characterization of the interbeds,such as theirthickness,dip,frequency,density and horizontal spacing.A typical point bar is first choosed for vertically projecting the interbeds interpreted fromsingle well data to its top surface along the surfaces of interbeds,and then a 3D spinning cube is calculated in combination with the depositional patterns of the point bar for strike characterization of the interbeds.The Snesim algorithm is used for 3D architectural modeling of lat-eral accretion,and the results are compared with those from sequential indicator simulation ( SIS) .It is concluded that SIS may represent the characteristics of interbeds,but it is inadequate to characterize their continuity.Besides,a lack of constraints from training images forbids satisfactory quantitative characterization.In contrast,the model based on Snesim is under the dual control of training images and 3D spinning cube,thus it can quantitatively characterize the distribution,ge-ometry and spatial structure of the interbeds.The data of watered-out reservoirs reveal that the architecture model built with multiple-point geostatistics method can guide the mapping of remaining oil.【期刊名称】《石油与天然气地质》【年(卷),期】2016(037)004【总页数】7页(P577-583)【关键词】训练图像;旋转数据体;密井网;多点地质统计学;点坝;剩余油分布【作者】刘可可;侯加根;刘钰铭;史燕青;柳琳;唐力;高兴军;周新茂【作者单位】中国石油大学北京地球科学学院,北京102249; 中国石油长城钻探工程有限公司测井公司,辽宁盘锦124011;中国石油大学北京地球科学学院,北京102249;中国石油大学北京地球科学学院,北京102249;中国石油大学北京地球科学学院,北京102249;中国石油大学北京地球科学学院,北京102249; 中国石油长城钻探工程有限公司地质研究院,辽宁盘锦124010;中国石油大学北京地球科学学院,北京102249; 中国石油华北油田公司勘探开发研究院,河北任丘062552;中国石油勘探开发研究院,北京100083;中国石油勘探开发研究院,北京100083【正文语种】中文【中图分类】TE122.2侧积层是曲流河储层中点坝内部砂体主要的渗流屏障,因此建立精确的点坝内部构型三维模型,对曲流河储层内部剩余油预测及挖潜具有重要实用价值[1-4]。
多点地质统计学python
多点地质统计学是一种统计方法,用于分析和解释地质数据中存在的空间相关性和变异性。
在Python中,有一些常用的库和工具可以用来实现多点地质统计学分析。
1. GeostatsPy:GeostatsPy是一个用于地质数据分析的Python 库,其中包含了多点地质统计学的功能。
它提供了各种常用的函数和方法,包括堆积图、变异函数、傅里叶分析等。
2. PyGSLIB:PyGSLIB是一个用于地质建模和矿业资源估计的Python库。
它包含了一些经典的多点地质统计学方法,如变异函数的建模、金属矿床建模等。
3. GS+:GS+是一款商业软件,用于地质数据分析和多点地质统计学。
它提供了丰富的功能和工具,包括数据可视化、空间模型拟合和预测等。
可以通过Python的相关库和GS+进行集成和数据交互。
4. SGeMS:SGeMS是一个面向地质建模的软件,具有强大的多点地质统计学功能。
它提供了一系列的算法和模型,如变异函数建模、空间预测等。
SGeMS支持Python接口,可以通过Python脚本进行批量处理和自定义分析。
以上是一些常用的Python库和工具,用于进行多点地质统计学分析。
根据具体的需求,选择合适的工具和方法进行分析和建模。
多点地质统计学在三维岩心孔隙分布建模中的应用张丽;孙建孟;孙志强;姜黎明;刘学锋【摘要】Multiple-point geostatistics and its principle and methods of building 3D distribution of core pore structure were presented. The 3D pore structure model was reconstructed using multiple-point geostatistical methods on the basis of 2D thin sections of Fontaineblean sandstone, and the accuracy of this method was evaluated by calculating local porosity distribution function and average percolation probability function in real 3D digital cores and the reconstructed digital cores. The results show that the characteristics of pore structures reconstructed by this method are similar to those of real 3D digital cores, and also have similar homogeneity and strong pore connectivity. The reconstruction of 3D pore structure model can be used to characterize the pore connectivity and homogeneity of the real core,but also can be used as the basis for simulation of other physical characteristics.%介绍多点地质统计学及其用于三维岩心孔隙结构建模的原理和方法.基于多点地质统计方法,以均质性较好的Fontaineblean砂岩的二维薄片为基础,重构三维孔隙结构模型,并且应用局部孔隙度理论和渗流概率函数进行准确性评价.结果表明:由此获取的孔隙结构模型与真实三维岩心孔隙分布十分相近,且具有相似的均质性和孔隙连通性;重构的三维孔隙结构模型不仅可以用来表征真实岩心的孔隙连通性和均质性,而且也可以作为其他物理特性模拟的基础.【期刊名称】《中国石油大学学报(自然科学版)》【年(卷),期】2012(036)002【总页数】5页(P105-109)【关键词】多点地质统计方法;二维薄片;三维数字岩心;孔隙分布;局部孔隙度理论【作者】张丽;孙建孟;孙志强;姜黎明;刘学锋【作者单位】中国石油大学地球科学与技术学院,山东青岛266580;中国石油大学地球科学与技术学院,山东青岛266580;山东科技大学研究生院,山东青岛266510;中国石油大学地球科学与技术学院,山东青岛266580;中国石油大学理学院,山东青岛266580【正文语种】中文【中图分类】TE19岩心是从地层获得的最直观、可靠的数据资料,其中包含的信息如粒径分布、沉积构造、孔隙结构[1]等特征,是进行数字岩石物理实验的基础。
多点地质统计学建模的发展趋势石书缘;尹艳树;冯文杰【摘要】从算法研究、训练图像处理和实际应用三个方面详细解剖了国内外多点地质统计学的发展历程,在此基础上,分析了多点地质统计学主流的几种算法的核心原理、适用范围及优缺点,以此来对储层建模的发展趋势作出展望.目前,多点地质统计学虽是随机建模的一种前沿研究热点,但由于其尚未成熟,仍需对建模算法进行研究.为此,在前人研究的基础上,重点分析了多点地质统计学的发展趋势:合理处理训练图像;合理利用软信息;选择合适的相似性方法;选择合适的标准化方法;合理利用平稳性;算法间的耦合;选择合适的过滤器;拓展缝洞型碳酸盐岩模拟.最后,提出多点地质统计学在储层建模方面,应从增加储层的模拟区域、提高模拟精度、扩大储层相的模拟范围和提高计算机模拟效率等方面进行改进.%Starting with algorithm designing, training image,and practical application,the authors analyzed multiple-point statistics research trends both in China and abroad. On such a basis, the core principles of the main four MPS algorithms, their applicable ranges and advantages as well as disadvantages were analyzed,so as to forecast the trend of reservoir modeling. Multiple-point geostatistics is an international forefront research tool in stochastic modeling; nevertheless, as the algorithm is not yet mature, it should be further improved. On the basis of previous researches, the existing problems of the multiple-point geostatistics that need to be modified and the direction of the processing of training images are proposed, such as suitable processing of training image, choice of similarity methods, choice of standardized data,smoothing,couplingamong algorithms,and expansion of the simulation extent for fractured-vuggy carbonate reservoir. In order to improve the usage of multiple-point geostatistics for reservoir stochastic modeling,we should spare no efforts to increase the simulation area,improve simulation accuracy,expand the scope of simulation of reservoir types,save simulation time and improve simulation efficiency.【期刊名称】《物探与化探》【年(卷),期】2012(036)004【总页数】6页(P655-660)【关键词】储层建模;多点地质统计学;模拟精度;缝洞型碳酸盐岩模拟;发展趋势【作者】石书缘;尹艳树;冯文杰【作者单位】中国石油勘探开发研究院,北京 100083;长江大学地球科学学院,湖北荆州 434023;中国石油大学地球科学学院,北京 102249【正文语种】中文【中图分类】TE319目前,中国东部大部分油田进入高含水开发阶段,在新区勘探难度加大的情况下,对老油田进行挖潜寻找剩余油及拓展深层碳酸盐岩勘探开发已成为当前油气勘探开发的几个主要发展方向。
多点地质统计学随机建模方法原理详细教程多点地质统计学随机建模是一种应用于地质领域的统计学建模方法,它主要用于处理地质参数在空间上的变化规律。
该方法的原理基于地质参数的随机性和空间相关性,通过构建具有地质属性的随机模型,可以模拟地质现象的空间分布。
具体而言,多点地质统计学随机建模方法主要包括以下几个步骤:1.数据准备:收集与地质参数相关的数据,例如岩性、厚度、含矿物质等。
要求数据具有一定的地质意义和空间分布规律。
2.变量描述:对收集到的数据进行统计分析,包括计算均值、方差、协方差等统计指标,以描述地质参数的分布特征。
3.变量变换:根据地质参数的实际特征,对数据进行变换,使其符合正态分布、对数正态分布或其他分布类型。
4.空间相关性建模:通过计算地质参数之间的空间相关性,可利用协方差函数、变差函数或半方差函数等,建立地质参数之间的空间相关模型。
5.随机模拟:根据变量的统计特征和空间相关模型,结合随机数生成算法,生成符合实际情况的具有随机性和空间相关性的地质参数数据。
6.模型验证:对生成的地质模型进行验证,比较随机模拟结果与实际数据的吻合程度。
可以使用统计指标如均值、方差、协方差等进行对比分析。
7.地质模型应用:根据随机模拟结果,可以进一步进行岩层插值、矿产资源评估和地质灾害风险评估等相关研究及应用。
总的来说,多点地质统计学随机建模方法是将统计学原理应用于地质参数的空间分布建模,通过对地质参数的统计特征和空间相关性的建模,生成具有随机性和空间相关性的地质模型。
这种方法可以提供地质领域研究的基础数据和分析手段,为地质灾害风险评估、资源勘探和环境评价等问题提供科学依据。
多点地质统计学Multiple-point geostatistic是相对于传统的两点地质统计学而言的,主要应用于储层表征与建模中.
传统的地质统计学在储层建模中主要应用于两大方面:其一,应用各种克里金方法建立确定性的模型,这类方法主要有简单克里金、普通克里金、泛克里金、协同克里金、贝叶斯克里金、指示克里金等;其二,应用各种随机建模的方法建立可选的、等可能的地质模型,这类方法主要有高斯模拟(如序贯高斯模拟)、截断高斯模拟、指示模拟(如序贯指示模拟)等。
上述方法的共同特点是空间赋值单元为象元(即网格),故在储层建模领域将其归属为基于象元的方法。
这些方法均以变差函数为工具,亦可将其归属为基于变差函数的方法。
变差函数局限性(传统地质统计学)变差函数只能把握空间上两点之间的相关性,亦即在二阶平稳或本征假设的前提下空间上任意两点之间的相关性,因而难于表征复杂的空间结构和再现复杂目标的几何形态(如弯曲河道)。
弯曲河道的3种不同的空间结构(图1a,b,c)在横向上(东西方向,图1d)和纵向上(南北方向,图1e)的变差函数十分相似,这说明应用变差函数不能区分这3种不同的空间结构及几何形态,因此,基于变差函数的传统地质统计学插值和模拟方法难于精确表征具有复杂空间结构和几何形态的地质体。
现有的储层随机建模的另一途径是基于目标的方法,它是以目标物体为基本模拟单元,进行离散物体的随机模拟(Haldorsen and Damsleth,1990;Holdenet al.,1998)。
主要方法为示性点过程(亦称标点过程),其根据先验地质知识、点过程理论及优化方法(如模拟退火)表征目标地质体的空间分布,因此这种方法可以较好地再现目标体几何形态。
但这种方法亦有其不足:1)每类具有不同几何形状的目标均需要有特定的一套参数(如长度、宽度、厚度等),而对于复杂几何形态,参数化较为困难;2)由于该方法属于迭代算法,因此当单一目标体内井数据较多时,井数据的条件化较为困难,而且要求大量机时
2多点地质统计学的基本概念
多点统计学着重表达多点之间的相关性。
建模方法
Guardiano and Srivasta-va(1993)提出了一种直接的(非迭代)算法,从训练图像中直接提取局部条件概率,并应用序贯指示模拟方法产生模拟实现。
由于该算法为非迭代算法,不存在收敛的问题,因而算法简单。
但由于在每模拟一个网格节点时均需重新扫描训练图像,以获取特定网格的局部条件概率,因此严重影响计算速度,难于进行实际应用。
Strebelle and Journel(2001)将算法加以改进,应用一种动态数据结构即“搜索树”一次性存储训练图像的条件概率分布,并保证在模拟过程中快速提取条件概率分布函数,从而大大减少了机时。
基于此,提出了多点统计随机模拟的Snesim算法(Strebelle and Journel,2001;Strebelle,2002)。
其建模基本步骤如下:
1)建立训练图像。
2)准备建模数据,将实测的井数据标注在最近的网格节点上。
3)应用自定义的与数据搜索邻域相联系的数据样板τn扫描训练图像,以构建搜索树。
4)确定一个访问未取样节点的随机路径。
在每一个未取样点u处,使得条件数据置于一个以u为中心的数据样板вn中。
令n′表示条件数据的个数,dn′为条件数据事件。
从搜索树中检索c(dn′)和ck(dn′)并求取u处的条件概率分布函数。
5)从u处的条件概率分布中提取一个值作为u处的随机模拟值。
该模拟值加入到原来的条件数据集中,作为后续模拟的条件数据。
6)沿随机路径访问下一个节点,并重复3)、4)步骤。
如此循环下去,直到所有节点都被模拟到为止,从而产生一个随机模拟实现。
7)改变随机路径,产生另一随机模拟实现。
多点地质统计学随机模拟方法(如Snesim算法)与传统的地质统计学随机模拟方法(如序贯指示模拟SIS)的本质差别在于未取样点处条件概率分布函数的求取方法不同。
前者应用多点数据样板扫描训练图像以构建搜索树并从搜索树中求取条件概率分布函数(上述第1步和第3步),而后者通过变差函数分析并应用克里金方法求取条件概率分布函数。
训练图像的建立
训练图像为研究区的定量地质模式。
一般地,单纯应用研究区井资料难于建立训练图像,而需要综合应用研究区资料及原型模型。
根据井资料及区域沉积相研究资料,研究区新近系明化镇组沉积相主要为中—低弯度河流相,沉积微相可分为河道砂岩相、河间溢岸砂岩相与泛滥平原泥岩相等3个微相,物源方向为北西—南东向。
河道与溢岸砂体呈迷宫状特征分布,河道砂体呈条带状,单河道宽度一般小于井距;溢岸砂体为透镜状。
定量的几何学特征需要结合原型模型来获取。
研究区沉积特征与黄骅坳陷港西开发区新近系明化镇组十分相似,两区同属一个含油气盆地(渤海湾盆地),层位及沉积相类型亦相同,具有较好的可比性。
港西开发区属于开发中后期的油田,其中一个区块具有井距50 m的密井网,对此曾作过系统的沉积微相及储层地质研究并建立了相应的地质模型(武军昌等,2002)。
在该原型模型中,河道砂体主要为条带状(可分叉合并),单河道砂体宽度约50~300 m,一般为100~200 m;河间漫溢砂体则呈透镜状。
单一漫溢砂体长、宽一般小于200 m。
这一密井网地质模型(微相几何学特征及组合模式)可作为研究区沉积微相研究与建模的原型模型。
综合两方面资料,应用数字化成图工具(Rc-tool),构建了研究区目的层的训练图像(如图3)。
从图中可以看出,训练图像反映了微相的定量分布模式(其中,棕色代表河道微相,浅兰色代表溢岸微相,黄色代表河道间微相),但并不要求忠实井数据,只要求反映储层变化的空间结构性,其作用相当于变差函数,但后者只反映空间两点间的结构性。
图4为根据训练图像及井资料手工编绘的研究区沉积微相分布图,该图反映了微相砂体的基本形态。
4.2多点统计随机建模应用
Snesim算法,对研究区沉积微相进行随机建模,包括准备数据、扫描训练图像以构建搜索树、选择随机路径、序贯求取各模拟点的条件概率分布函数并通过抽样获得模拟实现。
这一过程通过程序来实现。
图5为其中的一个模拟实现(其中,棕色代表河道微相,浅兰色代表溢岸微相,黄色代表河道间微相)。
从图中可以看出,模拟实现反映了训练图像的结构性,同时基本再现了微相砂体的几何形态(与图4相比较),并完全忠实于井信息。
为了对比,我们应用基于变差函数的序贯指示模拟方法对研究区目的层进行了随机模拟,如图6所示,该结果没有再现微相砂体的几何形态,而且部分微相砂体呈零散分布。
由此可知,多点统计随机建模比传统二点统计学建模方法具有明显的优越性。