本发明属于油气勘探、开发领域，尤其涉及一种钙屑砂岩储层孔隙类型识别方法。
背景技术
储层孔隙类型是认识、评价储层重点内容之一，孔隙类型在一定程度上可以反映储层的形成机制；不同的孔隙类型具有不同的孔隙结构，最终可导致受渗透率大小控制的开发效果存在巨大差异。目前，常见的利用钻井岩心样品，磨制铸体薄片，通过显微镜下观察来确定，或者直接利用岩心样品在扫描电镜下观察，进而确定某一目标储层的孔隙类型，描述其分布特征，这些方法对于相对均质的储层比较适用，但是对于储层非均质性较强的钙屑砂岩储层则不太适用，主要存在如下两个方面的问题：
(1)由于钻井取心成本等的限制，不可能在钻井过程中进行大规模的取心。钙屑砂岩是一种特殊类型的主要由碳酸盐岩屑构成的碎屑岩，孔隙类型复杂多样，钙屑砂岩与砾岩、常规砂岩等，在纵向上交互、横向上叠置分布，非均质性非常强。这使得实际取心与设计取心存在误差，导致实际取心样品做分析测试反映的储层孔隙类型，不能真正、全面的反映全井段储层的孔隙类型。
(2)岩屑作为完钻井录井的重要成果之一，目的层段每隔0.5m连续取样，通过精选岩屑样品进一步磨片观察也可识别孔隙类型，但在实际工作中利用率还是相对较低，且对于井数较多研究区，岩屑的选取工作量太大，导致储层孔隙类型在纵向和横向上的分布特征难以准确描述。
研究表明，利用常规测井曲线识别储集岩孔隙类型、描述其空间分布是解决上述问题有效办法之一。但现有技术中实际应用该类方法时，有的仅适用于纯灰岩或白云岩，且计算方法不明确，有的仅适用于火山碎屑岩，有的只能部分地识别孔隙类型，无法识别微孔隙，有的则判断标准过于笼统，难于实施。因此，对于孔隙类型多样的钙屑砂岩，有待总结一套新的方法对孔隙类型进行识别，描述其空间分布特征，从而为钙屑砂岩的进一步的储层评价提供基础。
技术实现要素：
为解决上述现有技术中在判断钙屑砂岩储层孔隙类型时，由于钻井取心成本等的限制，使得取心少，由于非均质性强，导致实际取心与设计取心存在误差，取心代表性差，所取样品不能真正、全面反映全井段储层孔隙类型；井数较多，岩屑的选取工作量太大，导致储层孔隙类型在纵向和横向上的分布特征难以准确描述；以及现有的利用常规测井曲线判断孔隙类型的方法难以适用钙屑砂岩储层的技术问题，本发明提供一种钙屑砂岩储层孔隙类型识别方法，具体方案如下：
一种钙屑砂岩储层孔隙类型识别方法，包括以下步骤：
选取钙屑砂岩发育段所在的目的层段；
选取所述目的层段中的钙屑砂岩干层段；
获取所述目的层段所属井的声波时差孔隙度曲线和密度孔隙度曲线，并将所述目的层段所属井的声波时差孔隙度曲线、密度孔隙度曲线与所述目的层段所属井的孔隙度较正曲线一起放置在同一测井曲线道内；
以所述目的层段所属井的孔隙度较正曲线为基准，深度坐标不变，在孔隙度坐标方向上将所述目的层段所属井的声波时差孔隙度曲线、密度孔隙度曲线平移，使所述目的层段所属井的密度孔隙度曲线、声波时差孔隙度曲线与所述目的层段所属井的孔隙度较正曲线在所述目的层段中的钙屑砂岩干层段重合；
设置孔隙类型的占比；
三条孔隙度曲线在所述目的层段中的钙屑砂岩干层段重合后，根据同一口井的测井曲线上，相同深度位置的密度孔隙度和声波时差孔隙度的相差程度与该深度位置处钙屑砂岩储层中满足预设占比条件的孔隙类型的对应关系识别所述目的层段所属井的钙屑砂岩储层中满足预设占比条件的孔隙类型。
优选的，所述孔隙类型的占比大于50％，所述同一口井的测井曲线上，相同深度位置的密度孔隙度和声波时差孔隙度的相差程度与该深度位置处钙屑砂岩储层中满足预设占比条件的孔隙类型的对应关系为：
当同一深度位置的密度孔隙度减去声波时差孔隙度的差值＞3％时，该深度位置的储层中占比大于50％的孔隙类型为粒间溶蚀孔隙；
当同一深度位置的密度孔隙度减去声波时差孔隙度的差值≥0.5％，且≤3％时，该深度位置的储层中占比大于50％的孔隙类型为粒间溶蚀充填孔隙；
当同一深度位置的密度孔隙度减去声波时差孔隙度的差值＜0.5％时，该深度位置的储层中占比大于50％的孔隙类型为微孔隙。
优选的，选取已知钙屑砂岩发育段所在的目的层段的岩屑包，从中采集岩屑样品，经磨制后形成岩屑薄片样品，连同岩心薄片样品，在显微镜下详细观察，分析并选取孔隙类型满足占比条件的样品；
选取所述已知目的层段中的钙屑砂岩干层段；
获取所述已知目的层段所属井的声波时差孔隙度曲线和密度孔隙度曲线，并将所述已知目的层段所属井的声波时差孔隙度曲线、密度孔隙度曲线与所述已知目的层段所属井的孔隙度较正曲线一起放置在同一测井曲线道内；
以所述已知目的层段所属井的孔隙度较正曲线为基准，深度坐标不变，在孔隙度坐标方向上将所述已知目的层段所属井的声波时差孔隙度曲线、密度孔隙度曲线平移，使已知目的层段所属井的密度孔隙度曲线、声波时差孔隙度曲线与已知目的层段所属井的孔隙度较正曲线在所述已知目的层段中的钙屑砂岩干层段重合；
读取所述已知目的层段所属井的密度孔隙度曲线和声波时差孔隙度曲线上，孔隙类型满足占比条件的样品所在深度相应的声波时差孔隙度和密度孔隙度，根据读取结果总结同一口井的测井曲线上，相同深度位置的密度孔隙度和声波时差孔隙度的相差程度与该深度位置处钙屑砂岩储层中满足预设占比条件的孔隙类型的对应关系。
优选的，根据钙屑砂岩的岩性识别标准选取所述钙屑砂岩发育段，所述钙屑砂岩的岩性识别标准为：自然伽玛测井的范围为40api＜gr＜55api，同时深侧向电阻率的范围为300欧姆·米<rd<4000欧姆·米。
优选的，选取已知钙屑砂岩发育段所在的目的层段的岩屑包，从中采集岩屑样品，经磨制后形成岩屑薄片样品，连同岩心薄片样品，在显微镜下详细观察，鉴定样品岩性；
读取所述已知目的层段所属井的常规测井曲线上不同岩性的样品所在深度对应的值，进而总结所述钙屑砂岩的岩性识别标准。
优选的，选取深侧向电阻率的范围为1000欧姆·米<rd<4000欧姆·米，密度大于2.65g/cm3，且补偿中子测井cnl小于5％的层段为所述钙屑砂岩干层段。
优选的，读取任一目的层段中的钙屑砂岩干层段各深度点对应的声波时差值和密度值，并将获得的所述声波时差值和密度值作为该目的层段的钙屑砂岩的声波时差骨架值和密度骨架值，利用声波时差孔隙度解释模型和密度孔隙度解释模型，分别计算该目的层段所属井的声波时差孔隙度和密度孔隙度，进而获得该目的层段所属井的的声波时差孔隙度曲线和密度孔隙度曲线。
优选的，所述孔隙度较正曲线为中子孔隙度曲线。
优选的，所述孔隙度较正曲线为岩心孔隙度曲线。
优选的，所述密度孔隙度和声波时差孔隙度的相差程度为密度孔隙度减去声波时差孔隙度的差值，或密度孔隙度与声波时差孔隙度的差值的绝对值，或密度孔隙度与声波时差孔隙度的差值的绝对值的倒数，或密度孔隙度与声波时差孔隙度的差值的平方。
与现有技术相比，本发明提供的一种钙屑砂岩储层孔隙类型识别方法，通过选取钙屑砂岩发育段所在的目的层段和所述目的层段中的钙屑砂岩干层段，并获取所述目的层段所属井的密度孔隙度曲线和声波时差孔隙度曲线，将两者与所述目的层段所属井的孔隙度较正曲线一起置于同一测井曲线道内，且以所述目的层段所属井的孔隙度较正曲线为基准，通过平移的方式使三条孔隙度曲线在所述目的层段中的钙屑砂岩干层段重合；设置孔隙度类型占比，三条孔隙度曲线在所述目的层段中的钙屑砂岩干层段重合后，根据同一口井的测井曲线上，相同深度位置的密度孔隙度和声波时差孔隙度的相差程度与该深度位置处钙屑砂岩储层中满足预设占比条件的孔隙类型的对应关系识别所述目的层段所属井的钙屑砂岩储层中满足预设占比条件的孔隙类型。该方法利用孔隙度测井曲线来识别钙屑砂岩储层中满足预设占比条件的孔隙类型，识别过程中不需要采集岩心样品和岩屑样品，因而没有由于取心成本限制而使得取心少；由于非均质性强致使实际取心与设计取心存在误差，取心代表性差而引起的所取样品不能真正、全面反映全井段储层孔隙类型的问题；以及，井数较多，岩屑的选取工作量太大，导致储层孔隙类型在纵向和横向上的分布特征难以准确描述的问题。该方法可操作性强，成本低，符合测井学原理，经中石化元坝须家河组须三段钙屑砂岩气藏应用证实，该方法可以方便、快捷并准确识别出不同的孔隙类型，判断标准清晰，易于实施，对钙屑砂岩储层孔隙类型的识别有很好的适应性，为钙屑砂岩储层的认识、综合评价奠定了良好的基础。
附图说明
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中：
图1为本发明钙屑砂岩储层孔隙类型识别方法流程图；
图2为本发明总结钙屑页岩岩性识别标准的实验方法流程图；
图3为本发明总结孔隙类型与孔隙度差值对应关系的实验方法流程图；
图4为本发明实施例中元陆6井的不同岩性的测井响应特征示意图；
图5为本发明实施例中元坝2井的钙屑砂岩薄片样品显示的占比大于50％的孔隙类型及其测井曲线响应特征；
图6为本发明实施例中元陆20井的钙屑砂岩薄片样品显示的占比大于50％的孔隙类型及其测井曲线响应特征；
图7为本发明实施例中元陆20井的钙屑砂岩薄片样品显示的占比大于50％的孔隙类型及其测井曲线响应特征。
在附图中，相同的部件采用相同的附图标记，附图并未按实际比例绘制。
具体实施方式
下面将结合附图对本发明作进一步的说明
钙屑砂岩发育段包括孔隙不发育的钙屑砂岩干层和发育粒间溶蚀孔隙、粒间溶蚀充填孔隙、微孔隙的钙屑砂岩储层。粒间溶蚀孔隙是指钙屑颗粒之间由于溶蚀作用而形成的孔隙，粒间溶蚀孔隙的孔径通常介于0.03mm到0.5mm之间。粒间溶蚀充填孔隙是指粒间溶蚀孔隙被粘土矿物全充填或部分充填所形成的孔隙，包括粒间溶蚀孔隙中充填的粘土矿物的晶间孔，及部分未被充填的残余粒间溶蚀孔隙，残余粒间溶蚀孔隙即粒间溶蚀孔隙被粘土矿物部分充填后余留下来的空间所形成的孔隙，例如某一粒间溶蚀孔隙被粘土矿物充填了70％，粒间溶蚀孔隙中余下的30％的空间所形成的孔隙即残余粒间溶蚀孔隙。上述比例只是示例性的表述，并不构成对本发明的限制。微孔隙是指呈侵染状分布于钙屑砂岩粒间或粒内的溶蚀微孔隙，其孔径通常小于0.03mm。此三种孔隙类型基本覆盖了钙屑砂岩中存在的所有孔隙类型，钙屑砂岩中即使存在其他的孔隙类型，其量也极少，在实际应用中可以忽略不计。
图4a、图5a、图6a、图7a的测井曲线图中，gr(gamma-ray)表示自然伽玛测井，单位为api。rd(deepresistivity)表示深侧向电阻率，单位为欧姆·米(ω·m)。rs(shortresistivity)表示浅侧向电阻率，单位为欧姆·米(ω·m)。den(density)表示密度，单位为g/cm3。cnl(compensatedneutronlogging)补偿中子测井，无因次量。ac(acoustic)声波时差，单位为微秒/英尺(μs/ft)。denpor(densityporosity)表示密度孔隙度，无因次量。acpor(acousticporosity)表示声波时差孔隙度，无因次量。cnlpor(compensatedneutronloggingporosity)表示补偿中子孔隙度，本发明中均简称为中子孔隙度，无因次量。denacpor(densityacousticporosity)是指密度孔隙度减去声波时差孔隙度的差值，无因次量。本发明中其他部分提及到的与上述缩写相同的缩写，其所表达的意思相同。图4a、图5a、图6a、图7a的测井曲线图中的rd和rs均采用对数坐标。图4c、图5b--图5e、图6b--图6e、图7b、图7c中右下角所示的数字表示该样品的放大比率。
图1为本发明识别钙屑砂岩储层孔隙类型的方法流程图，如图1所示，选取钙屑砂岩发育段所在的目的层段，所述目的层段为一个砂组或一个单砂体，一个砂组中至少包含一个单砂体，优选的，选取钙屑砂岩发育段所在的砂组为目的层段，选取目的层段中的钙屑砂岩干层段，获取所述目的层段所属井的声波时差孔隙度曲线和密度孔隙度曲线，并将所述目的层段所属井的声波时差孔隙度曲线、密度孔隙度曲线与所述目的层段所属井的孔隙度较正曲线一起放置在同一测井曲线道内；以所述目的层段所属井的孔隙度较正曲线为基准，深度坐标不变，在孔隙度坐标方向上将所述目的层段所属井的声波时差孔隙度曲线、密度孔隙度曲线平移，使所述目的层段所属井的密度孔隙度曲线、声波时差孔隙度曲线与所述目的层段所属井的孔隙度较正曲线在所述目的层段中的钙屑砂岩干层段重合；设置孔隙类型的占比；三条孔隙度曲线在所述目的层段中的钙屑砂岩干层段重合后，根据同一口井的测井曲线上，相同深度位置的密度孔隙度和声波时差孔隙度的相差程度与该深度位置处钙屑砂岩储层中满足预设占比条件的孔隙类型的对应关系识别所述目的层段所属井的钙屑砂岩储层中满足预设占比条件的孔隙类型。该方法利用孔隙度测井曲线来识别钙屑砂岩储层中满足预设占比条件的孔隙类型，识别过程中不需要采集岩心样品和岩屑样品，因而没有由于取心成本限制而使得取心少；由于非均质性强致使实际取心与设计取心存在误差，取心代表性差而引起的所取样品不能真正、全面反映储层孔隙类型的问题；以及，岩屑的选取工作量太大，导致储层孔隙类型在纵向和横向上的分布特征难以准确描述的问题。该方法可操作性强，成本低，符合测井学原理，经中石化元坝须家河组须三段钙屑砂岩气藏应用证实，该方法可以方便、快捷并准确识别出不同的孔隙类型，判断标准清晰，易于实施，对钙屑砂岩储层孔隙类型的识别有很好的适应性，为钙屑砂岩储层的认识、综合评价奠定了良好的基础。
优选的，利用钙屑砂岩岩性识别标准识别钙屑砂岩发育段，进而选取钙屑砂岩发育段所在的目的层段。该钙屑砂岩岩性识别标准通过实验总结而来，图2为总结钙屑砂岩岩性识别标准的方法流程图，图4为元坝须三段元坝6井中不同岩性的测井响应特征，如图2、图4所示，选取元坝须三段钙屑砂岩气藏中取心少且取心代表性差的已知钙屑砂岩发育段所在的目的层段的岩屑包，从中采集岩屑样品，经磨片制作后，连同岩心磨片样品，在显微镜下详细观察，鉴定样品岩性。优选的，为了提高钙屑砂岩岩性识别标准的准确性，还可以选取其他砂组一起实验。具体地，以元坝6井的三砂组为已知钙屑砂岩发育段所在的目的层段，结合元坝6井的二砂组一起实验，采集元坝6井二砂组、三砂组中的岩心、岩屑样品，经制样后在显微镜下详细观察，鉴定样品岩性。根据元坝6井二砂组和三砂组的样品岩性鉴定结果及样品所在深度，详细标定元坝6井的常规测井曲线以读取样品所在深度对应的常规测井曲线上的值，在深度方向上所取样品越密集，则钙屑砂岩岩性识别标准的准确度越高。对反映岩性的敏感测井曲线进行分析，由于研究区范围内，总体上岩性属于低渗致密储集岩，孔隙及流体信息在测井曲线上响应微弱，电阻率值受岩性的成分和粒度的影响较大，因此，电阻率曲线是识别钙屑砂岩等岩性的关键曲线之一，根据标定结果，总结钙屑砂岩岩性识别标准如表1所示。
表1
图4示例性地给出了其中几个样品的标定，如图4b、图4e所示分别为图4a中位置a4和位置d4处的样品，也就是元坝6井4267.2m和4301.2m处的取心样品，该样品为砾岩或砂砾岩。图4a为元坝6井测井曲线图，所示位置a4和位置d4所对应的常规测井曲线上，gr小于40api，rd大于4000欧姆·米。如图4c、图4d所示分别为图4a中位置b4和位置c4处的样品，也就是元坝6井4295.2m和4297m处的取心样品，该样品为钙屑砂岩。图4a所示位置b4和位置c4所对应的常规测井曲线上，gr小于55api，300欧姆·米≤rd≤4000欧姆·米。将钙屑砂岩岩性识别标准应用于元陆702井、元陆11井和元陆7井，所得结果与相应样品的观测鉴定结果吻合。利用上述钙屑砂岩岩性识别标准识别钙屑砂岩发育段，操作简单，识别准确度高，可利用常规测井曲线快速识别钙屑砂岩发育段。利用上述钙屑砂岩岩性识别标准识别钙屑砂岩发育段并不构成对本发明的限制，现有技术中只要能识别出钙屑砂岩发育段的技术手段均可应用到本发明中，如利用钙屑砂岩元素识别法识别钙屑砂岩发育段。
优选的，同一口井的测井曲线上，相同深度位置的密度孔隙度和声波时差孔隙度的相差程度与该深度位置处钙屑砂岩储层中满足预设占比条件的孔隙类型的对应关系通过实验总结获得，如图3所示为总结该对应关系的实验方法流程图。选取元坝须三段钙屑砂岩气藏中取心少且取心代表性差的已知钙屑砂岩发育段所在的目的层段的岩屑包，从中采集岩屑样品，经磨片制作后，连同岩心磨片样品，在显微镜下详细观察，选取孔隙类型符合预设占比条件的样品。用于观察孔隙类型的已知目的层段和用于岩性鉴定的已知目的层段可以是相同的目的层段，也可以是不同的目的层段。本实施例中，用于观察孔隙类型的已知目的层段和用于岩性鉴定的已知目的层段采用不同的目的层段，具体地，本实施例中用于观察孔隙类型的已知目的层段为元坝2井的三砂组，图5为元坝2井中孔隙类型符合预设占比条件的样品及其孔隙度测井曲线响应特征。具体地，设定孔隙类型占比大于50％。如图3、图5所示，采集元坝2井三砂组的岩心、岩屑样品，经制样后在显微镜下详细观察，选取孔隙类型占比大于50％的样品，具体地，选取粒间溶蚀孔隙大于50％的样品、粒间溶蚀充填孔隙大于50％的样品、微孔隙大于50％的样品。
如图5a所示，根据元坝2井的测井资料可知，位置a5附近层段的深测向电阻率偏高，大于4000欧姆·米，为砾岩段。如图3、图5所示，根据元坝2井的测井资料，选取元坝2井的三砂组中深侧向电阻率的范围为1000欧姆·米<rd<4000欧姆·米，密度大于2.65g/cm3，且补偿中子测井cnl小于5％的层段为钙屑砂岩干层段，该钙屑砂岩干层段为深度处于4371.9m--4372.2m的深度段。经薄片样品观察证实位于4371.9m--4372.2m的深度段为钙屑砂岩干层段，图5b示例性地展示了钙屑砂岩干层段中位于4372m深度处的样品，也即图5a中位置b5处的样品，该样品显示为钙屑砂岩干层，孔隙不发育。读取4371.9m--4372.2m的深度段各深度点对应的密度值和声波时差值，作为元坝2井钙屑砂岩的密度骨架值和声波时差骨架值，利用密度孔隙度解释模型和密度骨架值计算元坝2井整个深度段的密度孔隙度，利用声波时差孔隙度解释模型和声波时差骨架值计算元坝2井整个深度段的声波时差孔隙度，进而得到元坝2井整个深度段的密度孔隙度曲线和声波时差孔隙度曲线。将元坝2井的密度孔隙度曲线、声波时差孔隙度曲线和孔隙度较正曲线置于同一测井曲线道内，形成孔隙度测井曲线组合，优选的，孔隙度较正曲线为中子孔隙度曲线。以元坝2井的中子孔隙度曲线为基准，深度坐标不变，在孔隙度坐标方向上将元坝2井的密度孔隙度曲线和声波时差孔隙度曲线平移，使元坝2井的密度孔隙度曲线、声波时差孔隙度曲线和中子孔隙度曲线在选取的元坝2井的钙屑砂岩干层段，也就是位于4371.9m--4372.2m的深度段处重合，以完成密度孔隙度曲线和声波孔隙度曲线的较正。上述优选方案中，孔隙度较正曲线为中子孔隙度曲线并不构成对本发明的限制，本领域的技术人员还能选取岩心孔隙度曲线作为孔隙度较正曲线。
利用选取的孔隙类型占比大于50％的样品和该样品所在深度详细标定元坝2井的孔隙度测井曲线中经重合步骤后的密度孔隙度曲线和声波时差孔隙度曲线，读取该样品所在深度对应的密度孔隙度和声波时差孔隙度，深度方向上采集的样品越密集则总结的结果越准确。具体地，读取粒间溶蚀孔隙大于50％的样品所在深度对应的密度孔隙度和声波时差孔隙度，读取粒间溶蚀充填孔隙大于50％的样品所在深度对应的密度孔隙度和声波时差孔隙度，读取微孔隙大于50％的样品所在深度对应的密度孔隙度和声波时差孔隙度，根据读取结果分别总结当粒间溶蚀孔隙大于50％、粒间溶蚀充填孔隙大于50％和微孔隙大于50％时所对应的元坝2井中相同位置密度孔隙度与声波时差孔隙度的相差程度，优选的，密度孔隙度与声波时差孔隙度的相差程度是指密度孔隙度减去声波时差孔隙度的差值。图5a为元坝2井的测井曲线图，示出了元坝2井三砂组中的部分孔隙类型符合占比条件的层段，其中位于4372.2m-4372.5m的深度段占比大于50％的孔隙类型为微孔隙，图5c示例性地展示了该深度段中位于4372.3m处的样品，也就是图5a中位置c5处的样品，该样品中，微孔隙占比大于50％呈浸染状分布，该深度段同一深度位置密度孔隙度减去声波孔隙度的差值的最大值为0.48％。位于4376.5m-4377m的深度段占比大于50％的孔隙类型为粒间溶蚀充填孔隙，图5d示例性地展示了该深度段中位于4376.5m处的样品，也就是图5a中位置d5处的样品，该样品显示为粒间溶蚀充填孔隙占比大于50％，粒间溶蚀孔隙中有高岭石(一种粘土矿物)充填，该深度段同一深度位置密度孔隙度减去声波孔隙度的差值的最小值为0.5％，最大值为2.73％，且随着同一深度位置密度孔隙度减去声波孔隙度的差值的减小溶蚀充填孔的数量也相应减少。位于4379.3m-4380.3m深度段占比大于50％的孔隙类型为粒间溶蚀孔隙，图5e示例性地展示出了该深度段中位于4380m处的样品，也就是图5a中位置e5处的样品，该样品显示为粒间溶蚀孔隙占比大于50％，该深度段同一深度位置密度孔隙度减去声波孔隙度的差值的最小值>3％，最大值为5％。
根据孔隙类型占比>50％的样品所在深度对应的密度孔隙度和声波孔隙度的读取结果，总结得到同一口井的测井曲线上，相同深度位置的密度孔隙度和声波时差孔隙度的相差程度与该深度位置处钙屑砂岩储层中满足预设占比条件的孔隙类型的对应关系为：
当同一深度位置的密度孔隙度减去声波时差孔隙度的差值＞3％时，该深度位置的储层中占比大于50％的孔隙类型为粒间溶蚀孔隙；
当同一深度位置的密度孔隙度减去声波时差孔隙度的差值≥0.5％，且≤3％时，该深度位置的储层中占比大于50％的孔隙类型为粒间溶蚀充填孔隙；
当同一深度位置的密度孔隙度减去声波时差孔隙度的差值＜0.5％时，该深度位置的储层中占比大于50％的孔隙类型为微孔隙。
上述对应关系是基于如下测井学原理，储层中的粒间溶蚀孔隙、粒间溶蚀充填孔隙和微孔隙均为次生孔隙，密度孔隙度反映的是总孔隙度，声波时差孔隙度反映的是原生孔隙度，次生孔隙度＝密度孔隙度(总孔隙度)-声波时差孔隙度(原生孔隙度)，以次生孔隙度的大小来确定钙屑砂岩储层中符合占比的条件的孔隙类型。
在钙屑砂岩岩性识别标准，以及同一口井的测井曲线上，相同深度位置的密度孔隙度和声波时差孔隙度的相差程度与该深度位置处钙屑砂岩储层中满足预设占比条件的孔隙类型的对应关系确立后，将结果直接应用于钙屑砂岩储层孔隙类型的识别。如图1、图6所示，其中图6a为元陆20井的测井曲线图，选取元陆20井中gr小于55api，300欧姆·米≤rd≤4000欧姆·米的深度段所在的三砂组为钙屑砂岩发育段所在的目的层。元陆20井三砂组中位于4135.4m-4135.7m的深度段的深侧向电阻率的范围为1000欧姆·米<rd<4000欧姆·米，密度大于2.65g/cm3，且补偿中子测井cnl小于5％，故选取该深度段为钙屑砂岩干层段，也就是图6a中位置a6所在的层段。读取位置a6所在的层段中各深度点对应的密度值和声波时差值作为元陆20井钙屑砂岩的密度骨架值和声波时差骨架值，并通过密度孔隙度解释模型和密度骨架值计算元陆20井整个深度段的密度孔隙度，通过声波时差孔隙度解释模型和声波时差孔骨架值计算元陆20井整个深度段的声波时差孔隙度，进而获得元陆20井整个深度段的密度孔隙度曲线和声波时差孔隙度曲线，将元陆20井的密度孔隙度曲线、声波时差孔隙度曲线与孔隙度较正曲线置于同一测井曲线道内，优选的，孔隙度较正曲线为中子孔隙度曲线。以元陆20井中子孔隙度曲线为基准，将元陆20井的密度孔隙度曲线、声波时差孔隙度曲线平移，使元陆20井的密度孔隙度曲线、声波时差孔隙度曲线和中子孔隙度曲线在位于4135.4m-4135.7m的钙屑砂岩干层段，也就是位置a6所在层段重合，以完成密度孔隙度曲线和声波孔隙度曲线的较正。读取元陆20井同一深度点的密度孔隙度和声波孔隙度，并计算密度孔隙度减去声波孔隙度的差值，设置孔隙类型占比>50％，根据已确立的同一口井的测井曲线上，相同深度位置的密度孔隙度和声波时差孔隙度的相差程度与该深度位置处钙屑砂岩储层中满足预设占比条件的孔隙类型的对应关系识别占比>50％的孔隙类型。具体地如图6a所示，位于4141.5m-4142.2m的深度段，即位置d6所在深度段，该深度段内，同一深度点对应的密度孔隙度减去声波时差孔隙度的差值均在0.95％-1.95％之间，该差值处于0.5％-3％的范围内，故该深度段中占比大于50％的孔隙类型为粒间溶蚀充填孔隙。位于4137m-4137.5m的深度段，即位置b6所在的深度段，该深度段内，同一深度点对应的密度孔隙度减去声波时差孔隙度的差值均在0.5％-0.85％之间，该差值处于0.5％-3％的范围内，故该深度段中占比大于50％的孔隙类型为粒间溶蚀充填孔隙，且位置b6所在深度段内同一深度点对应的密度孔隙度减去声波时差孔隙度的差值比位置d6所在深度段内同一深度点对应的密度孔隙度减去声波时差孔隙度的差值要小，所以位置b6所在深度段内的粒间溶蚀充填孔隙的数量小于位置d6所深度段内的粒间溶蚀充填孔隙的数量。位于4140.1m-4141.3m的深度段，即位置c6所在的深度段，该深度段内，同一深度点对应的密度孔隙度减去声波时差孔隙度的差值均<0.5％，故该深度段中占比大于50％的孔隙类型为微孔隙。
将利用本发明所提供的钙屑砂岩储层孔隙类型识别方法对元陆20井进行评估所得到的结果与元陆20井相应深度的样品进行比对，上述评价结果与样品显示的孔隙类型相符，图6b-图6e给出了其中的部分对比样品，图6b为位置a6，也就是深度点4135.5m处的样品，该样品显示为钙屑砂岩干层，孔隙不发育，图6c为位置b6，也就是深度点4137.1m处的样品，该样品显示为粒间溶蚀充填孔隙占比大于50％，粒间溶蚀孔隙中有高岭石(一种粘土矿物)充填。图6d为位置c6，也就是深度点4140.25处的样品，该样品显示为微孔隙占比大于50％。图6e为位置d6，也就是深度点4141.8m处的样品，该样品显示为粒间溶蚀充填孔隙大于50％，粒间溶蚀孔隙中有高岭石充填。
采用与元陆20井完全相同的方法步骤对元陆7井进行评估，识别其钙屑砂岩储层的孔隙类型。图7a为元陆7井的测井曲线图，如图7a所示，位于3467.7m-3468.9m的深度段，也就是位置a7所在的深度段为钙屑砂岩干层段。位于3462.5m-3463m的深度段，也就是位置b7所在深度段，该深度段内，同一深度点的密度孔隙度减去声波时差孔隙度的差值介于0.53％-2.38％之间，故该深度段内占比大于50％的孔隙类型为粒间溶蚀充填孔隙。位于3465.5m-3466.5m的深度段，也就是位置c7所在深度段，该深度段内，同一深度点的密度孔隙度减去声波时差孔隙度的差值介于1.7％-2.87％之间，故该深度段内占比大于50％的孔隙类型为粒间溶蚀充填孔隙。位于3464.8m-3465.4m的深度段，也就是位置e7所在深度段，该深度段内，同一深度点的密度孔隙度减去声波时差孔隙度的差值大于3％，故该深度段内占比大于50％的孔隙类型为粒间溶蚀孔隙，同样也能判断出3469.4m-3470.2m的深度段，也就是位置d7所在深度段内占比大于50％的孔隙类型为粒间溶蚀孔隙，这与图7a中3461m-3471m测试段的测试结果能获得120万方/天高产工业气流相吻合，表明3461m-3471m的深度段内有相对集中的粒间溶蚀孔隙的存在，也就是3461m-3471m的深度段中，存在某一深度段其粒间溶蚀孔隙的占比大于50％。图7b、7c分别给出了位置b7，位置c7处的样品，样品显示占比大于50％的孔隙类型为粒间溶蚀充填孔隙，在粒间溶蚀孔隙中充填了高岭石。样品显示结果也与评价结果相符。
上述实施例是本发明优选的实施方式，本发明中密度孔隙度与声波时差孔隙度的相差程度并不局限于密度孔隙度减去声波孔隙度的差值，本领域的技术人员还能选择由密度孔隙度减去声波孔隙度的差值所衍生的其他值作为密度孔隙度与声波时差孔隙度的相差程度，例如，密度孔隙度与声波孔隙度的差值的绝对值，或者密度孔隙度减去声波孔隙度的差值的倒数，或者密度孔隙度与声波时差孔隙度的差值的平方，并同样利用实施例中样品观察和样品标定的方式总结出相应的符合占比条件的孔隙类型和密度孔隙度与声波孔隙度相差程度的对应关系。本发明并不局限于孔隙类型占比大于50％，也可以是其他比例，如75％，并同样利用实施例中样品观察和样品标定的方式总结出占比大于75％的孔隙类型和密度孔隙度与声波孔隙度相差程度的对应关系。
虽然已经参考优选实施例对本发明进行了描述，但在不脱离本发明的范围的情况下，可以对其进行各种改进并且可以对其中部分或者全部技术特征进行等同替换。尤其是，只要不存在逻辑或结构冲突，各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例，而是包括落入权利要求的范围内的所有技术方案。