CN104850732B - 一种基于砂体统计学的油藏小层划分方法及装置 - Google Patents
一种基于砂体统计学的油藏小层划分方法及装置 Download PDFInfo
- Publication number
- CN104850732B CN104850732B CN201510162800.8A CN201510162800A CN104850732B CN 104850732 B CN104850732 B CN 104850732B CN 201510162800 A CN201510162800 A CN 201510162800A CN 104850732 B CN104850732 B CN 104850732B
- Authority
- CN
- China
- Prior art keywords
- sand body
- sand
- well
- mrow
- msub
- 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.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及地质勘探技术领域，尤其涉及一种基于砂体统计学的油藏小层划分方法及装置。该小层划分方法包括：基于岩心资料及测井曲线确定隔夹层划分基准，划分出单井目的层段内的砂体；运用统计学方法统计油田所有单井中砂体的砂体厚度和砂体个数；当砂体厚度频率曲线及砂体个数频率曲线均呈现单一峰态，且偏度和峰度在预期范围内时，选取砂体个数频率曲线的峰值所对应的砂体个数作为小层数目，并选取砂体个数与小层数目相等的井，结合旋回特征建立标准井；将砂体个数多余标准井砂体个数的井，进行砂体合并；将砂体个数少于标准井砂体个数的井进行砂体劈分。本发明的油藏小层划分方法及装置，能够实现单砂体级别高分辨率层序地层的划分与对比。
Description
技术领域
本发明涉及地质勘探技术领域，具体而言，涉及一种基于砂体统计学的油藏小层划分方法及装置。
背景技术
自80年代起，我国开始引进和开展油藏描述，其在油田开发实践中的应用成效显著。目前，油藏描述与剩余油分布规律研究的结合，作为精细油藏描述的一项重要基础工作，已成为确定挖掘潜力和提高采收率等措施的重要途径。
在油藏描述与剩余油分布规律结合研究的过程中，小层的合理划分和对比是油藏描述中由单井评价到多井评价的重要程序，是描述储层形态和参数空间分布的基础。小层的划分与对比相铺相成，合理的层组划分是正确对比的基础，只有通过反复对比，才能在一定范围实现统一的分层。
伴随着我国石油工业的发展，国内地学工作者先后提出“旋回对比、分级控制”、“时间单元划分、等高程”对比、“等高程切片对比”，然而对于构造复杂和沉积复杂的油藏，由于相变频繁并缺少稳定标志层，导致小层的划分的对比困难重重，由此，发展复杂油藏条件下的小层对比技术意义重大。
发明内容
本发明的目的在于提供一种基于砂体统计学的油藏小层划分方法及装置，以实现单砂体级别高分辨率层序地层的划分与对比。
第一方面，本发明实施例提供了一种基于砂体统计学的油藏小层划分方法，包括：基于岩心资料及测井曲线确定隔夹层划分基准，根据所述隔夹层划分基准划分出单井目的层段内的砂体；运用统计学方法统计油田所有单井中划分出的所述砂体的砂体厚度和砂体个数，并根据统计的结果绘制砂体厚度频率曲线和砂体个数频率曲线；当所述砂体厚度频率曲线及所述砂体个数频率曲线均呈现单一峰态，且偏度和峰度在预期范围内时，选取砂体个数频率曲线的峰值所对应的砂体个数作为小层数目，并选取砂体个数与所述小层数目相等的井，结合旋回特征建立标准井；将砂体个数多于所述标准井砂体个数的井，进行砂体合并；将砂体个数少于所述标准井砂体个数的井以等厚原则进行砂体劈分。
结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，所述基于岩心资料及测井曲线确定隔夹层划分基准，根据所述隔夹层划分基准划分出单井目的层段内的砂体，包括：根据岩心资料分析确定隔夹层划分标准；将所述岩心资料分析的结论与测井曲线进行关联，并基于所述隔夹层划分标准建立测井曲线识别隔夹层的定量标准；依据所述定量标准划分出单井目的层段内的砂体。
结合第一方面，本发明实施例提供了第一方面的第二种可能的实施方式，其中，所述根据岩心资料分析确定隔夹层划分标准，包括：根据岩心资料进行砂质泥岩界面分析，泥岩、砂质泥岩及钙质砂层的厚度和孔隙度分析；根据所述岩心资料分析的结果确定隔层划分标准及夹层划分标准，其中所述隔层划分标准为：泥岩厚度大于0.4米，砂质泥岩厚度大于0.6米，钙质砂层厚度大于1米且孔隙度小于8md的层划分为隔层；所述夹层划分标准为：泥岩厚度大于0.2米，砂质泥岩厚度大于0.3米，钙质砂层厚度大于0.5米且孔隙度小于8md的层划分为夹层。
结合第一方面，本发明实施例提供了第一方面的第三种可能的实施方式，其中，所述将所述岩心资料分析的结论与测井曲线进行关联，包括：岩心归位；对比相同深度段内的岩心资料、录井资料、成像测井及测井资料的特征；根据所述对比的结果，建立岩性、岩性界面与测井曲线上自然伽马GR值的对应关系。
结合第一方面，本发明实施例提供了第一方面的第四种可能的实施方式，其中，该方法还包括：根据所述岩心资料分析的结论与测井曲线关联的结果，利用公式h＝0.125×N计算砂体的厚度；式中：0.125为测井曲线中测量点间距为0.125m；N为测井曲线自然伽马GR值中小于砂泥分界值GRlimit的个数；
式中：为测井序列GR平均值；MEgr为测井序列GR中位数；GEOgr为测井序列GR平均值；MADgr为测井序列GR平均绝对值偏差。
结合第一方面，本发明实施例提供了第一方面的第五种可能的实施方式，其中，所述将砂体个数多于所述标准井砂体个数的井，进行砂体合并，包括：选取平行物源方向和垂直物源方向的过标准井的连井剖面；基于所述连井剖面，在预设标志层限定的层段内，将砂体个数多于标准井砂体个数的井中的多个砂体统一编号为相同小层砂体。
结合第一方面，本发明实施例提供了第一方面的第六种可能的实施方式，其中，所述将砂体个数少于所述标准井砂体个数的井以等厚原则进行砂体劈分，包括：选择测井曲线中所显示的隔层位置或夹层位置作为劈分界面进行劈分，所述劈分后得到的各个小层之间的厚度差值符合预设差值范围。
第二方面，本发明实施例还提供了一种基于砂体统计学的油藏小层划分装置，包括：砂体划分模块，用于基于岩心资料及测井曲线确定隔夹层划分基准，根据所述隔夹层划分基准划分出单井目的层段内的砂体；统计模块，用于运用统计学方法统计油田所有单井中划分出的所述砂体的砂体厚度和砂体个数，并根据统计的结果绘制砂体厚度频率曲线和砂体个数频率曲线；标准井建立模块，用于当所述砂体厚度频率曲线及所述砂体个数频率曲线均呈现单一峰态，且偏度和峰度在预期范围内时，选取砂体个数频率曲线的峰值所对应的砂体个数作为小层数目，并选取砂体个数与所述小层数目相等的井，结合旋回特征建立标准井；合并与劈分模块，用于将砂体个数多于所述标准井砂体个数的井，进行砂体合并；将砂体个数少于所述标准井砂体个数的井以等厚原则进行砂体劈分。
结合第二方面，本发明实施例提供了第二方面的第一种可能的实施方式，其中，所述砂体划分模块，包括：划分标准确定单元，用于根据岩心资料分析确定隔夹层划分标准；关联单元，用于将所述岩心资料分析的结论与测井曲线进行关联，并基于所述隔夹层划分标准建立测井曲线识别隔夹层的定量标准；砂体划分单元，用于依据所述定量标准划分出单井目的层段内的砂体。
结合第二方面，本发明实施例提供了第二方面的第二种可能的实施方式，其中，该装置中还包括：砂体厚度计算模块，用于根据所述岩心资料分析的结论与测井曲线关联的结果，利用公式h＝0.125×N计算砂体的厚度；式中：0.125为测井曲线中测量点间距为0.125m；N为测井曲线自然伽马GR值中小于砂泥分界值GRlimit的个数；
式中：为测井序列GR平均值；MEgr为测井序列GR中位数；GEOgr为测井序列GR平均值；MADgr为测井序列GR平均绝对值偏差。
本发明实施例的油藏小层划分方法及装置，将砂体统计学特征与高分辨率层序地层学理论相结合，根据砂体厚度频率曲线和砂体个数频率曲线特征，确定小层数目，并选取砂体个数与小层数据相等的井结合旋回特征建立标准井，标准井建立后，根据各个井的砂体数目进行砂体合并及砂体劈分，从而实现目的层段内小层的划分与对比，克服了相关技术中基于基准面旋回分析难以确定目的层段内小层划分方案，以及难以对单砂体油水关系进行准确解释的缺点。本发明的基于砂体统计学的油藏小层划分方法有助于准确实现单砂体级别的小层划分方案的确定，有助于实现单砂体联通关系的分析，从而提高油藏开发效果。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1示出了本发明实施例中基于砂体统计学的油藏小层划分方法的流程图；
图2示出了本发明实施例中油藏小层划分技术流程图；
图3示出了本发明实施例中砂体个数频率曲线图；
图4示出了本发明实施例中基于砂体统计学的油藏小层划分装置的结构框图。
具体实施方式
下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
自上世纪80年代我国引进和开展油藏描述以来，其在油田开发实践中的应用成效显著。目前，油藏描述与剩余油分布规律研究的结合，作为精细油藏描述的一项重要基础工作，已成为确定挖掘潜力和提高采收效率等措施的重要途径。
油藏描述中，所谓“小层”通常是指单砂体或单砂层，属于油田最低级别的储层单元，为油气开发的基本单元。只有建立正确的等时对比，才能在油田范围内统一层组的划分，才能将砂体确定在同一个时间单元之内，其划分与对比的可靠程度直接关系到油藏描述的成败。
相关技术中，所采用的小层划分与对比的方案主要有“旋回对比、分级控制”、“时间单元划分、等高程”对比及“等高程切片对比”等，然而相关技术中的上述方式对于构造复杂和沉积复杂的油藏，由于相变频繁及缺少稳定标志层，导致小层的划分对比困难重重。
为了克服相关技术中，难以确定出小层划分方案的技术问题，本发明提供了一种基于砂体统计学的油藏小层划分方法，该方法将砂体统计学特征与高分辨率层序地层学理论相结合，并运用到地质分层方案及高分辨率层序地层划分与对比中去，有助于准确实现单砂体级别小层划分方案的确定，降低了油藏小层划分的难度，能够有效提高油藏开发的效果。
基于上述思想，本发明实施例提供了一种基于砂体统计学的油藏小层划分方法，如图1所示，该方法的主要处理步骤包括：
步骤S11：基于岩心资料及测井曲线确定隔夹层划分基准，根据隔夹层划分基准划分出单井目的层段内的砂体。
在矿产勘探和开发过程中，需要按地质设计的地层层位和深度，开展钻井工作，向井内下入取心工具，钻取出岩石样品。对岩心样品进行岩心分析是认识油气层地质特征的必要手段，岩心分析的主要方法有显微技术和分光技术等，通过对岩心样品进行岩心分析获得岩心资料。
测井是指利用岩层的电化学特性、导电特性、声学特性、放射性等地球物理特性，测量地球物理参数的方法。测井曲线是测井后形成的曲线。
本申请中利用岩心资料及测井曲线所反应的地层特征，确定出隔夹层的划分标准，其中隔夹层是隔层和夹层的统称，隔层是指在一定压差范围内能阻止流体在层组之间互相渗透的非渗透岩层；夹层是指单砂层(单油气层)之间或内部分布不稳定的不渗透或极低渗透的薄层。
油层内部夹有的隔夹层对提供油流不起作用，确定出隔夹层划分基准后，根据划分基准能够首先识别出目的层段内的隔层和夹层，确定隔层和夹层后即可划分出目的层段内具有储油功能的砂体。
在基于确定出的隔夹层划分基准，划分出单井目的段内的单砂体时，所采用的具体方式如下：
(1)根据岩心资料分析确定隔夹层划分标准；
(2)将岩心资料分析的结论与测井曲线进行关联，并基于隔夹层划分标准建立测井曲线识别隔夹层的定量标准；
(3)依据定量标准划分出单井目的层段内的砂体。
隔夹层划分中，关键是识别不同岩性的交接界面，岩心资料是识别各岩层分布的重要资料，因此在对岩心资料进行分析时，主要对砂质泥岩界面进行分析，同时还对泥岩、砂质泥岩及钙质砂层的厚度和孔隙度进行分析，由此获取目的段内各岩性地质的分布状况。
获取岩心资料分析结果后，可以结合实际所需的分层精度及分层工作难易程度，确定隔夹层的划分标准，其中确定出的隔夹层划分标准包括隔层划分标准及夹层划分标准。
考虑到隔夹层划分精度越高，实际划分时的难度越大，本申请中提供了一种隔夹层划分标准，利用该划分标准既能保证划分精度同时降低隔夹层实际划分时的工作难度，具体实施时，所述隔层划分标准为：泥岩厚度大于0.4米，砂质泥岩厚度大于0.6米，钙质砂层厚度大于1米且孔隙度小于8md的层划分为隔层；所述夹层划分标准为：泥岩厚度大于0.2米，砂质泥岩厚度大于0.3米，钙质砂层厚度大于0.5米且孔隙度小于8md的层划分为夹层。
本申请中对单井中的砂体划分是在测井曲线上进行的，为了实现对砂体的划分首先需要将岩心资料分析的结论与测井曲线进行关联，从而根据确定出的隔夹层划分标准建立测井曲线识别隔夹层的定量标准，根据该定量标准在测井曲线上划分出单井目的层段内的砂体。
在将岩心资料分析的结论与测井曲线进行关联时，所采用的方法包括：
岩心归位，即恢复岩心所在真实层位的工作。钻井取心时，岩心筒中可能有上次取心残留下来的岩心，而且岩心收获率一般达不到100％，以及钻具长度测量上产生的误差，使得岩心深度不准。因此在整理岩心时，必须对照电测资料，加上地质人员的判断，同时要校正钻具长度和测井深度上的系统误差，要将岩心的不同岩性与电测曲线解释的岩性仔细对应，恢复岩心所在真实深度。
对比相同深度段内的岩心资料、录井资料、成像测井及测井资料的特征；根据所述对比的结果，建立岩性、岩性界面与测井曲线上自然伽马GR值的对应关系，由此实现岩心分析结论与测井曲线的关联。
步骤S12：运用统计学方法统计油田所有单井中划分出的砂体的砂体厚度和砂体个数，并根据统计的结果绘制砂体厚度频率曲线和砂体个数频率曲线。
本步骤中，运用砂体统计学方法统计油田所有单井中划分出的砂体的砂体厚度和砂体个数，并根据结果绘制砂体厚度频率曲线和砂体个数频率曲线。
其中砂体厚度频率曲线是指横坐标为砂体厚度，纵坐标为对应厚度的砂体在砂体总样本中所占比重的曲线。例如砂体总数为100，砂体厚度为10-50的砂体个数占砂体总数的20％，50-100的砂体个数占砂体总数的30％……，根据该统计的结果绘制砂体厚度频率曲线。
砂体个数频率曲线是指，横坐标为砂体个数，纵坐标为具有相应砂体个数的单井在总井数量中所占的比重的曲线。例如，油田中的单井总量为100个，具有砂体个数为7的单井个数为2个，具有砂体个数为10的单井个数为11个，则在砂体个数频率曲线上，横坐标7所对应的纵坐标为2％，横坐标10所对应的纵坐标为11％。
本发明中，在统计砂体厚度时，首先需要计算单井中砂体的厚度，本方法中还进一步提供了砂体厚度的计算方法，具体计算方法如下：根据岩心资料分析的结论与测井曲线关联的结果，利用公式h＝0.125×N计算砂体的厚度；式中：0.125为测井曲线中测量点间距为0.125m；N为测井曲线自然伽马GR值中小于砂泥分界值GRlimit的个数；
式中：MEgr为测井序列GR中位数；GEOgr为测井序列GR平均值；MADgr为测井序列GR平均绝对值偏差。
根据上述公式计算得到划分出的每个砂体的厚度，根据统计得到的砂体厚度绘制砂体厚度频率曲线。
步骤S13：当砂体厚度频率曲线及砂体个数频率曲线均呈现单一峰态，且偏度和峰度在预期范围内时，选取砂体个数频率曲线的峰值所对应的砂体个数作为小层数目，并选取砂体个数与小层数目相等的井，结合旋回特征建立标准井。
步骤S14：将砂体个数多于标准井砂体个数的井，进行砂体合并；将砂体个数少于标准井砂体个数的井以等厚原则进行砂体劈分。
本发明中建立标准井，基于标准井实现小层的划分与对比。
本发明中，将砂体个数多于标准井砂体个数的井，进行砂体合并，包括：选取平行物源方向和垂直物源方向的过标准井的连井剖面；基于连井剖面，在预设标志层限定的层段内，将砂体个数多于标准井砂体个数的井中的多个砂体统一编号为相同小层砂体。
本发明中，将砂体个数少于标准井砂体个数的井以等厚原则进行砂体劈分，包括：选择测井曲线中所显示的隔层位置或夹层位置作为劈分界面进行劈分，劈分后得到的各个小层之间的厚度差值符合预设差值范围。
如图2示出了油藏小层划分技术流程图，图3示出了砂体个数频率曲线。
以下将结合图2及图3对本发明实施例的油藏小层划分方法进行详细说明，如图2所示，根据岩心及测井资料确定隔夹层划分标准，将岩心分析结论与测井曲线进行关联，建立测井曲线识别隔夹层的定量标准，依次划分出单井目的层段内的砂体；运用统计学方法统计油田所有井的单井砂体的单砂体厚度统计参数和砂体个数参数，并将所有单井参数作为样本统计砂体厚度和砂体个数以下参数：最大值，最小值、平均值、标志偏差、偏度和峰态；如果砂体厚度频率曲线和砂体个数频率曲线呈现单一峰态其偏度和分度在预期范围内，则选取砂体个数频率曲线的峰值所对应的砂体个数作为小层数目(如图3中，根据砂体个数频率曲线选取砂体个数11作为小层数目)，选取砂体个数与小层数目相同的井，结合旋回特征识别建立标准井；选取平行物源方向和垂直物源方向的过标准井连井剖面，在标志层确定的层段内，对砂体个数多于标准井砂体数目的井进行砂体合并，对少于标准井砂体数目的井，以等厚度原则，结合测井界面进行砂体劈分，实现小层的划分。基于标准井、砂体合并或劈分的结果以及成像测井资料建立对比剖面，从而基于对比剖面实现小层的对比。
对应于上述油藏小层划分方法，本发明实施例还提供了一种基于砂体统计学的油藏小层划分装置，参见图4所示的基于砂体统计学的油藏小层划分装置的结构框图，该装置包括以下模块：
砂体划分模块41，用于基于岩心资料及测井曲线确定隔夹层划分基准，根据隔夹层划分基准划分出单井目的层段内的砂体；
统计模块42，用于运用统计学方法统计油田所有单井中划分出的砂体的砂体厚度和砂体个数，并根据统计的结果绘制砂体厚度频率曲线和砂体个数频率曲线；
标准井建立模块43，用于当砂体厚度频率曲线及砂体个数频率曲线均呈现单一峰态，且偏度和峰度在预期范围内时，选取砂体个数频率曲线的峰值所对应的砂体个数作为小层数目，并选取砂体个数与小层数目相等的井，结合旋回特征建立标准井；
合并与劈分模块44，用于将砂体个数多于标准井砂体个数的井，进行砂体合并；将砂体个数少于标准井砂体个数的井以等厚原则进行砂体劈分。
本实施例中砂体划分模块，包括：划分标准确定单元，用于根据岩心资料分析确定隔夹层划分标准；关联单元，用于将岩心资料分析的结论与测井曲线进行关联，并基于隔夹层划分标准建立测井曲线识别隔夹层的定量标准；砂体划分单元，用于依据定量标准划分出单井目的层段内的砂体。
本实施例中砂体厚度计算模块，用于根据岩心资料分析的结论与测井曲线关联的结果，利用公式h＝0.125×N计算砂体的厚度；式中：0.125为测井曲线中测量点间距为0.125m；N为测井曲线自然伽马GR值中小于砂泥分界值GRlimit的个数；
式中：MEgr为测井序列GR中位数；GEOgr为测井序列GR平均值；MADgr为测井序列GR平均绝对值偏差。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
本发明实施例还提供了进行网络素材获取方法、界面控件显示方法的计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种基于砂体统计学的油藏小层划分方法，其特征在于，包括：
基于岩心资料及测井曲线确定隔夹层划分基准，根据所述隔夹层划分基准划分出单井目的层段内的砂体；
运用统计学方法统计油田所有单井中划分出的所述砂体的砂体厚度和砂体个数，并根据统计的结果绘制砂体厚度频率曲线和砂体个数频率曲线；
当所述砂体厚度频率曲线及所述砂体个数频率曲线均呈现单一峰态，且偏度和峰度在预期范围内时，选取砂体个数频率曲线的峰值所对应的砂体个数作为小层数目，并选取砂体个数与所述小层数目相等的井，结合旋回特征建立标准井；
将砂体个数多于所述标准井砂体个数的井，进行砂体合并；将砂体个数少于所述标准井砂体个数的井以等厚原则进行砂体劈分。
2.根据权利要求1所述的方法，其特征在于，所述基于岩心资料及测井曲线确定隔夹层划分基准，根据所述隔夹层划分基准划分出单井目的层段内的砂体，包括：
根据岩心资料分析确定隔夹层划分标准；
将所述岩心资料分析的结论与测井曲线进行关联，并基于所述隔夹层划分标准建立测井曲线识别隔夹层的定量标准；
依据所述定量标准划分出单井目的层段内的砂体。
3.根据权利要求2所述的方法，其特征在于，所述根据岩心资料分析确定隔夹层划分标准，包括：
根据岩心资料进行砂质泥岩界面分析，泥岩、砂质泥岩及钙质砂层的厚度和孔隙度分析；
根据所述岩心资料分析的结果确定隔层划分标准及夹层划分标准，其中所述隔层划分标准为：泥岩厚度大于0.4米，砂质泥岩厚度大于0.6米，钙质砂层厚度大于1米且孔隙度小于8md的层划分为隔层；所述夹层划分标准为：泥岩厚度大于0.2米，砂质泥岩厚度大于0.3米，钙质砂层厚度大于0.5米且孔隙度小于8md的层划分为夹层。
4.根据权利要求3所述的方法，其特征在于，所述将所述岩心资料分析的结论与测井曲线进行关联，包括：
岩心归位；
对比相同深度段内的岩心资料、录井资料、成像测井及测井资料的特征；
根据所述对比的结果，建立岩性、岩性界面与测井曲线上自然伽马GR值的对应关系。
5.根据权利要求2所述的方法，其特征在于，还包括：根据所述岩心资料分析的结论与测井曲线关联的结果，利用公式h＝0.125×N计算砂体的厚度；
式中：0.125为测井曲线中测量点间距为0.125m；
N为测井曲线自然伽马GR值中小于砂泥分界值GRlimit的个数；
<mrow>
<msub>
<mi>GR</mi>
<mrow>
<mi>lim</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>9000</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>ME</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>GEO</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>1000</mn>
</mfrac>
<msub>
<mi>MAD</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>ME</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>GEO</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中：
MEgr为测井序列GR中位数；
GEOgr为测井序列GR平均值；
MADgr为测井序列GR平均绝对值偏差。
6.根据权利要求1所述的方法，其特征在于，所述将砂体个数多于所述标准井砂体个数的井，进行砂体合并，包括：
选取平行物源方向和垂直物源方向的过标准井的连井剖面；
基于所述连井剖面，在预设标志层限定的层段内，将砂体个数多于标准井砂体个数的井中的多个砂体统一编号为相同小层砂体。
7.根据权利要求1所述的方法，其特征在于，所述将砂体个数少于所述标准井砂体个数的井以等厚原则进行砂体劈分，包括：
选择测井曲线中所显示的隔层位置或夹层位置作为劈分界面进行劈分，所述劈分后得到的各个小层之间的厚度差值符合预设差值范围。
8.一种基于砂体统计学的油藏小层划分装置，其特征在于，包括：
砂体划分模块，用于基于岩心资料及测井曲线确定隔夹层划分基准，根据所述隔夹层划分基准划分出单井目的层段内的砂体；
统计模块，用于运用统计学方法统计油田所有单井中划分出的所述砂体的砂体厚度和砂体个数，并根据统计的结果绘制砂体厚度频率曲线和砂体个数频率曲线；
标准井建立模块，用于当所述砂体厚度频率曲线及所述砂体个数频率曲线均呈现单一峰态，且偏度和峰度在预期范围内时，选取砂体个数频率曲线的峰值所对应的砂体个数作为小层数目，并选取砂体个数与所述小层数目相等的井，结合旋回特征建立标准井；
合并与劈分模块，用于将砂体个数多于所述标准井砂体个数的井，进行砂体合并；将砂体个数少于所述标准井砂体个数的井以等厚原则进行砂体劈分。
9.根据权利要求8所述的装置，其特征在于，所述砂体划分模块，包括：
划分标准确定单元，用于根据岩心资料分析确定隔夹层划分标准；
关联单元，用于将所述岩心资料分析的结论与测井曲线进行关联，并基于所述隔夹层划分标准建立测井曲线识别隔夹层的定量标准；
砂体划分单元，用于依据所述定量标准划分出单井目的层段内的砂体。
10.根据权利要求8所述的装置，其特征在于，还包括：砂体厚度计算模块，用于根据所述岩心资料分析的结论与测井曲线关联的结果，利用公式h＝0.125×N计算砂体的厚度；
式中：0.125为测井曲线中测量点间距为0.125m；
N为测井曲线自然伽马GR值中小于砂泥分界值GRlimit的个数；
<mrow>
<msub>
<mi>GR</mi>
<mrow>
<mi>lim</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>9000</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>ME</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>GEO</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>1000</mn>
</mfrac>
<msub>
<mi>MAD</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mover>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>ME</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>GEO</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中：
MEgr为测井序列GR中位数；
GEOgr为测井序列GR平均值；
MADgr为测井序列GR平均绝对值偏差。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510162800.8A CN104850732B (zh) | 2014-07-11 | 2015-04-08 | 一种基于砂体统计学的油藏小层划分方法及装置 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2014203829125 | 2014-07-11 | ||
| CN201420382912 | 2014-07-11 | ||
| CN201510162800.8A CN104850732B (zh) | 2014-07-11 | 2015-04-08 | 一种基于砂体统计学的油藏小层划分方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104850732A CN104850732A (zh) | 2015-08-19 |
| CN104850732B true CN104850732B (zh) | 2018-01-12 |
Family
ID=53850373
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510162800.8A Expired - Fee Related CN104850732B (zh) | 2014-07-11 | 2015-04-08 | 一种基于砂体统计学的油藏小层划分方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104850732B (zh) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105735978B (zh) * | 2016-02-19 | 2018-12-04 | 中国石油集团川庆钻探工程有限公司 | 碳酸盐层间岩溶型储层布井方法 |
| CN106772679A (zh) * | 2016-12-29 | 2017-05-31 | 中国石油天然气股份有限公司 | 一种地层对比的方法及装置 |
| CN109459791B (zh) * | 2017-09-06 | 2020-09-08 | 中国石油化工股份有限公司 | 一种利用测井曲线确定河道位置的方法和系统 |
| CN107730150A (zh) * | 2017-11-18 | 2018-02-23 | 东营汇烃石油技术有限公司 | 一种河流相储层小层自动对比的方法 |
| CN109403960B (zh) * | 2018-09-06 | 2022-03-29 | 中国石油集团川庆钻探工程有限公司 | 利用录井气测峰形态判断储层流体性质的方法 |
| CN109581509A (zh) * | 2018-11-15 | 2019-04-05 | 中国石油化工股份有限公司 | 旋回结构法小层划分对比的方法 |
| CN110134918B (zh) * | 2019-04-04 | 2021-05-28 | 中国石油大学(北京) | 基于滑动窗口法的岩心自动归位方法和系统 |
| CN112459766B (zh) * | 2019-09-09 | 2023-07-25 | 中国石油天然气股份有限公司 | 划分复式油气藏的方法及装置 |
| CN112668137B (zh) * | 2019-10-15 | 2022-11-01 | 中国石油天然气股份有限公司 | 重力流砂岩有效厚度的校正方法及校正装置 |
| CN113514879B (zh) * | 2020-04-09 | 2024-03-01 | 中国石油天然气股份有限公司 | “耳朵”层识别方法 |
| CN113534282B (zh) * | 2020-04-18 | 2024-03-29 | 中国石油化工股份有限公司 | 基于岩性对差值的断层活动频率计算方法 |
| CN113775337B (zh) * | 2020-06-10 | 2023-08-22 | 中国石油天然气股份有限公司 | 单砂体地质模型的建立方法、装置和计算机可读存储介质 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1873444A (zh) * | 2005-06-03 | 2006-12-06 | 中国石油大学(北京) | 用墨西哥帽小波对测井曲线进行小层对比的方法 |
| CN103184868A (zh) * | 2011-12-28 | 2013-07-03 | 陈美玲 | 一种砂岩油藏隔、夹层与剩余油分布识别方法 |
-
2015
- 2015-04-08 CN CN201510162800.8A patent/CN104850732B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1873444A (zh) * | 2005-06-03 | 2006-12-06 | 中国石油大学(北京) | 用墨西哥帽小波对测井曲线进行小层对比的方法 |
| CN103184868A (zh) * | 2011-12-28 | 2013-07-03 | 陈美玲 | 一种砂岩油藏隔、夹层与剩余油分布识别方法 |
Non-Patent Citations (5)
| Title |
|---|
| High Resolution sequence strtigraphic frameworks research in the 12 group of the member 3 of Qing Shankou formation in the Qianan oilfield;Cunlei Li等;《Applied mechanics and materials》;20130731;第339卷;728-731 * |
| Seismic sedimentology of conglomeratic sandbodies in lower third member of Shahejie Formation (Palaeogene) in Shengtuo area, East China;Yuan Yong等;《Journal of Central South University》;20141231(第12期);4630-4639 * |
| 基于统计学与专家系统的测井相自动识别;李存磊等;《地球物理学进展》;20110831;第26卷(第4期);1400-1408 * |
| 大情字井地区青山口组二段沉积相研究;李存磊;《中国优秀硕士学位论文全文数据库基础科学辑》;20120315;A011-118 * |
| 小层划分对比技术方法;温亮等;《内蒙古石油化工》;20101031(第19期);104-105 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104850732A (zh) | 2015-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104850732B (zh) | 一种基于砂体统计学的油藏小层划分方法及装置 | |
| US11704748B2 (en) | System and method for automatically correlating geologic tops | |
| Ringrose et al. | Reservoir model design | |
| CN106526693B (zh) | 裂缝识别方法和装置 | |
| EP2987003B1 (en) | System and method for automatically correlating geologic tops | |
| US9684084B2 (en) | Three-dimensional multi-modal core and geological modeling for optimal field development | |
| EP1779145B1 (en) | Method for creating facies probability cubes based upon geologic interpretation | |
| US20150088424A1 (en) | Identifying geological formation depth structure using well log data | |
| CN106951660A (zh) | 一种海相碎屑岩水平井储层测井解释方法及装置 | |
| US20140244227A1 (en) | Method for validating a training image for the multipoint geostatistical modeling of the subsoil | |
| CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
| CN108897066A (zh) | 碳酸盐岩裂缝密度定量预测方法及装置 | |
| CN107436452A (zh) | 基于概率神经网络算法的烃源岩预测方法及装置 | |
| CN106154322A (zh) | 测井曲线校正方法和装置 | |
| Kamali et al. | 3D geostatistical modeling and uncertainty analysis in a carbonate reservoir, SW Iran | |
| KR20210150917A (ko) | 웰 로그를 학습하여 암상을 추정하는 방법 및 장치 | |
| CN105986819B (zh) | 用于测井资料自动处理与综合解释的方法和装置 | |
| CN108663714A (zh) | 一种沉积微相刻画方法 | |
| CN116201535B (zh) | 油气藏目标井标志地层自动划分方法、装置及设备 | |
| Wylie Jr et al. | Well-log tomography and 3-D imaging of core and log-curve amplitudes in a Niagaran reef, Belle River Mills field, St. Clair County, Michigan, United States | |
| CN109165312A (zh) | 一种勘测数据管理系统 | |
| CN105989541B (zh) | 一种基于大数据的测井资料的处理方法 | |
| CN109598049B (zh) | 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法 | |
| Svanes et al. | Integration of subsurface applications to develop a dynamic stochastic modeling workflow | |
| CN114153002A (zh) | 储层天然裂缝三维地质建模方法、装置、电子设备及介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate 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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180112 Termination date: 20210408 |