一种基于乱序矩阵的半确定性裂缝建模方法及系统

摘要

本发明涉及裂缝建模技术领域，公开了一种基于乱序矩阵的半确定性裂缝建模方法及系统，所述方法包括：提取地震数据中裂缝敏感性高的属性数据；根据所述属性数据构建乱序权值矩阵和均匀随机变量矩阵；根据所述乱序权值矩阵、均匀随机变量矩阵和预设的裂缝片数量生成裂缝片的中心点；根据获取的裂缝片的走向和倾角、以及中心点构建旋转角和旋转轴；根据旋转角和旋转轴构建旋转矩阵，旋转矩阵用于表征裂缝片模型。本发明通过乱序权值矩阵和均匀随机变量矩阵，能够减少循环次数，进而能有效降低运算时间，提高建模的效率。且旋转矩阵能够实现对裂缝片产状进行准确描述，使所建立的裂缝模型更加接近真实情况，为油气藏的开采提供信息。

权利要求

1.一种基于乱序矩阵的半确定性裂缝建模方法，其特征在于，所述方法包括：提取地震数据中裂缝敏感性高的属性数据；

根据所述属性数据构建乱序权值矩阵和均匀随机变量矩阵；

构建乱序权值矩阵，包括：

对属性数据进行乱序排列，得到乱序权值数据；

对乱序权值数据依次进行极差正则化处理和归一化处理，得到乱序权值矩阵；

根据所述乱序权值矩阵、所述均匀随机变量矩阵和预设的裂缝片数量生成裂缝片的中心点，包括：判断乱序权值矩阵与均匀随机变量矩阵中相同位置的数值点大小；

确定为乱序权值矩阵的数值点大于均匀随机变量矩阵的数值点，则保留乱序权值矩阵的数值点作为裂缝种子点，根据裂缝种子点生成裂缝种子点矩阵；

根据预设的裂缝片数量从裂缝种子点矩阵提取对应数量裂缝种子点作为中心点；

根据获取的裂缝片的走向和倾角、以及所述中心点构建旋转角和旋转轴；

根据旋转角和旋转轴构建旋转矩阵，包括：根据旋转角和旋转轴构建四元数；

根据四元数构建旋转矩阵；

所述旋转矩阵用于表征裂缝片模型。

2.根据权利要求1所述的一种基于乱序矩阵的半确定性裂缝建模方法，其特征在于，在得到乱序权值矩阵时，同时构建与乱序权值矩阵相同大小的均匀随机变量矩阵。

3.根据权利要求1所述的一种基于乱序矩阵的半确定性裂缝建模方法，其特征在于，根据获取的裂缝片的走向和倾角、以及所述中心点构建旋转角和旋转轴，包括：根据中心点确定旋转轴的中心，并根据走向和倾角确定旋转轴的坐标位置；

根据坐标位置和中心构建旋转轴；

根据旋转轴的坐标位置构建旋转角。

4.根据权利要求3所述的一种基于乱序矩阵的半确定性裂缝建模方法，其特征在于，所述走向和倾角根据地质露头和/或测井裂缝成像获得。

5.一种基于乱序矩阵的半确定性裂缝建模系统，所述系统用于实现权利要求1‑4中任一项所述的基于乱序矩阵的半确定性裂缝建模方法，其特征在于，所述系统包括：提取模块，用于提取地震数据中裂缝敏感性高的属性数据；

第一构建模块，用于根据所述属性数据构建乱序权值矩阵和均匀随机变量矩阵，构建乱序权值矩阵为：对属性数据进行乱序排列，得到乱序权值数据；

对乱序权值数据依次进行极差正则化处理和归一化处理，得到乱序权值矩阵；

生成模块，用于根据所述乱序权值矩阵、所述均匀随机变量矩阵和预设的裂缝片数量生成裂缝片的中心点，具体为：判断乱序权值矩阵与均匀随机变量矩阵中相同位置的数值点大小；

确定为乱序权值矩阵的数值点大于均匀随机变量矩阵的数值点，则保留乱序权值矩阵的数值点作为裂缝种子点，根据裂缝种子点生成裂缝种子点矩阵；

