储层随机建模方法研究毕业论文

以二区六区块明3-5-1单层为例,通过二维相建模对设计的基于储层目标骨架的多点地质统计学储层随机建模方法以及Simpat方法进行对比分析。

1.建模前参数的准备

多点地质统计学储层随机建模方法最重要的输入参数就是训练图像。前文已经探讨了训练图像建立的过程,并通过地震、地质以及测井等资料建立了明3-5-1单层的训练图像,这里就不再赘述了。

图5-24 明3-5-1单层微相的训练图像(130*85网格)

另外,数据样板选择也比较重要。数据样板太小不能考虑太多的信息,而数据样板太大则增加了计算量。数据样板一般选择和目标体宽度接近为好。从明3-5-1单层沉积微相平面分布分析可以得之,河道宽度在100~240m之间,如果采取20m网格,则河道占据5~12个网格大小,因此采用7*7数据样板、两重网格搜索策略即可以满足要求。

在建模过程中还需要准备条件数据,对于研究区来说就是井数据和地震数据。对研究区明3-5-1进行优势相分析,判别准则为:如果井资料存在河道条件数据,则认为在此井处为河道微相;如果没有河道微相,而有溢岸相(天然堤、决口扇等),则定为溢岸相;如果仅存在泥岩,则定为泛滥平原。经过这样处理,获得了研究区56口井在明3-5-1单层的优势相。对优势相进行了编码,泛滥平原用0表示,河道用1表示,溢岸用2表示。

从地震振幅分布图看,由于在研究区北部断层较发育,影响了地震效果。只是在中部断层较少,构造较平稳,地震与砂体有较好的对应关系。而中部井资料也相对较丰富。因此,本次研究仅综合地震属性识别的主流线,指导河流骨架模型建立。

图5-25 结合地震及沉积微相预测的明3-5-1单层河道骨架模型

对于将地震属性作为软数据约束建模,由于断层影响了部分区域的地震效果,因此暂不做考虑。从地震沿层切片图上,可以较清楚看到中部河流分汊和汇合特征。在北上方,从北北西方向以及北北东方向流入两条河流,并在MQ4-8井附近汇合;随后继续向南方向流动至G1-64井附近产生分汊,分别向南和南东方向流出研究区。结合井资料获得砂体等厚图以及沉积微相分布图,采用数字化工具,获得了研究区河道流线(图5-25、图5-26)。在随后基于储层骨架的方法中,此河道主流线将作为河道骨架模型,指导不和约束随机建模。

图5-26 利用Simpa方t法和基于河道骨架的方法建立的明3-5-1单层微相模型比较

2.明3-5-1单层随机建模

利用Simpat方法和基于储层骨架的方法,分别建立起明3-5-1单层沉积微相储层地质模型(图5-26),以比较两种方法的建模效果。从模型效果来看,基于储层骨架的方法以及Simpat方法都较好地再现了河道形态及空间分布,模拟结果基本与实际勾绘的沉积微相分布图具有较高的相似性。但通过比较发现,基于储层骨架的方法模拟的河道在研究区连续分布,不存在断续的现象,河道分汊合并现象在模拟实中现得到了很好的反映;而Simpat方法模拟的河道整体特征也与实际微相分布基本相符合,但是河道出现了中断,而河道分汊合并现象再现效果也不是很好。这说明基于骨架的方法较Simpat方法在再现储层目标形态及分布上具有较明显的优势。Simpat方法以及基于储层骨架的方法都是通过模拟数据模式来再现储层目标的。数据模式反映的是储层形态特征,因此储层形态再现好坏可以通过数据模式选择正确与否来评价。前已述及,数据模式选择正确与否可以从模拟匹配度图来衡量,分别计算了各自的匹配度图(图5-27)。通过比较Simpat方法以及基于储层骨架的方法模拟实现的匹配度图可以看出,基于储层骨架的方法的匹配度图值较低,而Simpat方法的模拟匹配度图值较高,这表明基于储层骨架的方法相对于Simpat方法其选择的数据模式整体准确性高。此外,通过比较发现尤其是在条件数据较少的地方,在骨架模型指导下,数据模式选择具有较高的准确性,反映在匹配度图上就是较低的值,如图中红色方框标定的地方。这也从侧面说明骨架模型能够更好的约束随机建模,获得更加满意的地质模型。

图5-27 Simpat方法和基于河道骨架预测的方法的匹配度图比较1—河道;2—溢岸;3—泛滥平原

