CN110007344B - 一种断溶体储层连通性的地震识别方法及装置 - Google Patents
一种断溶体储层连通性的地震识别方法及装置 Download PDFInfo
- Publication number
- CN110007344B CN110007344B CN201910276171.XA CN201910276171A CN110007344B CN 110007344 B CN110007344 B CN 110007344B CN 201910276171 A CN201910276171 A CN 201910276171A CN 110007344 B CN110007344 B CN 110007344B
- Authority
- CN
- China
- Prior art keywords
- point
- seismic
- maximum likelihood
- earthquake
- formation pressure
- 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
- 238000000034 method Methods 0.000 title claims abstract description 66
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 123
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 118
- 239000002904 solvent Substances 0.000 claims abstract description 15
- 238000005553 drilling Methods 0.000 claims description 30
- 238000001914 filtration Methods 0.000 claims description 27
- 238000004891 communication Methods 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 19
- 230000002159 abnormal effect Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000010276 construction Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 5
- 239000000243 solution Substances 0.000 description 48
- 239000012530 fluid Substances 0.000 description 11
- 230000006870 function Effects 0.000 description 11
- 239000003921 oil Substances 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 7
- 230000018109 developmental process Effects 0.000 description 7
- 238000002347 injection Methods 0.000 description 7
- 239000007924 injection Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 239000002253 acid Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000000700 radioactive tracer Substances 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000012896 Statistical algorithm Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 230000033558 biomineral tissue development Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种断溶体储层连通性的地震识别方法及装置，该方法包括：基于三维地震数据体，构建地震最大似然体；利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓；在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地层压力以及第二点处的地层压力；比较第一点处的地层压力以及第二点处的地层压力，根据比较结果判断第一点以及第二点之间的储层是否连通。本发明利用三维地震数据体构建地震最大似然体，计算同一海拔高度的两点的地层压力，基于同一缝洞单元内同一深度地层压力相等的原理，进行储层连通性识别，解决了无井区域或井外区域储层连通性的有效识别问题。
Description
技术领域
本发明涉及碳酸盐岩储层储量预测技术领域，具体涉及一种断溶体储层连通性的地震识别方法及装置。
背景技术
断溶体是受多期构造挤压、拉张作用影响，沿断裂带发育的碳酸盐岩破碎带，多期地表水沿断裂下渗或局部热液上涌使破碎带内发生溶蚀而形成的缝洞系统，简单理解为断裂控制的次生溶蚀孔洞、裂缝和断裂自身的集合体。断溶体具有强烈的非均质性，横向沿断裂带分布，具有分段性，纵向成网状分布，具有不规则性。
在实现本发明的过程中，发现断溶体油气藏的复杂性给钻探带来巨大困难。由于钻头尺寸小，直接命中主体储集体的概率非常低，往往需要采用大型酸压沟通储集体来获得工业油气流，从而大大增加了勘探开发成本。因此，储层的连通性描述在钻探过程中就显得格外重要，通过描述储集体的连通性可以寻找规模储集体、描述缝洞单元、优化钻井设计、提高钻探的命中率，为油气勘探开发节约成本。
目前，储层的连通性常规描述方法是基于已钻井来实施的，主要的方法有注水见效法、干扰试井法、井间示踪剂法、压力系统分析法、流体性质分析法和生产特征相似法等。
注水见效法是指通过注水井向目的层注水，观察相邻井组是否受到影响。当油田进入注水开发阶段时期，注水井在邻井见效是井间连通非常直观的反映。但是注水见效也有其局限性，前提是油藏本身为非无限大的水体或者油藏没有与其他无限大水体沟通。
干扰试井法是一种多井试井判断方法，可解剖一个区块，确定井间连通的地层参数。具体实施是通过激动井改变工作制度(投产、换油嘴或关井等)，在另一口或数口观察井中通过高精度压力计接受干扰压力反应，进而可研究激动井与观测井之间地层的连通性。该方法的缺点是井中下入压力计的成功率不高，甚至有些井没有及时测压，获得的压力资料不充分，并且多井干扰试验会影响油田生产任务的完成。
井间示踪剂法是指通过注入井把示踪剂注入目的地层，一旦进入地层，示踪剂便与流体一起横向运移并能在相邻井中检测到。使用化学井间示踪剂不仅可以确定储层非均质性，而且可以确定流量的大小，既可以确定注入井与生产井之间的连通性，还可以了解连通范围，但是成本相对较高。
压力系统分析法可以分为三种：①油藏压降特征法是指一个独立的缝洞单元，在开发初期处于弹性阶段，其地下累积采出量与地层压力之间呈现明显的直线关系；②压力梯度法是指同一缝洞单元其原始地层压力随深度成线性关系；③压力趋势分析法是指同一缝洞单元，各井处于同一压力系统内，由于压力可以相互传递，当油气藏处于动态平衡状态后，各井的压降趋势一致。
流体性质分析法是指同一个缝洞单元，长期运移成藏过程(当达到平衡状态的时候)可以消除流体的成分差异，流体的粘度、密度、矿化度等参数具有同步的变化趋势。反之，如果油藏是分隔的，那么油藏的成分差异不会得到消除，流体的性质也不会有同步变化的趋势。
生产特征相似法是指进行近井组、投产时间相近的井分析，在工作制度基本一致的情况下，其生产特征的变化趋势一致的井组是连通的。
上述储层连通性常规描述方法的优点在于能准确描述井间储集体的连通性、确定地层参数，缺点在于无井区域或井之外区域储集体的连通性无法评价。目前对于无井或井之外储集体连通性的识别方法，国内外还没有一个有效的方法。
发明内容
鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的一种断溶体储层连通性的地震识别方法及装置。
根据本发明的一个方面，提供了一种断溶体储层连通性的地震识别方法，该方法包括：
基于三维地震数据体，构建地震最大似然体；
利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓；
在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地震均方根速度以及第二点处的地震均方根速度，将第一点处的地震均方根速度以及第二点处的地震均方根速度分别转化为第一点处的地层压力以及第二点处的地层压力；
比较第一点处的地层压力以及第二点处的地层压力，根据比较结果判断第一点以及第二点之间的储层是否连通。
可选的，基于三维地震数据体，构建地震最大似然体具体包括：对三维地震数据体进行滤波处理，根据滤波处理后的三维地震数据体构建地震最大似然体。
可选的，利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓具体包括：在有钻井区域，利用已钻井获得的储层信息对地震最大似然体的属性进行标定，确定地震最大似然体的属性的门槛值；
和/或在无钻井区域，借鉴相邻区域的地震最大似然体的属性的门槛值，结合地质规律认识，确定本区域的地震最大似然体的属性的门槛值。
可选的，利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓进一步包括：
确定由属性的门槛值限定的门槛范围区间；
将三维地震数据体中包含的不属于门槛范围区间的数据值以及不符合地质规律认识的数据值剔除；
将三维地震数据体中剩余的数据值作为样本点，用来描述断溶体的空间轮廓。
可选的，在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地震均方根速度以及第二点处的地震均方根速度，将第一点处的地震均方根速度以及第二点处的地震均方根速度分别转化为第一点处的地层压力以及第二点处的地层压力包括：
去除地震均方根速度的异常值，对地震均方根速度进行三维空间插值和平滑滤波处理，并将地震均方根速度转化为地震层速度；
根据地震层速度，将地震层速度转化为地层压力。
可选的，比较第一点处的地层压力以及第二点处的地层压力，根据比较结果判断第一点以及第二点之间的储层是否连通进一步包括：
比较第一点处的地层压力以及第二点处的地层压力，当第一点处和第二点处的地层压力的误差属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层连通；
当第一点处和第二点处的地层压力的误差不属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层不连通。
根据本发明另一个方面，提供了一种断溶体储层连通性的地震识别装置，包括：
地震最大似然体构建模块，基于三维地震数据体，构建地震最大似然体；
断溶体空间轮廓雕刻模块，利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓；
地层压力计算模块，在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地震均方根速度以及第二点处的地震均方根速度，将第一点处的地震均方根速度以及第二点处的地震均方根速度分别转化为第一点处的地层压力以及第二点处的地层压力；
储层连通性判断模块，比较第一点处的地层压力以及第二点处的地层压力，根据比较结果判断第一点以及第二点之间的储层是否连通。
可选的，地震最大似然体构建模块适于：对三维地震数据体进行滤波处理，根据滤波处理后的三维地震数据体构建地震最大似然体。
可选的，断溶体空间轮廓雕刻模块适于：在有钻井区域，利用已钻井获得的储层信息对地震最大似然体的属性进行标定，确定地震最大似然体的属性的门槛值；
和/或在无钻井区域，借鉴相邻区域的地震最大似然体的属性的门槛值，结合地质规律认识，确定本区域的地震最大似然体的属性的门槛值。
可选的，断溶体空间轮廓雕刻模块进一步适于：
确定由属性的门槛值限定的门槛范围区间；
将三维地震数据体中包含的不属于门槛范围区间的数据值以及不符合地质规律认识的数据值剔除；
将三维地震数据体中剩余的数据值作为样本点，用来描述断溶体的空间轮廓。
可选的，地层压力计算模块适于：
去除地震均方根速度的异常值，对地震均方根速度进行三维空间插值和平滑滤波处理，并将地震均方根速度转化为地震层速度；
根据地震层速度，将地震层速度转化为地层压力。
可选的，储层连通性判断模块适于：
比较第一点处的地层压力以及第二点处的地层压力，当第一点处和第二点处的地层压力的误差属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层连通；
当第一点处和第二点处的地层压力的误差不属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层不连通。
根据本发明的再一方面，提供了一种电子设备，包括：处理器、存储器、通信接口和通信总线，处理器、存储器和通信接口通过通信总线完成相互间的通信；
存储器用于存放至少一可执行指令，可执行指令使处理器执行上述一种断溶体储层连通性的地震识别方法对应的操作。
根据本发明的再一方面，提供了一种计算机存储介质，存储介质中存储有至少一可执行指令，可执行指令使处理器执行如上述一种断溶体储层连通性的地震识别方法对应的操作。
在本发明的一种断溶体储层连通性的地震识别方法及装置中，基于三维地震数据体，构建地震最大似然体；利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓；在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地震均方根速度以及第二点处的地震均方根速度，将第一点处的地震均方根速度以及第二点处的地震均方根速度分别转化为第一点处的地层压力以及第二点处的地层压力；比较第一点处的地层压力以及第二点处的地层压力，根据比较结果判断第一点以及第二点之间的储层是否连通。本发明利用三维地震数据体构建地震最大似然体，计算同一海拔高度的两点的地层压力，基于同一缝洞单元内同一深度地层压力相等的原理，进行储层连通性识别，解决了无井区域或井外区域储层连通性的有效识别问题。通过断溶体的储层连通性识别，为寻找规模储集体、描述缝洞单元、优化钻进设计、提高钻探命中率提供基础资料，减少酸压节约勘探开发成本。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
图1示出了依据实施例一的一种断溶体储层连通性的地震识别方法的流程图；
图2示出了依据实施例二的一种断溶体储层连通性的地震识别方法的流程图；
图3a示出了顺北油田5号断裂带原始地震剖面；
图3b示出了顺北油田5号断裂带倾角导向滤波后的地震剖面；
图3c示出了顺北油田5号断裂带最大似然体与地震剖面的叠合图；
图4a为顺北油田5号断裂带断溶体立体图；
图4b为顺北油田5号断裂带断溶体俯视图；
图4c为顺北油田5号断裂带断溶体A点地层压力曲线；
图4d为顺北油田5号断裂带断溶体B点地层压力曲线；
图5示出了依据实施例三的一种断溶体储层连通性的地震识别装置的装置结构图；
图6示出了根据本发明实施例的一种电子设备的结构示意图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。
实施例一
图1示出了依据实施例一的一种断溶体储层连通性的地震识别方法的流程图。如图1所示，该方法包括如下步骤：
步骤S110：基于三维地震数据体，构建地震最大似然体。
其中，三维地震数据体具体包括：地震勘探过程中收集到的用于解释地震剖面的数据集合。
根据地震勘探过程中收集到的三维地震数据体，使用最大似然估计算法，构建地震最大似然体。其中，最大似然估计算法是建立在最大似然原理基础上的一个统计算法。最大似然估计算法提供了一种给定观察数据来评估模型参数的方法，即：“模型已定，参数未知”。具体的，样本集中的各类样本都独立分布，只考虑一类样本集D，根据样本集D来估计参数向量θ。
记已知的样本集为：
D＝{x1，x1，...，xn}
式中，D为三维地震数据集合。
概率密度函数P(D|θ)称为相对于样本集D的参数向量θ的似然函数。
如果Θ是参数空间中能使似然函数P(D|θ)最大的θ值，则Θ是θ的极大似然估计量。参数集合Θ是样本集D的函数，记为：
Θ＝d(x1，x2,...,xn)＝d(D)
Θ(x1，x2,...,xn)称作极大似然函数估计值。根据参数集合Θ构建地震最大似然体。
步骤S120：利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓。
获取地震最大似然体的属性的门槛值，具体实施时，分两种情况。
情况一：在有钻井区域，利用已钻井获得的储层信息对所述地震最大似然体的属性进行标定，确定所述地震最大似然体的属性的门槛值。
情况二：在无钻井区域，借鉴相邻区域的地震最大似然体的属性的门槛值，结合地质规律认识，确定本区域的地震最大似然体的属性的门槛值。
其中，上述属性的门槛值具体包括：属性的上限门槛值和属性的下限门槛值。根据属性的上限门槛值和属性的下限门槛值，确定由属性的门槛值限定的门槛范围区间。将三维地震数据体中包含的不属于所述门槛范围区间的数据值以及不符合地质规律认识的数据值剔除，将三维地震数据体中剩余的数据值作为样本点，用来描述断溶体的空间轮廓。
步骤S130：在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地震均方根速度以及第二点处的地震均方根速度，将第一点处的地震均方根速度以及第二点处的地震均方根速度分别转化为第一点处的地层压力以及第二点处的地层压力。
具体的，选取断溶体的空间内位于同一海拔高度的第一点以及第二点，其中，第一点以及第二点的选取除了要求海拔高度相同外，其他不作限制。
获取第一点处的地震均方根速度，采用镂空法去除第一点处的地震均方根速度的异常值，对第一点处的地震均方根速度进行三维空间差值和平滑滤波处理，并将第一点处的地震均方根速度转化为第一点处的地震层速度。根据第一点处的地震层速度，把第一点处的地震层速度转化为第一点处的地层压力。其中，地震均方根速度的异常值具体包括：地震均方根速度的数据集合中超出预设的地震均方根速度范围区间的非正常值。
获取第二点处的地震均方根速度，采用镂空法去除第二点处的地震均方根速度的异常值，对第二点处的地震均方根速度进行三维空间差值和平滑滤波处理，并将第二点处的地震均方根速度转化为第二点处的地震层速度。根据第二点处的地震层速度，把第二点处的地震层速度转化为第二点处的地层压力。
步骤S140：比较第一点处的地层压力以及第二点处的地层压力，根据比较结果判断第一点以及第二点之间的储层是否连通。
分别获取到第一点处的地层压力以及第二点处的地层压力，比较第一点处和第二点处的地层压力的误差，具体计算公式如下：
式中，δAB为第一点处和第二点处的地层压力的误差，PAB为第一点处和第二点处的地层压力的差值，PA为第一点处的地层压力，PB为第二点处的地层压力。
当第一点处和第二点处的地层压力的误差属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层连通；当第一点处和第二点处的地层压力的误差不属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层不连通。
在本发明的一种断溶体储层连通性的地震识别方法中，基于三维地震数据体，构建地震最大似然体；利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓；在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地震均方根速度以及第二点处的地震均方根速度，将第一点处的地震均方根速度以及第二点处的地震均方根速度分别转化为第一点处的地层压力以及第二点处的地层压力；比较第一点处的地层压力以及第二点处的地层压力，根据比较结果判断第一点以及第二点之间的储层是否连通。该方式利用三维地震数据体构建地震最大似然体，计算同一海拔高度的两点的地层压力，基于同一缝洞单元内同一深度地层压力相等的原理，进行储层连通性识别，解决了无井区域或井外区域储层连通性的有效识别问题。通过断溶体的储层连通性识别，为寻找规模储集体、描述缝洞单元、优化钻进设计、提高钻探命中率提供基础资料，减少酸压节约勘探开发成本。
实施例二
图2示出了依据实施例二的一种断溶体储层连通性的地震识别方法的流程图。如图2所示，该方法包括如下步骤：
步骤S210：基于三维地震数据体，构建地震最大似然体。
其中，三维地震数据体具体包括：地震勘探过程中收集到的用于解释地震剖面的数据集合。
根据地震勘探过程中收集到的三维地震数据体，使用最大似然估计算法，构建地震最大似然体。其中，最大似然估计算法是建立在最大似然原理基础上的一个统计算法。最大似然估计算法提供了一种给定观察数据来评估模型参数的方法，即：“模型已定，参数未知”。具体的，样本集中的各类样本都独立分布，只考虑一类样本集D，根据样本集D来估计参数向量θ。
记已知的样本集为：
D＝{x1，x1，...，xn}
式中，D为三维地震数据集合。
概率密度函数P(D|θ)称为相对于样本集D的参数向量θ的似然函数。
如果Θ是参数空间中能使似然函数P(D|θ)最大的θ值，则Θ是θ的极大似然估计量。参数集合Θ是样本集D的函数，记为：
Θ＝d(x1，x2,...,xn)＝d(D)
Θ(x1，x2,...,xn)称作极大似然函数估计值。根据参数集合Θ构建地震最大似然体。
可选的，为了提高地震资料的品质，更易于断溶体特征识别，对三维地震数据体进行滤波处理，根据滤波处理后的三维地震数据体构建地震最大似然体。具体的，倾角导向滤波处理是沿地震同相轴的方向去除随机噪音，以达到增加同相轴连续性的目的。倾角导向滤波的关键是在去除叠置在反射层上的噪声的同时区分反射层的倾角和方位角，估算反射层的倾角和方位角后，采用一个滤波器增强沿反射层的信号，其中最重要的环节就是倾角和方位角的估算。具体实施时，
首先计算瞬时频率：
式中，Φ代表瞬时相位；μ(t,x,y)代表输入数据，μH(t，x，y)代表对应时间t的希尔伯特变换；tan-1表示反正切函数，输出的变化范围是(-π,π)。
其次计算瞬时波数kx和ky：
式中，kx、ky分别为x、y方向的瞬时波数，Φ代表瞬时相位；μ(t,x,y)代表输入数据，μH(t，x，y)代表对应时间t的希尔伯特变换。
其次计算瞬时时间倾角：
式中，θx、θy分别为x、y方向的倾角，kx、ky分别为x、y方向的瞬时波数。
最后计算地层倾角和方位角：
式中，θ为地层倾角，φ为方位角。
进一步，为了消除随机噪音，保证后续地层压力计算的精确，对地震最大似然体的数据集合进行平滑滤波处理。
步骤S220：获取地震最大似然体的属性的门槛值。
获取地震最大似然体的属性的门槛值，具体实施时，分两种情况。
情况一：在有钻井区域，利用已钻井获得的储层信息对所述地震最大似然体的属性进行标定，确定所述地震最大似然体的属性的门槛值。具体实施时，根据已钻井获得的储层深度、储层厚度、储层空间结构、储层物性等对地震最大似然体的属性进行井震标定，获得地震最大似然体标定后的属性的门槛值。
情况二：在无钻井区域，借鉴相邻区域的地震最大似然体的属性的门槛值，结合地质规律认识，确定本区域的地震最大似然体的属性的门槛值。具体实施时，借助相邻工区内井震标定的地震最大似然体的属性的门槛值，通过标准化数学手段，把本工区的地震最大似然体的属性的门槛值归一化到相邻工区的地震最大似然体的属性的值域范围。
其中，上述属性的门槛值具体包括：属性的上限门槛值和属性的下限门槛值。
步骤S230：利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓。
获取地震最大似然体的属性的门槛值，其中，属性的门槛值具体包括：属性的上限门槛值和属性的下限门槛值。根据属性的上限门槛值和属性的下限门槛值，确定由属性的门槛值限定的门槛范围区间。将三维地震数据体中包含的不属于所述门槛范围区间的数据值以及不符合地质规律认识的数据值剔除，将三维地震数据体中剩余的数据值作为样本点，用来描述断溶体的空间轮廓。
步骤S240：在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地震均方根速度以及第二点处的地震均方根速度，将第一点处的地震均方根速度以及第二点处的地震均方根速度分别转化为第一点处的地层压力以及第二点处的地层压力。
具体的，选取断溶体的空间内位于同一海拔高度的第一点以及第二点，其中，第一点以及第二点的选取除了要求海拔高度相同外，其他不作限制。
获取第一点处的地震均方根速度，采用镂空法去除第一点处的地震均方根速度的异常值，对第一点处的地震均方根速度进行三维空间差值和平滑滤波处理，并将第一点处的地震均方根速度转化为第一点处的地震层速度。根据第一点处的地震层速度，把第一点处的地震层速度转化为第一点处的地层压力。其中，地震均方根速度的异常值具体包括：地震均方根速度的数据集合中超出预设的地震均方根速度范围区间的非正常值。
获取第二点处的地震均方根速度，采用镂空法去除第二点处的地震均方根速度的异常值，对第二点处的地震均方根速度进行三维空间差值和平滑滤波处理，并将第二点处的地震均方根速度转化为第二点处的地震层速度。根据第二点处的地震层速度，把第二点处的地震层速度转化为第二点处的地层压力。
将地震均方根速度转化为地震层速度的具体公式如下：
式中，Vi为预测层段的地震层速度，Vri、Vr(i-1)分别为预测层段与其上覆层段的地震均方根速度，Ti、T(i-1)分别为预测层段与其上覆地层的地震反射时间。
利用Fillippone法将地震层速度转化为地层压力的具体公式如下：
式中，Pf为地层压力，Vmax、Vmin分别为孔隙率接近于零和刚性接近于零时的地震波速度，Vmax近似于基质速度，Vmin近似于孔隙流体速度；Vi为预测层段的地震层速度；Pov为上覆地层的压力。其中，
Vmax＝1.4Vr0+3KT
Vmin＝0.7Vr0+0.5KT
式中，Vr0为T0时间为零时的均方根速度，K为变化率，T为预测层段的地震反射时间，G为单位换算系数，G＝9.8*10-3，ρ为上覆地层平均密度，单位为g/cm3，H为地层埋深，单位为m。其中，
式中，Vr为预测层段地震均方根速度，T为预测层段地震反射时间，Vr0为T0时间为零时的地震均方根速度，T0为常数0。
进一步，为了确保后续地层压力计算精确，提取第一点处的地震层速度体的平面属性和第二点处的地震层速度体的平面属性，根据地震层速度体的平面属性判断地震层速度的合理性。若地震层速度体的平面属性符合地质认识及规律，则地震层速度可用于后续地层压力计算；若地震层速度体的平面属性不符合地质认识及规律，则重新对地震层速度体进行优化处理直到满足后续计算要求。其中，地震层速度体具体包括：描述地震层速度的地震数据集合；地震层速度体的平面属性具体包括：根据描述地震层速度的地震数据，经数学方法处理后获得的有关地震层速度的属性特征。
步骤S250：比较所述第一点处的地层压力以及所述第二点处的地层压力，根据比较结果判断所述第一点以及第二点之间的储层是否连通。
分别获取到第一点处的地层压力以及第二点处的地层压力，比较第一点处和第二点处的地层压力的误差，具体计算公式如下：
式中，δAB为第一点处和第二点处的地层压力的误差，PAB为第一点处和第二点处的地层压力的差值，PA为第一点处的地层压力，PB为第二点处的地层压力。
当第一点处和第二点处的地层压力的误差属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层连通；当第一点处和第二点处的地层压力的误差不属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层不连通。其中，预设的地层压力差误差范围参照石油地震勘探技术规范，预设的地层压力差误差取值5％。具体实施时，若第一点处和第二点处的地层压力的误差δAB∈[-5％,5％]，则第一点和第二点两点之间的储层连通；若第一点处和第二点处的地层压力的误差
则第一点和第二点两点之间的储层不连通。
为了更透彻理解本发明，以顺北油田5号断裂带为例，对本发明提供的一种断溶体储层连通性识别方法进一步地进行说明。
图3a示出了顺北油田5号断裂带原始地震剖面；图3b示出了顺北油田5号断裂带倾角导向滤波后的地震剖面；图3c示出了顺北油田5号断裂带最大似然体与地震剖面的叠合图。对比图3a和图3b，可以清晰看出经过倾角导向滤波后地震剖面信噪比提高、断溶体边界更清晰，为后续地震最大似然体计算奠定良好基础。图3c为地震最大似然体与地震剖面的叠合图，地震最大似然体详细刻画了断溶体边界，展示了断溶体纵向结构特征。
图4a为顺北油田5号断裂带断溶体立体图；图4b为顺北油田5号断裂带断溶体俯视图；图4c为顺北油田5号断裂带断溶体A点地层压力曲线；图4d为顺北油田5号断裂带断溶体B点地层压力曲线。从图4a和图4b分析得到断溶体轮廓和空间展布符合地质认识及规律，可以作为断溶体单元进行油藏描述和储层连通性分析。图4b展示的是四组断溶体单元，在第二组断溶体单元内选择任意相同海拔高度的两点A和B，根据地震均方根速度计算A和B两点的地层压力曲线，比较A和B两点的地层压力，通过A和B两点的地层压力的相对误差判断储层的连通性。
表1为A和B两点的目的层段的地层压力及相对误差。从图4c和图4d可以看出，A和B两点地层压力趋势一致，符合该地区地层压实规律。从表1可以看出，目的层段A和B两点压力误差绝对值小于2％，包含于系统误差内δAB∈[-5％,5％]，因此判断断溶体内A和B两点间储层是连通的。
| 深度(m) | A点压力(Mpa) | B点压力(Mpa) | A-B压力差(Mpa) | 相对压力误差 |
| 7010 | 80.006 | 80.953 | -0.948 | -1.171％ |
| 7060 | 80.695 | 81.650 | -0.954 | -1.169％ |
| 7115 | 81.386 | 82.348 | -0.962 | -1.168％ |
| 7174 | 82.079 | 83.048 | -0.969 | -1.167％ |
| 7234 | 82.773 | 83.751 | -0.978 | -1.168％ |
| 7293 | 83.469 | 84.456 | -0.987 | -1.169％ |
| 7353 | 84.170 | 85.165 | -0.995 | -1.168％ |
| 7412 | 84.874 | 85.877 | -1.003 | -1.168％ |
| 7471 | 85.581 | 86.591 | -1.010 | -1.166％ |
| 7529 | 86.292 | 87.309 | -1.017 | -1.165％ |
| 7588 | 87.007 | 88.030 | -1.024 | -1.163％ |
| 7646 | 87.724 | 88.754 | -1.030 | -1.161％ |
| 7704 | 88.446 | 89.481 | -1.036 | -1.158％ |
| 7762 | 89.170 | 90.199 | -1.029 | -1.141％ |
| 7820 | 89.887 | 90.995 | -1.108 | -1.218％ |
| 7877 | 90.649 | 91.931 | -1.282 | -1.394％ |
| 7935 | 91.585 | 92.874 | -1.289 | -1.388％ |
| 7992 | 92.532 | 93.824 | -1.292 | -1.377％ |
表1
在本发明的一种断溶体储层连通性的地震识别方法及装置中，基于三维地震数据体，构建地震最大似然体；利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓；在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地震均方根速度以及第二点处的地震均方根速度，将第一点处的地震均方根速度以及第二点处的地震均方根速度分别转化为第一点处的地层压力以及第二点处的地层压力；比较第一点处的地层压力以及第二点处的地层压力，根据比较结果判断第一点以及第二点之间的储层是否连通。该方式利用三维地震数据体构建地震最大似然体，计算同一海拔高度的两点的地层压力，基于同一缝洞单元内同一深度地层压力相等的原理，进行储层连通性识别，解决了无井区域或井外区域储层连通性的有效识别问题。在构建地震最大似然体之前通过滤波处理，消除随机噪音，提高地震资料品质，更易于断溶体特征识别。通过断溶体的储层连通性识别，为寻找规模储集体、描述缝洞单元、优化钻进设计、提高钻探命中率提供基础资料，减少酸压节约勘探开发成本。
实施例三
图5示出了依据实施例三的一种断溶体储层连通性的地震识别装置的装置结构图。如图5所示，本装置包括：
地震最大似然体构建模块51，基于三维地震数据体，构建地震最大似然体；
断溶体空间轮廓雕刻模块52，利用地震最大似然体的属性的门槛值对地震最大似然体进行空间雕刻，获得断溶体的空间轮廓；
地层压力计算模块53，在断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取第一点处的地震均方根速度以及第二点处的地震均方根速度，将第一点处的地震均方根速度以及第二点处的地震均方根速度分别转化为第一点处的地层压力以及第二点处的地层压力；
储层连通性判断模块54，比较第一点处的地层压力以及第二点处的地层压力，根据比较结果判断第一点以及第二点之间的储层是否连通。
可选的，地震最大似然体构建模块51适于：对三维地震数据体进行滤波处理，根据滤波处理后的三维地震数据体构建地震最大似然体。
可选的，断溶体空间轮廓雕刻模块52适于：在有钻井区域，利用已钻井获得的储层信息对地震最大似然体的属性进行标定，确定地震最大似然体的属性的门槛值；
和/或在无钻井区域，借鉴相邻区域的地震最大似然体的属性的门槛值，结合地质规律认识，确定本区域的地震最大似然体的属性的门槛值。
可选的，断溶体空间轮廓雕刻模块52进一步适于：
确定由属性的门槛值限定的门槛范围区间；
将三维地震数据体中包含的不属于门槛范围区间的数据值以及不符合地质规律认识的数据值剔除；
将三维地震数据体中剩余的数据值作为样本点，用来描述断溶体的空间轮廓。
可选的，地层压力计算模块53适于：
去除地震均方根速度的异常值，对地震均方根速度进行三维空间插值和平滑滤波处理，并将地震均方根速度转化为地震层速度；
根据地震层速度，将地震层速度转化为地层压力。
可选的，储层连通性判断模块54适于：
比较第一点处的地层压力以及第二点处的地层压力，当第一点处和第二点处的地层压力的误差属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层连通；
当第一点处和第二点处的地层压力的误差不属于预设的地层压力差误差范围时，判断第一点和第二点两点之间的储层不连通。
上述各个模块的具体结构和工作原理可参照方法实施例中相应步骤的描述，此处不再赘述。
本申请实施例提供了一种非易失性计算机存储介质，计算机存储介质存储有至少一可执行指令，该计算机可执行指令可执行上述任意方法实施例中的基于多级网络节点的故障定位方法。
图6示出了根据本发明实施例的一种电子设备的结构示意图，本发明具体实施例并不对电子设备的具体实现做限定。
如图6所示，该电子设备可以包括：处理器(processor)602、通信接口(Communications Interface)604、存储器(memory)606、以及通信总线608。
其中：
处理器602、通信接口604、以及存储器606通过通信总线608完成相互间的通信。
通信接口604，用于与其它设备比如客户端或其它服务器等的网元通信。
处理器602，用于执行程序610，具体可以执行上述基于多级网络节点的故障定位方法实施例中的相关步骤。
具体地，程序610可以包括程序代码，该程序代码包括计算机操作指令。
处理器602可能是中央处理器CPU，或者是特定集成电路ASIC(ApplicationSpecific Integrated Circuit)，或者是被配置成实施本发明实施例的一个或多个集成电路。电子设备包括的一个或多个处理器，可以是同一类型的处理器，如一个或多个CPU；也可以是不同类型的处理器，如一个或多个CPU以及一个或多个ASIC。
存储器606，用于存放程序610。存储器606可能包含高速RAM存储器，也可能还包括非易失性存储器(non-volatile memory)，例如至少一个磁盘存储器。
程序610具体可以用于使得处理器602执行上述方法实施例中的各项操作。
在此提供的算法和显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述，构造这类系统所要求的结构是显而易见的。此外，本发明也不针对任何特定编程语言。应当明白，可以利用各种编程语言实现在此描述的本发明的内容，并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
在此处所提供的说明书中，说明了大量具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。
类似地，应当理解，为了精简本公开并帮助理解各个发明方面中的一个或多个，在上面对本发明的示例性实施例的描述中，本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而，并不应将该公开的方法解释成反映如下意图：即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说，如下面的权利要求书所反映的那样，发明方面在于少于前面公开的单个实施例的所有特征。因此，遵循具体实施方式的权利要求书由此明确地并入该具体实施方式，其中每个权利要求本身都作为本发明的单独实施例。
Claims (9)
1.一种断溶体储层连通性的地震识别方法，其特征在于，包括：
基于三维地震数据体，构建地震最大似然体；
利用所述地震最大似然体的属性的门槛值对所述地震最大似然体进行空间雕刻，获得断溶体的空间轮廓；
在所述断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取所述第一点处的地震均方根速度以及所述第二点处的地震均方根速度，将所述第一点处的地震均方根速度以及所述第二点处的地震均方根速度分别转化为所述第一点处的地层压力以及所述第二点处的地层压力；
比较所述第一点处的地层压力以及所述第二点处的地层压力，根据比较结果判断所述第一点以及第二点之间的储层是否连通；
其中，所述利用所述地震最大似然体的属性的门槛值对所述地震最大似然体进行空间雕刻，获得断溶体的空间轮廓具体包括：在有钻井区域，利用已钻井获得的储层信息对所述地震最大似然体的属性进行标定，确定所述地震最大似然体的属性的门槛值；
在无钻井区域，借鉴相邻区域的地震最大似然体的属性的门槛值，结合地质规律认识，确定本区域的所述地震最大似然体的属性的门槛值。
2.根据权利要求1所述的方法，其特征在于，所述基于三维地震数据体，构建地震最大似然体具体包括：对所述三维地震数据体进行滤波处理，根据滤波处理后的三维地震数据体构建地震最大似然体。
3.根据权利要求1所述的方法，其特征在于，所述利用所述地震最大似然体的属性的门槛值对所述地震最大似然体进行空间雕刻，获得断溶体的空间轮廓进一步包括：
确定由所述属性的门槛值限定的门槛范围区间；
将所述三维地震数据体中包含的不属于所述门槛范围区间的数据值以及不符合地质规律认识的数据值剔除；
将三维地震数据体中剩余的数据值作为样本点，用来描述断溶体的空间轮廓。
4.根据权利要求1所述的方法，其特征在于，所述在所述断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取所述第一点处的地震均方根速度以及所述第二点处的地震均方根速度，将所述第一点处的地震均方根速度以及所述第二点处的地震均方根速度分别转化为所述第一点处的地层压力以及所述第二点处的地层压力包括：
去除地震均方根速度的异常值，对所述地震均方根速度进行三维空间插值和平滑滤波处理，并将所述地震均方根速度转化为地震层速度；
根据所述地震层速度，将所述地震层速度转化为地层压力。
5.根据权利要求1所述的方法，其特征在于，所述比较所述第一点处的地层压力以及所述第二点处的地层压力，根据比较结果判断所述第一点以及第二点之间的储层是否连通进一步包括：
比较所述第一点处的地层压力以及所述第二点处的地层压力，当所述第一点处和第二点处的地层压力的误差属于预设的地层压力差误差范围时，判断所述第一点和第二点两点之间的储层连通；
当所述第一点处和第二点处的地层压力的误差不属于预设的地层压力差误差范围时，判断所述第一点和第二点两点之间的储层不连通。
6.一种断溶体储层连通性的地震识别装置，其特征在于，包括：
地震最大似然体构建模块，基于三维地震数据体，构建地震最大似然体；
断溶体空间轮廓雕刻模块，利用所述地震最大似然体的属性的门槛值对所述地震最大似然体进行空间雕刻，获得断溶体的空间轮廓；
地层压力计算模块，在所述断溶体的空间内选择位于同一海拔高度的第一点以及第二点，获取所述第一点处的地震均方根速度以及所述第二点处的地震均方根速度，将所述第一点处的地震均方根速度以及所述第二点处的地震均方根速度分别转化为所述第一点处的地层压力以及所述第二点处的地层压力；
储层连通性判断模块，比较所述第一点处的地层压力以及所述第二点处的地层压力，根据比较结果判断所述第一点以及第二点之间的储层是否连通；
其中，所述断溶体空间轮廓雕刻模块具体适于：在有钻井区域，利用已钻井获得的储层信息对所述地震最大似然体的属性进行标定，确定所述地震最大似然体的属性的门槛值；
在无钻井区域，借鉴相邻区域的地震最大似然体的属性的门槛值，结合地质规律认识，确定本区域的所述地震最大似然体的属性的门槛值。
7.根据权利要求6所述的装置，其特征在于，所述地震最大似然体构建模块适于：对所述三维地震数据体进行滤波处理，根据滤波处理后的三维地震数据体构建地震最大似然体。
8.一种电子设备，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；
所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如权利要求1-5中任一项所述的一种断溶体储层连通性的地震识别方法对应的操作。
9.一种计算机存储介质，所述存储介质中存储有至少一可执行指令，所述可执行指令使处理器执行如权利要求1-5中任一项所述的一种断溶体储层连通性的地震识别方法对应的操作。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910276171.XA CN110007344B (zh) | 2019-04-08 | 2019-04-08 | 一种断溶体储层连通性的地震识别方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910276171.XA CN110007344B (zh) | 2019-04-08 | 2019-04-08 | 一种断溶体储层连通性的地震识别方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110007344A CN110007344A (zh) | 2019-07-12 |
| CN110007344B true CN110007344B (zh) | 2021-02-09 |
Family
ID=67170234
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910276171.XA Active CN110007344B (zh) | 2019-04-08 | 2019-04-08 | 一种断溶体储层连通性的地震识别方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110007344B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113313610B (zh) * | 2020-02-26 | 2024-02-13 | 中国石油化工股份有限公司 | 一种断溶体油藏注水井流道调整方法 |
| CN113495295B (zh) * | 2020-04-07 | 2024-01-30 | 中国石油天然气股份有限公司 | 弱连通缝洞型储层的识别方法、装置、设备及存储介质 |
| CN112394416B (zh) * | 2020-10-20 | 2024-06-25 | 中国石油天然气集团有限公司 | 非均质断裂控制储层预测方法及装置 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| BRPI0611500A2 (pt) * | 2005-05-26 | 2011-02-22 | Exxonmobil Upstream Res Co | método para analisar a qualidade conectada de um reservatório de hidrocarbonetos |
| CN103902827B (zh) * | 2014-04-03 | 2017-04-12 | 刘鑫 | 一种碳酸盐岩水平井流动单元划分方法 |
| US10132147B2 (en) * | 2014-07-02 | 2018-11-20 | Weatherford Technology Holdings, Llc | System and method for modeling and design of pulse fracturing networks |
| CN106353806B (zh) * | 2015-07-17 | 2018-04-10 | 中国石油化工股份有限公司 | 一种检测生物礁储层连通性的方法 |
| CN107515425A (zh) * | 2016-06-16 | 2017-12-26 | 中国石油化工股份有限公司 | 适用于陡坡带砂砾岩沉积连通体的地震预测方法 |
| CN107083939A (zh) * | 2017-06-26 | 2017-08-22 | 中国石油化工股份有限公司 | 用于碳酸盐岩断溶体油藏的立体开发方法 |
| CN107272065A (zh) * | 2017-08-04 | 2017-10-20 | 中国石油化工股份有限公司 | 碳酸盐岩断溶体轮廓检测方法 |
-
2019
- 2019-04-08 CN CN201910276171.XA patent/CN110007344B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110007344A (zh) | 2019-07-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110007344B (zh) | 一种断溶体储层连通性的地震识别方法及装置 | |
| CN108301825B (zh) | 一种高温高压储层孔隙压力分析方法 | |
| US10495782B2 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
| US20100191514A1 (en) | Method, Program and Computer System for Scaling Hydrocarbon Reservoir Model Data | |
| CN107895092B (zh) | 一种基于复杂非线性注采建模的井间连通定量评价方法 | |
| CN109798108B (zh) | 深层低渗-致密储层复杂流体性质的识别方法和装置 | |
| CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
| CN113175322B (zh) | 地层漏失压力剖面的建立方法、计算机装置及存储介质 | |
| CN113486539B (zh) | 一种深水钻井全过程中地层压力不确定性分析方法 | |
| CN113361111A (zh) | 一种低渗透油藏储层特征模型的划分方法 | |
| Saafan et al. | Inferring capillary pressure curve from 2D rock images based on fractal theory in low-permeability sandstone: a new integrated approach | |
| CN110007345B (zh) | 一种断溶体油气藏洞穴型储层储量计算方法及装置 | |
| CN110320562B (zh) | 页岩气储层水平井中声波时差的校正方法 | |
| CN108457646B (zh) | 确定储层流体性质的方法 | |
| CN115577616A (zh) | 基于深度学习的碳酸岩缝洞地震刻画方法及装置 | |
| CN112394416A (zh) | 非均质断裂控制储层预测方法及装置 | |
| CN113803054B (zh) | 油水界面深度确定方法及防止生产井水窜的预警方法 | |
| CN113822489B (zh) | 浅水流钻井风险的确定方法、装置和服务器 | |
| WO2014092713A2 (en) | System, method and computer program product for evaluating and ranking geobodies using a euler characteristic | |
| CN114139386B (zh) | 一种基于多层合采的连通性评估方法及装置 | |
| CN118346267A (zh) | 基于随钻电阻率的地层压力预测方法、装置、设备及介质 | |
| CN118349801A (zh) | 一种岩溶作用影响深度的确定方法、装置、设备及介质 | |
| CN114645709A (zh) | 一种聚合物驱地层系数的试井方法和装置 | |
| CN118363077A (zh) | 天然气水合物储层类型现场识别方法、装置及设备 | |
| CN114776277A (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 |