基岩岩性识别和物性计算方法及装置制造方法及图纸

本发明专利技术提供一种基岩岩性识别和物性计算方法及装置，该方法包括：对基岩岩样的铸体薄片进行分析确定基岩岩样中的岩石类型；在获取基岩岩样对应深度的多个测井曲线数据的基础上，将多个测井曲线数据进行两两交会作散点图，从多个散点图中优选对基岩岩石类型敏感的测井曲线；根据岩心分析结果和测井敏感曲线，建立测井敏感曲线与基岩岩性和基岩孔隙度的关系，得出基岩二元岩性分类与孔隙度具有很好的相关关系，依据上述两者的关系确定基岩未取心段的岩石类型和孔隙度。该方案为新区没有岩心资料标定情况下进行基岩岩性的识别和物性的计算提供了较为可靠的方法。

【技术实现步骤摘要】
基岩岩性识别和物性计算方法及装置
本专利技术涉及基岩潜山二元岩性和孔隙度计算
，特别涉及一种基岩岩性识别和物性计算方法及装置。
技术介绍
基岩是不整合覆盖在沉积层序之上的变质岩或火成岩组合(不论年代)，是盆地形成和演化的基础，也是盆地勘探的重要目标之一。基岩潜山虽然在油气勘探和开发中所占比重不大，是一种非常规储层，但也可以形成世界级大油田。从岩石学上看，基岩的岩性十分复杂，依据岩心和岩屑的实验室铸体薄片(铸体薄片是将有色液态胶在真空加压下注入岩石孔隙空间，待液态胶固化后磨制成的岩石薄片)分析等手段，已经识别出岩浆岩、正变质岩和副变质岩等上百种岩石类型。一些学者在大量岩心分析对比的基础上提出了“优势岩性”的概念，并根据常规测井曲线的特征划分出了9大类岩性。有人为了简化起见，有时会将混合花岗岩等岩石也当作花岗岩来研究，并根据敏感测井曲线视深电阻率和中子交会划分出独立1区、2区和3区等。这一分区实际可能和岩性有密切关系。但实际在勘探中很难发育纯的花岗岩潜山。变质岩是最为常见的基岩岩石类型。基岩的电阻率变化差异巨大，而且与钻井用的泥浆矿物度有密切的关系。且在一些地区基岩的视深电阻率变化在10-300万欧姆·米之间。因此，用电阻-中子交会图来进行储层的划分可能无法推广应用。由于基岩岩性复杂，经历了复杂的构造演化和风化淋滤等后生成岩作用和埋藏后的热液作用等过程，储层的发育和岩性之间的关系变得非常复杂，岩石学角度的岩石分类不能较好地为测井基岩岩性判识和储层物性预测提供可靠的方法，分别通过声波、密度和中子计算的孔隙度或两两交互得到的孔隙度之间差异大，关键是对岩石的基质密度、基质声波和基质中子值无法准确获得。而孔隙度的计算对储层评价和储量计算非常重要。
技术实现思路
本专利技术实施例提供了一种基岩岩性识别和物性计算方法及装置，使用岩石的矿物组成来对基岩岩性进行区分，使得岩性的界限特征更加明显；基于上述对基岩岩性的区分方法获得的基岩岩性和常规的测井数据计算得到的孔隙度更加可信。该基岩岩性识别和物性计算方法包括：对基岩岩样的铸体薄片进行分析，确定基岩岩样中的矿物组成及矿物含量，根据所述矿物组成及矿物含量确定基岩岩样中的岩石类型；获取所述基岩岩样对应深度的多个测井曲线数据；以所述基岩岩样中的岩石类型为基础，将所述多个测井曲线数据进行两两交会作散点图，获得多个散点图，从所述多个散点图中选择符合预设要求的散点图对应的两个测井曲线作为测井敏感曲线；基于所述测井敏感曲线，对基岩未取心段进行测井岩性和物性分析，确定基岩未取心段的岩石类型和孔隙度。该基岩岩性识别和物性计算装置包括：岩石类型确定模块，用于对基岩岩样的铸体薄片进行分析，确定基岩岩样中的矿物组成及矿物含量，根据所述矿物组成及矿物含量确定基岩岩样中的岩石类型；测井曲线数据获取模块，用于获取所述基岩岩样对应深度的多个测井曲线数据；散点图绘制模块，用于以所述基岩岩样中的岩石类型为基础，将所述多个测井曲线数据进行两两交会作散点图，获得多个散点图，从所述多个散点图中选择符合预设要求的散点图对应的两个测井曲线作为测井敏感曲线；孔隙度确定模块，用于基于所述测井敏感曲线，对基岩未取心段进行测井岩性和物性分析，确定基岩未取心段的岩石类型和孔隙度。本专利技术还提出一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述所述基岩岩性识别和物性计算方法。本专利技术实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述所述基岩岩性识别和物性计算方法的计算机程序。在本专利技术实施例中，根据分析后得到的基岩岩样的矿物组成及矿物含量确定基岩岩样中的岩石类型，相比较使用岩石学角度的岩石分类，本专利技术可以较好地为测井基岩岩性判识提供可靠的方法；基于上述对基岩岩性的区分方法获得的基岩岩性和常规的测井数据计算基岩岩样的孔隙度，相比较通过声波、密度和中子计算的孔隙度或两两交互得到的孔隙度，本专利技术得到的孔隙度更加可信。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1是本专利技术实施例提供的一种基岩岩性识别和物性计算方法流程图(一)；图2是本专利技术实施例提供的一种密度-自然伽马(API)交会图；图3是本专利技术实施例提供的一种中子-密度交会图；图4是本专利技术实施例提供的一种岩心标定后的长英质和铁镁质岩石的测井孔隙度和密度交会图；图5是本专利技术实施例提供的一种双矿物系岩石骨架交会三角图；图6是本专利技术实施例提供的一种中子-密度交会图(二)；图7是本专利技术实施例提供的一种中子-声波交会图；图8是本专利技术实施例提供的一种基岩岩性识别和物性计算方法流程图(二)；图9是本专利技术实施例提供的一种BC-2井536-562米井段常规测井长英质和岩心长英质含量对比图；图10是本专利技术实施例提供的一种BC-2井922-983米井段常规测井铁镁质和岩心铁镁质含量对比图；图11是本专利技术实施例提供的一种基岩岩性识别和物性计算装置结构框图(一)；图12是本专利技术实施例提供的一种基岩岩性识别和物性计算装置结构框图(二)。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。针对上述
技术介绍
中存在的问题，本专利技术的目的是提供一种应用测井进行基岩岩性识别和物性计算方法和装置。如图1所示，该基岩岩性识别和物性计算方法包括：步骤101：对基岩岩样的铸体薄片进行分析，确定基岩岩样中的矿物组成及矿物含量，根据所述矿物组成及矿物含量确定基岩岩样中的岩石类型。具体的，按照如下方式执行步骤101：在以岩心柱样、井壁取心的铸体薄片分析的基础上，统计岩石的主要矿物成分，如石英、斜长石、钾长石、角闪石、黑云母、粘土矿物含量和类型等；二元岩性归集，将浅色矿物(石英和长石)和深色矿物(角闪石、黑云母等)进行含量的归集和归一化为100％，当暗色矿物的含量大于预设比例含量(本专利技术设定为15％)时，所述暗色矿物为深色矿物，所述深色矿物为铁镁质岩石；当暗色矿物的含量小于预设比例含量时，所述暗色矿物为浅色矿物，所述浅色矿物为长英质岩石。长英质矿物系：含有硅、钾、钠等元素的无色或颜色较浅的矿物，由于硅、钾、钠等元素原子序数较小，因此在常规测井响应上表现为低密度值；因长石类矿物不含结晶水和结构水，因此中子值较低；含碱性长石类的矿物中富含高放射性元素K40，具有较高的自然伽马值，含斜长石类的矿物由于不含高放射性元素K40，自然伽马数值较低。铁镁质矿物系：含有铁、镁、钙、钛等元素，具有深浅不同颜色的矿物，主要有黑云母、辉石、普通角闪石、橄榄石、绿泥石等，该类矿物因含有原子序数及密度较大的铁、镁、钙、钛等元素，且富含结晶水或结构水，不含高放射性元素K40，因此，其测井响应特征为高密度、高中子、低自然伽马。步骤102：获取本文档来自技高网...

