4. Алгоритм интерпретации квд модифицированным методом детерминированных моментов

Алгоритм получения оценок (1) (сценарий работы комплекса программ) состоит из следующей последовательности действий:

1. Решается задача идентификации КВД, а именно, определяются параметры модели КВД (4) выбранным методом оптимизации и управляющий параметр(5).

Оптимизационная задача идентификации (4), при использовании метода оптимизации Гаусса-Ньютона, сводится к последовательному решению систем линейных уравнений (СЛУ) вида :

(12)

где - матрица частных производных (в программе производные заменяются их конечными разностями);- вектор невязок; _{(доказательство приведено в приложении В).}

3. На основе приближений пластового давления определяются оценки детерминированных моментов, параметры нефтяного пласта (гидропроводность, пъезопроводность) и диагностического коэффициент (1), (2).

5. Результаты решения тестовых заданий

Значения забойного давления, динамического уровня и дебиты со скважин представлены в Таблица 1. (Приложение А.), исходные данные и дополнительные априорные данные в Таблица 2. (Приложение А.).

1. Тестовый пример №1

Таблица 1. Результаты интерпретации КВД скважины № 121

Результаты интерпретации КВД скважины №121
Перечень результатов интерпретации скважины	Значение
Гидропроводность (см/сП)	1.06
Пъезопроводность (см2/сек)	0.002448273
Пластовое давление (ат)	146,5
Детерменированный момент нулевого порядка ()	89815389
Детерменированный момент первого порядка ()	2,34785*1013
Детерменированный момент второго порядка ()	8.00496*1018
Диагностический коэффициент (d*)	1,304271528
Тип нефтяного пласта (залежи)	Трещиновато-пористый

При выполнении расчетов учитывается экспертная оценка пластового давления.

На рисунке 1: модель 1 – фактичекое значение забойного давления, модель 2 – значения забойного давления, расчитанные по методу деформированного многогранника, при значениях параметров модели

МН -

146,7096

0,298216

0,08629

модель 3 – значение забойного давления в режиме ручного подбора параметров при

ALo(1)	ALo(2)	ALo(3)
155	0,45	0,1

Рис. 1 – Графики моделей КВД, по методу деформированного многогранника.

Рис. 2 – Ошибка оценок пластового давления по методу деформированного многогранника.

Рис. 3 – Точность прогноза забойного давления по методу деформированного многогранника.

На рисунке 4: модель 1 построена при параметрах модели, заданных вручную

ALo(1)	ALo(2)	ALo(3)
180,00	0,18	0,69

Модель 2 построена по методу Гаусса-Ньютона при параметрах модели

ALo(1)	ALo(2)	ALo(3)
179,63	0,18	0,70

Рис. 4 – График моделей КВД по методу Гаусса-Ньютона.

Вид матрицы :

0,285229	51,23572	216,482
0,435029	156,2885	179,8839
0,544264	293,2987	163,0451
0,628562	451,6347	153,0631
0,69534	624,5193	146,4309
0,749038	807,2979	141,7291
0,792636	996,6686	138,2545
0,82827	1190,257	135,6127
0,857533	1386,349	133,5626
0,881651	1583,709	131,9472
0,901581	1781,46	130,6595
0,918085	1978,987	129,6237
0,931775	2175,871	128,7847
0,943146	2371,842	128,1012
0,9526	2566,733	127,5419
0,960467	2760,46	127,0825
0,967019	2952,996	126,704
0,972478	3144,353	126,3915
0,977029	3334,57	126,1328
0,980824	3523,708	125,9185
0,98399	3711,835	125,7406

Вид матрицы :

0,285229	0,435029	0,544264	0,628562	0,69534	0,749038	0,792636	0,82827	0,857533	0,881651
51,23572	156,2885	293,2987	451,6347	624,5193	807,2979	996,6686	1190,257	1386,349	1583,709
216,482	179,8839	163,0451	153,0631	146,4309	141,7291	138,2545	135,6127	133,5626	131,9472

Продолжение:

0,901581	0,918085	0,931775	0,943146	0,9526	0,960467	0,967019	0,972478	0,977029	0,980824	0,98399
1781,46	1978,987	2175,871	2371,842	2566,733	2760,46	2952,996	3144,353	3334,57	3523,708	3711,835
130,6595	129,6237	128,7847	128,1012	127,5419	127,0825	126,704	126,3915	126,1328	125,9185	125,7406

Вид матрицы :

46,11

224987,58

6567,88

Вид матрицы :

46,11

224987,58

6567,77