CN112149312B - 沉积微相的确定方法及装置、电子设备和存储介质 - Google Patents
沉积微相的确定方法及装置、电子设备和存储介质 Download PDFInfo
- Publication number
- CN112149312B CN112149312B CN202011087264.7A CN202011087264A CN112149312B CN 112149312 B CN112149312 B CN 112149312B CN 202011087264 A CN202011087264 A CN 202011087264A CN 112149312 B CN112149312 B CN 112149312B
- Authority
- CN
- China
- Prior art keywords
- grid
- determining
- updated
- river
- type
- 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.)
- Active
Links
- 230000008021 deposition Effects 0.000 title claims abstract description 143
- 238000000034 method Methods 0.000 title claims abstract description 102
- 238000003860 storage Methods 0.000 title claims abstract description 29
- 238000009826 distribution Methods 0.000 claims abstract description 104
- 238000010586 diagram Methods 0.000 claims abstract description 72
- 238000011160 research Methods 0.000 claims abstract description 24
- 239000004576 sand Substances 0.000 claims description 145
- 238000004590 computer program Methods 0.000 claims description 10
- 238000011835 investigation Methods 0.000 claims description 9
- 239000013049 sediment Substances 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 8
- 238000012545 processing Methods 0.000 description 22
- 238000005516 engineering process Methods 0.000 description 12
- 238000004891 communication Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 230000003287 optical effect Effects 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 230000005236 sound signal Effects 0.000 description 4
- 238000003491 array Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 2
- 238000011166 aliquoting Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 235000017166 Bambusa arundinacea Nutrition 0.000 description 1
- 235000017491 Bambusa tulda Nutrition 0.000 description 1
- 241001330002 Bambuseae Species 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 235000015334 Phyllostachys viridis Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000011425 bamboo Substances 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000010422 painting Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000000819 phase cycle Methods 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本公开涉及一种沉积微相的确定方法及装置、电子设备和存储介质，涉及油气田开发储层技术领域，所述的沉积微相的确定方法，包括：获取研究区的平面网格分布图；根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型；根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系；根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型。本公开实施例可解决了目前存在沉积微相不够准确问题。不能了满足储层精细描述高效、准确、快速的需求；本公开符合各沉积微相接触关系绘图特点，能实现准确地确定沉积微相。
Description
技术领域
本公开涉及油气田开发储层技术领域，尤其涉及一种沉积微相的确定方法及装置、电子设备和存储介质。
背景技术
大庆长垣油田经过50多年的开发，已进入特高含水期，剩余油呈现高度零散分布，调整挖潜对象逐渐由河道向河间砂转变，经过长期的注水开发，注采不完善的薄差油层动用状况仍然较差，剩余油相对富集，存在着较大精细开发调整潜力。为改善薄差油层开发效果，加快攻关薄差油层提高采收率技术，探索改善薄差油层开发效果新途径，对薄差油层密井网精细刻画迫在眉睫。
有关文献报道沉积微相成图技术，参见（1）李建雄，谷跃民，党虎强等“高密度井网开发区井震联合储层预测方法”（石油地球物理勘探，2011年第3期），（2）夏竹，李中超，贾瑞忠等“井震联合薄储层沉积微相表征实例研究”（石油地球物理勘探，2016年第5期），（3）李伟，岳大力，胡光义等“分频段地震属性优选及砂体预测方法——秦皇岛32-6油田北区实例”（石油地球物理勘探，2017年第2期），（4）吴景超，黄江波，王思权等“滚动勘探中薄层砂体厚度的预测研究及应用”（地球物理学进展，2017年第10期），（5）肖佃师，卢双舫，王海生等“三角洲外前缘薄砂体地震综合预测方法”（中国石油大学学报，2015年第8期）。上述（1）针对该区特定地震地质条件和实际资料条件，利用随机协同模拟技术为核心的反演技术以及井震结合预测砂岩厚度的方法；上述（2）以地质统计学及平面插值算法为数学工具，建立了储层目标砂比预测模型，并依据模型相关度，选择合适的平面插值算法进行砂地比成图，利用砂地比与薄储层沉积相的转换模式，结合单井相及沉积背景，获得了相对合理的河道砂厚度分布趋势；上述（3）综合测井与地震资料，针对薄层砂体与厚层砂体不同的地震属性响应特征，采用先优选地震资料频段，再提取并优选地震属性的方法进行砂体预测；上述（4）从经典的楔状模型出发，通过正演得出砂体厚度与峰值频率的相关规律，进而研究峰值频率和砂体不同厚度间的关系，从而拟合出适合目标区的线性关系式，从而实现目标区的厚度预测；上述（5）结合曲线重构、地震正演等技术，分别对适用于薄储层的地震反演及属性分析方法进行研究，依据地震属性对不同规模砂体的分辨能力指导属性优选及阈值选取，定性刻画砂层组内薄砂体的空间展布。
可见，上述方法均为利用地震属性、反演技术预测砂体厚度，不能自动绘制相带图，该方法存在预测河间砂体厚度时不够准确，且河间微相成图速度较慢的问题。为了满足储层精细描述高效、准确、快速的需求，发明了一种符合相控模式控制，符合各沉积微相接触关系绘图特点的相序模式控制的河间砂体厚度分布趋势的预测方法。
发明内容
本公开提出了一种沉积微相的确定方法及装置、电子设备和存储介质技术方案，解决了目前存在沉积微相不够准确问题，不能了满足储层精细描述高效、准确、快速的需求。
根据本公开的一方面，提供了一种沉积微相的确定方法，包括：获取研究区的平面网格分布图；
根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型；
根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系；
根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型。
优选地，在所述获取研究区的平面网格分布图之前，确定所述平面网格分布图；所述确定所述平面网格分布图的方法，包括：
获取研究区；
根据所述研究区的沉积相绘制的河道砂的分布特征；
对所述分布特征进行网格化，得到分布特征网格图；
根据所述研究区内已有井的数据对所述分布特征网格图的网格进行赋值，得到所述平面网格分布图。
优选地，所述研究区内已有井的数据，至少包括：二类砂岩厚度、一类砂岩厚度、有效厚度的一种或几种；
以及/或，
所述根据所述研究区内已有井的数据对所述分布特征网格图的网格进行赋值，得到所述平面网格分布图的方法，包括：
若所述网格内包含井，则根据所述井的数据对所述分布特征网格图的网格进行赋值；
若所述网格内不包含井，则根据其他网格内所述井的数据进行插值，得到插值数据；根据所述插值数据对不包含井的所述网格进行赋值。
优选地，所述根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型的方法，包括：
确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格；
基于所述河道砂网格以及/或所述非河道砂网格确定所述平面网格分布图内每个网格的待更新沉积微相类型。
优选地，所述确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格的方法，包括：
获取设定的地层厚度；
根据所述地层厚度及所述每个网格内的赋值确定河道砂最优门槛值；
根据所述河道砂最优门槛值确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格。
优选地，所述根据所述地层厚度及所述每个网格内的赋值确定河道砂最优门槛值的方法，包括：
根据所述每个网格内的赋值以及所述设定的地层厚度，确定所述每个网格内的砂地比；
根据所述砂地比确定对应的河道概率点及河间概率点；
分别根据所述河道概率点及河间概率点确定对应的河道概率线及河间概率线；
基于所述的河道概率线及河间概率线确定河道砂最优门槛值。
优选地，所述根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系的方法，包括：
根据所述积微相类型，对每个网格的所述沉积微相类型进行赋值，得到类型值；
将相邻网格的所述类型值作差，得到相邻网格的多个接触关系；
以及/或，
所述根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型的方法，包括：
根据所述平面网格分布图中垂直主河道方向的相邻接触关系确定所述相邻接触关系对应的网格内待更新沉积微相类型是否可能更新；
若不可能更新，则将待更新沉积微相类型确定为更新沉积微相类型；
若可能更新，进一步确定是否更新；
若不更新，则将待更新沉积微相类型确定为更新沉积微相类型；
若更新，确定需要更新的网格，并对所述需要更新的网格进一步细分得到细分网格，以及确定所述细分网格内的更新沉积微相类型。
优选地，所述根据所述平面网格分布图中垂直主河道方向的相邻接触关系确定所述相邻接触关系对应的网格内待更新沉积微相类型是否可能更新的方法，包括：
若所述相邻接触关系为0，说明河道与河道接触或非河道与非河道接触，此时不可能更新；
若所述相邻接触关系不为0，说明河道与非河道接触，此时可能更新。
优选地，所述若可能更新，进一步确定是否更新的方法，包括：
获取所述相邻接触关系沿主河道走向的相邻上侧的相邻接触关系；
若相邻上侧的相邻接触关系为0，则不更新；否则，更新。
优选地，所述若更新，确定需要更新的网格的方法，包括：
根据所述相邻接触关系确定所述河道与所述非河道接触的位置关系；
根据所述位置关系确定需要更新的网格。
优选地，所述根据所述相邻接触关系确定所述河道与所述非河道接触的位置关系的方法，包括：
若所述相邻接触关系大于0，所述河道在所述非河道左侧；
若所述相邻接触关系小于0，所述河道在所述非河道右侧。
优选地，所述根据所述位置关系确定需要更新的网格的方法，包括：
若所述河道在所述非河道左侧，则将相邻上侧的相邻接触关系中最大值对应的网格确定为需要更新的网格；
若所述河道在所述非河道右侧，则将相邻上侧的相邻接触关系中最小值对应的网格确定为需要更新的网格。
优选地，所述对所述需要更新的网格进一步细分得到细分网格，以及确定所述细分网格内的更新沉积微相类型的方法，包括：
根据沉积微相原理的预设规则，确定所述细分网格内的更新沉积微相类型。
根据本公开的一方面，提供了沉积微相的确定装置，包括：
获取单元，用于获取研究区的平面网格分布图；
第一确定单元，用于根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型；
第二确定单元，用于根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系；
更新确定单元，用于根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型。
根据本公开的一方面，提供了一种电子设备，包括：
处理器；
用于存储处理器可执行指令的存储器；
其中，所述处理器被配置为：执行上述沉积微相的确定方法。
根据本公开的一方面，提供了一种计算机可读存储介质，其上存储有计算机程序指令，所述计算机程序指令被处理器执行时实现上述沉积微相的确定方法。
本公开提出了一种沉积微相的确定方法及装置、电子设备和存储介质技术方案，解决了目前存在沉积微相不够准确问题，不能了满足储层精细描述高效、准确、快速的需求。本公开符合各沉积微相接触关系绘图特点，能实现准确地确定沉积微相，即可以确定河间砂体厚度分布趋势，更符合沉积特征。本公开为储层精细描述提供了可靠的技术支撑，大大提高了工作效率。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，而非限制本公开。
根据下面参考附图对示例性实施例的详细说明，本公开的其它特征及方面将变得清楚。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，这些附图示出了符合本公开的实施例，并与说明书一起用于说明本公开的技术方案。
图1示出根据本公开实施例的沉积微相的确定方法的流程图；
图2示出根据本公开实施例确定河道砂最优门槛值的原理图；
图3示出根据本公开实施例的单井相成图对应的分布特征网格图；
图4示出根据本公开实施例的分布特征网格图网格化及赋值得到的平面网格分布图；
图5示出根据本公开实施例确定沉积微相后的示意图；
图6是根据一示例性实施例示出的一种电子设备800的框图；
图7是根据一示例性实施例示出的一种电子设备1900的框图。
具体实施方式
以下将参考附图详细说明本公开的各种示例性实施例、特征和方面。附图中相同的附图标记表示功能相同或相似的元件。尽管在附图中示出了实施例的各种方面，但是除非特别指出，不必按比例绘制附图。
在这里专用的词“示例性”意为“用作例子、实施例或说明性”。这里作为“示例性”所说明的任何实施例不必解释为优于或好于其它实施例。
本文中术语“和/或”，仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。另外，本文中术语“至少一种”表示多种中的任意一种或多种中的至少两种的任意组合，例如，包括A、B、C中的至少一种，可以表示包括从A、B和C构成的集合中选择的任意一个或多个元素。
另外，为了更好地说明本公开，在下文的具体实施方式中给出了众多的具体细节。本领域技术人员应当理解，没有某些具体细节，本公开同样可以实施。在一些实例中，对于本领域技术人员熟知的方法、手段、元件和电路未作详细描述，以便于凸显本公开的主旨。
可以理解，本公开提及的上述各个方法实施例，在不违背原理逻辑的情况下，均可以彼此相互结合形成结合后的实施例，限于篇幅，本公开不再赘述。
此外，本公开还提供了沉积微相的确定装置、电子设备、计算机可读存储介质、程序，上述均可用来实现本公开提供的任一种沉积微相的确定方法，相应技术方案和描述和参见方法部分的相应记载，不再赘述。
图1示出根据本公开实施例的沉积微相的确定方法的流程图，如图1所示，所述沉积微相的确定方法，包括：步骤S101：获取研究区的平面网格分布图；步骤S102：根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型；步骤S103：根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系；步骤S104：根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型。解决了目前存在沉积微相不够准确问题，不能了满足储层精细描述高效、准确、快速的需求。本公开符合各沉积微相接触关系绘图特点，能实现准确地确定沉积微相。
骤S101：获取研究区的平面网格分布图。
在本公开中，在所述获取研究区的平面网格分布图之前，确定所述平面网格分布图；所述确定所述平面网格分布图的方法，包括：获取研究区；根据所述研究区的沉积相绘制的河道砂的分布特征；对所述分布特征进行网格化，得到分布特征网格图；根据所述研究区内已有井的数据对所述分布特征网格图的网格进行赋值，得到所述平面网格分布图。其中，所述的分布特征网格图的网格的数值一般沿着主河道走向以及垂直主河道方向逐渐增大。
在本公开的实施例中，首先选取研究区，根据所述研究区单井沉积相绘制的河道砂的分布特征，将研究区平面分布区域按照x及y两个方向进行网格化（按照需求10m、15m及20m均可），得到分布特征网格图。通常沿主河道走向为y方向，自北向南为网格的数值逐渐增大；垂直主河道方向为x方向，自西向东为网格的数值逐渐增大。
在本公开中，所述研究区内已有井的数据，至少包括：二类砂岩厚度、一类砂岩厚度、有效厚度的一种或几种。
在本公开中，所述根据所述研究区内已有井的数据对所述分布特征网格图的网格进行赋值，得到所述平面网格分布图的方法，包括：若所述网格内包含井，则根据所述井的数据对所述分布特征网格图的网格进行赋值；若所述网格内不包含井，则根据其他网格内所述井的数据进行插值，得到插值数据；根据所述插值数据对不包含井的所述网格进行赋值。
在本公开的实施例中，统计所选研究区内已有井的数据：二类砂岩厚度、一类砂岩厚度、井的有效厚度数据，对所述分布特征网格图的每个网格分别用二类砂岩厚度、一类砂岩厚度、有效厚度赋值；过井点的网格用井点测井解释数据赋值，非过井点网格应用克里金方法进行平面插值，分别得到其二类砂岩厚度、一类砂岩厚度、有效厚度的平面网格分布图。
步骤S102：根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型。
在本公开中，所述根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型的方法，包括：确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格；基于所述河道砂网格以及/或所述非河道砂网格确定所述平面网格分布图内每个网格的待更新沉积微相类型。
在本公开中，仅仅需要确定非河道砂对应的所述平面网格分布图内每个网格的待更新沉积微相类型。具体地说，获取多个设定值，根据所述多个设定值以及所述非河道砂网格的已有井的数据（砂岩厚度）确定所述平面网格分布图内每个非河道砂网格的待更新沉积微相类型。
在本公开的实施例中，河道相（河道砂）网格中的待更新沉积微相类型为河道砂；非河道相（非河道砂）网格中的待更新沉积微相类型根据二类砂岩厚度、一类砂岩厚度、有效厚度以及多个设定值确定主体砂、非主体砂、表外砂及尖灭泥对应的待更新沉积微相类型，多个设定值分别为0及0.5。即有效厚度大于0.5米为主体砂，有效厚度大于0米小于0.5米为非主体砂；有效厚度等于0、二类砂岩厚度大于0为表外砂；二类砂岩厚度小于0为尖灭泥。最后，对网格进行赋值（沉积微相）。其中，河道砂代表1，主体砂代表2，非主体砂代表3，表外砂代表4，尖灭泥代表5。也就是说，每个非河道砂网格的数值Z配置为2-5任一值，河道砂网格的数值Z配置为1。
在本公开中，所述确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格的方法，包括：获取设定的地层厚度；根据所述地层厚度及所述每个网格内的赋值确定河道砂最优门槛值；根据所述河道砂最优门槛值确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格。
在本公开的实施例中，统计研究区内已有的井点的砂地比（砂岩厚度与地层厚度的比值），通过概率统计分析，明确河道砂最优门槛值，这样就可以应用门槛值区分每个网格是河道砂（河道砂网格），还是非河道砂（非河道砂网格）；砂地比大于或等于河道砂最优门槛值为河道砂，砂地比小于河道砂最优门槛值为非河道砂。其中，砂岩厚度至少包括二类砂岩厚度、一类砂岩厚度、有效厚度的一种或几种；每口井的地层厚度是定值。
其中，所述根据所述地层厚度及所述每个网格内的赋值确定河道砂最优门槛值的方法，包括：根据所述每个网格内的赋值以及所述设定的地层厚度，确定所述每个网格内的砂地比；根据所述砂地比确定对应的河道概率点及河间概率点；分别根据所述河道概率点及河间概率点确定对应的河道概率线及河间概率线；基于所述的河道概率线及河间概率线确定河道砂最优门槛值。
图2示出根据本公开实施例确定河道砂最优门槛值的原理图。在本公开的具体实施例的图2中，基于所述河道概率线及河间概率线确定河道砂最优门槛值的方法，包括：基于所述砂地比，将所述河道概率线及河间概率线绘制到同一坐标系，将所述的河道概率线及河间概率线的交点确定河道砂最优门槛值。
在本公开的实施例的图2中，不同的砂地比对应不同的河道概率点、河间概率点；分别绘制河道概率点、河间概率点对应的河道概率曲线、河间概率曲线；河道概率曲线与河间概率曲线的交点为河道砂最优门槛值。砂地比的取值为0-1之间。
在本公开的实施例的图2中，砂地比大于或等于河道砂最优门槛值对应的河间概率曲线为非河道砂，所述非河道砂所在的平面网格分布图内的网格为非河道砂网格。砂地比小于或等于河道砂最优门槛值对应的河道概率曲线为河道砂，所述河道砂所在的平面网格分布图内的网格为河道砂网格。同理，我们可以认为砂地比小于或等于河道砂最优门槛值对应的河间概率曲线为河道砂，所述河道砂所在的平面网格分布图内的网格为河道砂网格。砂地比大于或等于河道砂最优门槛值对应的河道概率曲线为非河道砂，所述非河道砂所在的平面网格分布图内的网格为非河道砂网格。
步骤S103：根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系。
在本公开中，所述根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系的方法，包括：根据所述积微相类型，对每个网格的所述沉积微相类型进行赋值，得到类型值；将相邻网格的所述类型值作差，得到相邻网格的多个接触关系。
在本公开的实施例中，河道砂对应的网格配置（赋值）为1，主体砂对应的网格配置（赋值）2，非主体砂对应的网格配置（赋值）3，表外砂对应的网格配置（赋值）4，尖灭泥对应的网格配置（赋值）5。其中，河道砂中每个网格的数值Z配置为1，非河道砂中每个网格的数值Z为2-5任一值。
在本公开的实施例中，沉积微相接触关系值（相邻网格的多个接触关系）为N，N（xnyn，xnyn-1）代表x方向第n个，y方向第n个网格的沉积微相与x方向第n个，y方向第n-1个网格的沉积微相的接触关系，即N（xnyn，xnyn-1）=Z差=Z（xnyn）-Z（xnyn-1）。
在本公开的实施例中，沉积微相接触关系值（相邻网格的多个接触关系），包括5种接触关系，其中同一种沉积微相赋值为N（0），渐变相接触（如河道砂与主体砂）赋值为N（±1），过渡相接触（如河道砂与非主体砂）赋值为N（±2），突变相接触（如河道砂与表外砂）赋值为N（±3），极速突变相接触（如河道砂与尖灭泥）赋值为N（±4）。N（0）表示，N（0）为N=0；N（±1）表示，河道砂Z（xnyn）与主体砂接触Z（xnyn-1）为N=1，主体砂与河道砂接触为N=-1；N（±2）表示，河道与主非主体砂接触为N=2，非主体砂与河道砂接触为N=-2；；N（±3）表示，河道与表外砂接触为N=3，表外砂与河道砂接触为N=-3；N（±4）表示，河道砂与尖灭泥接触为N=4，尖灭泥与河道接触为N=-4。
步骤S104：根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型。
在本公开中，所述根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型的方法，包括：根据所述平面网格分布图中垂直主河道方向的相邻接触关系确定所述相邻接触关系对应的网格内待更新沉积微相类型是否可能更新；若不可能更新，则将待更新沉积微相类型确定为更新沉积微相类型；若可能更新，进一步确定是否更新；若不更新，则将待更新沉积微相类型确定为更新沉积微相类型；若更新，确定需要更新的网格，并对所述需要更新的网格进一步细分得到细分网格，以及确定所述细分网格内的更新沉积微相类型。
在本公开的实施例中，按照n=1，2，3，4…方向，其中，n≥1。先分析x方向（垂直主河道方向），再分析y（沿主河道走向）方向，计算所有网格的沉积微相的接触关系值N，分以下几种情况处理。
先分析得到河道与河道接触或非河道与非河道接触或河道与非河道接触。所述根据所述平面网格分布图中垂直主河道方向的相邻接触关系确定所述相邻接触关系对应的网格内待更新沉积微相类型是否可能更新的方法，包括：若所述相邻接触关系为0，说明河道与河道接触或非河道与非河道接触，此时不可能更新；此时，则将待更新沉积微相类型确定为更新沉积微相类型，待更新沉积微相类型没有改变。若所述相邻接触关系不为0，说明河道与非河道接触，此时可能更新。
例如，当N（xnyn，xn-1yn）=Z（xnyn）-Z（xn-1yn）= 0时，则沉积微相值不变，河间砂体厚度分布相近，河道与河道接触或非河道与非河道接触。说明河道与河道接触或非河道与非河道接触，此时不可能更新。
又例如，当N（xnyn，xn-1yn）=Z（xnyn）-Z（xn-1yn）不等于 0时，说明河道与非河道接触，此时可能更新。
在本公开的实施例中，所述若可能更新，进一步确定是否更新的方法，包括：获取所述相邻接触关系沿主河道走向的相邻上侧的相邻接触关系；若相邻上侧的相邻接触关系为0，则不更新；否则，更新。即，若相邻上侧的相邻接触关系不为0，则更新。
例如，当N（xnyn，xn-1yn）=Z（xnyn）-Z（xn-1yn）不等于 0时，分析与Z（xnyn）相邻上侧的网格Z（xnyn+1）的关系。
如果Z（xnyn）=Z（xnyn+1），则沉积微相值不变，即待更新沉积微相类型不更新；将待更新沉积微相类型确定为更新沉积微相类型。
如果Z（xnyn）与Z（xnyn+1）不相等，则更新。所述若更新，确定需要更新的网格的方法，包括：根据所述相邻接触关系确定所述河道与所述非河道接触的位置关系；根据所述位置关系确定需要更新的网格。
在本公开的实施例中，所述根据所述相邻接触关系确定所述河道与所述非河道接触的位置关系的方法，包括：若所述相邻接触关系大于0，所述河道在所述非河道左侧；若所述相邻接触关系小于0，所述河道在所述非河道右侧。
例如，当N（xnyn，xn-1yn）= Z（xnyn）-Z（xn-1yn）>0时，河道与非河道接触，河道在非河道左侧。当N（xnyn，xn-1yn）= Z（xnyn）-Z（xn-1yn）<0时，非河道与河道接触，河道在非河道右侧。
在本公开的实施例中，所述根据所述位置关系确定需要更新的网格的方法，包括：若所述河道在所述非河道左侧，则将相邻上侧的相邻接触关系中最大值对应的网格确定为需要更新的网格；若所述河道在所述非河道右侧，则将相邻上侧的相邻接触关系中最小值对应的网格确定为需要更新的网格。
所述对所述需要更新的网格进一步细分得到细分网格，以及确定所述细分网格内的更新沉积微相类型的方法，包括：根据沉积微相原理的预设规则，确定所述细分网格内的更新沉积微相类型。
其中，所述预设规则，包括：若河道在非河道左侧，则将则将相邻上侧的相邻接触关系中最大值对应的网格确定为需要更新的网格，对所述需要更新的网格进一步细分得到细分网格，将所述细分网格的右上角网格更新为在相同横坐标下的相邻上侧的网格值，所述细分网格的其他网格的待更新沉积微相类型保持不变，即与未进一步细分的更新的网格内的待更新沉积微相类型保持一致。若河道在非河道右侧，则将则将相邻上侧的相邻接触关系中最小值对应的网格确定为需要更新的网格，对所述需要更新的网格进一步细分得到细分网格，将所述细分网格的右上角网格更新为在相同横坐标下的相邻上侧的网格值，所述细分网格的其他网格的待更新沉积微相类型保持不变，即与未进一步细分的更新的网格内的待更新沉积微相类型保持一致。
例如，当N（xnyn，xn-1yn）= Z（xnyn）-Z（xn-1yn）>0时，河道与非河道接触，河道在非河道左侧。
如果Z（xnyn）<Z（xnyn+1），则将两个网格xnyn，xnyn+1中Z值最大的网格xnyn+1进行四等分，四等分的网格：左上角网格xnyn+1（11），右上角网格xnyn+1（12），左下角网格xnyn+1（21），右下角网格xnyn+1（22）；将Z（xnyn+1（21））=Z（xnyn）。
又如果Z（xnyn）>Z（xnyn+1），则将两个网格xnyn，xnyn+1中Z值最大的网格xnyn进行四等分，四等分的网格：左上角网格xnyn（11），右上角网格xnyn（12），左下角网格xnyn（21），右下角网格xnyn（22）；左上角网格Z（xnyn（11））=Z（xnyn+1）。
例如，当N= Z（xnyn）-Z（xn-1yn）<0时，非河道与河道接触，河道在非河道右侧。
如果Z（xn-1yn）<Z（xn-1yn+1），则将两个网格xn-1yn，xn-1yn+1中Z值最小的网格xn-1yn进行四等分，左上角网格即为xn-1yn（11），右上角网格即为xn-1yn（12），左下角网格即为xn-1yn（21），右下角网格即为xn-1yn（22），将Z（xn-1yn（12））=Z（xn-1yn+1）。
又如果Z（xn-1yn）>Z（xn-1yn+1），则将两个网格xn-1yn，xn-1yn+1中Z值最小的网格xn-1yn+1进行四等分，左上角网格即为xn-1yn（11），右上角网格即为xn-1yn（12），左下角网格即为xn-1yn（21），右下角网格即为xn-1yn（22），右下角网格Z（xn-1yn+1（22））=Z（xn-1yn）。
沉积微相的确定方法的执行主体可以是沉积微相的确定装置，例如，沉积微相的确定方法可以由终端设备或服务器或其它处理设备执行，其中，终端设备可以为用户设备（User Equipment，UE）、移动设备、用户终端、终端、蜂窝电话、无绳电话、个人数字处理（Personal Digital Assistant，PDA）、手持设备、计算设备、车载设备、可穿戴设备等。在一些可能的实现方式中，该沉积微相的确定方法可以通过处理器调用存储器中存储的计算机可读指令的方式来实现。
本领域技术人员可以理解，在具体实施方式的上述方法中，各步骤的撰写顺序并不意味着严格的执行顺序而对实施过程构成任何限定，各步骤的具体执行顺序应当以其功能和可能的内在逻辑确定。
本公开还提出了一种沉积微相的确定装置，所述沉积微相的确定装置，包括：获取单元，用于获取研究区的平面网格分布图；第一确定单元，用于根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型；第二确定单元，用于根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系；更新确定单元，用于根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型。解决了目前存在沉积微相不够准确问题，不能了满足储层精细描述高效、准确、快速的需求。本公开符合各沉积微相接触关系绘图特点，能实现准确地确定沉积微相。
图3示出根据本公开实施例的单井相成图对应的分布特征网格图。是通过井点的沉积微相类型进行储层精细刻画，识别出河道、主体、非主体、表外、尖灭几种沉积微相类型，并在图上绘制出河道的走向及趋势。
图4示出根据本公开实施例的分布特征网格图网格化及赋值得到的平面网格分布图。对河道区域内网格赋值为1，对主体区域内网格赋值为2，对非主体区域内网格赋值为3，对表外区域内网格赋值为4，对尖灭区域内网格赋值为5。
图5示出根据本公开实施例确定沉积微相后的示意图。如图5所示，能够较好地描述沉积微相的接触关系，后验井统计沉积微相预测符合率达到85%以上，可以确定河间砂体厚度分布趋势，更符合沉积特征。
在一些实施例中，本公开实施例提供的装置具有的功能或包含的模块可以用于执行上文方法实施例描述的方法，其具体实现可以参照上文方法实施例的描述，为了简洁，这里不再赘述。
本公开实施例还提出一种计算机可读存储介质，其上存储有计算机程序指令，所述计算机程序指令被处理器执行时实现上述沉积微相的确定方法。计算机可读存储介质可以是非易失性计算机可读存储介质。
本公开实施例还提出一种电子设备，包括：处理器；用于存储处理器可执行指令的存储器；其中，所述处理器被配置为上述沉积微相的确定方法。电子设备可以被提供为终端、服务器或其它形态的设备。
图6是根据一示例性实施例示出的一种电子设备800的框图。例如，电子设备800可以是移动电话，计算机，数字广播终端，消息收发设备，游戏控制台，平板设备，医疗设备，健身设备，个人数字助理等终端。
参照图6，电子设备800可以包括以下一个或多个组件：处理组件802，存储器804，电源组件806，多媒体组件808，音频组件810，输入/输出（I/ O）的接口812，传感器组件814，以及通信组件816。
处理组件802通常控制电子设备800的整体操作，诸如与显示，电话呼叫，数据通信，相机操作和记录操作相关联的操作。处理组件802可以包括一个或多个处理器820来执行指令，以完成上述的方法的全部或部分步骤。此外，处理组件802可以包括一个或多个模块，便于处理组件802和其他组件之间的交互。例如，处理组件802可以包括多媒体模块，以方便多媒体组件808和处理组件802之间的交互。
存储器804被配置为存储各种类型的数据以支持在电子设备800的操作。这些数据的示例包括用于在电子设备800上操作的任何应用程序或方法的指令，联系人数据，电话簿数据，消息，图片，视频等。存储器804可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器（SRAM），电可擦除可编程只读存储器（EEPROM），可擦除可编程只读存储器（EPROM），可编程只读存储器（PROM），只读存储器（ROM），磁存储器，快闪存储器，磁盘或光盘。
电源组件806为电子设备800的各种组件提供电力。电源组件806可以包括电源管理***，一个或多个电源，及其他与为电子设备800生成、管理和分配电力相关联的组件。
多媒体组件808包括在所述电子设备800和用户之间的提供一个输出接口的屏幕。在一些实施例中，屏幕可以包括液晶显示器（LCD）和触摸面板（TP）。如果屏幕包括触摸面板，屏幕可以被实现为触摸屏，以接收来自用户的输入信号。触摸面板包括一个或多个触摸传感器以感测触摸、滑动和触摸面板上的手势。所述触摸传感器可以不仅感测触摸或滑动动作的边界，而且还检测与所述触摸或滑动操作相关的持续时间和压力。在一些实施例中，多媒体组件808包括一个前置摄像头和/或后置摄像头。当电子设备800处于操作模式，如拍摄模式或视频模式时，前置摄像头和/或后置摄像头可以接收外部的多媒体数据。每个前置摄像头和后置摄像头可以是一个固定的光学透镜***或具有焦距和光学变焦能力。
音频组件810被配置为输出和/或输入音频信号。例如，音频组件810包括一个麦克风（MIC），当电子设备800处于操作模式，如呼叫模式、记录模式和语音识别模式时，麦克风被配置为接收外部音频信号。所接收的音频信号可以被进一步存储在存储器804或经由通信组件816发送。在一些实施例中，音频组件810还包括一个扬声器，用于输出音频信号。
I/ O接口812为处理组件802和***接口模块之间提供接口，上述***接口模块可以是键盘，点击轮，按钮等。这些按钮可包括但不限于：主页按钮、音量按钮、启动按钮和锁定按钮。
传感器组件814包括一个或多个传感器，用于为电子设备800提供各个方面的状态评估。例如，传感器组件814可以检测到电子设备800的打开/关闭状态，组件的相对定位，例如所述组件为电子设备800的显示器和小键盘，传感器组件814还可以检测电子设备800或电子设备800一个组件的位置改变，用户与电子设备800接触的存在或不存在，电子设备800方位或加速/减速和电子设备800的温度变化。传感器组件814可以包括接近传感器，被配置用来在没有任何的物理接触时检测附近物体的存在。传感器组件814还可以包括光传感器，如CMOS或CCD图像传感器，用于在成像应用中使用。在一些实施例中，该传感器组件814还可以包括加速度传感器，陀螺仪传感器，磁传感器，压力传感器或温度传感器。
通信组件816被配置为便于电子设备800和其他设备之间有线或无线方式的通信。电子设备800可以接入基于通信标准的无线网络，如WiFi，2G或3G，或它们的组合。在一个示例性实施例中，通信组件816经由广播信道接收来自外部广播管理***的广播信号或广播相关信息。在一个示例性实施例中，所述通信组件816还包括近场通信（NFC）模块，以促进短程通信。例如，在NFC模块可基于射频识别（RFID）技术，红外数据协会（IrDA）技术，超宽带（UWB）技术，蓝牙（BT）技术和其他技术来实现。
在示例性实施例中，电子设备800可以被一个或多个应用专用集成电路（ASIC）、数字信号处理器（DSP）、数字信号处理设备（DSPD）、可编程逻辑器件（PLD）、现场可编程门阵列（FPGA）、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述方法。
在示例性实施例中，还提供了一种非易失性计算机可读存储介质，例如包括计算机程序指令的存储器804，上述计算机程序指令可由电子设备800的处理器820执行以完成上述方法。
图7是根据一示例性实施例示出的一种电子设备1900的框图。例如，电子设备1900可以被提供为一服务器。参照图7，电子设备1900包括处理组件1922，其进一步包括一个或多个处理器，以及由存储器1932所代表的存储器资源，用于存储可由处理组件1922的执行的指令，例如应用程序。存储器1932中存储的应用程序可以包括一个或一个以上的每一个对应于一组指令的模块。此外，处理组件1922被配置为执行指令，以执行上述方法。
电子设备1900还可以包括一个电源组件1926被配置为执行电子设备1900的电源管理，一个有线或无线网络接口1950被配置为将电子设备1900连接到网络，和一个输入输出（I/O）接口1958。电子设备1900可以操作基于存储在存储器1932的操作***，例如Windows ServerTM，Mac OS XTM，UnixTM, LinuxTM，FreeBSDTM或类似。
在示例性实施例中，还提供了一种非易失性计算机可读存储介质，例如包括计算机程序指令的存储器1932，上述计算机程序指令可由电子设备1900的处理组件1922执行以完成上述方法。
本公开可以是***、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本公开的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子（非穷举的列表）包括：便携式计算机盘、硬盘、随机存取存储器（RAM）、只读存储器（ROM）、可擦式可编程只读存储器（EPROM或闪存）、静态随机存取存储器（SRAM）、便携式压缩盘只读存储器（CD-ROM）、数字多功能盘（DVD）、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波（例如，通过光纤电缆的光脉冲）、或者通过电线传输的电信号。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构（ISA）指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机（例如利用因特网服务提供商来通过因特网连接）。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列（FPGA）或可编程逻辑阵列（PLA），该电子电路可以执行计算机可读程序指令，从而实现本公开的各个方面。
这里参照根据本公开实施例的方法、装置（***）和计算机程序产品的流程图和/或框图描述了本公开的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。
这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。
也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。
附图中的流程图和框图显示了根据本公开的多个实施例的***、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的***来实现，或者可以用专用硬件与计算机指令的组合来实现。
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中技术的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (11)
1.一种沉积微相的确定方法，其特征在于，包括：
获取研究区的平面网格分布图；在所述获取研究区的平面网格分布图之前，确定所述平面网格分布图；其中，所述确定所述平面网格分布图的方法，包括：根据所述研究区内已有井的数据对分布特征网格图的网格进行赋值，得到所述平面网格分布图；其中，所述研究区内已有井的数据，至少包括：二类砂岩厚度、一类砂岩厚度、有效厚度的一种或几种；
根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型；
根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系；其中，所述根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系的方法，包括：根据所述积微相类型，对每个网格的所述沉积微相类型进行赋值，得到类型值；将相邻网格的所述类型值作差，得到相邻网格的多个接触关系；
根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型；其中，所述根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型的方法，包括：根据所述平面网格分布图中垂直主河道方向的相邻接触关系确定所述相邻接触关系对应的网格内待更新沉积微相类型是否可能更新；若不可能更新，则将待更新沉积微相类型确定为更新沉积微相类型；若可能更新，进一步确定是否更新；若不更新，则将待更新沉积微相类型确定为更新沉积微相类型；若更新，确定需要更新的网格，并对所述需要更新的网格进一步细分得到细分网格，以及确定所述细分网格内的更新沉积微相类型；其中，所述根据所述平面网格分布图中垂直主河道方向的相邻接触关系确定所述相邻接触关系对应的网格内待更新沉积微相类型是否可能更新的方法，包括：若所述相邻接触关系为0，此时不可能更新；若所述相邻接触关系不为0，此时可能更新。
2.根据权利要求1所述的确定方法，其特征在于，所述确定所述平面网格分布图的方法，还包括：
获取研究区；
根据所述研究区的沉积相绘制的河道砂对应的分布特征；
对所述分布特征进行网格化，得到分布特征网格图。
3.根据权利要求1-2任一项所述的确定方法，其特征在于，所述根据所述研究区内已有井的数据对所述分布特征网格图的网格进行赋值，得到所述平面网格分布图的方法，包括：
若所述网格内包含井，则根据所述井的数据对所述分布特征网格图的网格进行赋值；
若所述网格内不包含井，则根据其他网格内所述井的数据进行插值，得到插值数据；根据所述插值数据对不包含井的所述网格进行赋值。
4.根据权利要求1-2任一项所述的确定方法，其特征在于，所述根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型的方法，包括：
确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格；
基于所述河道砂网格以及/或所述非河道砂网格确定所述平面网格分布图内每个网格的待更新沉积微相类型。
5.根据权利要求3所述的确定方法，其特征在于，所述根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型的方法，包括：
确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格；
基于所述河道砂网格以及/或所述非河道砂网格确定所述平面网格分布图内每个网格的待更新沉积微相类型。
6.根据权利要求4所述的确定方法，其特征在于，所述确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格的方法，包括：
获取设定的地层厚度；
根据所述地层厚度及所述每个网格内的赋值确定河道砂最优门槛值；
根据所述河道砂最优门槛值确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格。
7.根据权利要求5所述的确定方法，其特征在于，所述确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格的方法，包括：
获取设定的地层厚度；
根据所述地层厚度及所述每个网格内的赋值确定河道砂最优门槛值；
根据所述河道砂最优门槛值确定所述平面网格分布图内的每个网格为河道砂网格或非河道砂网格。
8.根据权利要求6-7任一项所述的确定方法，其特征在于，所述根据所述地层厚度及所述每个网格内的赋值确定河道砂最优门槛值的方法，包括：
根据所述每个网格内的赋值以及所述设定的地层厚度，确定所述每个网格内的砂地比；
根据所述砂地比确定对应的河道概率点及河间概率点；
分别根据所述河道概率点及河间概率点确定对应的河道概率线及河间概率线；
基于所述的河道概率线及河间概率线确定河道砂最优门槛值。
9.一种沉积微相的确定装置，其特征在于，包括：获取单元，用于获取研究区的平面网格分布图；在所述获取研究区的平面网格分布图之前，确定所述平面网格分布图；其中，所述确定所述平面网格分布图，包括：根据所述研究区内已有井的数据对分布特征网格图的网格进行赋值，得到所述平面网格分布图；其中，所述研究区内已有井的数据，至少包括：二类砂岩厚度、一类砂岩厚度、有效厚度的一种或几种；
第一确定单元，用于根据所述平面网格分布图确定所述平面网格分布图内每个网格的待更新沉积微相类型；
第二确定单元，用于根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系；其中，所述根据每个网格的所述沉积微相类型确定相邻网格的多个接触关系，包括：根据所述积微相类型，对每个网格的所述沉积微相类型进行赋值，得到类型值；将相邻网格的所述类型值作差，得到相邻网格的多个接触关系；
更新确定单元，用于根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型；其中，所述根据所述多个接触关系对所述待更新沉积微相类型进行更新，确定所述每个网格的更新沉积微相类型，包括：根据所述平面网格分布图中垂直主河道方向的相邻接触关系确定所述相邻接触关系对应的网格内待更新沉积微相类型是否可能更新；若不可能更新，则将待更新沉积微相类型确定为更新沉积微相类型；若可能更新，进一步确定是否更新；若不更新，则将待更新沉积微相类型确定为更新沉积微相类型；若更新，确定需要更新的网格，并对所述需要更新的网格进一步细分得到细分网格，以及确定所述细分网格内的更新沉积微相类型；其中，所述根据所述平面网格分布图中垂直主河道方向的相邻接触关系确定所述相邻接触关系对应的网格内待更新沉积微相类型是否可能更新，包括：若所述相邻接触关系为0，此时不可能更新；若所述相邻接触关系不为0，此时可能更新。
10.一种电子设备，其特征在于，包括：
处理器；
用于存储处理器可执行指令的存储器；
其中，所述处理器被配置为调用所述存储器存储的指令，以执行权利要求1至8中任意一项所述的沉积微相的确定方法。
11.一种计算机可读存储介质，其上存储有计算机程序指令，其特征在于，所述计算机程序指令被处理器执行时实现权利要求1至8中任意一项所述的沉积微相的确定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011087264.7A CN112149312B (zh) | 2020-10-13 | 2020-10-13 | 沉积微相的确定方法及装置、电子设备和存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011087264.7A CN112149312B (zh) | 2020-10-13 | 2020-10-13 | 沉积微相的确定方法及装置、电子设备和存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112149312A CN112149312A (zh) | 2020-12-29 |
| CN112149312B true CN112149312B (zh) | 2024-01-30 |
Family
ID=73953128
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011087264.7A Active CN112149312B (zh) | 2020-10-13 | 2020-10-13 | 沉积微相的确定方法及装置、电子设备和存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112149312B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104809277A (zh) * | 2015-04-14 | 2015-07-29 | 中国石油天然气股份有限公司 | 一种超低渗致密储层水平井地质建模方法 |
| CN106960076A (zh) * | 2017-02-28 | 2017-07-18 | 长江大学 | 一种基于水动力学的沉积相预测方法 |
| CN107861171A (zh) * | 2017-11-08 | 2018-03-30 | 中国石油大港油田勘探开发研究院 | 一种陆相湖盆沉积微相平面分布图编制方法 |
| CN109783944A (zh) * | 2019-01-21 | 2019-05-21 | 中国石油大学(北京) | 三维流动单元模型建立方法、装置和设备 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140136171A1 (en) * | 2012-11-13 | 2014-05-15 | Chevron U.S.A. Inc. | Unstructured Grids For Modeling Reservoirs |
-
2020
- 2020-10-13 CN CN202011087264.7A patent/CN112149312B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104809277A (zh) * | 2015-04-14 | 2015-07-29 | 中国石油天然气股份有限公司 | 一种超低渗致密储层水平井地质建模方法 |
| CN106960076A (zh) * | 2017-02-28 | 2017-07-18 | 长江大学 | 一种基于水动力学的沉积相预测方法 |
| CN107861171A (zh) * | 2017-11-08 | 2018-03-30 | 中国石油大港油田勘探开发研究院 | 一种陆相湖盆沉积微相平面分布图编制方法 |
| CN109783944A (zh) * | 2019-01-21 | 2019-05-21 | 中国石油大学(北京) | 三维流动单元模型建立方法、装置和设备 |
Non-Patent Citations (4)
| Title |
|---|
| 基于ArcGIS地质统计学模块的沉积微相平面图绘制;高峰;陈恭洋;;勘探地球物理进展;第33卷(第01期);64-67 * |
| 沉积微相的定量化实现过程及研究方法;虞云岩 等;测井技术;第33卷(第4期);388-393 * |
| 油藏地质建模技巧及质量控制方法;于金彪;;新疆石油地质;第38卷(第02期);188-192 * |
| 用ArcGIS提取古沉积面恢复栅格点的方法;高明星;刘少峰;张铁强;;地理空间信息;第07卷(第04期);35-37 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112149312A (zh) | 2020-12-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113762623B (zh) | 一种滑坡方向与趋势的识别和预测方法、装置及电子设备 | |
| CN112948445B (zh) | 用于预测煤中稀土矿产资源靶区的方法及电子设备 | |
| CN113781518B (zh) | 神经网络结构搜索方法及装置、电子设备和存储介质 | |
| CN112149312B (zh) | 沉积微相的确定方法及装置、电子设备和存储介质 | |
| CN111489010B (zh) | 钻爆法隧道掌子面前方围岩级别智能预测方法及装置 | |
| CN113139484B (zh) | 人群定位方法及装置、电子设备和存储介质 | |
| CN112382346A (zh) | 含油饱和度的确定方法及装置、电子设备和存储介质 | |
| CN110456957B (zh) | 显示交互方法、装置、设备、存储介质 | |
| CN115935816B (zh) | 钻井参数确定方法、装置、设备及存储介质 | |
| CN110059260B (zh) | 一种推荐方法、装置、设备和介质 | |
| CN116796503A (zh) | 一种页岩有机物碳含量计算模型建立方法和装置 | |
| CN115495865A (zh) | 黏土矿物含量的确定方法及装置、电子设备和存储介质 | |
| CN111949924B (zh) | 渗透率的确定方法及装置、电子设备和存储介质 | |
| CN112211624A (zh) | 储层的破裂压力确定方法及装置、电子设备和存储介质 | |
| CN114329878A (zh) | 页岩气井小层划分方法、装置、设备和存储介质 | |
| CN112901134B (zh) | 距离确定、速度预测的方法及装置、电子设备和存储介质 | |
| CN114415238B (zh) | 同期异相海相碳酸岩地震解释方法及装置、设备和介质 | |
| CN115496311A (zh) | 致密油储层有利区的确定方法及装置、电子设备和介质 | |
| CN112922592B (zh) | 窜流通道、速度的预测方法及装置、电子设备和存储介质 | |
| CN114035227B (zh) | 基于随钻xrd全岩录井的变质岩潜山储层孔隙度预测方法 | |
| CN112462423A (zh) | 生物礁厚度预测方法及装置 | |
| CN113406697B (zh) | 储层的横波品质因子获取方法和装置 | |
| CN117248879A (zh) | 水平井压裂链路优化方法、装置、电子设备及存储介质 | |
| CN115704309A (zh) | 一种复杂砂砾岩储层的产能确定方法、装置及电子设备 | |
| CN118068419A (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 |