利用基于储层骨架的多点地质统计学随机建模方法建立了明3-5-1单层微相的10个实现(图5-28),求取了河道出现的概率(图5-29)。从10个实现结果来看,模拟结果具有较高的相似性,其分布与勾绘的沉积微相平面分布吻合程度较高。不过在河道边部以及井资料较少区域,模拟实现之间还是有一定的差异的。表明模拟实现能够真实反映河道分布及其不确定性。从河道概率图上可以看出,河道出现概率高的区域与实际勾绘的沉积微相分布极为相似,表明基于储层骨架的多点地质统计学随机建模方法正确预测了河道分布。

以上建模结果分析表明,与Simpat方法相比,基于储层骨架的多点地质统计学随机建模方法能够更好再现储层的目标形态及分布,因此,更适用于储层目标体的地质建模。

图5-28 利用基于储层骨架的方法获得的10个实现

图5-29 利用基于河道骨架预测的方法建立的10个实现获得的河道概率

随机建模（Stochastic Modeling）方法承认地质参数的分布有一定的随机性，而人们对它的认识总会存在一些不确定的因素，因此建立地质模型时考虑了这些随机性引起的多种可能出现的实际情况，供地质人员选择。

随机建模方法认为，在现有技术情况下，对地下储层的认识存在一定的不确定性，一是已知资料控制点有限，以300m井距井网为例，井孔提示的储层体积所占整个储层体积，以百万至千万分之一数量级计，绝大部分储层性质是依靠这些少数已知点去推测的；二是描述这些控制点储层性质的技术本身还存在一定的误差，如测井解释渗透率，经常可达数倍的误差。随机建模方法同时又认为，作为地质体的储层，其各项属性的非均质分布，由于其有一定的地质成因，应存在一定的地质统计特征，用这一地质统计特征去表征储层非均质性的总体面貌，而不追求每一个预测点的确定的数值，仍然在一定时间、一定条件下可以为油气田开发提供合理的地质模型，保证流体流动模拟的可信和开发决策的正确。

8.1.3.1 随机建模的类别

储层随机建模通常又分为条件模拟和非条件模拟。其根本区别在于条件模拟较非条件模拟不仅要求模拟产生的储层随机图像（包含储层分布和物性等方面信息的图像）符合实际资料所观测到的储层属性空间分布的相关结构（地质统计特征），而且要求在井位处（或资料点处）的模拟结果与实际资料一致。通常讲的随机模拟一般指条件模拟。随机模拟方法分为以下两类：

1）离散性模拟方法：离散性模拟主要建立储层岩相的分布模型，用来描述离散性的地质特征，包括确定储层、隔层、砂体（储渗体）的空间分布边界和空间几何形态等。实际上就是实现气藏描述中的储层分布预测。所采用的模拟方法包括：示点性过程模拟、马尔可夫-贝叶斯指示模拟、序贯指示模拟、镶嵌过程模拟、截断高斯模拟等方法；对于非条件模拟，则可采用布尔模拟。

2）连续性模拟方法：连续性模拟主要建立岩相边界控制下的储层参数（孔隙度、渗透率、含水饱和度、泥质含量、碳酸盐含量等）的分布模型，即油气藏描述中的储层参数预测。所采用的模拟方法包括：退火模拟、序贯指示模拟、分形随机函数法、高斯随机函数法以及马尔可夫随机域法等。

对于非条件模拟，则可采用转带法。

8.1.3.2 储层随机建模的基本步骤

储层建模工作的实施主要包括以下三个基本步骤。

（1）建立储层原型模型

建立储层原型模型是随机建模的基础，所谓原型模型就是储层的实体地质模型，任何油藏（储层）描述方法都是只由零散信息对储层实体所进行的一种推断，这种推断可以是确定性的（如地震储层横向预测），也可以是不确定性的（如统计推断）。在不了解地质实体的前提下，任何一次研究结果，只能看作一次对地质实体的随机抽样，抽样结果的准确性依赖于统计的概率把握程度，这种把握程度只能来自于原型模型的建立。

储层原型模型的建立就是为了构筑一个与实际储层尽可能接近的储层信息标准答案库，从可见的实体模型描述入手，来建立各种地质知识库（这其中包括了各种储层的边界和储层参数的空间分布），建立相应的先验概率知识，如参数分布的范围、均值、方差、分布函数等。只有这样，储层随机建模才有依据。

原型模型的建立方法较多，目前主要采用的有：

1）物理模拟——以水槽模拟为主要代表；

2）野外露头精细描述——国内外已广泛开展；

3）现代沉积研究——在沉积学领域已有大量实例，是构筑沉积学理论的重要基础；

4）密井网精细对比与描述——主要在老开发区进行；