【技术保护点】
1.一种基岩岩性识别和物性计算方法，其特征在于，包括：对基岩岩样的铸体薄片进行分析，确定基岩岩样中的矿物组成及矿物含量，根据所述矿物组成及矿物含量确定基岩岩样中的岩石类型；获取所述基岩岩样对应深度的多个测井曲线数据；以所述基岩岩样中的岩石类型为基础，将所述多个测井曲线数据进行两两交会作散点图，获得多个散点图，从所述多个散点图中选择符合预设要求的散点图对应的两个测井曲线作为测井敏感曲线；基于所述测井敏感曲线，对基岩未取心段进行测井岩性和物性分析，确定基岩未取心段的岩石类型和孔隙度。

【技术特征摘要】
1.一种基岩岩性识别和物性计算方法，其特征在于，包括：对基岩岩样的铸体薄片进行分析，确定基岩岩样中的矿物组成及矿物含量，根据所述矿物组成及矿物含量确定基岩岩样中的岩石类型；获取所述基岩岩样对应深度的多个测井曲线数据；以所述基岩岩样中的岩石类型为基础，将所述多个测井曲线数据进行两两交会作散点图，获得多个散点图，从所述多个散点图中选择符合预设要求的散点图对应的两个测井曲线作为测井敏感曲线；基于所述测井敏感曲线，对基岩未取心段进行测井岩性和物性分析，确定基岩未取心段的岩石类型和孔隙度。2.如权利要求1所述的基岩岩性识别和物性计算方法，其特征在于，对基岩岩样的铸体薄片进行分析，确定基岩岩样中的矿物组成及矿物含量，根据所述矿物组成及矿物含量确定基岩岩样中的岩石类型，包括：所述矿物组成包括暗色矿物，暗色矿物包括浅色矿物和/或深色矿物，当暗色矿物的含量大于预设比例含量时，所述暗色矿物为深色矿物，所述深色矿物为铁镁质岩石；当暗色矿物的含量小于预设比例含量时，所述暗色矿物为浅色矿物，所述浅色矿物为长英质岩石。3.如权利要求1所述的基岩岩性识别和物性计算方法，其特征在于，所述多个测井曲线数据包括密度测井值、自然伽马测井值、自然电位测井值、补偿中子测井值、声波测井值和电阻率测井值。4.如权利要求1所述的基岩岩性识别和物性计算方法，其特征在于，所述预设要求为散点图中的预设比例的散点区分的基岩岩性与基岩岩样的岩性一致。5.如权利要求2所述的基岩岩性识别和物性计算方法，其特征在于，基于所述测井敏感曲线，对基岩未取心段进行测井岩性和物性分析，确定基岩未取心段的岩石类型和孔隙度，包括：在测井敏感曲线对应的交会图中绘制交会三角形；根据绘制的交会三角形确定计算长英质矿物含量和铁镁质矿物含量的测井解释所需参数，根据所述测井解释所需参数计算基岩未取心段孔隙度。6.如权利要求5所述的基岩岩性识别和物性计算方法，其特征在于，根据所述测井解释所需参数按照如下公式计算基岩未取心段孔隙度：DEN校＝DEN-V泥质×DEN泥质；CN校＝CN-V泥质×CN泥质；DEN校＝V1×DEN1+V2×DEN2+φ×DENf；CN校＝V1×CN1+V2×CN2+φ×CNf；1＝V1+V2+φ；其中，SH表示自然伽马相对变化值；GR表示自然伽马测井值；GRmax表示自然伽马极大值；GRmin表示自然伽马极小值；V泥质表示泥质含量；GCUR表示计算泥质含量所用的经验参数；DEN校表示校正后的密度测井值；DEN表示密度测井值；DEN泥质表示泥岩密度值；CN校表示校正后的中子测井值；CN表示中子测井值；CN泥质表示泥岩中子值；V1表示长英质矿物含量；DEN1表示长英质矿物密度骨架值；V2表示铁镁质矿物含量；DEN2表示铁镁质矿物密度骨架值；φ表示孔隙度；DENf表示孔隙流体密度值；CN1表示长英质矿物中子骨架值；CN2表示铁镁质矿物中子骨架值；CNf表示孔隙流体中子值。7.如权利要求1或5所述的基岩岩性识别和物性计算方法，其特征在于，还包括：将所述基岩未取心段孔隙度与实验得到的基岩岩样的孔隙度进行比较，根据比较结果确定是否对所述基岩未取心段孔隙度进行校正。8.如权利要求1所述的基岩岩性识别和物性计算方法，其特征在于，还包括：根据基岩未取心段测井对应的深度、测井参数、岩石类型和孔隙度绘制剖面图；所述基岩未取心段的测井参数包括密度测井值、自然伽马测井值、自然电位测井值、补偿中子测井值、声波测井值、电阻率测井值其中之一或多个。9.一种基岩岩性识别和物性计算装置，其特征在于，包括：岩石类型确定模块，用于对基岩岩样的铸体薄片进行分析，确定基岩岩样中的矿物组成及矿物含量，根据所述矿物组成及矿物含量确定基岩岩样中的岩石类型；...

【专利技术属性】
技术研发人员：窦立荣，肖玉峰，李国平，
申请(专利权)人：中国石油天然气集团有限公司，中国石油国际勘探开发有限公司，
类型：发明
国别省市：北京,11