技术领域
[0001] 本发明涉及石油天然气勘探与开发领域，特别是一种砾岩岩心颗粒结构中粒度的分析方法。
相关背景技术
[0002] 碎屑岩结构是指碎屑岩内各结构组分的特点和相互关系，包括颗粒结构、杂基和胶结物结构、孔隙结构、碎屑颗粒与杂基和胶结物之间关系。碎屑岩颗粒结构包括粒度、圆度、球度、形状、分选性以及颗粒表面结构。碎屑岩颗粒结构中粒度特征（大小、分布）是碎屑岩分类命名的依据，也是碎屑岩颗粒结构中分选性求取的基础。碎屑岩颗粒结构中粒度特征及分选性是分析碎屑岩沉积时自然地理环境以及水动力条件的重要标志。同时，碎屑岩颗粒结构中粒度特征及分选性对碎屑岩埋藏演化过程中储集物性演化也具有重要的控制作用。因此，碎屑岩颗粒结构的粒度分析是油气勘探开发领域中的重要研究内容之一。
[0003] 目前，碎屑岩粒度分析方法主要有直接测量法、筛析法、沉降法、激光粒度分析法及图像分析法。直接测量法是利用直尺或测规直接测量松散沉积物中直径大于2mm的砾石的直径，并分析其特征的粒度分析方法。筛析法是分析细砾级和砂级松散沉积物粒度特征的主要方法，通常是取50g以上样品，在震筛机上筛约10min，然后分级称重。沉降法的基本原理是利用颗粒沉降速度来划分粒度较细的粉砂和粘土的粒级分布。激光粒度分析法是根据激光照射到颗粒后，颗粒能使激光产生衍射或散射的现象来测试粒径范围在0.02-2mm的粒度分布。图像分析法包括薄片图像分析和岩心扫描图像分析，薄片图像分析是指将显微镜下的图像采集到计算机中，在计算机上对颗粒的二维图像进行测量统计处理，分析碎屑岩中粒度直径小于2mm的粒度分布特征；岩心扫描图像分析是通过选取一定长度（10-50cm）的岩心样品，运用岩心剖切、抛光和图像扫描技术，形成高分辨率的岩心彩色扫描图像，利用先进的图像分析方法提取图像中的砾石（直径大于2mm），获取单个砾石的粒径、形态、磨圆等数据，统计砾石的分选、分布、定向等参数。
[0004] 砾岩是指直径大于2mm的碎屑颗粒含量大于50%的碎屑岩。砾岩储层是石油、天然气重要的储层之一。但是，目前已有的碎屑岩粒度分析方法不能够准确表征砾岩的粒度特征。直接测量法和筛析法仅适用于松散沉积物，而对于已固结成岩的砾岩，若采用机械破碎，容易使颗粒发生破碎或破碎不完全，导致分析结果不准确；若采用化学试剂溶解填隙物的方法，会使颗粒同时发生溶解，导致分析结果不准确；激光粒度分析法、薄片图像分析法仅适用于直径小于2mm的颗粒；岩心扫描图像分析法仅适用于分析砾岩中直径大于2mm的砾石。因此，为了准确的表征砾岩颗粒结构中粒度的特征，用以分析碎屑岩沉积时自然地理环境以及水动力条件，以及控制碎屑岩埋藏演化过程中储集物性演化，需要建立适用于砾岩颗粒结构中粒度的分析方法。
具体实施方式
[0022] 以胜利油田某地区砾岩颗粒结构求取方法来说明该发明的具体技术方案。
[0023] 一种砾岩岩心颗粒结构中粒度的分析方法，具体步骤如下：
[0024] （1）对于砾岩岩心中直径大于2mm的砾石颗粒，借助透明纸，采用人工素描的方法对岩心1：1精描，然后利用CorelDRAW X4SP2绘图软件进行清绘，获取砾级颗粒图像（如图5所示），再通过CorelDRAWX4SP2图像分析软件求取每个砾石的面积SAi及其等面积圆直径Di、岩心总面积S，将每个砾石的面积求和得砾石的总面积SA。表1所示为某砾岩岩心砾石颗粒1:1精描分析粒度数据表。
[0025] 砾石的总面积计算公式为：
[0026]
[0027] 式中，SAi-指岩心中每个砾石颗粒的面积；n-指岩心中砾石颗粒的个数。
[0028] 表1
[0029]
[0030]
[0031] （2）对于砾岩中直径0.03-2mm的颗粒，从岩心中选择代表性区域钻取样品，磨制岩石薄片，利用蔡司AxioscopeA1APOL.数字透反偏光显微镜及摄像系统获取0.03-2mm的颗粒图像，然后通过Axio Vision Software Rel.图像分析软件求取薄片中每个颗粒的面积Sbi及其等面积圆直径di、薄片总面积Sc（如图6所示），将每个颗粒的面积求和得薄片中颗粒的总面积Sb；然后，计算砾岩岩心中每个直径-等面积圆直径为0.03-2mm的颗粒所对应的面积SBi及岩心中0.03-2mm颗粒的总面积SB。表2所示为某砾岩岩心0.03-2mm颗粒薄片图像分析粒度数据表。
[0032] 岩心中0.03-2mm颗粒的总面积SB计算公式为：
[0033]
[0034] 式中，Sbi-指薄片中每个0.03-2mm颗粒的面积；n-指薄片中砂级颗粒的个数。
[0035] SBi=Sbi/Sc×（S-SA），i=1，2，3……n（3）
[0036] 式中，SBi-砾岩岩心中每个等面积圆直径为0.03-2mm的颗粒所对应的面积。
[0037]
[0038] 式中，n-指岩心中0.03-2mm颗粒的个数；
[0039] 表2
[0040]
[0041]
[0042] （3）将砾岩岩心宏观精描图像分析与典型区域微观薄片图像分析的分析结果相结合，计算砾岩岩心中颗粒总面积Sg以及每个等面积圆直径对应颗粒面积占颗粒总面积的百分含量（如表3所示）。
[0043] 计算公式如下：
[0044] 砾岩岩心中颗粒总面积：
[0045] Sg=SA+SB （5）
[0046] 砾岩岩心中每个等面积圆直径的砾石占颗粒总面积的百分含量：
[0047] PAi=SAi/Sg×100% （6）
[0048] 砾岩岩心中每个等面积圆直径2-0.03mm的颗粒占颗粒总面积的百分含量：
[0049] PBi=SBi/Sg×100% （7）
[0050] 表3所示为某砾岩岩心各等面积圆直径对应颗粒面积占颗粒总面积百分含量计算表。
[0051] 表3
[0052]
[0053]
[0054] 表3中：φ=-log2D，D为砾岩中颗粒等面积圆直径（单位mm）。
[0055] （4）利用砾岩岩心中每个等面积圆直径对应颗粒面积占颗粒总面积的百分含量数据（如表3所示），通过Grapher绘图软件绘制砾岩岩心粒度概率累积曲线图（如图7所示），求取砾岩颗粒的结构参数：Mz=-2.57φ，So=3.08，其中Mz为平均粒径，So为分选系数。
[0056] 表4所示为本发明求取的颗粒结构参数与现有分析法求取的颗粒结构参数对比表。
[0057] 表4
[0058]粒度结构参数 岩心扫描图像分析法 薄片图像分析法 本发明方法
Mz/φ -3.83 1.3 -2.85
So 1.8 1.6 2.83
[0059] 对于同一块砾岩岩心，图2为利用岩心扫描图像分析法获得的砾岩粒度概率累积曲线，图3为薄片图像分析法获取的砾岩粒度概率累积曲线，图4是本发明方法获取的砾岩粒度概率累积曲线，表4是利用三种方法获取的砾岩粒度概率累积曲线求取的分选系数（So）、平均粒径（Mz）。对比分析可知，利用三种方法所获取的砾岩粒度直径范围、粒度概率累积曲线形态、分选系数（So）、平均粒径（Mz）存在明显差异，并且本发明方法获取的砾岩粒度直径范围及粒度概率累积曲线形态、分选系数（So）、平均粒径（Mz）更符合理论上砾岩的粒度特征。
