[发明专利]基于统计火山矿物含量的深时火山活动强度重建方法

权利要求书

1.基于统计火山矿物含量的深时火山活动强度重建方法，其特征在于，包括以下步骤：

S1：利用某一地质历史时期单位面积火山矿物的含量百分比来代表该时刻火山活动的强度；其中，采用Glagolev-Chayes统计方法，通过网格进行多点自动计数，获得火山矿物颗粒的点数来代表其相对含量；

S2：通过建立某一段地质历史时期火山矿物的含量变化序列，来表征火山活动强度的演化。

2.根据权利要求1所述的基于统计火山矿物含量的深时火山活动强度重建方法，其特征在于，步骤S1中还包括样品采集与制片，具体如下，

按照需求的采样密度对钻孔岩石样本进行采样，并完成普通岩石薄片制作；

对准噶尔盆地石炭系火山型沉积岩按照0.3-0.5m采样间距进行高分辨率样品采集，并标记好样本的地层顶底；按照平行于顶底面的方式进行切片，制作岩石微观薄片3-5片，为后续显微镜下的定量统计做准备；西北缘研究样品采集于新疆太勒古拉剖面、201剖面、贝拉舍剖面、恰尔奈也剖面以及柳树沟剖面，南缘样品采集于新疆庙尔沟剖面；共采集新鲜样品3500个，其中西北缘1000个，南缘2500个，且样品表面清洗干净，将以上样品磨样并进行制片。

3.根据权利要求2所述的基于统计火山矿物含量的深时火山活动强度重建方法，其特征在于，步骤S1中还包括单样本图像采集数量参数设定，具体如下，

确定单个样本采集微观组分图像的数量，将统计结果受颗粒分选影响程度降到最低，保障统计结果的置信度达到95％以上；

将碎屑统计认为是抽样统计，对于随机分布的总体，按照绝对精度决定样本量n的公式为：

其中α为显著性水平；是N(0，1)分布的分位数；S²为总体方差，S²≈P(1-p)；d为绝对精度(绝对误差)；N为样品总数，P为颗粒百分比；当α取0.05时，达到置信度为1-α＝95％，此时μ_1-α/2＝1.96；当P取0.5时，有最大总体方差0.25，d为绝对精度即绝对误差。

4.根据权利要求3所述的基于统计火山矿物含量的深时火山活动强度重建方法，其特征在于，步骤S1中还包括岩石微观组分图像采集，具体如下，

通过显微镜对制作的岩石微观薄片进行图像采集，每张薄片采集不同放大倍数下图像合计100张，同一岩石岩本合计400张。

5.根据权利要求4所述的基于统计火山矿物含量的深时火山活动强度重建方法，其特征在于，步骤S1中还包括基于灰度火山矿物自动识别参数设定，具体如下，

在电子显微镜下，通过对火山矿物的鉴定、提取火山矿物的灰度参数信息，为火山矿物自动识别做好参数准备；

微观薄片观察表明，研究区准噶尔盆地石炭系岩浆岩以火山岩为主，基于大量样本，通过电子显微镜灰度识别功能，确定火山矿物的灰度范围值。

6.根据权利要求5所述的基于统计火山矿物含量的深时火山活动强度重建方法，其特征在于，步骤S1中还包括基于粒度分析的网格参数设定，具体如下，

通过火山矿物组分粒径分布特征，来确定显微镜中图像识别系统的网格参数；在显微镜下，通过统计一定数量样本的火山矿物的粒度参数，取样本中累计频率分布在10％的最小粒度的平均值作为网格长度，依据矿物粒度的平均粒径大小，可分类进行参数设定。

7.根据权利要求6所述的基于统计火山矿物含量的深时火山活动强度重建方法，其特征在于，通过显微镜对制作的岩石微观薄片进行图像采集时，还包括以下步骤：

识别图像数据的摄像角度信息，记为实时摄像角度；

获取采集图像数据时岩石微观薄片所在环境的照明信息，记为实时照明强度；

根据图像采集过程得到最优图像数据，设定标准摄像角度和标准照明强度；

若所述实时摄像角度与所述标准摄像角度不匹配、或所述实时照明强度与所述标准照明强度不匹配，则自动舍弃对应的图像数据。

8.计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有一个或多个计算机程序，所述一个或多个计算机程序被一个或多个处理器执行时实现如权利要求1-7任一项所述的基于统计火山矿物含量的深时火山活动强度重建方法。