1.一种砂岩型铀矿地震勘探方法，其特征在于，所述方法包括： 获取工区的叠后地震数据集； 根据所述叠后地震数据集反演得到地震波阻抗数据集； 利用反演得到的地震波阻抗数据集构建砂岩判别因子模型； 基于所述砂岩判别因子模型进行砂岩地层预测； 其中，根据所述叠后地震数据集反演得到地震波阻抗数据集包括： 基于稀疏脉冲反演方法，利用所述叠后地震数据集反演得到地震波阻抗数据集； 所述稀疏脉冲反演是基于地震数据利用波阻抗趋势进行约束，计算过程如下： 叠后地震数据的计算公式为： S＝R*W+N (1) 其中，S为叠后地震数据，R为反射系数，W为地震子波，N为噪音； 地震波阻抗的计算公式为： 其中，AIi为第i层的波阻抗值，Ri为第i层的反射系数； 建立目标函数： F＝Lp(R)+λLq(S-Syn)+αL1(ΔAI) (3) 其中，F为目标函数值，L为模运算，p、q为模因子，S为地震数据，Syn为合成地震道，λ为残差权重因子，α为趋势权重因子，ΔAI为阻抗趋势差； 求解目标函数最终得到地震波阻抗数据集； 其中，利用反演得到的地震波阻抗数据集构建砂岩判别因子模型包括： 利用反演得到的地震波阻抗数据集构建砂岩判别因子模型，计算公式为： 其中，D表示砂岩判别因子，AI表示地震波阻抗，f表示频率，单位为Hz，表示地震波阻抗一阶导数，k0和k1表示加权系数； 在计算过程中，地震波阻抗频域的一阶导数根据泰勒展开式进行计算，计算公式为： 其中，f1表示地震低频端频率，单位为Hz，f2表示地震主频，单位为Hz，表示频率为f1的单频波阻抗数据，/>表示频率为f2的单频波阻抗数据； 根据研究工区地震资料频谱分析结果确定f1＝10Hz，f2＝40Hz，并且根据研究工区经验参数确定k0＝1.2，k1＝156.7， 单频波阻抗数据集通过短时傅里叶变换由地震波阻抗数据获得。

2.根据权利要求1所述的方法，其特征在于，所述基于所述砂岩判别因子模型进行砂岩地层预测包括： 利用直方图确定砂岩地层对应的判别因子值域范围； 在所述砂岩判别因子模型中依据判别因子值域范围圈定砂岩地层分布范围。

3.一种砂岩型铀矿地震勘探装置，其特征在于，所述装置包括： 地震数据输入模块，用于获取工区的叠后地震数据集； 反演模块，用于根据所述叠后地震数据集反演得到地震波阻抗数据集； 砂岩判别因子构建模块，用于利用反演得到的地震波阻抗数据集构建砂岩判别因子模型； 预测模块，用于基于所述砂岩判别因子模型进行砂岩地层预测； 其中，根据所述叠后地震数据集反演得到地震波阻抗数据集包括： 基于稀疏脉冲反演方法，利用所述叠后地震数据集反演得到地震波阻抗数据集； 所述稀疏脉冲反演是基于地震数据利用波阻抗趋势进行约束，计算过程如下： 叠后地震数据的计算公式为： S＝R*W+N (1) 其中，S为叠后地震数据，R为反射系数，W为地震子波，N为噪音； 地震波阻抗的计算公式为： 其中，AIi为第i层的波阻抗值，Ri为第i层的反射系数； 建立目标函数： F＝Lp(R)+λLq(S-Syn)+αL1(ΔAI) (3) 其中，F为目标函数值，L为模运算，p、q为模因子，S为地震数据，Syn为合成地震道，λ为残差权重因子，α为趋势权重因子，ΔAI为阻抗趋势差； 求解目标函数最终得到地震波阻抗数据集； 其中，利用反演得到的地震波阻抗数据集构建砂岩判别因子模型包括： 利用反演得到的地震波阻抗数据集构建砂岩判别因子模型，计算公式为： 其中，D表示砂岩判别因子，AI表示地震波阻抗，f表示频率，单位为Hz，表示地震波阻抗一阶导数，k0和k1表示加权系数； 在计算过程中，地震波阻抗频域的一阶导数根据泰勒展开式进行计算，计算公式为： 其中，f1表示地震低频端频率，单位为Hz，f2表示地震主频，单位为Hz，表示频率为f1的单频波阻抗数据，/>表示频率为f2的单频波阻抗数据； 根据研究工区地震资料频谱分析结果确定f1＝10Hz，f2＝40Hz，并且根据研究工区经验参数确定k0＝1.2，k1＝156.7， 单频波阻抗数据集通过短时傅里叶变换由地震波阻抗数据获得。

4.根据权利要求3所述的装置，其特征在于，所述预测模块包括： 值域确定单元，用于利用直方图确定砂岩地层对应的判别因子值域范围； 选择单元，用于在所述砂岩判别因子模型中依据判别因子值域范围圈定砂岩地层分布范围。

请联系平台

请联系平台