上篇文章，我们已经生成过一个采空区网格。
储量计算时，我们需要用石膏矿层的网格减去采空区网格，然后再计算储量。
一、合并B1，B2，B3
将三个石膏矿层合并到一起，可以先复制一个B1，然后把B2，B3都merge进来。
B1复制为B
将B2，B3合并进B：
B2合并进入B
B3同样操作，合并进入B。
二、B减去采空区
从石膏矿层中减去采空区
然后B改名为"B-CaiKongqu"，方便自己查看。
三、参数准备
GOCAD中储量计算公式是以体积来表达的：
Vpore= Vnet × porosity
- Vnet 指有效网格，本实例中指石膏品位大于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 × 107吨，下面详细解释一下。
疑问：
- 数据单位问题
- 数据中使用的的密度单位为g/cm3，而软件中默认的单位为Kg/m3；1g/cm3 = 1000 Kg/m3，二者相差1000倍。
所以，计算的最终结果应该是9.9 × 107 × 103 kg，也就是9.9 × 107吨。 - 软件中显示的体积单位，但因为我们在处理孔隙度时，把密度乘进去了，所以我们实际得到的数据为质量单位（kg）而非体积单位(m3)。
- 结果中几个数字的意思：
- 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%的石膏的质量。
