/// <summary>
        /// Функция моделирования работы обводненной скважины с погружным насосом, возвращает дебит газа (тыс.м3/сут)
        /// </summary>
        /// <param name="eps">Точность с которой будет расчитываться забойное давление (безразмерная)</param>
        /// <param name="dynamicH">Динамический уровень пластовой жидкости, отсчитывается от входа в насос (м)</param>
        /// <param name="bottomHolePressure">Забойное давление (МПа)</param>
        /// <param name="waterRate">Расход жидкости (м3/сут)</param>
        /// <param name="nglRate">Дебит газового конденсата (т/сут)</param>
        /// <returns>Дебит газа в стандартных условиях (тыс.куб.м/сут)</returns>
        public double ModelingPump(double eps, double dynamicH, out double bottomHolePressure, out double waterRate, out double nglRate)
        {
            //Скважина делится на 3 части:
            //topPartWell - часть затрубного пространства, выше динамического уровня
            //bottomPartWell - часть затрубного пространства, ниже динамического уровня
            //innerPartWell - НКТ
            double topPartLength       = Tubing.Length - dynamicH;
            double dynamicHTemperature = BottomholeTemperature + (WellheadTemperature - BottomholeTemperature) * (dynamicH / Tubing.Length);


            AnnularTubing = new Tubing(Tubing);
            AnnularTubing.TubingDiameter = Tubing.CalcEquivalentAnnularDiameter();

            //Для затрубного пространства свободного от жидкости создаем эквивалент сухой скважины
            DryWell topPartWell = new DryWell(this);

            topPartWell.Tubing             = AnnularTubing;
            topPartWell.BottomholePressure = dynamicHTemperature;


            double Ph;             //давление в точке динамического уровня
            double gasQ = topPartWell.Modeling(eps, out Ph, out nglRate);

            bottomHolePressure = Ph + WaterFluid.CalcColumnPressure(dynamicH);
            //Давление для притока воды берется как давление на середине динамического уровня
            double waterPress = Ph + WaterFluid.CalcColumnPressure(dynamicH / 2);

            waterRate = Layer.CalcWaterRate(bottomHolePressure);

            return(gasQ);
        }
 /// <summary>
 /// Конструктор копирования DryWell
 /// </summary>
 /// <param name="dryWell"></param>
 public DryWell(DryWell dryWell)
 {
     Layer                 = dryWell.Layer;
     GasFluid              = dryWell.GasFluid;
     Tubing                = new Tubing(dryWell.Tubing);
     WellheadPressure      = dryWell.WellheadPressure;
     WellheadTemperature   = dryWell.WellheadTemperature;
     BottomholeTemperature = dryWell.BottomholeTemperature;
     NglFluid              = dryWell.NglFluid;
 }
        /// <summary>
        /// Функция высчитывания забойного давления остановленной скважины по затрубному давлению (МПа)
        /// </summary>
        public double CalcStaticBottomholePressureByAnnularTubing(double annularWellheadPressure)
        {
            AnnularTubing = new Tubing(Tubing);
            AnnularTubing.TubingDiameter = Tubing.CalcEquivalentAnnularDiameter();

            DryWell well = new DryWell(this);

            well.Tubing           = AnnularTubing;
            well.WellheadPressure = annularWellheadPressure;

            GasFlow gasFlow = new GasFlow(GasFluid, 0, annularWellheadPressure, WellheadTemperature);
            double  Pb      = well.Tubing.CalcStaticBottomholePressure(gasFlow, BottomholeTemperature);

            return(Pb);
        }
        /// <summary>
        /// Функция высчитывания затрубного давления на устье остановленной скважины по забойному давлению (МПа)
        /// </summary>
        public double CalcStaticAnnularTubingByBottomholePressure(double bottomholePressure)
        {
            AnnularTubing = new Tubing(Tubing);
            AnnularTubing.TubingDiameter = Tubing.CalcEquivalentAnnularDiameter();

            DryWell well = new DryWell(this);

            well.Tubing             = AnnularTubing;
            well.WellheadPressure   = 0;
            well.BottomholePressure = bottomholePressure;

            GasFlow gasFlow = new GasFlow(GasFluid, 0, bottomholePressure, BottomholeTemperature, true);

            double annularWellheadPressure = well.Tubing.CalcStaticWellheadPressure(gasFlow, WellheadTemperature);

            return(annularWellheadPressure);
        }
        /// <summary>
        /// Функция получения точек кривых притока P2 и оттока P1 по затрубному пространству для узлового анализа
        /// обводненной скважины (с откачкой насосом)
        /// </summary>
        /// <param name="dynamicH">Динамический уровень пластовой жидкости, отсчитывается от входа в насос (м)</param>
        /// <param name="num">Число точек</param>
        /// <param name="P1">Массив значений забойного давления для точек кривой оттока с забоя (МПа)</param>
        /// <param name="P2">Массив значений забойного давления для точек кривой притока из пласта к забою (МПа)</param>
        /// <param name="gasDischarge">Массив значений дебита газа (тыс.куб.м/сут)</param>
        /// <param name="waterDischarge">Массив значений дебит воды (м3/сут)</param>
        public void GetPointsAnnularTubingNodeAnalize(double dynamicH, int num, out double[] P1, out double[] P2, out double[] gasDischarge, out double[] waterDischarge)
        {
            double topPartLength       = Tubing.Length - dynamicH;
            double dynamicHTemperature = BottomholeTemperature + (WellheadTemperature - BottomholeTemperature) * (dynamicH / Tubing.Length);


            P1             = new double[num];
            P2             = new double[num];
            gasDischarge   = new double[num];
            waterDischarge = new double[num];

            double maxQ = Layer.CalcMaxGasRate();
            double step = maxQ / num;
            double gasQ = 0;

            //создается новый Tubing чтобы не перезаписывать SegmentFlows
            Tubing annularTubing = new Tubing(AnnularTubing);

            //Для затрубного пространства свободного от жидкости создаем эквивалент сухой скважины
            DryWell topPartWell = new DryWell(this);

            topPartWell.Tubing             = annularTubing;
            topPartWell.BottomholePressure = dynamicHTemperature;

            for (int i = 0; i < num; i++)
            {
                gasDischarge[i] = gasQ;

                GasFlow gasFlow = new GasFlow(GasFluid, gasQ, WellheadPressure, WellheadTemperature);

                P2[i] = topPartWell.Layer.CalcBottomholePressure(gasQ);
                P1[i] = topPartWell.Tubing.CalcBottomholePressure(gasFlow, dynamicHTemperature);
                //Давление для притока воды берется как давление на середине динамического уровня
                double waterPres = P2[i] + WaterFluid.CalcColumnPressure(dynamicH / 2);
                waterDischarge[i] = Layer.CalcWaterRate(waterPres);

                gasQ += step;
            }
        }