[发明专利]基于石榴子石热红外特征的光谱地质变质压力分析方法、系统以及存储介质

权利要求书

1.一种基于石榴子石热红外特征的光谱地质变质压力分析方法,其特征在于，包括：

获取发育石榴子石变质带岩石样本中石榴子石的热红外波谱信息；

提取10-13μm间石榴子石第一主吸收谷波长λ；

根据λ与石榴子石晶胞参数a轴值的经验关系计算所测样本中石榴子石的晶胞参数a₀值；

基于晶胞参数a₀值作变质带的相对变质压力分析。

2.根据权利要求1所述的一种基于石榴子石热红外特征的光谱地质变质压力分析方法，其特征在于，获取发育石榴子石变质带岩石样本中石榴子石的热红外波谱包括：

定位到发育石榴子石变质带岩石样本中的石榴子石单晶，并检测获取波谱数据；

判断上述波谱数据在10-13μm光谱区间是否具有诊断性双峰式吸收特征，如果是则保留波谱数据，否则重新定位到岩石样本中石榴子石单晶，重新检测获取波谱数据，直到所获取的波谱数据在10-13μm光谱区间具有诊断性双峰式吸收特征，则判断所获取的石榴子石波谱数据满足要求；

重复上述步骤预设次数，且每个样本保留预设条数的波谱数据。

3.根据权利要求2所述的一种基于石榴子石热红外特征的光谱地质变质压力分析方法，其特征在于，定位到发育石榴子石变质带岩石样本中的石榴子石单晶，并检测获取波谱数据的设备为傅里叶变换红外光谱仪，光谱范围包含10-13μm，波谱分辨率优于25nm，采样探头的光斑优于2mm。

4.根据权利要求1所述的一种基于石榴子石热红外特征的光谱地质变质压力分析方法，其特征在于，提取10-13μm间石榴子石第一主吸收谷波长λ包括：

截取10-13μm区间的波谱数据；

通过基于小波变换的峰值检测算法分析计算获取10-13μm区间的波谱数据第一主吸收谷波长。

5.根据权利要求4所述的一种基于石榴子石热红外特征的光谱地质变质压力分析方法，其特征在于，通过基于小波变换的峰值检测算法分析计算获取10-13μm区间的波谱数据第一主吸收谷波长所提及的基于小波变换的峰值检测算法系由scipy模块的scipy.signal.find_peaks_cw函数提供，小波选用scipy.signal.morlet2。

6.根据权利要求1所述的一种基于石榴子石热红外特征的光谱地质变质压力分析方法，其特征在于，根据λ与石榴子石晶胞参数a轴值的经验关系计算所测样本中石榴子石的晶胞参数a₀值包括：

从RRUFF波谱库筛选获取同时具有红外波谱数据和粉末衍射数据的石榴子石族矿物样本，分别读取其在10-13μm光谱区间第一主吸收谷波长λ及其晶胞参数a₀；

从USGS波谱库中筛选获取由Nicolet测量且附有X射线衍射分析数据的石榴子石族矿物样本，分别读取其在10-13μm光谱区间第一主吸收谷波长λ及其晶胞参数a₀；

对在10-13μm光谱区间第一主吸收谷波长λ及晶胞参数a₀进行回归运算得到λ与a₀值的经验关系；

利用上述经验公式对每个变质带岩石样本的多条石榴子石波谱数据分别计算对应的石榴子石晶胞参数值，并取平均值作为样本石榴子石的代表晶胞参数值。

7.根据权利要求6所述的一种基于石榴子石热红外特征的光谱地质变质压力分析方法，其特征在于，λ与a₀值的经验关系如下：

a₀=0.5874·λ+5.2397；

a₀单位为埃米；

λ单位为微米。