Да, для каждой фазы отдельные карты.

Можно сделать как вы описали - внутри внутренного контура взять карту эффективных толщин и проинтерполировать к нулю на внешнем контуре газо/нефтеносности.

Есть еще такой варинт:

1. Из кровли вычитаете ваш контакт, получаете объем под кровлей (V1)

2. Из подошвы также вычитаете контакт, получаете объем под подошвой (V2)

3. Нужный объем залежи находите как V1-V2 (отрицательный объем удаляете)

4. Далее это объем умножаете на коэфф. песчанистости и получаете карту эфф.нн. толщин

бонус: коэфф. песчанистости можно посчитать путем деления карты эффективных на общих толщин.

исходя из задачи, как я понял, речь о 2Д моделировании. в этом случае алгоритм следующий:

1. строится карта кровли коллектора (не пласта, методом схождения от карты кровли моделируемого пласта или ближайшего пласта, если имеется структурка по сейсмике, либо от карты подошвы коллектора вышележащего пласта)

2. строится карта подошвы коллектора (не пласта, методом схождения от карты кровли коллектора)

    на основании пп.1-2 определяются внешний контур нефте- и газоносности по карте кровли коллектора, внутренний контур нефте- и газоносности по карте подошвы коллектора. НЕ ПЛАСТА, А КОЛЛЕКТОРА.

3. строится карта песчанистости. здесь или интерполяцией по скважинным данным, либо с использованием атрибутного тренда, если таковой имеется.

  первые 3 карты - исходные данные для построения 2Д модели. 

4. кровля - подошва = карта общих толщин

5. Нобщ*Кпесч = Нэфф

6. с учетом контуров продуктивности (внешний, внутренний ВНК, ГВК, ГНК) строится карта продуктивных толщин. при этом внутри внутреннего контура Нэфф = Нпрод. 

между внешним и внутренним контуром Нэфф доумножается на распределение от 1 (на внутреннем контуре) до 0 на внешнем контуре, с учетом имеющихся скважинных данных в ВНЗ, ГВЗ, ГНЗ.

если залежь не имеет внутреннего контура нефте- и газоносности (массивная залежь), то распределение зависит от скважинных данных.

ПО, разработанное для 2Д моделирования (изолайн, геоглоб, NGT) имеют "кнопку" построить карту продуктивных толщин с учетом выполненных пп. 1-5. ПО для 3Д моделирования аналогичных "кнопок" не имеют, в них придется самостоятельно построить "карту-костыль" по алгоритму коэффициентов, описанных в п.6.

при таком подходе будут иметься неточности в зонах роста или падении структуры в зонах раздела (ВНЗ, ГВЗ, ГНЗ), но обычно ими можно пренебречь.

при считывании карты продуктивных толщин с 3Д модели данных неточностей быть не должно...

p.s. ну и надо следить, чтобы Нэфф не превышала Нобщ, Нпрод не превышала Нэфф... при 2Д интерполяции бывают такие ошибки.