Методы повышения эффективности процессов добычи и транспорта газа, страница 22

Исходной информацией при числовом моделировании слу­жили результаты решения системы дифференциальных уравне­ний, описывающих нестационарное изотермическое течение га­за. К полученным в результате числового решения величинам давления добавляли центрированную адаптивную помеху с дис­персией, соответствующей погрешности измерительного прибо­ра класса 1. Аналогично имитировались погрешности измере­ния расхода, которые соответствовали погрешностям измери­тельного прибора класса 2.

В табл. II.6 приведены массивы исходной информации и ре­зультаты моделирования.

Массив р — двумерный массив относительных величин р, 5 = 37 (37 временных слоев); массив Qo — относительных вели­чин расхода; D—диаметр газопровода, м; Ах — протяженность элементарного__участка, х=15км; |0 — начальное значение ве­личины £— 1/УА, для элементарного участка (выбирается произ­вольно); т, п — постоянные величины; B — sign{...}; Аг — на­чальное значение параметра А (выбирается произвольно); А/— интервал дискретизации по времени; Д£=1800 с; с — скорость

47


Таб.

пи

ца Ц.6

Ряехоп

Давление по элементарным участкам

Время

Q.10-7,

тыс. м*/сут

1

2

3

4

0

222

1,203283

,210403

1,206970

1,204463

1

222

1,168380

1

[,158418

1,165866

1,166906

2

222

1,125972

1

,117540

1,141580

1,148640

3

222

1,094797

1

[,094040

1,078330

1,098059

4

223

1,068330

1

,057070

1,056600

1,063440

5

224

1,007820

1

,013815

1,014646

t ,026064

6

226

0,9845051

0,965390

0,981010

0,969186

7

228

0,9387153

0,938252

0,915440

0,900970

8

233

0,896400

0,862295

0,841510

0,824728

9

233

0,808133

0,778330

0,737000

0,710309

Продолжение табл.  II.6

Время

Расход Q-10—?, тыс. м*/сут

Давление

; По элементарным

участкам

5

6

7

0

222

1,208233

,196553

1,198063

1

222

1,164050

,162250

1,175530

2

222

1,130810

,132637

1,119620

3

222

1,099190

,091330

1,090800

4

223

1,056715

,043909

1,075590

5

224

1,017264

,996140

0,999590

6

226

0,947717

0,955930

0,939980

7

228

0,893700

0,884560

0,865727

8

233

0,810580

0,798950

0,765944

9

233

0,687550

0,635520

0,622060

звука в газе, с = 375 м/с; хо = 0,000084 с/м; ро = 0,0000307 с/м. Расчет 1/А приведен в табл. II.9, 11.10, 11.11. Эти таблицы от­личаются выбором коэффициента ys.

Размерности и диапазоны изменения констант, входных и выходных величин алгоритма приведены ниже.


Параметр

Я М

At, с N D

Ц

Диапазон измере­ния

0—200

0—200

0—3600

0—5

0,5-1,4

0—200

ps, ps-!, кгс/см2 20 • 10«—76 • 104
Qs_i, кг/с                50—150

Т, К                      220—350


Параметр

Ах, м

тп, п

/V кгс/см2


Диапазон измере­ния

10 000—50 000

3—5 0,1—0,3