5）地震资料的确定性建模方法——主要依靠地震资料空间大信息量的优势，依靠资料处理，确定储层分布的宏观模型，重点是砂体的分布，同时也能对储层孔隙度、渗透率参数进行趋势性的估计。

（2）建立储层的随机模型

取得了储层原型模型以后，就可以建立储层的随机模型，它是以反映储层各项特征的参数统计为手段，建立相应的概率模型，如储层厚度、孔隙度、渗透率、含油饱和度等参数的分布规律和空间结构。对参数分布规律的认识主要以传统概率统计为基础，确定参数分布的大小范围、均值、方差、分布函数类型等，进而对空间结构进行分析（变差函数的计算）。

通过对储层特征建立随机模型，可以把各种地质认识（定性描述）和观测数据有机地结合起来，并可以反映由于信息缺乏而引起的不确定性。在已经建立的随机模型的基础上，再进行随机模拟，产生出反映储层非均质性的一系列等概率实现。每个实现就是一种可能的储层参数的空间分布，它们之间的差异反映了随机模型中所包含的不确定性，也就是我们常常谈到的研究中的多解性问题。

（3）储层的随机模拟

建立了储层随机模型后，就可以进行储层的随机模拟，随机模拟分为条件模拟和非条件模拟两种。非条件模拟只是要求再现地质特征的空间分布规律及相关性，而条件模拟不仅要求再现地质特征的空间分布和相关性，而且还要求在抽样位置上与实测数据一致或在指定位置上具有指定的特性。

对于不同的储层属性，具有不同的随机模型，应采用不同的模拟方法。由于大型计算机的出现，使细网格和高维空间的模拟得以实现，在实际应用中，寻求一种快速有效的模拟算法成为众多的研究者所探求的目标。

8.1.3.3 储层随机建模的基本流程

储层随机建模一般分为两个阶段进行，即先采用离散型模拟方法，建立储层的骨架模型；然后在储层骨架模型边界的控制下，应用针对连续性变量（如储层物性）的模拟方法建立储层参数模型。这就是目前大多数研究者使用的两阶段建模的基本流程。

陈恭洋［4］根据两阶段建模的思路，提出了一个基本的随机建模流程（图8.1），该流程图中包括了9个方面的研究内容。

图8.1 储层随机建模总体设计流程框图［4］

1）地层模型：以克里格插值技术为基本手段，主要研究储层顶、底界面的空间展布特征，并通过地质统计对比确定小断层带的空间分布。大的断层可由地震资料解释予以确定。该项研究主要提供后续储层和油气藏模拟的大的边界信息。

2）沉积相分析：包括大相和微相分析两部分研究内容，并以后者为研究重点。大相分析以区域沉积背景知识为指南，结合地震相的分析，明确研究工区较大范围内的沉积体系及空间展布特征。最后确定出油气藏范围内储层所处的相带沉积部位，为微相研究奠定坚实的基础。

微相分析重点研究沉积成因单元的结构要素及其组合型式以及它们的空间展布规律，为储层随机建模提供必要的地质先验知识，主要依据沉积学的研究手段进行。

3）高分辨率层序地层分析：主要应用于油气藏规模的储层对比技术，依靠岩心和测井资料，进行开发阶段的储层表征中储层的精细对比。因为储层岩性、几何形态、连续性及岩石物理特征等是在沉积物堆积过程中产生的，精确的地层对比可以在四维空间中对这些特征有更清楚的认识，高分辨地层对比是识别非均质性的有效方法。另外，具时间意义的地层界面通常与流体流动单元的岩石物理面相一致，可通过精细地层对比，划分流动单元。随着时间分辨率的提高，对地层形态和规模、相的位置和岩石物理特征的预测也就更加精确。与沉积相的分析相结合，是目前油田覆盖区建立储层原型地质模型最有效的方法。

4）储层岩相分布的离散型随机模拟：这是储层随机建模的核心内容之一，一般作为储层随机建模的第一步，为储层参数空间分布的连续性模拟提供边界控制信息。序贯指示模拟（SISIM）和示性点过程模拟（MPPS）被认为是两种有效的研究方法。序贯指示模拟以指示理论为基础，将各种沉积微相带视为空间分布的离散性随机变量，进行地质统计学的条件模拟，其缺点是难以描述储层的形态特征。而示性点过程模拟是一种面向对象的方法，十分符合沉积学的思想和推理过程，将沉积学研究所认定的储层砂体几何形态、位置、大小、连通方式等储层参数作为服从一定分布的离散型随机变量，建立相应的随机模型进行随机模拟，其缺点是难以实现条件模拟。将两者有机地结合起来可能是一种好的途径。

