上篇文章，我们已经生成过一个采空区网格。
储量计算时，我们需要用石膏矿层的网格减去采空区网格，然后再计算储量。

一、合并B1，B2，B3

将三个石膏矿层合并到一起，可以先复制一个B1，然后把B2，B3都merge进来。

B1复制为B
将B2，B3合并进B：
B2合并进入B
B3同样操作，合并进入B。

二、B减去采空区

从石膏矿层中减去采空区

然后B改名为"B-CaiKongqu"，方便自己查看。

三、参数准备

GOCAD中储量计算公式是以体积来表达的：
V_pore= V_net × porosity

• V_net 指有效网格，本实例中指石膏品位大于55%的网格
• porosity ：孔隙度，这个值在0到1之间。
  孔隙度这里用品位×密度来代替 ，但要保证这两个值相乘之后要位于0到1之间。我们把品位数据乘以0.01，这个是正常的；同时也把密度数据乘以0.1，这个是作了处理的。所以最终储量算出来以后要乘以10.

定义孔隙度

这个所谓的孔隙度=品位×0.01×密度×0.1

创建孔隙度

然后计算得到我们设置的这个孔隙度值：

孔隙度值计算

定义NTG

NTG用来判断网格的有效性或无效性。

创建NTG

如果品位>55,值为1，否则值为0；

NTG判别函数

操作如下：

计算NTG
计算完成后，白色为NTG为0的区域，这些地方品位不足55%，将不参与储量计算。
NTG分布图

四、储量计算

储量计算
计算结果
前面说过，对密度乘以0.1后计算结果要乘以10。
所以石膏矿的储量就直接为9.9 × 10⁷吨，下面详细解释一下。

疑问：

1. 数据单位问题

• 数据中使用的的密度单位为g/cm³，而软件中默认的单位为Kg/m³；1g/cm³ = 1000 Kg/m³，二者相差1000倍。
  所以，计算的最终结果应该是9.9 × 10⁷ × 10³ kg，也就是9.9 × 10⁷吨。
• 软件中显示的体积单位，但因为我们在处理孔隙度时，把密度乘进去了，所以我们实际得到的数据为质量单位（kg）而非体积单位(m³)。

2. 结果中几个数字的意思：

• geology_98×173×460，为工区总网格数
• Net Number of cells for oil Zone，为矿体的网格数量
• Gross Volume for Oil Zone，为矿体的总体积
• Net volume for Oil Zone，为品位大于指定值（55%）的矿石总体积
• Porous Volume for Oil Zone，为乘以孔隙度之后的体积，即孔隙体积（石油、天然气）。这里因为对孔隙度作了处理，计算结果为品位大于55%的石膏的质量。