根据预设的裂缝片数量从裂缝种子点矩阵提取对应数量裂缝种子点作为中心点；

第二构建模块，用于根据获取的裂缝片的走向和倾角、以及所述中心点构建旋转角和旋转轴；

第三构建模块，用于根据旋转角和旋转轴构建旋转矩阵，具体为：根据旋转角和旋转轴构建四元数；

根据四元数构建旋转矩阵；

所述旋转矩阵用于表征裂缝片模型。

说明书

一种基于乱序矩阵的半确定性裂缝建模方法及系统

技术领域

[0001] 本发明涉及裂缝建模技术领域，具体涉及一种基于乱序矩阵的半确定性裂缝建模方法及系统。

背景技术

[0002] 裂缝是致密储层的有效储集空间和主要渗流通道，影响致密油气藏的井网部署与注水开发效果，裂缝分布规律研究以及三维模型建立是致密油气藏开发的重要地质依据。

[0003] 利用随机方法模拟裂隙网络最早于二十世纪八十年代末由Neuman提出，该方法早期多应用于模拟等效连续介质模型，九十年代后期逐步应用于裂隙网络的随机模拟。目前，随机模拟被认为是刻画介质非均质性的一种高效可行的方法，已广泛应用于采矿工程、水资源工程、石油工程和核废料地质处置等多个领域。尤其在模拟复杂地质结构且数据资料较少时，随机模拟方法有很好的效果。其中离散裂缝网络建模最为常用，1983年Baecher提出以自己名字命名的模型，1995年Ivanova提出了GEOFRAC模型。DFN多以Baecher模型为主，该模型基于示性点随机过程的布尔模拟是一种常用的随机模拟方法。点过程确定裂缝的中心位置，示性过程确定点的属性，如裂缝形状、倾角、倾向、开度等属性。2018年，国内学者董少群提出基于裂缝密度约束的离散裂缝网络建模方法，通过该方法可以实现对裂缝发育位置的约束。

[0004] 目前，就该方法在国内外的应用情况而言，这种建模方法有着两个明显的缺陷，致使该方法在使用时存在着一定的局限性。第一个缺点是裂缝密度情况通常需要考察野外露头情况或地震属性再次反演获得，步骤繁复且不够准确。另外一点则是随机约束方法裂缝片建模方法存在生成耗时长、速度慢、准确性低的问题。

发明内容

[0005] 本发明的目的在于提供一种基于乱序矩阵的半确定性裂缝建模方法及系统，解决现有裂缝片建模方法存在生成耗时长、速度慢、准确性低的问题。

[0006] 为实现上述发明目的，本发明所采用的技术方案是：一种基于乱序矩阵的半确定性裂缝建模方法，所述方法包括：

[0007] 提取地震数据中裂缝敏感性高的属性数据；

[0008] 根据所述属性数据构建乱序权值矩阵和均匀随机变量矩阵；

[0009] 根据所述乱序权值矩阵、所述均匀随机变量矩阵和预设的裂缝片数量生成裂缝片的中心点；

[0010] 根据获取的裂缝片的走向和倾角、以及所述中心点构建旋转角和旋转轴；

[0011] 根据旋转角和旋转轴构建旋转矩阵，所述旋转矩阵用于表征裂缝片模型。

[0012] 优选地，所述方法还包括：构建乱序权值矩阵，包括：

[0013] 对属性数据进行乱序排列，得到乱序权值数据；

[0014] 对乱序权值数据依次进行极差正则化处理和归一化处理，得到乱序权值矩阵。

[0015] 优选地，在得到乱序权值矩阵时，同时构建与乱序权值矩阵相同大小的均匀随机变量矩阵。

[0016] 优选地，根据所述乱序权值矩阵、所述均匀随机变量矩阵和预设的裂缝片数量生成裂缝片的中心点，包括：

[0017] 判断乱序权值矩阵与均匀随机变量矩阵中相同位置的数值点大小；

[0018] 确定为乱序权值矩阵的数值点大于均匀随机变量矩阵的数值点，则保留乱序权值矩阵的数值点作为裂缝种子点，根据裂缝种子点生成裂缝种子点矩阵；

