一种沉积期微幅度古隆起的推算方法及装置的制造方法

一种沉积期微幅度古隆起的推算方法及装置的制造方法
【专利摘要】本发明提供了一种沉积期微幅度古隆起的推算方法及装置，涉及地质勘探技术领域。方法包括：根据现今构造图数据确定现今构造趋势面图数据，进而确定各井间网格点的拟构造幅度；拟构造幅度为现今构造图上的点相对于现今构造趋势面的竖直方向幅度；根据现今构造图的剥蚀线数据确定待推算地区地层的沉积期构造趋势面图数据；根据沉积期构造趋势面图数据和拟构造幅度确定沉积期构造的相对高低数据，并根据沉积期构造的相对高低数据生成沉积期微幅度古隆起图数据，进而确定沉积微相图数据。本发明可以解决现有技术仅能推算隆起幅度大约在50米及其以上的古隆起，对于隆起幅度仅有20米甚至更小的沉积期微幅度古隆起推算精度较差的问题。
【专利说明】
一种沉积期微幅度古隆起的推算方法及装置
技术领域
[0001] 本发明涉及地质勘探技术领域，尤其涉及一种沉积期微幅度古隆起的推算方法及 装置。
【背景技术】
[0002] 当前，在世界范围内，与古地貌直接相关的油气储量大于5 X 108bbl油当量的大油 气田众多。古地貌是指地质历史时期中由构造变形、沉积充填、差异压实和风化剥蚀等因素 形成并残留至今或埋在地下的某些地貌类型。与古地貌有关的油气圈闭分成两大类:一类 为与沉积期古地貌有关的油气圈闭，如:礁体、颗粒滩、砂坝、砂丘、侵蚀河谷等;另一类为与 剥蚀期古地貌有关的油气圈闭，如：因古地貌或表生淋滤作用产生的溶洞储层，潜山储层 等。而当前的古地貌反演方法也主要分为剥蚀古地貌和沉积古地貌反演两大方面。其中，沉 积古地貌反演方法主要有印模法、沉积学法、层序地层学法和层拉平法。印模法优点是简单 易操作，而缺点是上覆地层参考层不易确定，压实校正不易确定，只能半定量表征沉积期古 地貌。沉积学法优点是理论基础充实，缺点是仅能定性展示沉积期古地貌的时空格局。层序 地层学和层拉平方法优点是借助层序地层学理论和物探技术实现古地貌的刻画，而缺点是 要求地震资料的品质要高，同时物探纵向分辨率低，对于物探资料品质差的深层较薄地层 中的沉积期微幅度起伏刻画不清。盆地模拟法技术用于古地貌恢复的优点是考虑因素全， 资料类型多，假设条件较为合理，但其缺点是参数太多，受多种因素控制，很难得到可靠的 参数。
[0003] 沉积期水下沉积期微幅度古隆起是控制潮坪相碳酸盐岩颗粒滩分布的关键因素。 颗粒滩又是后期白云化并发育为优质储层最有利的微相。因此，刻画沉积期水下古隆起是 寻找潮坪相优质白云岩储层的先决条件。现有技术一运用层序地层学与印模法相结合的方 法，根据层序地层学理论选取目标层，通过残余厚度求取、真厚度校正、压实恢复、潜山剥蚀 区确定、古水深恢复等步骤，恢复曙北地区沙四段SQ2主要滩坝沉积时期的古地貌。根据恢 复古地貌的特征，将研究区分为沿岸带、潜山带、北部斜坡带、南部斜坡带和洼陷带5个地貌 单元。滩坝砂体主要分布在沿岸带和南、北部斜坡带，以南部斜坡带最为发育。古地貌通过 沉积地貌形态及其所控制的水体分布与相对水深、水动力、物源搬运通道和卸载场所等因 素控制滩坝砂体的发育。研究区砂质滩坝主要发育于SQ2低位体系域沉积时期，至湖侵体系 域沉积时期水体加深，在潜山带发育碳酸盐岩滩坝。现有技术二利用盆地模拟反演古地貌 技术和依据2D和3D地震资料解释及钻井分层数据得到的不同时代地层的残余厚度图和现 今埋深图，在对主要构造不整合面剥蚀量恢复和去压实校正的基础上，进行古构造、古地貌 恢复;计算了不同地层的在不同构造时期的古埋深和古厚度，系统深入研究了巴楚-麦盖提 地区古构造演化历史，阐明了巴楚隆起与和田古隆起构造演化关系及迀移规律，结合奥陶 系岩溶储层发育特征，探讨了古隆起演化与岩溶储层发育之间关系，建立了岩溶储层发育 的构造演化模式，推算了奥陶系有利岩溶储层分布地区。
[0004] 发明人在研究本发明以及上述现有技术时发现，现有技术一和现有技术二均是为 了推算隆起幅度大约在50米及其以上的古隆起。对于隆起幅度仅有20米甚至更小的沉积期 微幅度古隆起，现有技术的推算精度较差。

【发明内容】

