新型杂卤石钾盐矿数据累积体积计算方法、系统及存储介质

本发明涉及含钾钾盐矿层勘探，尤其涉及一种新型杂卤石钾盐矿数据累积体积计算方法、系统及存储介质。

背景技术：

1、当前，钾盐是钾肥的主要原料，作为农业大国，钾盐属于中国当前紧缺的战略性矿产资源之一。含钾岩矿藏勘探过程中，如何有效预测矿产资源潜力，实现钾盐资源评价是一个重要的问题。近年来，伴随着“油钾兼探”工作的逐步开展，应用地球物理勘探技术对含钾钾盐矿层的地下分布范围进行预测和描述，成为实现地下目的层钾盐矿资源发现和富集区描述的重要途径。

2、钾盐资源按其物态可分为固体和液体两大类型，其中固体钾盐矿层资源量估算方法与其他固体矿估算类似。常用方法有两种：1、体积法；2、成矿地质体体积法。

3、方法一，体积法是传统固体矿产资源估算方法，主要利用矿体体积、品位、比重来估算矿产资源储量，操作简单易行。该方法基于地质条件均一的简化模型基础上，通过相似类比法进行外推，应用控制区内含矿地质建造的体积与资源总量计算含矿系数，通过计算预测区的体积估算出预测区潜在的资源量。

4、方法二，成矿地质体体积法是传统体积法的改进方法，其基本思想与传统体积法类似，是基于矿床模型综合地质信息进行资源量预测的方法。估算过程首先研究预测区内资料，确定模型区，求其含矿率；然后利用体积法计算公式估算其他预测区的矿产资源量，并依据预测区的综合地质信息情况，确定其与模型区的相似性，乘以模型区含矿系数和相似系数，最终完成各工作区资源量预测评价。

5、以上两种方法都是建立在地下钾盐矿层地质模型化分析基础上，可准确描述钾盐矿层体积的估算技术。对于地下构造明确、矿床分布稳定的矿体比较适合。

6、然而，当矿床分布不稳定时，现有的计算方法难以准确估算，例如在川东北地区嘉陵江地层中发育的新型杂卤石钾盐矿，是一种不同于国外典型固体钾盐矿藏的固体钾盐矿，由于其赋存地层构造复杂，地质成因特殊，矿床分布极不稳定，新型杂卤石钾盐矿在地层中往往呈团块状分布，横向连续性差，厚度差异大，采用稳定地质模型对其体积进行估算误差相对较大。

技术实现思路

1、为克服以上所述技术缺陷，本发明实施例提供一种新型杂卤石钾盐矿数据累积体积计算方法、系统及存储介质，能够提高不稳定分布钾盐矿体体积计算的准确性。

2、本发明一实施例提供了一种新型杂卤石钾盐矿数据累积体积计算方法，包括：

3、获取目标层段中有效钾盐矿层的钾参数和密度对应的有效门槛值；

4、通过地震反演计算得到所述目标层段的所述钾参数和密度的参数值，根据所述有效门槛值对所述钾参数和密度进行过滤，得到有效钾参数的时间域空间分布；其中，所述有效钾参数即所述目标层段中有效钾盐矿层的钾参数，所述有效钾盐矿层即符合品位要求的新型杂卤石钾盐矿层；

5、对所述有效钾参数的时间域的有效钾参数进行时间-深度转换后，得到深度域有效钾参数，通过累加计算得到所述有效钾盐矿层的累积厚度，对所述有效钾盐矿层的累积厚度进行校正，得到所述有效钾盐矿层有效厚度；

6、根据所述有效钾盐矿层的所述深度域有效钾参数空间分布计算所述有效钾盐矿层的平均有效钾参数，并根据数学公式转换得到所述有效钾盐矿层的钾含量系数平面分布值；

7、根据所述有效钾盐矿层有效厚度和所述有效钾盐矿层的钾含量系数平面分布值得到所述目标层段中新型杂卤石钾盐矿的累积体积估算结果。

8、作为上述方案的改进，所述获取目标层段中有效钾盐矿层的钾参数和密度对应的有效门槛值，包括：

9、获取符合期望条件的有效钾盐矿层特征，有效钾盐矿层即符合品位要求的新型杂卤石钾盐矿层，其中，所述有效钾盐矿层特征包括所述钾参数和密度；

10、获取所述钾参数和密度的参考测井曲线，根据所述参考测井曲线得到符合期望条件的所述钾参数和密度对应的有效门槛值。

11、作为上述方案的改进，通过地震反演计算得到所述目标层段的所述钾参数和密度的值，根据所述有效门槛值对所述钾参数和密度进行过滤，得到有效钾参数的时间域空间分布，有效钾参数即目标层段中有效钾盐矿层的钾参数，包括：

12、通过地震反演，得到所述目标层段的所述钾参数和密度的值；

13、对所述参数值不满足所述有效门槛值的所述钾参数和密度进行过滤；

14、保留所述参数值满足所述有效门槛值的所述钾参数和密度，得到有效钾参数的时间域空间分布，其中所述有效钾参数即所述目标层段中有效钾盐矿层的钾参数。

