Приведенная формула предназначена для случая, когда среднее пластовое давление Рпл выше Рнас, а Рзаб, при котором определен дебит жидкости Q, ниже давления насыщения. Индикаторную кривую в этом случае называют композитной. На участке от Рзаб=Рпл до Рзаб=Рнас индикаторная линия представляет собой прямую, а на участке от Рзаб=Рнас до Рзаб = 0 индикаторная линия загибается к оси Х. Физический смысл этого следующий: при Рзаб >= Рнас Кпрод постоянен, при Рзаб <= Рнас Кпрод снижается пропорционально снижению Рзаб.

Когда Рзаб >= Рнас для расчетов используется формула Кпрод=Q/(Рпл-Рнас).

Когда Рзаб <= Рнас используется уравнение Вогеля, одна из модификаций которого приведена выше.

Литературные ссылки прилагаю. Это отрывок из книги H. Dale Beggs "Production optimization using Nodal analysis".IPR_models.pdf

Для случая, когда Рпл <= Рнас приведенную тобой формулу композитной индикаторной линии использовать не нужно. Для этого случая есть другая формула:

Qmax=Qж/[1-0.2*Pзаб/Рпл-0.8*(Рзаб/Рпл)^2]

где Qmax = AOFP/1.8 = Кпрод*Рпл/1.8

Qmax - максимально-возможный дебит для двухфазного потока (нефть+газ) при Рзаб = 0 атм
AOFP - Absolute Open Flow Potential = Кпр*Рпл - абсолютно-свободный дебит нефти при Рзаб=0 атм
Рзаб - забойное давление, при котором замерили дебит жидкости Qж
Рпл - среднее пластовое давление (оно должно быть ниже Рнас)
Кпрод - коэффициент продуктивности

Пример использования этой формулы приведен в приложенном выше файле IPR models