低孔、致密灰岩岩石结构组分测井参数的计算方法及装置的制造方法
【技术领域】
[0001] 本发明涉及石油地质勘探和测井技术领域，尤其涉及一种低孔、致密灰岩岩石结 构组分测井参数的计算方法及装置。
【背景技术】
[0002] 目前，利用测井资料识别灰岩岩石结构组分一直是一个世界难题，而如何获得有 效测井参数是提高识别精度的关键。当前获得有效测井参数的研究方法和思路主要分为以 下几个步骤：1)岩-电关系分析:首先对取心井的岩心钻孔，获取样品，对样品进行磨片，获 取薄片，利用显微镜进行薄片鉴定，获取岩石结构组分信息，并将薄片深度归位到测井深 度，然后进行岩石结构组分与测井相关分析，确定与岩石结构组分相关性最好的测井参数 (常规测井曲线、合成参数曲线等），建立岩石结构组分的测井特征参数样本库。2)技术识 另IJ:主要采用神经网络、聚类分析、判别分析、模糊数学等数学算法，神经网络法对建立的样 本测井参数库进行训练，获取人工智能;聚类分析法建立测井相与岩石结构组分对应关系； 判别分析方法建立不同岩石结构组分的测井判别公式。模糊数学法建立不同岩石结构组分 测井分类模式的隶属度函数。3)技术应用:通过对研究区取心井的研究，确定最佳的测井参 数，选取最合适的识别技术，对未取心井开展岩石结构组分识别，并用工区内其他取心井验 证识别效果。
[0003] 现有技术中关于岩石结构组分测井识别方法在国内研究较少。而在国外专门针对 岩石结构组分测井识别研究一般是根据灰岩岩颗粒或者晶粒尺寸进行岩石物理三分类，第 一类将尺寸大于lOOym定义为颗粒结构，第二类将尺寸在20μηι-100μηι定义为泥晶颗粒结构， 第三类将尺寸小于20μπι定义为泥晶为主结构(包括部分泥晶颗粒结构、颗粒泥晶结构、泥晶 结构），通过岩心孔隙度、渗透率、含水饱和度参数建立了岩石结构数（Rock Fabric Number)计算公式，通过测井计算孔隙度、渗透率及含水饱和度代替岩心参数，实现岩石结 构组分测井识别。该方法建立在以下两个地质条件下：（1)地层主要以粒间、晶间孔或溶孔 为主，裂缝不发育；（2)地层在气水界面之上，不存在纯水层，满足这两个条件下地层岩石结 构组分测井识别能够取得很好识别效果。
[0004] 现有技术的该方法有效解决岩石结构组分测井识别的技术瓶颈，但是其所应用的 岩石结构数公式也存在不足:其建立基础是地层岩心孔隙度范围在4%_40%之间，而对于 国内低孔、致密灰岩岩石储层，岩心孔隙度范围在0.5%-10%之间，通过现有技术的岩石结 构数公式已经无法准确确定低孔、致密灰岩岩石的结构数参数。
【发明内容】
[0005] 本发明实施例提供一种低孔、致密灰岩岩石结构组分测井参数的计算方法及装 置，以解决现有技术中的岩石结构数公式已经不适用于岩心孔隙度范围在〇.5%-10%之间 的低孔、致密灰岩岩石储层，造成无法准确确定低孔、致密灰岩岩石的结构数参数的问题。
[0006] 为了达到上述目的，本发明采用如下技术方案：
[0007] -种低孔、致密灰岩岩石结构组分测井参数的计算方法，包括：
[0008] 获取待测灰岩岩石储层的测井数据;所述测井数据包括:深电阻率测井曲线、冲洗 带电阻率测井曲线、泥浆滤液电阻率、自然伽马曲线以及声波时差曲线；
[0009] 根据所述深电阻率测井曲线和冲洗带电阻率测井曲线确定地层裂缝孔隙度；
[0010] 根据深电阻率测井曲线中的深电阻率大小、自然伽马曲线中的自然伽马值大小、 地层裂缝孔隙度的大小以及声波时差曲线中的声波时差值的大小，确定第一岩石结构数参 数、第二岩石结构数参数、第三岩石结构数参数以及第四岩石结构数参数；
[0011] 将所述第一岩石结构数参数、第二岩石结构数参数、第三岩石结构数参数以及第 四岩石结构数参数叠加，确定岩石结构数参数总和；
[0012] 根据所述岩石结构数参数总和和深电阻率测井曲线中的深电阻率大小，确定岩石 结构数参数最终值。
[0013] 具体的，根据所述深电阻率测井曲线和冲洗带电阻率测井曲线确定地层裂缝孔隙 度，包括：
[0014] 判断同一时刻下深电阻率测井曲线中的深电阻率和冲洗带电阻率测井曲线中的 冲洗带电阻率的大小；
[0015] 若深电阻率大于冲洗带电阻率，则根据公式：
[0017] 计算地层裂缝孔隙度?frc;其中RX0为冲洗带电阻率;RT为深电阻率;Rmf为泥浆滤 液电阻率；
[0018] 若深电阻率小于冲洗带电阻率，则根据公式：
[0020]计算地层裂缝孔隙度〇frc。
[0021 ]具体的，根据深电阻率测井曲线中的深电阻率大小、自然伽马曲线中的自然伽马 值大小、地层裂缝孔隙度的大小以及声波时差曲线中的声波时差值的大小，确定第一岩石 结构数参数、第二岩石结构数参数、第三岩石结构数参数以及第四岩石结构数参数，包括：
[0022] 根据所述自然伽马值大小和地层裂缝孔隙度的大小，确定所述第一岩石结构数参 数；
[0023] 所述根据所述自然伽马值大小和地层裂缝孔隙度的大小，确定所述第一岩石结构 数参数，包括：
[0024] 若自然伽马值大于50API，且地层裂缝孔隙度小于0.08%，根据公式：
[0026]计算所述第一岩石结构数参数ARFN1;其中，RT为深电阻率;GR为自然伽马值；
[0027]若自然伽马值大于50API，且地层裂缝孔隙度大于0.08%，根据公式：
[0029]计算所述第一岩石结构数参数ARFNl
[0030] 此外，根据深电阻率测井曲线中的深电阻率大小、自然伽马曲线中的自然伽马值 大小、地层裂缝孔隙度的大小以及声波时差曲线中的声波时差值的大小，确定第一岩石结 构数参数、第二岩石结构数参数、第三岩石结构数参数以及第四岩石结构数参数，包括：
[0031] 根据所述声波时差值的大小和自然伽马值大小，确定所述第二岩石结构数参数；
[0032] 所述根据所述声波时差值的大小和自然伽马值大小，确定所述第二岩石结构数参 数，包括：
[0033]若声波时差值小于等于51us/f，且自然伽马值小于等于50API，根据公式：
[0035]计算所述第二岩石结构数参数ARFN2;其中，e为自然底数;AC为声波时差值;RT为 深电阻率。
[0036] 此外，根据深电阻率测井曲线中的深电阻率大小、自然伽马曲线中的自然伽马值 大小、地层裂缝孔隙度的大小以及声波时差曲线中的声波时差值的大小，确定第一岩石结 构数参数、第二岩石结构数参数、第三岩石结构数参数以及第四岩石结构数参数，包括：
[0037] 根据所述声波时差值的大小、自然伽马值大小和深电阻率大小，确定所述第三岩 石结构数参数；
[0038] 所述根据所述声波时差值的大小、自然伽马值大小和深电阻率大小，确定所述第 三岩石结构数参数，包括：
[0039] 若声波时差值大于51us/f，自然伽马值小于等于50API，且深电阻率大于200ohmm， 根据公式：
[0041 ]计算所述第三岩石结构数参数ARFN3;其中，e为自然底数;AC为声波时差值;RT为 深电阻率。
[0042] 此外，根据深电阻率测井曲线中的深电阻率大小、自然伽马曲线中的自然伽马值 大小、地层裂缝孔隙度的大小以及声波时差曲线中的声波时差值的大小，确定第一岩石结 构数参数、第二岩石结构数参数、第三岩石结构数参数以及第四岩石结构数参数，包括：
[0043] 根据所述声波时差值的大小、自然伽马值大小、深电阻率和地层裂缝孔隙度的大 小，确定所述第四岩石结构数参数；
[0044] 所述根据所述声波时差值的大小、自然伽马值大小、深电阻率和地层裂缝孔隙度 的大小，确定所述第四岩石结构数参数，包括：
[0045] 若声波时差值大于51us/f，自然伽马值小于等于40API，深电阻率小于200ohmm且 地层裂缝孔隙度小于0.08%，根据公式：
[0047]计算所述第四岩石结构数参数ARFN4;其中，e为自然底数;AC为声波时差值;RT为 深电阻率;GR为自然伽马值；
[0048] 若声波时差值大于51us/f，自然伽马值大于40API且小于等于50API，深电阻率小 于200ohmm且地层裂缝孔隙度大于0.08 %，根据公式：
[0050]计算所述第四岩石结构数参数ARFN4。
[0051]具体的，将所述第一岩石结构数参数、第二岩石结构数参数、第三岩石结构数参数 以及第四岩石结构数参数叠加，确定岩石结构数参数总和，包括：
[0052]根据公式：
[0053] ARFN=ARFN1+ARFN2+ARFN3+ARFN4 [0054]确定岩石结构数参数总和ARFN。
[0055]具体的，根据所述岩石结构数参数总和和深电阻率测井曲线中的深电阻率大小， 确定岩石结构数参数最终值，包括：
[0056]若ARFN小于等于第一预设值，确定所述岩石结构数参数最终值ARFf等于所述 ARFN；
[0057] 若ARFN大于所述第一预设值，且深电阻率大于9000ohmm，确定所述岩石结构数参 数最终值ARFN'等于第二预设值；
[0058] 若ARFN大于所述第一预设值，且深电阻率小于等于9000ohmm，确定所述岩石结构 数参数最终值ARFN'等于所述第一预设值；
[0059]若ARFN小于等于第三预设值，确定所述岩石结构数参数最终值ARFf等于所述第 三预设值。
[0060] 一种低孔、致密灰岩岩石结构组分测井参数的计算装置，包括：
[0061] 测井数据获取单元，用于获取待测灰岩岩石储层的测井数据;所述测井数据包括： 深电阻率测井曲线、冲洗带电阻率测井曲线、泥浆滤液电阻率、自然伽马曲线以及声波时差 曲线；
[0062] 地层裂缝孔隙度确定单元，用于根据所述深电阻率测井曲线和冲洗带电阻率测井 曲线确定地层裂缝孔隙度；
[0063]岩石结构数参数确定单元，用于根据深电阻率测井曲线中的深电阻率大小、自然 伽马曲线中的自然伽马值大小、地层裂缝孔隙度的大小以及声波时差曲线中的声波时差值 的大小，确定第一岩石结构数参数、第二