5）测井和地震资料处理：这方面的技术已在现代油气藏描述中被大量采用。更重要的是补充建模时仅依靠井点信息的不足，使储层建模不仅在油气藏开发阶段发挥重大作用，而且在勘探的各个时期也能充分发挥作用，提供新的储层预测方法。

6）分形和地质统计学条件模拟：这是解决储层参数空间分布的关键性模拟方法。地质统计学模型可以很好地刻画储层参数分布的空间结构和变异性。而分形方法则能精确地表征储层的非均质性，并能克服由克里金方法所带来的光滑效应。两者的结合已被大量的研究实例证明是一种有效的储层预测途径。

7）网格粗化：储层建模阶段的细网格模拟可以尽可能精细地提示储层的非均质特征。但遗憾的是，在油藏动态模拟器中，由于受到目前计算能力的限制，难以接受这种细网格的参数输入。因此，必须进行网格的粗化，粗化的准则一般需要考虑到储层孔隙容积和储层的渗流能力（即孔隙度和渗透率），其中尤以储层对流体传导能力（渗透率）的近似最为关键。

8）油气藏数值模拟动态拟合与静态资料约束决策：这是对前述储层随机建模所产生的多幅等概率实现的图像进行优选决策的过程。研究的重点并不在于动态模拟，因此无需考虑复杂条件下的数值模拟问题。主要是对油气藏压力、产油气量和含水率三项参数进行历史拟合，并结合静态地质资料的各项条件约束（包括储层参数的统计规律和地质认识等），选取一个最符合动态和静态条件的随机图像作为所建立的储层地质模型。这一模型是以各种参数场的形式所表示的。

9）三维可视化：即将前面所建立的反映储层地质模型的各种参数场通过计算机进行三维成像或制图。目前，三维可视化的研究与设计已经成为计算机成像领域中的一项热门课题，它使所取得的成果大大地增强了油气藏的研究与管理的可操作性和直观性。

综上所述，储层建模实际上是对油田各类数据资料通过计算机技术进行有效的综合。因此，从地质角度上讲，要形成一套比较先进而有效的建模方法，更大程度上还是要依赖于先进的地质、地球物理和分析测试资料处理技术来获取可靠的输入参数。

8.1.3.4 储层随机建模的软件系统

在随机模型方法和理论发展的同时，模拟软件也得到了一定的发展，美国斯坦福大学、墨西哥矿业技术学院、荷兰皇家/壳牌公司、雪飞龙公司、GeoQuest公司等都开发和研制了自己的地质统计学和储层模拟软件。加拿大GeoStat系统公司和McGill大学联合推出了智能模拟或专家系统软件GeoStat，法国石油研究院和地质统计中心联合开发的HERESIM软件包也取得了较大的影响。这些软件的主要功能如下：

1）以转带法和指示克里格法相结合，用于储层的横向和垂向对比，其数学基础是Bessel函数和指示相关函数（美国墨西哥矿业技术学院开发TUBA软件）；

2）用于SGI图形工作站的地质模型软件，其特色是可以采取任意切片的方法来展现储层孔隙度、渗透率和砂体在连续断面或切片上的分布特征，其数理基础是随机模拟（美国Strata-Model公司研制SGM软件）；

3）以条件概率法为基础设计，主要用于模拟砂岩油藏中的三维储层的连通性和构形（荷兰皇家/壳牌集团公司推出MONARCH软件）；

4）以BP神经网络技术为主、依据地质统计学和地震特征进行随机建模的软件，其关键方法是分析并拟合储层物理特性和岩石属性的直方图和变差函数分布，求出它的特征值，以建立数学模型（荷兰Jason公司推出Stat Mod软件）；

5）将地质统计和智能模拟技术相结合，不仅包括各种数值运算、多元统计，还包含可引导、承担、评价和推断地质统计运行的知识和专家经验。因此，该软件具有两大特色：一是储层地质特性模拟及立体化定量显示；二是具有地质解释中的专家知识和经验（加拿大GeoStat系统公司和McGill大学联合推出GeoStat系统）。

上述软件都在各自的使用中发挥了很大的效益，也取得了不少有意义的成果。尽管每套软件各有侧重，但考察它们的共同之处，主要体现在三个方面：①强调储层描述的高度定量化，体现了油气储层研究已从定性发展到了定量的水平；②均从储层骨架分布和储层参数特征两个方面进行建模，把握了储层特征的关键要素；③体现了多学科、多信息的综合研究趋势。因此，从储层建模软件的发展，也显示出了储层随机建模在当前油气勘探开发研究中的重要意义和良好前景。