[0019] 根据预设的裂缝片数量从裂缝种子点矩阵提取对应数量裂缝种子点作为中心点。

[0020] 优选地，根据旋转角和旋转轴构建旋转矩阵，包括：

[0021] 根据旋转角和旋转轴构建四元数；

[0022] 根据四元数构建旋转矩阵。

[0023] 优选地，根据获取的裂缝片的走向和倾角、以及所述中心点构建旋转角和旋转轴，包括：

[0024] 根据中心点确定旋转轴的中心，并根据走向和倾角确定旋转轴的坐标位置；

[0025] 根据坐标位置和中心构建旋转轴；

[0026] 根据旋转轴的坐标位置构建旋转角。

[0027] 优选地，所述走向和倾角根据地质露头和/或测井裂缝成像获得。

[0028] 本发明还提供一种基于乱序矩阵的半确定性裂缝建模系统，所述系统用于实现上述的基于乱序矩阵的半确定性裂缝建模方法，所述系统包括：

[0029] 提取模块，用于提取地震数据中裂缝敏感性高的属性数据；

[0030] 第一构建模块，用于根据所述属性数据构建乱序权值矩阵和均匀随机变量矩阵；

[0031] 生成模块，用于根据所述乱序权值矩阵、所述均匀随机变量矩阵和预设的裂缝片数量生成裂缝片的中心点；

[0032] 第二构建模块，用于根据获取的裂缝片的走向和倾角、以及所述中心点构建旋转角和旋转轴；

[0033] 第三构建模块，用于根据旋转角和旋转轴构建旋转矩阵，所述旋转矩阵用于表征裂缝片模型。

[0034] 本发明还提供一种电子设备，包括存储器、处理器以及存储在所述存储中并可在所述处理器上运行的计算机程序，所述处理器执行该计算机程序时实现上述的基于乱序矩阵的半确定性裂缝建模方法。

[0035] 本发明还提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被所述处理器执行时实现上述的基于乱序矩阵的半确定性裂缝建模方法。

[0036] 本发明的有益效果集中体现在：

[0037] 1、本发明通过乱序权值矩阵和均匀随机变量矩阵，能够减少循环次数，进而能有效降低运算时间，提高建模的效率。

[0038] 2、本发明的旋转矩阵能够实现对裂缝片产状进行准确描述，使所建立的裂缝模型更加接近真实情况，为油气藏的开采提供信息。

附图说明

[0039] 图1是本发明一种实施方式提供的基于乱序矩阵的半确定性裂缝建模方法的流程图；

[0040] 图2是本发明一种实施方式提供的基于乱序矩阵的半确定性裂缝建模系统的框图；

[0041] 图3是本发明100裂缝片数量的裂缝模型图；

[0042] 图4是本发明1000裂缝片数量的裂缝模型图；

[0043] 图5是本发明10000裂缝片数量的裂缝模型图；

[0044] 图6是本发明50000裂缝片数量的裂缝模型图；

[0045] 图7是本发明建模方法与现有建模方法的建模时间对比图。

具体实施方式

[0046] 为了使本领域的技术人员更好地理解本发明的技术方案，下面结合附图和具体实施例对本发明作进一步的详细说明。

[0047] 图1是本发明一种实施方式提供的基于乱序矩阵的半确定性裂缝建模方法的流程图，如图1所示，一种基于乱序矩阵的半确定性裂缝建模方法，所述方法包括：

[0048] 步骤S101：提取地震数据中裂缝敏感性高的属性数据；

[0049] 在本实施例，根据地震属性数据结果，选取曲率、曲度、方差进行进一步优选。曲率属性可以有效的刻画地层的挠曲、褶皱、断裂等地质构造。曲度属性通过拉普拉斯算子将振幅的变化参与到曲率的计算过程中，通过观察振幅值的大小变化识别断裂或沉积体。地震方差属性体使用层位标定建立地震层序，确定地层地质构造描述范围、平面网格及模拟层的划分、地质模型的建立。地震方差体属性可以很好的描述储集层砂体展布和地震道特征等方面，能够表达出地质构造间不连续的断层与褶皱关系。