[0005] 本发明的实施例提供一种沉积期微幅度古隆起的推算方法及装置，以解决现有技 术仅能推算隆起幅度大约在50米及其以上的古隆起，对于隆起幅度仅有20米甚至更小的沉 积期微幅度古隆起现有技术的推算精度较差的问题。
[0006] 为达到上述目的，本发明采用如下技术方案：
[0007] -种沉积期微幅度古隆起的推算方法，包括：
[0008] 获取待推算地区地层的现今构造图数据；
[0009] 根据所述现今构造图数据确定现今构造趋势面图数据；
[0010] 根据所述现今构造图数据和所述现今构造趋势面图数据确定各井间网格点的拟 构造幅度;所述拟构造幅度为现今构造图上的点相对于现今构造趋势面的竖直方向幅度；
[0011] 获取现今构造图的剥蚀线数据，并根据所述现今构造图的剥蚀线数据确定待推算 地区地层的沉积期构造趋势面图数据；
[0012] 根据所述沉积期构造趋势面图数据和所述拟构造幅度确定沉积期构造的相对高 低数据，并根据所述沉积期构造的相对高低数据生成沉积期微幅度古隆起图数据；
[0013] 根据所述沉积期微幅度古隆起图数据确定沉积微相图数据。
[0014] 具体的，根据所述现今构造图数据确定现今构造趋势面图数据，包括：
[0015] 在现今构造图数据对应的现今构造图中获取第一数据点和第二数据点；其中，所 述第一数据点的坐标为(1〇,7〇,2())，所述第二数据点的坐标为(&，7 11，21〇;所述第一数据点 为现今构造图中的最低点；所述第二数据点为现今构造图中的以最高点为中心，以预设距 离为半径的范围内的最低点；
[0016] 在第一数据点和第二数据点连接的直线上每隔一预设高差Δ z确定一第三数据 点;所述第三数据点的坐标为(Xk，yk，Zk)

[0018]确定空间内各点坐标(x，y，z)，并根据克里格插值算法对空间内各点坐标进行插 值，生成现今构造趋势面图数据；
1上的点的高程值均 为Zko
[0019] 具体的，根据所述现今构造图数据和所述现今构造趋势面图数据确定各井间网格 点的拟构造幅度，包括：
[0020] 将现今构造图数据中的各井间网格点的坐标与现今构造趋势面图数据中的各点 坐标在竖直方向上 对应；
[0021] 确定现今构造图数据中的各井间网格点与其对应的现今构造趋势面图数据中的 点在竖直方向上的幅度差，并根据各幅度差确定各井间网格点的拟构造幅度。
[0022] 具体的，获取现今构造图的剥蚀线数据，并根据所述现今构造图的剥蚀线数据确 定待推算地区地层的沉积期构造趋势面图数据，包括：
[0023] 在现今构造图的剥蚀线数据对应的剥蚀线上选取任意两点，并以其中一点为坐标 原点建立方位坐标系;所述方位坐标系的正北方向为〇°，正东方向为90° ;其中，在所述方位 坐标系中原点坐标为(1'0,7'0)，所述任意两点中的另一点坐标为(1'1，7'1);
[0024] 将从(x'o，y'o)指向（x'i，y'i)的向量的方向作为剥蚀线走向和沉积期构造趋势面 走向；将从(X'Q，y'Q)指向（1'1，7'1)的向量的垂直方向作为沉积期构造趋势面倾向；
[0026] 确定沉积期构造趋势面上各点坐标(X'，y'，z'）；其中，（1'：1，7'：1，2'：1)为过沉积期 构造趋势面上的某点坐标；（x'hy'hZ'J为沿沉积期构造趋势面倾向预设距离为Γ的某点 的坐标;α为预先设置的沉积期构造趋势面倾角；Θ为剥蚀线走向与方位坐标系的正北方向
，且沿沉积期构造趋势面走向方向的 直线方程;zx",y" = z ' j为过(x'j，y'j，z'j)点且沿沉积期构造趋势面走向方向的直线的高程 值；
[0027] 根据沉积期构造趋势面上各点坐标(X'，y'，z'）构成沉积期构造趋势面图数据。
[0028] 具体的，根据所述沉积期构造趋势面图数据和所述拟构造幅度确定沉积期构造的 相对高低数据，包括：
[0029] 将各网格点的沉积期构造趋势面图数据与其拟构造幅度相加，确定沉积期构造的 相对高低数据。
[0030] -种沉积期微幅度古隆起的推算装置，包括：
[0031] 现今构造图数据获取单元，用于获取待推算地区地层的现今构造图数据；
[0032] 现今构造趋势面图数据确定单元，用于根据所述现今构造图数据确定现今构造趋 势面图数据；
[0033] 拟构造幅度确定单元，用于根据所述现今构造图数据和所述现今构造趋势面图数 据确定各井间网格点的拟构造幅度;所述拟构造幅度为现今构造图上的点相对于现今构造 趋势面的竖直方向幅度；
[0034] 沉积期构造趋势面图数据确定单元，用于获取现今构造图的剥蚀线数据，并根据 所述现今构造图的剥蚀线数据确定待推算地区地层的沉积期构造趋势面图数据；
[0035] 沉积期微幅度古隆起图数据生成单元，用于根据所述沉积期构造趋势面图数据和 所述拟构造幅度确定沉积期构造的相对高低数据，并根据所述沉积期构造的相对高低数据 生成沉积期微幅度古隆起图数据；
[0036]沉积微相图数据确定单元，用于根据所述沉积期微幅度古隆起图数据确定沉积微 相图数据。
[0037] 此外，所述现今构造趋势面图数据确定单元，具体用于：
[0038] 在现今构造图数据对应的现今构造图中获取第一数据点和第二数据点；其中，所 述第一数据点的坐标为(1〇,7〇,2())，所述第二数据点的坐标为(&，7 11，21〇;所述第一数据点 为现今构造图中的最低点；所述第二数据点为现今构造图中的以最高点为中心，以预设距 离为半径的范围内的最低点；
[0039]在第一数据点和第二数据点连接的直线上每隔一预设高差Δ Z确定一第三数据 点;所述第三数据点的坐标为(xk，yk，zk):

[0041] 确定空间内各点坐标(x，y，z)，并根据克里格插值算法对空间内各点坐标进行插 值，生成现今构造趋势面图数据;：上的点的高程值均 为Zko
[0042] 此外，所述拟构造幅度确定单元，具体用于：
[0043] 将现今构造图数据中的各井间网格点的坐标与现今构造趋势面图数据中的各点 坐标在竖直方向上 对应；
[0044]确定现今构造图数据中的各井间网格点与其对应的现今构造趋势面图数据中的 点在竖直方向上的幅度差，并根据各幅度差确定各井间网格点的拟构造幅度。
[0045] 此外，沉积期构造趋势面图数据确定单元，具体用于：
[0046] 在现今构造图的剥蚀线数据对应的剥蚀线上选取任意两点，并以其中一点为坐标 原点建立方位坐标系;所述方位坐标系的正北方向为〇°，正东方向为90° ;其中，在所述方位 坐标系中原点坐标为(1'0,7'0)，所述任意两点中的另一点坐标为(1'1，7'1);
[0047] 将从(x'o，y'o)指向（x'i，y'i)的向量的方向作为剥蚀线走向和沉积期构造趋势面 走向；将从(X'Q，y'Q)指向（1'1，7'1)的向量的垂直方向作为沉积期构造趋势面倾向；
[0049] 确定沉积期构造趋势面上各点坐标(X'，y'，z'）；其中，（1'：1，7'：1，2'：1)为过沉积期 构造趋势面上的某点坐标；（x'hy'hZ'J为沿沉积期构造趋势面倾向预设距离为Γ的某点 的坐标;α为预先设置的沉积期构造趋势面倾角；Θ为剥蚀线走向与方位坐标系的正北方向
为过(X、，y ' j，Z '」）点，且沿沉积期构造趋势面走向方向的 直线方程;zx",y" = z ' j为过(x'j，y'j，z'j)点且沿沉积期构造趋势面走向方向的直线的高程 值；
[0050] 根据沉积期构造趋势面上各点坐标(X'，y'，z'）构成沉积期构造趋势面图数据。
[0051] 此外，所述沉积期微幅度古隆起图数据生成单元，具体用于：
[0052]将各网格点的沉积期构造趋势面图数据与其拟构造幅度相加，确定沉积期构造的 相对高低数据。
[0053]本发明实施例提供的一种沉积期微幅度古隆起的推算方法及装置，其所针对的环 境为弱构造变形区中的深埋藏潮坪相碳酸盐岩的沉积期。本发明实施例根据现今构造图数 据确定现今构造趋势面图数据;进而根据现今构造图数据和现今构造趋势面图数据确定各 井间网格点的拟构造幅度;其中，该拟构造幅度为现今构造图上的点相对于现今构造趋势 面的竖直方向幅度;之后获取现今构造图的剥蚀线数据，并根据现今构造图的剥蚀线数据 确定待推算地区地层的沉积期构造趋势面图数据;根据沉积期构造趋势面图数据和拟构造 幅度确定沉积期构造的相对高低数据，并根据沉积期构造的相对高低数据生成沉积期微幅 度古隆起图数据;根据沉积期微幅度古隆起图数据确定沉积微相图数据。通过本发明可以 对隆起幅度仅有20米甚至更小的沉积期微幅度古隆起进行推算，所形成的沉积期微幅度古 隆起图精度更高。可以解决现有技术仅能推算隆起幅度大约在50米及其以上的古隆起，对 于隆起幅度仅有20米甚至更小的沉积期微幅度古隆起，现有技术的推算精度较差的问题。
【附图说明】
[0054]为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可 以根据这些附图获得其他的附图。
[0055] 图1为本发明实施例提供一种沉积期微幅度古隆起的推算方法的流程图一；
[0056] 图2为本发明实施例提供一种沉积期微幅度古隆起的推算方法的流程图二；
[0057]图3为本发明实施例中苏东区块奥陶纪古沉积期地势格局与苏东区块三叠纪构造 反转后的地势格局的示意图；
[0058]图4为本发明实施例中海相碳酸盐岩沉积环境的示意图；
[0059] 图5为本发明实施例中海相碳酸盐岩沉积环境简化示意图；
[0060] 图6为本发明实施例中现今构造图的剥蚀线示意图；
[0061] 图7为本发明实施例中的一方位坐标系以及剥蚀线示意图；
[0062] 图8为本发明实施例中一应用实例区域的示意图；
[0063] 图9为本发明实施例中的无阻流量和颗粒滩储层厚度的关系示意图；
[0064] 图10为本发明实施例提供的一种沉积期微幅度古隆起的推算装置的结构示意图。
【具体实施方式】
[0065] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0066] 在实现本发明实施例时，发明人发现，现有的古地形恢复方法主要是针对动辄几 万平方公里的范围二级构造单元或地下地形差异幅度非常大的三级构造单元。对于面积只 有几百平方公里、仅在单斜构造上发育幅度20米甚至更小鼻状构造、井距已达到1000米甚 至更小的气田来说，仍然需要了解优质储层分布的非均质性，以期提高加密井的成功率。而 上述古地貌反演方法并不适用面积和构造幅度如此之小的气田。其原因是，工区目的层埋 深在3000米左右，地震波传播距离远、中间反射层多，导致反映目的层的信号弱，分辨率降 低;地表为疏松沙漠地带，引起地震采集过程中高频信号严重衰减，导致地震分辨率降低； 目的层之上又有煤层引起地震波传播速度低、形成强振幅、高连续的地震反射、导致其下部 地层信号弱、分辨率降低。地表疏松、埋藏深、煤层干扰导致目的层地震资料品质很差，无法 有效反映厚度仅有50米左右的地层层序关系。因此，基于地震资料的古地貌反演方法无法 应用。
[0067] 如图1所示，本发明实施例提供一种沉积期微幅度古隆起的推算方法，包括：
[0068] 步骤101、获取待推算地区地层的现今构造图数据。
[0069] 步骤102、根据现今构造图数据确定现今构造趋势面图数据。
[0070] 步骤103、根据现今构造图数据和现今构造趋势面图数据确定各井间网格点的拟 构造幅度。
[0071] 其中，拟构造幅度为现今构造图上的点相对于现今构造趋势面的竖直方向幅度。
[0072] 步骤104、获取现今构造图的剥蚀线数据，并根据现今构造图的剥蚀线数据确定待 推算地区地层的沉积期构造趋势面图数据。
[0073] 步骤105、根据沉积期构造趋势面图数据和拟构造幅度确定沉积期构造的相对高 低数据，并根据沉积期构造的相对高低数据生成沉积期微幅度古隆起图数据。
[0074]步骤106、根据沉积期微幅度古隆起图数据确定沉积微相图数据。
[0075]本发明实施例提供的一种沉积期微幅度古隆起的推算方法，其所针对的环境为弱 构造变形区中的深埋藏潮坪相碳酸盐岩的沉积期。本发明实施例根据现今构造图数据确定 现今构造趋势面图数据;进而根据现今构造图数据和现今构造趋势面图数据确定各井间网 格点的拟构造幅度;其中，该拟构造幅度为现今构造图上的点相对于现今构造趋势面的竖 直方向幅度;之后获取现今构造图的剥蚀线数据，并根据现今构造图的剥蚀线数据确定待 推算地区地层的沉积期构造趋势面图数据;根据沉积期构造趋势面图数据和拟构造幅度确 定沉积期构造的相对高低数据，并根据沉积期构造的相对高低数据生成沉积期微幅度古隆 起图数据;根据沉积期微幅度古隆起图数据确定沉积微相图数据。通过本发明可以对隆起 幅度仅有20米甚至更小的沉积期微幅度古隆起进行推算，所形成的沉积期微幅度古隆起图 精度更高。可以解决现有技术仅能推算隆起幅度大约在50米及其以上的古隆起，对于隆起 幅度仅有20米甚至更小的沉积期微幅度古隆起，现有技术的推算精度较差的问题。
[0076] 为了使本领域的技术人员更好的了解本发明，下面列举一个更为详细的实施例， 如图2所示，本发明实施例提供一种沉积期微幅度古隆起的推算方法，其所针对的环境为弱 构造变形区中的深埋藏潮坪相碳酸盐岩的沉积期，其中弱构造变形区是指地层未产生明显 的弯曲、褶皱或断裂的区域。为了实现如图2所示的沉积期微幅度古隆起的推算方法，此处 本发明存在如下假设条件：如图3所示，在弱构造变形区，构造趋势面的变动并没有改变局 部点的拟构造幅度，如苏东区块奥陶纪古沉积期地势格局在经过构造反转后变为苏东区块 三叠纪构造反转后的地势格局。或者构造趋势面的变动对于局部井点的拟构造幅度改动程 度相同。也就是说，现今的各点拟构造幅度等于或正比于沉积期各点的拟构造幅度。
[0077] 该沉积期微幅度古隆起的推算方法包括：
[0078]步骤201、获取待推算地区地层的现今构造图数据。
[0079]此处的现今构造图的形成技术已经较为成熟，例如：利用标志层法、沉积旋回法、 等厚法、特征曲线法等多种方法，先选择地层的标志层，再建立标准对比剖面和建立地层对 比网络，最终确定气田范围内所有钻井合理的地层顶界和底界的深度等步骤进行地层对 比，并换算成海拔深度。利用克里格插值算法推算井间网格的海拔深度，将相同海拔连线得 到构造等高线，从而形成现今构造图。
[0080] 步骤202、在现今构造图数据对应的现今构造图中获取第一数据点和第二数据点。
[0081] 其中，第一数据点的坐标为(xo，yo，Z())，第二数据点的坐标为(xn，yn，z n)。
[0082] 第一数据点为现今构造图中的最低点;所述第二数据点为现今构造图中的以最高 点为中心，以预设距离(例如5000米)为半径的范围内的最低点；为了使最终得到的拟构造 幅度为正值，ZQ-般要小于现今构造图中的最小等高线值，Zn-般要小于现今构造图中的最 高点局部区域的最小值。
[0083]步骤203、在第一数据点和第二数据点连接的直线上每隔一预设高差Δ Z确定一第 三数据点。
[0084]其中，该第三数据点的坐标为(xk，yk，zk)。
[0085]
，表示在第一数据点和第二数据点之间所划分的区域份数。
[0087] 确定空间内各点坐标(x，y，z)，并根据克里格插值算法对空间内各点坐标进行插 值，生成现今构造趋势面图数据
上的点的高程值均 为Zko
[0088] 其中，kE [ 1，n_l ]。
[0089]
表示过第一数据点和第二数据点的直线方程。
[0090]
第三数据点， 的直线方程。
[0091] 步骤205、将现今构造图数据中的各井间网格点的坐标与现今构造趋势面图数据 中的各点坐标在竖直方向上 对应。
[0092]步骤206、确定现今构造图数据中的各井间网格点与其对应的现今构造趋势面图 数据中的点在竖直方向上的幅度差，并根据各幅度差确定各井间网格点的拟构造幅度。
[0093] 值得说明的是，如图4所示，海相碳酸盐岩沉积环境总体是由陆地到海洋的斜坡， 在此斜坡之上发育局部的隆起，这些隆起或沉没于水下或露出水面。这些隆起的顶面相对 于斜坡的趋势面具有一定的正向幅度。
[0094] 海相碳酸盐岩的这种沉积地貌可以简化如图5所示，地面的起伏线如图5中的曲线 所示，斜坡的趋势面如图5中的直线所示。在本发明实施例中，拟构造幅度为现今构造图上 的点相对于现今构造趋势面的竖直方向幅度，即图5中的曲线上的点相对于直线上的点的 竖直方向距离301。这区别于常规的现有技术中的构造幅度302，即以图5中的曲线上的点向 直线做垂线段，该垂线段的长度即为构造幅度302。
[0095]步骤207、在现今构造图的剥蚀线数据对应的剥蚀线上选取任意两点，并以其中一 点为坐标原点建立方位坐标系。
[0096]所述方位坐标系的正北方向为0°，正东方向为90° ;其中，在所述方位坐标系中原 点坐标为(X ' 〇，y ' 〇)，所述任意两点中的另一点坐标为(X ' 1，y ' 1)。
[0097]值得说明的是，如图6所示，现今构造图中的剥蚀线401是指在剥蚀线401的一侧地 层厚度大于0，而在其另一侧地层厚度为0。
[0098]另外，该方位坐标系可以如图7所示，其正北方向为0°，正东方向为90°，则(x'〇，y
[0099]步骤208、将从(x'o，y'o)指向（x'i，y'i)的向量的方向作为剥蚀线走向和沉积期构 造趋势面走向；将从(x'o，y'o)指向（1'1，7'1)的向量的垂直方向作为沉积期构造趋势面倾 向。
[0101] 确定沉积期构造趋势面上各点坐标(X'，y'，Z'），并根据沉积期构造趋势面上各点 坐标(X'，y'，Z'）构成沉积期构造趋势面图数据。
[0102] 其中，（1'：1，7'：1，2'：1)为过沉积期构造趋势面上的某点坐标；（1'」，7'」，2'」)为沿沉 积期构造趋势面倾向预设距离为Γ的某点的坐标;α为预先设置的沉积期构造趋势面倾角； θ为剥蚀线走向与方位坐标系的正北方向的夹角；
为过(x'j，y'j，z 'J)点，且沿沉积期构造趋势面走向方向的直线方程;2^/=2'伪过(4,7、，4)点且沿 沉积期构造趋势面走向方向的直线的高程值。
[0103] 步骤210、将各网格点的沉积期构造趋势面图数据与其拟构造幅度相加，确定沉积 期构造的相对高低数据，并根据沉积期构造的相对高低数据生成沉积期微幅度古隆起图数 据。
[0104] 步骤211、根据沉积期微幅度古隆起图数据确定沉积微相图数据。
[0105] 例如沉积期微幅度古隆起图数据对应的古隆起图呈西北方向高，东南方向低的总 体格局。若在钻井时取出的岩石样品，也称为岩芯的岩石颜色表现为灰黑色则其为水下环 境，可划归潮下带;若岩石类型为颗粒灰岩或者具有残余颗粒结构的白云岩，则其沉积微相 为颗粒滩;若白云岩的白云化类型为蒸发栗，则其为水上环境，可划归到潮上带。通过此种 方式可以得到准确的沉积微相图。
[0106] 为了使本领域的技术人员更好的了解本发明的技术效果，下面列举一个本发明的 应用实例：
[0107] 以沉积学理论为指导、利用测井地质学方法追踪各井马五5顶面和底面深度数据， 并结合井口补心海拔将钻井深度转化为海拔深度，在考虑区域构造格局情况下，利用克里 格方法计算各点的海拔深度，并形成现今构造图（图8的A部分）。研究区的井距以达到800米 左右，利用井资料得到的构造图的精度已能满足刻画微幅度鼻状构造。从现今构造面特征 出发，计算地层的平均倾角和倾向，利用平均倾角和倾向，计算现今构造趋势面（图8的B部 分）。利用现今构造减去现今构造趋势，即得到现今的拟构造幅度（图8的C部分）。研究区在 马五段沉积后即整体抬升剥蚀，马五5地层剥蚀线的走向即为古构造的平均走向，则古构造 倾向就与剥蚀线平均走向垂直。在确定了古构造的倾向和倾角之后，计算沉积期构造趋势 面(图8的D部分）。计算沉积期构造趋势面，现今潮坪相沉积地貌的小于3°的地形坡度，采用 小于3°的某个角度作为沉积期构造趋势面倾角。
[0108] 基于本发明的假设，由沉积期构造趋势面变化到现今构造趋势面并没有改变局部 点的拟构造幅度，或者构造趋势面的变动对于局部点的拟构造幅度改变程度相同。也就是 说，现今的各点拟构造幅度等于或正比于沉积期各点的拟构造幅度。利用现今构造得到拟 构造幅度加上沉积期构造趋势面，即可得到沉积期古地形（图8的E部分）。结合反演得到古 地形、取芯井岩石学，白云化类型等资料，可以进一步反演弱变形构造区沉积微相图（图8的 F部分）。研究区共有266 口井钻遇马五5地层，微相类型有潮上带云坪、颗粒滩、灰坪、灰泥丘 四种，该工区内物性和储层发育概率的表现为下表1:
[0109] 表1:
[0110]
[0111] 从表1中可以看出，颗粒滩物性比其他微相物性要好;研究区共有266 口井钻遇马 五5地层，其中76 口井钻遇颗粒滩微相，颗粒滩钻遇率仅为28.6 %，但在钻遇颗粒滩的76 口 井中，有68 口井的颗粒滩演化为储层，钻遇颗粒滩后找到储层的概率为89.5%。
[0112] 通过由颗粒滩演化而来的白云岩储层厚度(简称颗粒滩储层厚度)与马五5单层测 试高产量井的无阻流量关系图（图9)来看，颗粒滩储层厚度与无阻流量呈非常强烈的正相 关关系，其相关系数可达到:R = 〇.95。因此，在马五5地层中寻找高产量气井的关键是找到 颗粒滩白云岩储层，颗粒滩储层又赋存于颗粒滩微相中。由油气勘探开发的实践知道，颗粒 滩微相主要存在于水下低隆起顶部水动力能量较强的部位。所以通过本发明的方法，寻找 到沉积期水下低隆起分布范围，则可以寻找到颗粒滩的分布范围，进而寻找到高产量气井 分布地点。以本实施例结果为例，找到颗粒滩并有高产量气井的概率为89.5% X 95 % = 85%〇
[0113] 本发明实施例提供的一种沉积期微幅度古隆起的推算方法，其所针对的环境为弱 构造变形区中的深埋藏潮坪相碳酸盐岩的沉积期。本发明实施例根据现今构造图数据确定 现今构造趋势面图数据;进而根据现今构造图数据和现今构造趋势面图数据确定各井间网 格点的拟构造幅度;其中，该拟构造幅度为现今构造图上的点相对于现今构造趋势面的竖 直方向幅度;之后获取现今构造图的剥蚀线数据，并根据现今构造图的剥蚀线数据确定待 推算地区地层的沉积期构造趋势面图数据;根据沉积期构造趋势面图数据和拟构造幅度确 定沉积期构造的相对高低数据，并根据沉积期构造的相对高低数据生成沉积期微幅度古隆 起图数据;根据沉积期微幅度古隆起图数据确定沉积微相图数据。通过本发明可以对隆起 幅度仅有20米甚至更小的沉积期微幅度古隆起进行推算，所形成的沉积期微幅度古隆起图 精度更高。可以解决现有技术仅能推算隆起幅度大约在50米及其以上的古隆起，对于隆起 幅度仅有20米甚至更小的沉积期微幅度古隆起，现有技术的推算精度较差的问题。
[0114] 对应于上述图1和图2所述的方法实施例，如图10所示，本发明实施例还提供一种 沉积期微幅度古隆起的推算装置，包括：
[0115] 现今构造图数据获取单元51，用于获取待推算地区地层的现今构造图数据。
[0116] 现今构造趋势面图数据确定单元52,用于根据所述现今构造图数据确定现今构造 趋势面图数据。
[0117] 拟构造幅度确定单元53,用于根据所述现今构造图数据和所述现今构造趋势面图 数据确定各井间网格点的拟构造幅度;所述拟构造幅度为现今构造图上的点相对于现今构 造趋势面的竖直方向幅度。
[0118]沉积期构造趋势面图数据确定单元54，用于获取现今构造图的剥蚀线数据，并根 据所述现今构造图的剥蚀线数据确定待推算地区地层的沉积期构造趋势面图数据。
[0119] 沉积期微幅度古隆起图数据生成单元55,用于根据所述沉积期构造趋势面图数据 和所述拟构造幅度确定沉积期构造的相对高低数据，并根据所述沉积期构造的相对高低数 据生成沉积期微幅度古隆起图数据。
[0120] 沉积微相图数据确定单元56,用于根据所述沉积期微幅度古隆起图数据确定沉积 微相图数据。
[0121] 此外，所述现今构造趋势面图数据确定单元52,具体可以：
[0122] 在现今构造图数据对应的现今构造图中获取第一数据点和第二数据点；其中，所 述第一数据点的坐标为(1〇,7〇,2())，所述第二数据点的坐标为(&，7 11，21〇;所述第一数据点 为现今构造图中的最低点；所述第二数据点为现今构造图中的以最高点为中心，以预设距 离为半径的范围内的最低点。
[0123] 在第一数据点和第二数据点连接的直线上每隔一预设高差Δ Z确定一第三数据 点;所述第三数据点的坐标为(xk，yk，zk);

[0125] 确定空间内各点坐标(x，y，z)，并根据克里格插值算法对空间内各点坐标进行插 值，生成现今构造趋势面图数据：
上的点的高程值均 为Zko
[0126] 此外，该拟构造幅度确定单元53，具体可以：
[0127] 将现今构造图数据中的各井间网格点的坐标与现今构造趋势面图数据中的各点 坐标在竖直方向上 对应。
[0128] 确定现今构造图数据中的各井间网格点与其对应的现今构造趋势面图数据中的 点在竖直方向上的幅度差，并根据各幅度差确定各井间网格点的拟构造幅度。
[0129] 此外，该沉积期构造趋势面图数据确定单元54,具体可以：
[0130] 在现今构造图的剥蚀线数据对应的剥蚀线上选取任意两点，并以其中一点为坐标 原点建立方位坐标系;所述方位坐标系的正北方向为〇°，正东方向为90° ;其中，在所述方位 坐标系中原点坐标为(1'0,7'0)，所述任意两点中的另一点坐标为(1'1，7'1);
[0131] 将从(x'o，y'o)指向（x'i，y'i)的向量的方向作为剥蚀线走向和沉积期构造趋势面 走向；将从(X'Q，y'Q)指向（1'1，7'1)的向量的垂直方向作为沉积期构造趋势面倾向；
[0133]确定沉积期构造趋势面上各点坐标(X'，y'，z'）；其中，（乂'：1，7'：1，2'：1)为过沉积期 构造趋势面上的某点坐标；（x'hy'hZ'J为沿沉积期构造趋势面倾向预设距离为Γ的某点 的坐标;α为预先设置的沉积期构造趋势面倾角；Θ为剥蚀线走向与方位坐标系的正北方向
过(X ?，y ?，ζ、）点，且沿沉积期构造趋势面走向方向的 直线方程;zx",y" = z ' j为过(x'j，y'j，z'j)点且沿沉积期构造趋势面走向方向的直线的高程 值；
[0134] 根据沉积期构造趋势面上各点坐标(x'，y'，z'）构成沉积期构造趋势面图数据。
[0135] 此外，该沉积期微幅度古隆起图数据生成单元55,具体可以将各网格点的沉积期 构造趋势面图数据与其拟构造幅度相加，确定沉积期构造的相对高低数据。
[0136] 值得说明的是，本发明实施例提供的一种沉积期微幅度古隆起的推算装置的具体 实现方式可以参见上述的方法实施例，此处不再赘述。
[0137] 本发明实施例提供的一种沉积期微幅度古隆起的推算装置，其所针对的环境为弱 构造变形区中的深埋藏潮坪相碳酸盐岩的沉积期。本发明实施例根据现今构造图数据确定 现今构造趋势面图数据;进而根据现今构造图数据和现今构造趋势面图数据确定各井间网 格点的拟构造幅度;其中，该拟构造幅度为现今构造图上的点相对于现今构造趋势面的竖 直方向幅度;之后获取现今构造图的剥蚀线数据，并根据现今构造图的剥蚀线数据确定待 推算地区地层的沉积期构造趋势面图数据;根据沉积期构造趋势面图数据和拟构造幅度确 定沉积期构造的相对高低数据，并根据沉积期构造的相对高低数据生成沉积期微幅度古隆 起图数据;根据沉积期微幅度古隆起图数据确定沉积微相图数据。通过本发明可以对隆起 幅度仅有20米甚至更小的沉积期微幅度古隆起进行推算，所形成的沉积期微幅度古隆起图 精度更高。可以解决现有技术仅能推算隆起幅度大约在50米及其以上的古隆起，对于隆起 幅度仅有20米甚至更小的沉积期微幅度古隆起，现有技术的推算精度较差的问题。
[0138] 本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序 产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实 施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机 可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产 品的形式。
[0139] 本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程 图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流 程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序 指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产 生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实 现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。 [0140]这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特 定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指 令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或 多个方框中指定的功能。
[0141]这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计 算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或 其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一 个方框或多个方框中指定的功能的步骤。
[0142]本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例 的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员， 依据本发明的思想，在【具体实施方式】及应用范围上均会有改变之处，综上所述，本说明书内 容不应理解为对本发明的限制。
【主权项】
1. 一种沉积期微幅度古隆起的推算方法，其特征在于，包括： 获取待推算地区地层的现今构造图数据； 根据所述现今构造图数据确定现今构造趋势面图数据； 根据所述现今构造图数据和所述现今构造趋势面图数据确定各井间网格点的拟构造 幅度;所述拟构造幅度为现今构造图上的点相对于现今构造趋势面的竖直方向幅度； 获取现今构造图的剥蚀线数据，并根据所述现今构造图的剥蚀线数据确定待推算地区 地层的沉积期构造趋势面图数据； 根据所述沉积期构造趋势面图数据和所述拟构造幅度确定沉积期构造的相对高低数 据，并根据所述沉积期构造的相对高低数据生成沉积期微幅度古隆起图数据； 根据所述沉积期微幅度古隆起图数据确定沉积微相图数据。2. 根据权利要求1所述的沉积期微幅度古隆起的推算方法，其特征在于，根据所述现今 构造图数据确定现今构造趋势面图数据，包括： 在现今构造图数据对应的现今构造图中获取第一数据点和第二数据点；其中，所述第 一数据点的坐标为(10,70,20)，所述第二数据点的坐标为(&，711，^);所述第一数据点为现 今构造图中的最低点；所述第二数据点为现今构造图中的以最高点为中心，以预设距离为 半径的范围内的最低点； 在第一数据点和第二数据点连接的直线上每隔一预设高差A Z确定一第三数据点；所 述第三数据点的坐标为(xk，yk，zk);? -1 = ; 根据公式确定空间内各点坐标(X，y，z )，并根据克里格插值算法对空间内各点坐标进行插值，生 成现今构造趋势面图数据;在直线^ = + 上的点的高程值均为zk。 :凡―λ y" -y〇3. 根据权利要求2所述的沉积期微幅度古隆起的推算方法，其特征在于，根据所述现今 构造图数据和所述现今构造趋势面图数据确定各井间网格点的拟构造幅度，包括： 将现今构造图数据中的各井间网格点的坐标与现今构造趋势面图数据中的各点坐标 在竖直方向上一一对应； 确定现今构造图数据中的各井间网格点与其对应的现今构造趋势面图数据中的点在 竖直方向上的幅度差，并根据各幅度差确定各井间网格点的拟构造幅度。4. 根据权利要求3所述的沉积期微幅度古隆起的推算方法，其特征在于，获取现今构造 图的剥蚀线数据，并根据所述现今构造图的剥蚀线数据确定待推算地区地层的沉积期构造 趋势面图数据，包括： 在现今构造图的剥蚀线数据对应的剥蚀线上选取任意两点，并以其中一点为坐标原点 建立方位坐标系;所述方位坐标系的正北方向为0°，正东方向为90° ;其中，在所述方位坐标 系中原点坐标为(1'〇,7'〇)，所述任意两点中的另一点坐标为&'1，7'1); 将从(X'Q，y'Q)指向（1'1，7'1)的向量的方向作为剥蚀线走向和沉积期构造趋势面走向； 将从(X'Q，y'Q)指向（1'1，7'1)的向量的垂直方向作为沉积期构造趋势面倾向； 根据公式确定沉积期构造趋势面上各点坐标(X'，y'，z'）；其中，（x'i，y'i，z'i)为过沉积期构造 趋势面上的某点坐标；（1^，7^，2^)为沿沉积期构造趋势面倾向预设距离为1'的某点的坐 标;α为预先设置的沉积期构造趋势面倾角；Θ为剥蚀线走向与方位坐标系的正北方向的夹 角；" =?Χ (.Υ"―1';.) + 7为过(X 'j，7 'j，z 'j)点，且'沿沉积其月构造趋势面走向方向的直线 方程;zx", / = z ' j为过(X ' j，y ' j，z ' j)点且沿沉积期构造趋势面走向方向的直线的高程值； 根据沉积期构造趋势面上各点坐标(X'，y'，z'）构成沉积期构造趋势面图数据。5. 根据权利要求4所述的沉积期微幅度古隆起的推算方法，其特征在于，根据所述沉积 期构造趋势面图数据和所述拟构造幅度确定沉积期构造的相对高低数据，包括:将各网格 点的沉积期构造趋势面图数据与其拟构造幅度相加，确定沉积期构造的相对高低数据。6. -种沉积期微幅度古隆起的推算装置，其特征在于，包括： 现今构造图数据获取单元，用于获取待推算地区地层的现今构造图数据； 现今构造趋势面图数据确定单元，用于根据所述现今构造图数据确定现今构造趋势面 图数据； 拟构造幅度确定单元，用于根据所述现今构造图数据和所述现今构造趋势面图数据确 定各井间网格点的拟构造幅度;所述拟构造幅度为现今构造图上的点相对于现今构造趋势 面的竖直方向幅度； 沉积期构造趋势面图数据确定单元，用于获取现今构造图的剥蚀线数据，并根据所述 现今构造图的剥蚀线数据确定待推算地区地层的沉积期构造趋势面图数据； 沉积期微幅度古隆起图数据生成单元，用于根据所述沉积期构造趋势面图数据和所述 拟构造幅度确定沉积期构造的相对高低数据，并根据所述沉积期构造的相对高低数据生成 沉积期微幅度古隆起图数据； 沉积微相图数据确定单元，用于根据所述沉积期微幅度古隆起图数据确定沉积微相图 数据。7. 根据权利要求6所述的沉积期微幅度古隆起的推算装置，其特征在于，所述现今构造 趋势面图数据确定单元，具体用于:在现今构造图数据对应的现今构造图中获取第一数据 点和第二数据点；其中，所述第一数据点的坐标为（ XQ，yQ，Z())，所述第二数据点的坐标为 (Xn，y n，Zn);所述第一数据点为现今构造图中的最低点;所述第二数据点为现今构造图中的 以最1?点为中心，以预设距尚为半径的范围内的最低点； 在第一数据点和第二数据点连接的直线上每隔一预设高差A Z确定一第三数据点；所 述第三数据点的坐标为(Xk，yk，Zk);0 -1 = "Δζ。?，: 根据公另确定空间内各点坐标(X，y，ζ )，并根据克里格插值算法对空间内各点坐标进行插值，生 成现今构造趋势面图数据;在直自t的点的高程值均为zk。8. 根据权利要求7所述的沉积期微幅度古隆起的推算装置，其特征在于，所述拟构造幅 度确定单元，具体用于： 将现今构造图数据中的各井间网格点的坐标与现今构造趋势面图数据中的各点坐标 在竖直方向上一一对应； 确定现今构造图数据中的各井间网格点与其对应的现今构造趋势面图数据中的点在 竖直方向上的幅度差，并根据各幅度差确定各井间网格点的拟构造幅度。9. 根据权利要求8所述的沉积期微幅度古隆起的推算装置，其特征在于，沉积期构造趋 势面图数据确定单元，具体用于： 在现今构造图的剥蚀线数据对应的剥蚀线上选取任意两点，并以其中一点为坐标原点 建立方位坐标系;所述方位坐标系的正北方向为0°，正东方向为90° ;其中，在所述方位坐标 系中原点坐标为(1'〇,7'〇)，所述任意两点中的另一点坐标为&'1，7'1); 将从(X'Q，y'Q)指向（1'1，7'1)的向量的方向作为剥蚀线走向和沉积期构造趋势面走向； 将从(X'Q，y'Q)指向（1'1，7'1)的向量的垂直方向作为沉积期构造趋势面倾向； 根据公式：确定沉积期构造趋势面上各点坐标(X'，y'，Z'）；其中，（x'i，y'i，z'i)为过沉积期构造 趋势面上的某点坐标；（1^，7^，2^)为沿沉积期构造趋势面倾向预设距离为1'的某点的坐 标;α为预先设置的沉积期构造趋势面倾角；Θ为剥蚀线走向与方位坐标系的正北方向的夹 角；少"=? χ (·γ "一.Y ) + 为过(X ' j，y ' j，z ' j)点，且沿沉积期构造趋势面走向方向的直线 方程;ZX", y" = z ' j为过(X ' j，y ' j，z ' j)点且沿沉积期构造趋势面走向方向的直线的高程值； 根据沉积期构造趋势面上各点坐标(X'，y'，z'）构成沉积期构造趋势面图数据。10.根据权利要求9所述的沉积期微幅度古隆起的推算装置，其特征在于，所述沉积期 微幅度古隆起图数据生成单元，具体用于:将各网格点的沉积期构造趋势面图数据与其拟 构造幅度相加，确定沉积期构造的相对高低数据。
【文档编号】G06Q50/02GK105868482SQ201610213582
【公开日】2016年8月17日
【申请日】2016年4月7日
【发明人】王建国, 常宏, 杜支文, 顾岱鸿, 孙振, 成育红, 王彦良, 纪彦博, 杨璐菁
【申请人】中国石油大学(北京), 北京海贝新源石油科技有限公司