15、作为上述方案的改进，对所述有效钾参数的时间域的有效钾参数进行时间-深度转换后，得到深度域有效钾参数，通过累加计算得到所述有效钾盐矿层的累积厚度，对所述有效钾盐矿层的累积厚度进行校正，得到有效钾盐矿层有效厚度，包括：

16、将所述时间域的有效钾参数数据转换为所述深度域有效钾参数数据；

17、对所述深度域有效参数数据中的非零值赋1值，在所述目标层段进行数值累加计算，并乘以深度域采样间隔，得到所述目标层段中的有效钾盐矿层的累积厚度；

18、应用所述目标层段根据测井解释得到的所述有效钾盐矿层的厚度结果，对所述目标层段中的有效钾盐矿层累积厚度进行标定和校正，得到所述有效钾盐矿层的有效厚度。

19、作为上述方案的改进，所述根据所述有效钾参数计算所述平均有效钾参数平面分布值，并根据数学公式转换得到所述有效钾盐矿层的钾含量系数平面分布值，包括：

20、根据以下数学公式计算所述有效钾盐矿层的钾含量系数平面分布值：

21、vk＝(k-kmin)/(kmax-kmin)

22、式中，vk为所述有效钾盐矿层的钾含量系数平面分布值，k为平均有效钾参数平面分布值，kmin为非钾盐矿层的钾参数，kmax为纯钾盐矿层的钾参数。

23、作为上述方案的改进，所述根据所述新型杂卤石钾盐矿层有效厚度和所述有效钾盐矿层钾含量系数平面分布值得到所述目标层段的新型杂卤石钾盐矿的累积体积估算结果，包括：

24、根据下式计算所述目标层段的新型杂卤石钾盐矿的累积体积估算结果：

25、a＝σihki*vki*si

26、式中，a为所述目标层段的新型杂卤石钾盐矿的累积体积估算结果，hki为所述有效钾盐矿层的有效厚度，vki为所述有效钾盐矿层钾含量系数，si为地震网格单元面积。

27、本发明另一实施例对应提供了一种新型杂卤石钾盐矿数据累积体积计算系统，包括：

28、门槛值获取模块，用于获取目标层段中有效钾盐矿层的钾参数和密度对应的有效门槛值；

29、有效参数过滤模块，用于通过地震反演计算得到所述目标层段的所述钾参数和密度的值，根据所述有效门槛值对所述钾参数和密度进行过滤，得到有效钾参数的时间域空间分布；其中，所述有效钾参数即所述目标层段中有效钾盐矿层的钾参数，所述有效钾盐矿层即符合品位要求的新型杂卤石钾盐矿层；

30、有效厚度计算模块，用于对所述有效钾参数的时间域的有效钾参数进行时间-深度转换后，得到深度域有效钾参数，通过累加计算得到所述有效钾盐矿层的累积厚度，对所述有效钾盐矿层的累积厚度进行校正，得到所述有效钾盐矿层的有效厚度；

31、钾含量系数平面分布值计算模块，用于根据所述有效钾参数计算所述平均有效钾参数平面分布值，并根据数学公式转换得到所述有效钾盐矿层的钾含量系数平面分布值；

32、累积体积估算模块，用于根据所述有效钾盐矿层有效厚度和所述有效钾盐矿层的钾含量系数平面分布值得到所述目标层段新型杂卤石钾盐矿的累积体积估算结果。

33、作为上述方案的改进，所述根据所述有效钾盐矿层有效厚度和所述有效钾盐矿层钾含量系数平面分布值得到所述目标层段中新型杂卤石钾盐矿的累积体积估算结果，包括：

34、根据下式计算所述目标层段的新型杂卤石钾盐矿的累积体积估算结果：

35、a＝σihki*vki*si

36、式中，a为所述目标层段新型杂卤石钾盐矿的累积体积估算结果，hki为所述有效钾盐矿层的有效厚度，vki为所述有效钾盐矿层的钾含量系数平面分布值，si为地震网格单元面积。

37、本发明另一实施例对应提供了一种新型杂卤石钾盐矿数据累积体积计算装置，包括：处理器、存储器以及存储在所述存储器中且被配置为由所述处理器执行的计算机程序，所述处理器执行所述计算机程序时实现如本发明实施例一所述的新型杂卤石钾盐矿数据累积体积计算方法。

38、本发明另一实施例对应提供了一种计算机可读存储介质，其特征在于，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行如本发明实施例一所述的新型杂卤石钾盐矿数据累积体积计算方法。

39、相比于现有技术，在本发明实施例中，通过有效门槛值对有效参数进行筛选，提高了起始数据的准确性；进一步通过地震反演计算获得新型杂卤石钾盐矿的参数数据，并采用采用体积精确累积方法计算不稳定分布新型杂卤石钾盐矿体的累积体积，提高了新型杂卤石钾盐矿不稳定分布时体积计算的准确性和可靠性，使得计算结果更贴近真实数据。