[0050] 根据实际地震属性体计算，可确定砂岩裂缝以断裂活动成因的构造缝和构造形变成因的构造缝两大组成为主，多属于斜缝和高角度缝；沉积成因缝以低角度缝和水平缝为主；倾角大于30度的构造缝与低于30度的沉积成因缝可组成网状缝，利于高产；从而选定方差和曲度两个作为致密砂岩裂缝敏感属性数据。

[0051] 步骤S102：根据所述属性数据构建乱序权值矩阵和均匀随机变量矩阵；

[0052] 作为本实施例的进一步优化，所述方法还包括：构建乱序权值矩阵，包括：

[0053] 步骤A1：对属性数据进行乱序排列，得到乱序权值数据；即将敏感的属性数据乱序排列，保留属性数据原有分布模式，则保证整个数据的概率分布不变。

[0054] 步骤A2：对乱序权值数据依次进行极差正则化处理和归一化处理，得到乱序权值矩阵；在得到乱序权值矩阵时，同时构建与乱序权值矩阵相同大小的均匀随机变量矩阵。

[0055] 在本实施例中，归一化处理的公式为：

[0056] ；

[0057] 式中， 表示归一化后的乱序权值数据， 表示原始的乱序权值数据，i、j分别代表坐标位置，min(A)表示原始的乱序权值数据中的最小值，max(A)表示原始的乱序权值数据中的最大值。

[0058] 在本实施例中，采用构建乱序权值矩阵和均匀随机变量矩阵的方式，相比现有的随机投点的方式，本方法能够减少运算时间，提高数据处理效率，进而提高建模的效率。

[0059] 步骤S103：根据所述乱序权值矩阵、所述均匀随机变量矩阵和预设的裂缝片数量生成裂缝片的中心点。

[0060] 作为本实施例的进一步优化，根据所述乱序权值矩阵、所述均匀随机变量矩阵和预设的裂缝片数量生成裂缝片的中心点，包括：

[0061] 步骤B1：判断乱序权值矩阵与均匀随机变量矩阵中相同位置的数值点大小；

[0062] 步骤B1：确定为乱序权值矩阵的数值点大于均匀随机变量矩阵的数值点，则保留乱序权值矩阵的数值点作为裂缝种子点，根据裂缝种子点生成裂缝种子点矩阵；

[0063] 步骤B1：根据预设的裂缝片数量从裂缝种子点矩阵提取对应数量裂缝种子点作为中心点，即可以确定裂缝片的中心位置。

[0064] 在本实施例中，由于地震属性数据量巨大，批量生成的裂缝种子点数一般情况满足预设裂缝数量。但当预设裂缝片数量过大时，则会重新建立乱序权值数据以及随机变量数据，再次进行判断并提取相应裂缝种子点，直至满足预设裂缝片数量。

[0065] 步骤S104：根据获取的裂缝片的走向和倾角、以及所述中心点构建旋转角和旋转轴；

[0066] 作为本实施例中的进一步优化，根据获取的裂缝片的走向和倾角、以及所述中心点构建旋转角和旋转轴，包括：

[0067] 步骤C1：根据中心点确定旋转轴的中心，并根据走向和倾角确定旋转轴的坐标位置，在本实施例中，所述走向和倾角根据地质露头和/或测井裂缝成像获得。

[0068] 步骤C2：根据坐标位置和中心构建旋转轴；

[0069] 在本实施例中，旋转轴 的坐标位置计算方法为：

[0070] ；

[0071] ；

[0072] ；

[0073] 上式中，Inclination为裂缝片的倾角，Trend为裂缝片的走向。

[0074] 步骤C3：根据旋转轴的坐标位置构建旋转角。

[0075] 在本实施例中，旋转角 的计算公式为：

[0076] 。

[0077] 步骤S105：根据旋转角和旋转轴构建旋转矩阵，所述旋转矩阵用于表征裂缝片模型。

[0078] 在本实施例中，由于裂缝面的法向量N具有唯一性，同时裂缝片可以通过法向量和过裂缝片的一个点进行确定；例如采用下式来征裂缝片所在面：

[0079] ；其中， ；

[0080] 则裂缝片所在面的单位法向量N为 ，任意点 到裂缝片的距离为。

