CN114059999A - 一种重力流沉积成因测井识别方法 - Google Patents

一种重力流沉积成因测井识别方法 Download PDF

Info

Publication number
CN114059999A
CN114059999A CN202111151678.6A CN202111151678A CN114059999A CN 114059999 A CN114059999 A CN 114059999A CN 202111151678 A CN202111151678 A CN 202111151678A CN 114059999 A CN114059999 A CN 114059999A
Authority
CN
China
Prior art keywords
gravity flow
facies
sedimentary
logging
lithofacies
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111151678.6A
Other languages
English (en)
Other versions
CN114059999B (zh
Inventor
杨田
田景春
余文强
李晓芳
何青
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN202111151678.6A priority Critical patent/CN114059999B/zh
Publication of CN114059999A publication Critical patent/CN114059999A/zh
Application granted granted Critical
Publication of CN114059999B publication Critical patent/CN114059999B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B25/00—Apparatus for obtaining or removing undisturbed cores, e.g. core barrels or core extractors
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开一种重力流沉积成因测井识别方法，包括以下步骤：先进行深水重力流沉积岩芯的岩相组合划分及厚度统计，再进行深水重力流沉积岩相类型测井识别，最后进行测井岩相组合类型划分与成因识别；本发明通过多测井曲线组合，有效解决了单一测井曲线识别重力流沉积岩相误差大，识别不准确的难题，使得重力流沉积岩相测井识别更加精确，以岩芯分析不同类型岩相组合厚度分布为约束条件，实现了重力流沉积岩相组合类型的测井识别，使得重力流沉积成因的测井识别成为了可能，整套方法有效解决了无取芯井段重力流沉积成因识别的难题，对指导无取芯或少取芯井区重力流沉积成因识别与砂体分布预测的研究具有重要意义。

Description

一种重力流沉积成因测井识别方法
技术领域
本发明涉及石油天然气勘探与开发技术领域，尤其涉及一种重力 流沉积成因测井识别方法。
背景技术
我国深水重力流沉积在海相地层、陆相地层、不同地质时代、不 同埋深均广泛发育，其中蕴含丰富的油气资源，勘探潜力巨大，已成 为油气勘探开发的重点领域，但由于深水重力流沉积特征多样、成因 机制复杂，导致对其分布规律认识不清，这是制约我国深水重力流沉 积油气高效勘探开发的主要原因，重力流沉积成因多样，一次重力流 事件可形成具有多种沉积成因的沉积产物，单一沉积成因的沉积物多 具有单层厚度薄且分布广的特征，精细的岩芯岩相和岩相组合分析是 识别其沉积成因的最可靠方法，但钻井取芯成本高昂，岩芯数量有限， 如何有效识别大部分未取芯的钻井中钻遇的重力流沉积的成因，是开 展重力流沉积分布规律总结进而有效指导重力流沉积油气高效勘探 开发的前提。
测井曲线纵向连续性好，且未取芯的钻井同样具有丰富的测井数 据，大量研究均通过开展沉积物的测井识别来明确未取芯井段沉积特 征，开展沉积物分布规律预测，即通过岩芯观察，建立不同岩相特征 岩芯与测井响应之间的联系，从而达到通过测井响应值来识别未取芯 段沉积特征的目的，但该方法在识别重力流沉积特征及其成因上还存 在以下问题：(1)单一测井响应值不能准确的识别重力流沉积类型， 由于重力流沉积具有杂基含量高，砂泥频繁互层的特征，导致其岩性 区分较为困难，存在大量的伽玛测井值相对较高的砂岩或者泥质砂 岩，所以单纯的采用以往对岩性识别较为准确的自然伽玛测井值并不能有效的区分重力流沉积的不同岩相；(2)单一的岩相类型测井识 别不能准确的判断重力流沉积成因，现已知重力流流体类型多样导致 其沉积成因多样，包含滑动滑塌沉积、砂质碎屑流沉积、泥质碎屑流 沉积、混合重力流沉积、低密度浊流沉积、异重流沉积等多种类型， 相同的岩相类型在不同成因中均有发育，如块状砂岩在滑动沉积、砂 质碎屑流沉积、混合重力流沉积、异重流沉积中均有发育，泥砂混杂 的块状沉积在滑塌沉积、泥质碎屑流沉积、混合重力流沉积中均匀发 育，砂泥频繁互层沉积在混合重力流沉积、低密度浊流沉积、异重流 沉积中均有发育，所以单一的岩相分析无法准确的判断重力流沉积成 因，进而无法开展重力流沉积分布预测。
综上所述，针对深水重力流沉积，需要结合实际地质条件开展重 力流沉积岩芯的岩相与岩相组合划分研究，在此基础上，通过基础岩 相类型的测井识别，测井岩相组合类型的划分，从而达到测井识别重 力流沉积成因的目的，因此，本发明提出一种重力流沉积成因测井识 别方法以解决现有技术中存在的问题。
发明内容
针对上述问题，本发明的目的在于提出一种重力流沉积成因测井 识别方法，该方法在充分考虑深水重力流沉积动力多样、重力流沉积 特征复杂等基础上，在重力流沉积岩芯岩相、岩相组合特征精细划分 和量化分析的基础上，通过多类型测井曲线组合贝叶斯函数判别准确 识别重力流沉积岩相类型，以岩芯岩相组合厚度分布为约束，开展测 井岩相组合划分，从而实现通过测井响应识别无取芯井重力流沉积成 因。
为了实现本发明的目的，本发明通过以下技术方案实现：一种重 力流沉积成因测井识别方法，包括以下步骤：
步骤一：深水重力流沉积岩芯的岩相组合划分及厚度统计
先对重力流沉积岩芯的岩相类型进行划分，接着在重力流沉积岩 相类型划分的基础上对重力流沉积岩芯的岩相组合类型进行划分，然 后在重力流沉积岩相和岩相组合类型划分的基础上对岩相组合类型 厚度进行量化统计；
步骤二：深水重力流沉积岩相类型测井识别
先对测井曲线数据进行规范化预处理，再利用岩相类型和测井曲 线关系确定判别函数，接着通过测井识别岩相类型；
步骤三：测井岩相组合类型划分与成因识别
先对测井岩相组合类型进行划分，再明确重力流沉积成因。
进一步改进在于：所述步骤一中，对重力流沉积岩芯的岩相类型 进行划分的具体步骤为：从重力流沉积近端到重力流沉积远端，选择 覆盖重力流沉积范围内不同分布位置重力流沉积岩芯，开展物质组 成、岩性、粒度、颜色和沉积构造方面的详细研究，以沉积物粒度为 基础，结合沉积构造和物质组成，划分重力流沉积岩相类型。
进一步改进在于：所述步骤一中，对重力流沉积岩芯的岩相组合 类型进行划分的具体步骤为：在重力流沉积岩相类型划分的基础上， 以稳定泥岩层或底冲刷层为分界标准，划分一次重力流事件在垂向上 形成的沉积物岩相组合类型，根据岩相组合类型与沉积动力之间的对 应关系，总结岩相组合模式和重力流沉积成因。
进一步改进在于：所述步骤一中，对岩相组合类型厚度进行量化 统计的具体步骤为：在重力流沉积岩相和岩相组合类型划分的基础 上，对取芯井开展系统的岩相、岩相组合划分，利用CoreLDRAW2018 软件，绘制岩芯观察分析得到的岩相、岩相组合在垂向上的叠置的沉 积综合柱状图，统计不同岩相组合类型的厚度并形成数据库，明确不 同成因的组合类型厚度分布范围。
进一步改进在于：所述步骤二中，对测井曲线数据进行规范化预 处理包括测井曲线标准化、曲线拼接、深度校正和岩芯校正，以保证 岩芯深度与测井曲线深度的准确对应。
进一步改进在于：所述步骤二中，识别判别函数的具体步骤为： 将划分出的岩相类型与测井曲线对应，从而明确不同岩相类型的测井 响应值，选取多类型的测井响应值，建立不同岩相类型测井响应值的 数据库，使用SPSS软件对数据库进行处理分析，建立不同岩相类型 基于多类型测井响应值综合的贝叶斯判别函数。
进一步改进在于：所述步骤二中，测井识别岩相类型的具体步骤 为：使用Matlab软件编写程序，按照每个测点间距为0.125m对测井 数据采样贝叶斯函数进行批量计算，根据贝叶斯判别函数后验概率值 最大这一判别规则，选取每个测点中计算结果的最大值，这一判别函 数代表的岩相类型即为该测点对应的岩相类型，将计算结果导入石文 或ResForm软件绘制图件，与岩芯观察岩性综合柱状图进行对比，检 验测井识别岩相类型的准确度。
进一步改进在于：所述步骤三中，测井岩相组合类型划分与成因 识别具体步骤为：根据测井识别获得的不同岩相类型垂向叠置关系， 结合岩芯统计获得的不同成因岩相组合厚度分布范围，对测井岩相组 合类型进行划分，综合岩相组合类型与沉积动力之间的关系，明确重 力流沉积成因。
本发明的有益效果为：本发明通过多测井曲线组合，有效解决了 单一测井曲线识别重力流沉积岩相误差大，识别不准确的难题，使得 重力流沉积岩相测井识别更加精确，以岩芯分析不同类型岩相组合厚 度分布为约束条件，实现了重力流沉积岩相组合类型的测井识别，使 得重力流沉积成因的测井识别成为了可能，整套方法有效解决了无取 芯井段重力流沉积成因识别的难题，对指导无取芯或少取芯井区重力 流沉积成因识别与砂体分布预测的研究具有重要意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面 将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而 易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域 普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这 些附图获得其他的附图。
图1是本发明的方法流程图；
图2是本发明实施例中的深水重力流沉积岩性柱状图；
图3是本发明实施例中的深水重力流沉积岩相组合类型图；
图4是本发明实施例中的深水重力流沉积测井数据预处理图；
图5a是本发明实施例中的砂岩相测井曲线特征图；
图5b是本发明实施例中的砂泥互层岩相测井曲线特征图；
图5c是本发明实施例中的砂质泥岩相测井曲线特征图；
图5d是本发明实施例中的暗色泥岩相测井曲线特征图；
图5e是本发明实施例中的黑色页岩相测井曲线特征图；
图6是本发明实施例中的深水重力流沉积岩相组合类型划分图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方 案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部 分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普 通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例， 都属于本发明保护的范围。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、 “左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方 位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述 本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特 定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限 制。此外，术语“第一”、“第二”、“第三”、“第四”等仅用于 描述目的，而不能理解为指示或暗示相对重要性。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定， 术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固 定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也 可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可 以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具 体情况理解上述术语在本发明中的具体含义。
参见图1，本实施例提供了一种重力流沉积成因测井识别方法， 包括以下步骤：
步骤一：深水重力流沉积岩芯的岩相组合划分及厚度统计
S1.先对重力流沉积岩芯的岩相类型进行划分
以长庆油田某地区深水重力流沉积岩为研究对象，结合前人研究 成果以及油田已有资料，研究区西南、南部至中心部位，主要发育半 深湖-深湖相沉积，西缘、西南缘、南缘及研究区东北部发育大范围 的三角洲相沉积，根据重力流发育规律，从重力流沉积近端到重力流 沉积远端，选择覆盖重力流沉积范围内不同分布位置重力流沉积岩 芯，开展物质组成、岩性、粒度、颜色、沉积构造等方面的详细研究， 以沉积物粒度为基础，结合沉积构造和物质组成，划分重力流沉积岩 相类型，研究区岩性以细砂岩、粉砂岩和泥岩为主，局部发育少量中 粗砂岩，部分井段可见薄层凝灰岩发育，主要发育的沉积构造包括块 状层理、平行层理、正粒序层理、逆粒序层理、韵律层理、沙纹层理、 水平层理和变形构造，部分岩芯中可见变形的泥岩撕裂屑或砂质团 块，在此基础上，以单一岩性及沉积构造的组合划分研究区内重力流 沉积岩相类型，分别是：块状层理细砂岩相、平行层理细砂岩相、正 粒序层理细砂岩相、逆粒序层理细砂岩相、块状层理粉砂岩相、平行 层理粉砂岩相、沙纹层理粉砂岩相、韵律层理粉砂岩相、泥质碎屑砂 岩相、泥质粉砂岩相、变形构造砂岩相、块状层理泥岩相、水平层理 泥岩相；
S2.对重力流沉积岩芯的岩相组合类型进行划分
在重力流沉积岩相类型划分的基础上，以稳定泥岩层或底冲刷层 为分界标准，划分一次重力流事件在垂向上形成的沉积物岩相组合类 型，根据岩相组合类型与沉积动力之间的对应关系，总结岩相组合模 式(参见图2)和重力流沉积成因，结合重力流沉积动力类型，划分 岩相组合类型，分别是：滑动滑塌沉积岩相组合、砂质碎屑流沉积岩 相组合、泥质碎屑流沉积岩相组合、低密度浊流沉积岩相组合、混合 重力流沉积岩相组合和异重流沉积岩相组合(参见图3)，其中，滑 动滑塌沉积岩相组合主要由变形构造粉砂岩相、泥质碎屑泥质粉砂岩 相、泥质粉砂岩相或粉砂质泥岩相组成；砂质碎屑流沉积岩相组合主 要由块状层理砂岩相组成，内部也可含有一定量的泥质碎屑；泥质碎 屑流沉积岩相组合主要由泥质碎屑砂岩相、泥质粉砂岩相或粉砂质泥 岩相组成；低密度浊流沉积岩相组合主要由平行层理细(粉)砂岩相、 正粒序细(粉)砂岩相、沙纹层理粉砂岩相、韵律层理粉砂岩相、块 状或水平层理泥岩相组成；混合重力流沉积岩相组合自下而上由H1 段和H3段成对组合出现，H1段主要是块状细砂岩相，H3段主要有泥 质碎屑砂岩相、泥质粉砂岩相或粉砂质泥岩相；异重流沉积岩相组合 由逆-正粒序细(粉)砂岩相、平行层理细(粉)砂岩相、沙纹层理 粉细砂岩相组成；
S2.对岩相组合类型厚度进行量化统计
在重力流沉积岩相和岩相组合类型划分的基础上，对取芯井开展 系统的岩相、岩相组合划分，利用CoreLDRAW2018软件，绘制岩芯 观察分析得到的岩相、岩相组合在垂向上的叠置的沉积综合柱状图， 统计不同岩相组合类型的厚度并形成数据库，明确不同成因的组合类 型厚度分布范围，形成岩芯厚度数据库，统计岩相组合的厚度分布范 围(参见表1)，滑动滑塌沉积岩相组合厚度一般在0.06-4.18m之 间，平均厚度1.29m，厚度最大可达7.22m；砂质碎屑流沉积岩相组 合厚度一般在0.07-2.30m之间，平均厚度0.90m，厚度最大可达 6.20m；泥质碎屑流沉积岩相组合厚度在0.21m至2.79m之间，平均 厚度为1.08m；混合重力流沉积岩相组合厚度一般在0.10m-2.20m之 间，平均厚度0.99m，最大厚度可达3.50m，其中下部H1段厚度一般 0.10m-1.60m，平均厚度0.56m，上部H3段厚度一般0.01m-1.20m，平均厚度0.43m；低密度浊流沉积岩相组合厚度一般在0.02m至1.21m 之间，平均厚度0.50m，最大可达6.16m；异重流沉积岩相组合单层 厚度差异较大，从0.10m-0.70m之间都有分布，平均厚度0.30m；
表1长庆油田某地区深水重力流沉积岩相类型厚度表
Figure BDA0003287356170000091
步骤二：深水重力流沉积岩相类型测井识别
S1.对测井曲线数据进行规范化预处理
在明确不同岩相类型的测井响应值前，需要对测井曲线进行规范 化处理，主要涉及测井曲线标准化、曲线拼接、深度校正和岩芯校正 等，以保证岩芯深度与测井曲线深度的准确对应(参见图4)，将测 井、录井数据导入ResForm软件中生成综合录井图，厚层泥岩的测井 曲线特征相对较为明显，以此为依据与符合特征的测井曲线对应，消 除岩电差；将厘米级的精细岩芯沉积学描述综合柱状图与录井图进行 比对，根据岩芯沉积与测井曲线的实际对应关系，调整岩芯采样深度 使之与测井曲线形态对应，实现岩芯深度校正，明确取芯段的实际深 度；
S2.利用岩相类型和测井曲线识别判别函数
在测井数据预处理的基础上，将划分出的岩相类型与测井曲线对 应，从而明确不同岩相类型的测井响应值大小，选取多类型的测井响 应值，建立不同岩相类型测井响应的数据库，根据沉积物粒度和黏土 矿物含量将岩相类型简化为五种，分别是Y1砂岩相(参见图5a)、 Y2砂泥互层岩相(参见图5b)、Y3砂质泥岩相(参见图5c)、Y4 暗色泥岩相(参见图5d)及Y5黑色页岩相(参见图5e)；结合测井 资料，选取自然电位(SP)、自然伽玛(GR)、密度(DEN)、声波 时差(AC)及补偿中子(CNL)这5种测井数据类型形成数据库并建 立测井识别函数；Y1具有低伽马、中低电位、中等密度、中低中子、 低声波时差的特征；Y2具有中等伽马、高电位、高密度、中低中子、 中低声波时差的特征；Y3具有中低伽马、中等电位、中高密度、中 等中子、中低声波时差的特征；Y4具有中伽马、中电位、中等密度、 中高中子、中声波时差的特征；Y5具有高伽马、中高电位、中低密 度、中高中子、中等声波时差的特征(参见表2)；使用SPSS软件 对数据库进行分类判别，得出Y1至Y5的识别准确度(参见表3)及 判别函数各常数(参见表4)，其中Y1、Y2及Y5识别准确度高，Y3 和Y4识别准确度较低，通过代表性测井响应的交会图分析，可进一 步区分Y3和Y4这2种岩相类型，在明确区分5种岩相类型的基础上 建立不同岩相类型基于多类型测井响应值综合的贝叶斯判别函数(参 见表5)；
表2长庆油田某地区深水重力流沉积岩相类型测井值特征表
Figure BDA0003287356170000111
表3长庆油田某地区深水重力流沉积岩相类型测井识别准确度
Figure BDA0003287356170000121
表4长庆油田某地区深水重力流沉积岩相类型贝叶斯判别函数常数 表
Figure BDA0003287356170000122
表5长庆油田某地区深水重力流沉积岩相类型贝叶斯判别函数
Figure BDA0003287356170000123
S3.通过测井识别岩相类型
使用Matlab软件编写程序，按照每个测点间距为0.125m对测井 数据进行采样，同时运行贝叶斯函数对每口井的测井数据进行计算， 根据贝叶斯判别函数后验概率值最大这一判别规则，选取每个测点中 计算结果的最大值，这一判别函数代表的岩相类型即为该测点对应的 岩相类型，将计算结果导入石文或ResForm软件生成岩相的测井识别 结果柱状图，与岩芯实际观察柱状综合图进行对比，检验测井识别岩 相类型的准确度，识别准确度的计算方式为所观察岩芯中与测井识别 对应准确的累积厚度与所观察岩芯的总厚度的比值，以研究区内3口 全取心井为例进行测井识别，并与钻井岩芯实际观察柱状综合图对比 岩相类型，结果表明对深水重力流沉积可以有效识别，识别准确度在 80％以上；
步骤三：测井岩相组合类型划分与成因识别
根据测井识别获得的不同岩相类型垂向叠置关系，结合岩芯统计 获得的不同成因岩相组合厚度分布范围，对测井岩相组合类型进行划 分；以里231井(参见图6a)及白522井(参见图6b)长73段地层 的岩相测井识别结果为例，结合厚度数据划分岩相组合。
根据岩相组合厚度分布范围，里231井2096m及2098m发育厚度 接近的块状砂岩相，但2098m处块状砂岩相上覆砂质泥岩相，根据厚 度划分为混合重力流沉积岩相组合，而前者划分为砂质碎屑流沉积岩 相组合；里231井2092m处发育厚度约1.2m的砂质泥岩相，白522 井1940m处发育厚度约1.25m的砂质泥岩相，根据岩芯厚度量化统计， 二者均可划归滑动滑塌沉积岩相组合，但后者顶底均为暗色泥岩相， 根据岩相叠置组合规律，白522井1940m处的砂质泥岩相应划归泥质 碎屑流沉积岩相组合；白522井1961m处发育厚度约0.38m的砂质泥 岩相与0.62m的块状砂岩相，顶底均为湖相泥岩，同时根据其叠置情 况，累计厚度符合混合事件层厚度分布，划分为混合重力流沉积岩相 组合。
综上所述，里231井及白522井长73段地层由底至顶表现为砂 质碎屑流沉积、滑动滑塌沉积、黑色页岩、混合重力流沉积、暗色泥 岩、低密度浊流沉积这六种岩相组合的有序叠置；综合岩相组合类型 与沉积动力之间的关系明确了深水重力流沉积动力机制，暗色泥岩及 黑色页岩这两种岩相组合对应正常湖相悬浮泥质沉积。
因此，里231井长73亚段从底至顶由砂质碎屑流、滑动滑塌过渡 为正常湖相泥质沉积，之后发育多期混合重力流，被正常湖相泥岩分 隔，顶部过渡为滑动滑塌，并向低密度浊流转化；白522井长73亚段 底部发育低密度浊流及砂质碎屑流沉积，向上过渡为正常湖相黑色页 岩及暗色泥岩，其间发育有多期泥质碎屑流及低密度浊流沉积。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明， 凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等， 均应包含在本发明的保护范围之内。

Claims (8)

1.一种重力流沉积成因测井识别方法，其特征在于：包括以下步骤：
步骤一：深水重力流沉积岩芯的岩相组合划分及厚度统计
先对重力流沉积岩芯的岩相类型进行划分，接着在重力流沉积岩相类型划分的基础上对重力流沉积岩芯的岩相组合类型进行划分，然后在重力流沉积岩相和岩相组合类型划分的基础上对岩相组合类型厚度进行量化统计；
步骤二：深水重力流沉积岩相类型测井识别
先对测井曲线数据进行规范化预处理，再利用岩相类型和测井曲线关系确定判别函数，接着通过测井识别岩相类型；
步骤三：测井岩相组合类型划分与成因识别
先对测井岩相组合类型进行划分，再明确重力流沉积成因。
2.根据权利要求1所述的一种重力流沉积成因测井识别方法，其特征在于：所述步骤一中，对重力流沉积岩芯的岩相类型进行划分的具体步骤为：从重力流沉积近端到重力流沉积远端，选择覆盖重力流沉积范围内不同分布位置重力流沉积岩芯，开展物质组成、岩性、粒度、颜色和沉积构造方面的详细研究，以沉积物粒度为基础，结合沉积构造和物质组成，划分重力流沉积岩相类型。
3.根据权利要求1所述的一种重力流沉积成因测井识别方法，其特征在于：所述步骤一中，对重力流沉积岩芯的岩相组合类型进行划分的具体步骤为：在重力流沉积岩相类型划分的基础上，以稳定泥岩层或底冲刷层为分界标准，划分一次重力流事件在垂向上形成的沉积物岩相组合类型，根据岩相组合类型与沉积动力之间的对应关系，总结岩相组合模式和重力流沉积成因。
4.根据权利要求1所述的一种重力流沉积成因测井识别方法，其特征在于：所述步骤一中，对岩相组合类型厚度进行量化统计的具体步骤为：在重力流沉积岩相和岩相组合类型划分的基础上，对取芯井开展系统的岩相、岩相组合划分，利用CoreLDRAW 2018软件，绘制岩芯观察分析得到的岩相、岩相组合在垂向上的叠置的沉积综合柱状图，统计不同岩相组合类型的厚度并形成数据库，明确不同成因的组合类型厚度分布范围。
5.根据权利要求1所述的一种重力流沉积成因测井识别方法，其特征在于：所述步骤二中，对测井曲线数据进行规范化预处理包括测井曲线标准化、曲线拼接、深度校正和岩芯校正，以保证岩芯深度与测井曲线深度的准确对应。
6.根据权利要求1所述的一种重力流沉积成因测井识别方法，其特征在于：所述步骤二中，识别判别函数的具体步骤为：将划分出的岩相类型与测井曲线对应，从而明确不同岩相类型的测井响应值，选取多类型的测井响应值，建立不同岩相类型测井响应值的数据库，使用SPSS软件对数据库进行处理分析，建立不同岩相类型基于多类型测井响应值综合的贝叶斯判别函数。
7.根据权利要求1所述的一种重力流沉积成因测井识别方法，其特征在于：所述步骤二中，测井识别岩相类型的具体步骤为：使用Matlab软件编写程序，按照每个测点间距为0.125m对测井数据采样贝叶斯函数进行批量计算，根据贝叶斯判别函数后验概率值最大这一判别规则，选取每个测点中计算结果的最大值，这一判别函数代表的岩相类型即为该测点对应的岩相类型，将计算结果导入石文或ResForm软件绘制图件，与岩芯观察岩性综合柱状图进行对比，检验测井识别岩相类型的准确度。
8.根据权利要求1所述的一种重力流沉积成因测井识别方法，其特征在于：所述步骤三中，测井岩相组合类型划分与成因识别具体步骤为：根据测井识别获得的不同岩相类型垂向叠置关系，结合岩芯统计获得的不同成因岩相组合厚度分布范围，对测井岩相组合类型进行划分，综合岩相组合类型与沉积动力之间的关系，明确重力流沉积成因。
CN202111151678.6A 2021-09-29 2021-09-29 一种重力流沉积成因测井识别方法 Expired - Fee Related CN114059999B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111151678.6A CN114059999B (zh) 2021-09-29 2021-09-29 一种重力流沉积成因测井识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111151678.6A CN114059999B (zh) 2021-09-29 2021-09-29 一种重力流沉积成因测井识别方法

Publications (2)

Publication Number Publication Date
CN114059999A true CN114059999A (zh) 2022-02-18
CN114059999B CN114059999B (zh) 2023-06-06

Family

ID=80233879

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111151678.6A Expired - Fee Related CN114059999B (zh) 2021-09-29 2021-09-29 一种重力流沉积成因测井识别方法

Country Status (1)

Country Link
CN (1) CN114059999B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117607965A (zh) * 2023-12-27 2024-02-27 成都理工大学 基于地震属性与水平井轨迹分析的河道砂体内幕结构解析方法

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4260020A (en) * 1979-09-04 1981-04-07 The Dow Chemical Company Method and tool for controlling fluid flow from a tubing string into a low pressure earth formation
US5844799A (en) * 1996-01-26 1998-12-01 Institut Francais Du Petrole Method for simulating the filling of a sedimentary basin
WO2006116122A2 (en) * 2005-04-22 2006-11-02 Shell Internationale Research Maatschappij B.V. Systems and processes for use in treating subsurface formations
US20140328454A1 (en) * 2013-05-06 2014-11-06 Chevron U.S.A. Inc. System and method for determining an orientation of reservoir geobodies from unoriented conventional cores
CN204087631U (zh) * 2014-07-11 2015-01-07 山东科技大学 一种重力流运动与沉积模拟装置
CN104564041A (zh) * 2014-07-24 2015-04-29 中国石油大学(华东) 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法
CN104989392A (zh) * 2015-07-10 2015-10-21 中国石油天然气股份有限公司 一种岩性识别方法
CN105116467A (zh) * 2015-08-14 2015-12-02 中国石油大学(华东) 近岸水下扇砂砾岩沉积单元体划分方法
CN105158802A (zh) * 2015-08-21 2015-12-16 西南石油大学 坳陷湖盆重力流沉积物测井定量识别方法
CN105277976A (zh) * 2015-09-14 2016-01-27 山东科技大学 基于岩石露头雷达探测的地震正演模拟方法
CN105803994A (zh) * 2016-03-10 2016-07-27 成都理工大学 一种水下浊流形成河道的预测方法及应用
CN106988737A (zh) * 2017-04-28 2017-07-28 中国石油大港油田勘探开发研究院 一种利用岩性组合识别沉积相的方法
CN109444379A (zh) * 2018-12-24 2019-03-08 中海石油(中国)有限公司 深水重力流砂岩储层的定量分类识别图版构建方法及系统
US20200283324A1 (en) * 2019-03-08 2020-09-10 Anglo American Services (Uk) Ltd Tailings deposition
CN111896674A (zh) * 2020-07-15 2020-11-06 成都理工大学 一种页岩油气储集体沉积微相处理系统及鉴别方法
CN112034526A (zh) * 2020-08-13 2020-12-04 中国石油大学(华东) 基于岩相组合的灰质泥岩发育区薄层浊积砂体的地震识别方法

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4260020A (en) * 1979-09-04 1981-04-07 The Dow Chemical Company Method and tool for controlling fluid flow from a tubing string into a low pressure earth formation
US5844799A (en) * 1996-01-26 1998-12-01 Institut Francais Du Petrole Method for simulating the filling of a sedimentary basin
WO2006116122A2 (en) * 2005-04-22 2006-11-02 Shell Internationale Research Maatschappij B.V. Systems and processes for use in treating subsurface formations
US20140328454A1 (en) * 2013-05-06 2014-11-06 Chevron U.S.A. Inc. System and method for determining an orientation of reservoir geobodies from unoriented conventional cores
CN204087631U (zh) * 2014-07-11 2015-01-07 山东科技大学 一种重力流运动与沉积模拟装置
CN104564041A (zh) * 2014-07-24 2015-04-29 中国石油大学(华东) 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法
CN104989392A (zh) * 2015-07-10 2015-10-21 中国石油天然气股份有限公司 一种岩性识别方法
CN105116467A (zh) * 2015-08-14 2015-12-02 中国石油大学(华东) 近岸水下扇砂砾岩沉积单元体划分方法
CN105158802A (zh) * 2015-08-21 2015-12-16 西南石油大学 坳陷湖盆重力流沉积物测井定量识别方法
CN105277976A (zh) * 2015-09-14 2016-01-27 山东科技大学 基于岩石露头雷达探测的地震正演模拟方法
CN105803994A (zh) * 2016-03-10 2016-07-27 成都理工大学 一种水下浊流形成河道的预测方法及应用
CN106988737A (zh) * 2017-04-28 2017-07-28 中国石油大港油田勘探开发研究院 一种利用岩性组合识别沉积相的方法
CN109444379A (zh) * 2018-12-24 2019-03-08 中海石油(中国)有限公司 深水重力流砂岩储层的定量分类识别图版构建方法及系统
US20200283324A1 (en) * 2019-03-08 2020-09-10 Anglo American Services (Uk) Ltd Tailings deposition
CN111896674A (zh) * 2020-07-15 2020-11-06 成都理工大学 一种页岩油气储集体沉积微相处理系统及鉴别方法
CN112034526A (zh) * 2020-08-13 2020-12-04 中国石油大学(华东) 基于岩相组合的灰质泥岩发育区薄层浊积砂体的地震识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
TIAN YANG: "Gravity-flow deposits caused by different initiation processes in a deeplake system", 《RESEARCHGATE》 *
任小军等: "坳陷湖盆缓坡重力流成因储集层沉积特征及发育条件-以准噶尔盆地腹部石南31井区为例", 《新疆石油地质》 *
操应长等: "东营凹陷沙三中亚段三角洲前缘滑塌型重力流岩相类型及其分布特征", 《沉积与特提斯地质》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117607965A (zh) * 2023-12-27 2024-02-27 成都理工大学 基于地震属性与水平井轨迹分析的河道砂体内幕结构解析方法

Also Published As

Publication number Publication date
CN114059999B (zh) 2023-06-06

Similar Documents

Publication Publication Date Title
CN101930082B (zh) 采用电阻率资料进行储层流体类型判别的方法
Klaja et al. Geological interpretation of spectral gamma ray (SGR) logging in selected boreholes
CN108303752B (zh) 砂砾岩有效储层常规测井定量识别方法
CN104636588B (zh) 计算烃源岩中有机碳含量的方法及装置
CN102012526A (zh) 采用电阻率资料进行储层流体类型判别的方法
CN109653725A (zh) 基于沉积微相和岩石相的混积储层水淹程度测井解释方法
Harris The role of geology in reservoir simulation studies
CN110727034A (zh) 一种砂岩型铀矿编图方法
Wimmers et al. Integration of sedimentology, petrophysics and rock typing as key to understanding a tight gas reservoir
CN109814174B (zh) 一种碎屑岩不整合结构体测井综合识别方法
CN113792499B (zh) 基于沉积特征的松散承压含水层富水性动态确定方法
Alfy et al. Quantitative hydro-geophysical analysis of a complex structural karst aquifer in Eastern Saudi Arabia
CN114059999B (zh) 一种重力流沉积成因测井识别方法
CN112784404A (zh) 一种基于常规测井资料的砂砾岩束缚水饱和度计算方法
Hasoon et al. A Comprehensive Review for Integrating Petrophysical Properties, Rock Typing, and Geological Modeling for Enhanced Reservoir Characterization
CN109917489B (zh) 一种地下承压水位确定的新方法
Martyns-Yellowe Integrated Lithostratigraphic, Formation and Structural Evaluation of the Migrant Structure, Sable Subbasin, Offshore Nova Scotia.
Jageler et al. Use of well logs and dipmeters in stratigraphic-trap exploration: geologic exploration methods
Sokhal et al. An example from the Algerian Sahara illustrates the use of the hydraulic flow unit technique to discriminate fluid flow routes in confined sand reservoirs
Pradhanitasari Petrophysical characterizations and shale distribution analysis of Kissenda reservoir sands, Gabon
Behera et al. Geostatistical Modelling and Vertical Effect Analysis of a Ferruginous East Coast Bauxite Deposit, India
Savitri et al. A Novel Estimation Approach of Unconfined Compressive Strength (UCS) and Improved Rock Mechanical Properties Correlations by Considering the Effect of Overpressure Condition, Mineralogy, and Shale Composition in Montney Shale Gas, British Columbia, Canada
Crowley et al. The stag oilfield
Ayodele An integrated study of the early cretaceous (Valanginian) reservoir from the Gamtoos Basin, offshore South Africa with special reference to seismic cacies, formation evaluation and static reservoir modeling
CN117552780A (zh) 对源内不同类型夹层的测井识别方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20230606

CF01 Termination of patent right due to non-payment of annual fee