[0081] 式中，a表示单位法向量的x轴坐标；b表示单位法向量的y轴坐标；c表示单位法向量的z轴坐标；x表示平面某点的x轴坐标；y表示平面某点的y轴坐标；z表示平面某点的z轴坐标。

[0082] 由于裂缝片所在面的法向量N具有唯一性，因而可以通过构造目标法向量N来产生目标裂缝片所在的面。假定法向量N由Z轴旋转而来，此时需要确定旋转轴U和旋转角 ，即Z轴单位向量e(0,0,1)绕旋转轴 ，旋转 ，得到目标法向量N。

[0083] 旋转轴 与单位向量e(0,0,1)和目标法向量N垂直，因而 ，单位向量e与目标法向量N的夹角为旋转角 ， 。

[0084] 因此，作为本实例的进一步优化，根据旋转角和旋转轴构建旋转矩阵，包括：

[0085] 步骤D1：根据旋转角和旋转轴构建四元数；

[0086] 即四元数q为：

[0087] 其中 表示裂缝片定义的旋转角， 表示裂缝片定义的旋转轴；

[0088] 步骤D2：根据四元数构建旋转矩阵；

[0089] 即旋转矩阵R为： 。

[0090] 因此，可构建符合目标倾角和走向的目标旋转矩阵R。为方便输入，也可以直接输入目标单位向量N的各个分量。采用以上方法可以获得准确获得满足目标倾角、走向的旋转矩阵，进而大量产生目标倾角和走向，以及目标迹长的裂缝片。

[0091] 图2是本发明一种实施方式提供的基于乱序矩阵的半确定性裂缝建模系统的框图，如图2所示，一种基于乱序矩阵的半确定性裂缝建模系统，所述系统用于实现上述的基于乱序矩阵的半确定性裂缝建模方法，所述系统包括：

[0092] 提取模块，用于提取地震数据中裂缝敏感性高的属性数据；

[0093] 第一构建模块，用于根据所述属性数据构建乱序权值矩阵和均匀随机变量矩阵；

[0094] 生成模块，用于根据所述乱序权值矩阵、所述均匀随机变量矩阵和预设的裂缝片数量生成裂缝片的中心点；

[0095] 第二构建模块，用于根据获取的裂缝片的走向和倾角、以及所述中心点构建旋转角和旋转轴；

[0096] 第三构建模块，用于根据旋转角和旋转轴构建旋转矩阵，所述旋转矩阵用于表征裂缝片模型。

[0097] 本发明通过乱序权值矩阵和均匀随机变量矩阵，能够减少循环次数，进而能有效降低运算时间，提高建模的效率。其次，本发明的旋转矩阵能够实现对裂缝片产状进行准确描述，使所建立的裂缝模型更加接近真实情况，为油气藏的开采提供信息。

[0098] 本发明还提供一种电子设备，包括存储器、处理器以及存储在所述存储中并可在所述处理器上运行的计算机程序，所述处理器执行该计算机程序时实现上述的基于乱序矩阵的半确定性裂缝建模方法。

[0099] 本发明还提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被所述处理器执行时实现上述的基于乱序矩阵的半确定性裂缝建模方法。

[0100] 实施例二：

[0101] 在本实施例中，分别100、1000、1000、50000数量的裂缝模型，裂缝模型图像如图3‑6所示。本发明的裂缝模型能满足不同倾角走向的模型，相比较于欧拉角方法更容易理解且计算效率高；更重要的是四元数旋转矩阵可以避免欧拉角出现万向节死锁的现象。

[0102] 图7为现有建模方法与本发明建模方法的建模时间对比图，从图中可以看出，本发明的建模方法相比现有方法提高80倍，因为现有方法中应用了投点法，不断循环迭代会消耗大量时长，本发明所采用的矩阵运算法则通过提前建立随机矩阵，减少循环次数从而实现高效裂缝建模。

[0103] 需要说明的是，对于前述的各个方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本申请并不受所描述的动作顺序的限制，因为依据本申请，某一些步骤可以采用其他顺序或者同时进行。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于优选实施例，所涉及的动作和单元并不一定是本申请所必须的。