Отправлено: 03.11.08 11:02. Заголовок: Модели Красильникова

Красным цветом в левой части уравнения, насколько я понимаю, отмечена дисперсная характеристика метода каким образом значение дисперсии оценивается по результатам исследования керна (откуда эти выделенные мной красным значения)?
Насколько универсальны дисперсии предложенные в примерах моделей, определение ФЕС в смешанных коллекторах(по аналогии с красным цветом)?




parameters <Константы>

//пластовый флюид
1 530 //время пробега в пластовом флюиде
2 1.1 //плотность пластового флюида
3 1 //водородосодержание пластового флюида

//известняк
4 155 //время пробега в известняке
5 2.71 //плотность известняка
6 0 //водородосодержание известняка

//доломит
7 138 //время пробега в доломите
8 2.87 //плотность доломита
9 0.03 //водородосодержание доломита

//ангидрит
10 164 //время пробега в ангидрите
11 2.98 //плотность ангидрита
12 -0.007//водородосодержание ангидрита

//кварц
13 168 //время пробега в кварце
14 2.64 //плотность кварца
15 -0.04 //водородосодержание кварца

//газ
16 700 //время пробега в газе
17 0.14 //плотность газа
18 0.09 //водородосодержание газа

//битум
19 350 //время пробега в битуме
20 1.3 //плотность битума
21 0.5 //водородосодержание битума

//глина
22 330 //время пробега в глине
23 2.71 //плотность глины
24 0.3 //водородосодержание глины

//кавернозность
25 156 //время пробега
26 1.16 //плотность
27 1 //водородосодержание

//нефть
28 620 //время пробега в нефти
29 0.8 //плотность нефти
30 0.8 //водородосодержание нефти

//аргиллит
31 240 //время пробега в аргиллите
32 2.71 //плотность аргиллита
33 -0.019 //водородосодержание аргиллита
end

//***************************************************************************************************
Program <Чистка результатов>
ColumnCreate(Кп,'0')
ColumnCreate(Кг,0)
ColumnCreate(Кп_кав,0)
ColumnCreate(Изв,0)
ColumnCreate(Дол,0)
ColumnCreate(Анг,0)
ColumnCreate(Арг,0)
ColumnCreate(Кварц,0)
ColumnCreate(Бит,0)
ColumnCreate(lit,'--')
SetDescription(Кп, 'POR', 0, 'Depth')
SetDescription(Кг, 'SGAS', 0, 'Depth')
SetDescription(Кп_кав, 'CPOR', 0, 'Depth')
SetDescription(Изв, 'COMP', 0, 'Depth')
SetDescription(Дол, 'COMP', 0, 'Depth')
SetDescription(Анг, 'COMP', 0, 'Depth')
SetDescription(Арг, 'COMP', 0, 'Depth')
SetDescription(Кварц, 'COMP', 0, 'Depth')
SetDescription(lit, 'LIT', 0, 'Depth')
End

//***************************************************************************************************
model <Кп+Изв+Дол>
equations
1/9.596 : АК = Кп*par[1] + Изв*par[4] + Дол*par[7]
1/0.062 : ГГК = Кп*par[2] + Изв*par[5] + Дол*par[8]
1/0.036 : W = Кп*par[3] + Изв*par[6] + Дол*par[9]
300 : Кп + Изв + Дол= 1
in
АК, ГГК, W
out
Кп : [0.0, 0.3]
Изв : [0.0, 1.0]
Дол : [0.0, 1.0]
end
program <Кп+Изв+Дол> parameters <Константы>
Optimization('Кп+Изв+Дол', 500, 1.0e-5)
АК_теор = Кп*par[1] + Изв*par[4] + Дол*par[7]
ГГК_теор = Кп*par[2] + Изв*par[5] + Дол*par[8]
W_теор = Кп*par[3] + Изв*par[6] + Дол*par[9]
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

//***************************************************************************************************
model <Кп+Изв+Дол+Анг>
equations
1/9.596 : АК = Кп*par[1] + Изв*par[4] + Дол*par[7] + Анг*par[10]
1/0.062 : ГГК = Кп*par[2] + Изв*par[5] + Дол*par[8] + Анг*par[11]
1/0.036 : W = Кп*par[3] + Изв*par[6] + Дол*par[9] + Анг*par[12]
400 : Кп + Изв + Дол + Анг = 1
in
АК, ГГК, W
out
Кп : [0.0, 0.3]
Изв : [0.0, 1.0]
Дол : [0.0, 1.0]
Анг : [0.0, 1.0]
end
program <Кп+Изв+Дол+Анг> parameters <Константы>
Optimization('Кп+Изв+Дол+Анг', 500, 1.0e-5)
АК_теор = Кп*par[1] + Изв*par[4] + Дол*par[7] + Анг*par[10]
ГГК_теор = Кп*par[2] + Изв*par[5] + Дол*par[8] + Анг*par[11]
W_теор = Кп*par[3] + Изв*par[6] + Дол*par[9] + Анг*par[12]
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

//***************************************************************************************************
model <Кп+Изв+Дол+Бит>
equations
1/15 : АК = Кп*par[1] + Изв*par[4] + Дол*par[7] + Бит*par[19]
1/0.1 : ГГК = Кп*par[2] + Изв*par[5] + Дол*par[8] + Бит*par[20]
1/0.05 : W = Кп*par[3] + Изв*par[6] + Дол*par[9] + Бит*par[21]
400 : Кп + Изв + Дол + Бит = 1
in
АК, ГГК, W
out
Кп : [0.0, 0.3]
Изв : [0.0, 1.0]
Дол : [0.0, 1.0]
Бит : [0.0, 1.0]
end
program <Кп+Изв+Дол+Бит> parameters <Константы>
Optimization('Кп+Изв+Дол+Бит', 500, 1.0e-5)
АК_теор = Кп*par[1] + Изв*par[4] + Дол*par[7] + Бит*par[19]
ГГК_теор = Кп*par[2] + Изв*par[5] + Дол*par[8] + Бит*par[20]
W_теор = Кп*par[3] + Изв*par[6] + Дол*par[9] + Бит*par[21]
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

//***************************************************************************************************
model <Кп+Дол+Арг+Кварц>
equations
1/15 : АК = Кп*par[1] + Дол*par[7] + Арг*par[31] + Кварц*par[13]
1/0.1 : ГГК = Кп*par[2] + Дол*par[8] + Арг*par[32] + Кварц*par[14]
1/0.05 : W = Кп*par[3] + Дол*par[9] + Арг*par[33] + Кварц*par[15]
400 : Кп + Дол + Арг + Кварц= 1
in
АК, ГГК, W
out
Кп : [0.0, 0.3]
Дол : [0.0, 1.0]
Арг : [0.0, 1.0]
Кварц : [0.0, 1.0]
end
program <Кп+Дол+Арг+Кварц> parameters <Константы>
Optimization('Кп+Дол+Арг+Кварц', 500, 1.0e-5)
АК_теор = Кп*par[1] + Дол*par[7] + Арг*par[31] + Кварц*par[13]
ГГК_теор = Кп*par[2] + Дол*par[8] + Арг*par[32] + Кварц*par[14]
W_теор = Кп*par[3] + Дол*par[9] + Арг*par[33] + Кварц*par[15]
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

//***************************************************************************************************
model <Кп+Изв+Дол+Анг+Кварц>
equations
1/9.596 : АК = Кп*par[1] + Изв*par[4] + Дол*par[7] + Анг*par[10] + Кварц*par[13]
1/0.062 : ГГК = Кп*par[2] + Изв*par[5] + Дол*par[8] + Анг*par[11] + Кварц*par[14]
1/0.036 : W = Кп*par[3] + Изв*par[6] + Дол*par[9] + Анг*par[12] + Кварц*par[15]
500 : Кп + Изв + Дол + Анг + Кварц = 1
in
АК, ГГК, W
out
Кп : [0.0, 0.3]
Изв : [0.0, 1.0]
Дол : [0.0, 1.0]
Анг : [0.0, 1.0]
Кварц : [0.0, 1.0]
end
program <Кп+Изв+Дол+Анг+Кварц> parameters <Константы>
Optimization('Кп+Изв+Дол+Анг+Кварц', 500, 1.0e-5)
АК_теор = Кп*par[1] + Изв*par[4] + Дол*par[7] + Анг*par[10] + Кварц*par[13]
ГГК_теор = Кп*par[2] + Изв*par[5] + Дол*par[8] + Анг*par[11] + Кварц*par[14]
W_теор = Кп*par[3] + Изв*par[6] + Дол*par[9] + Анг*par[12] + Кварц*par[15]
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

//***************************************************************************************************
model <Кп+Кп_кав+Изв+Дол>
equations
1/15 : АК = Кп*par[1] + Кп_кав*par[25] + Изв*par[4] + Дол*par[7]
1/0.1 : ГГК = (Кп + Кп_кав)*par[2] + Изв*par[5] + Дол*par[8]
1/0.05 : W = (Кп + Кп_кав)*par[3] + Изв*par[6] + Дол*par[9]
400 : Кп + Кп_кав + Изв + Дол = 1
in
АК, ГГК, W
out
Кп : [0.0, 0.3]
Кп_кав : [0.0, 0.3]
Изв : [0.0, 1.0]
Дол : [0.0, 1.0]
end
program <Кп+Кп_кав+Изв+Дол> parameters <Константы>
Optimization('Кп+Кп_кав+Изв+Дол', 500, 1.0e-5)
АК_теор = Кп*par[1] + Кп_кав*par[25] + Изв*par[4] + Дол*par[7]
ГГК_теор = (Кп + Кп_кав)*par[2] + Изв*par[5] + Дол*par[8]
W_теор = (Кп + Кп_кав)*par[3] + Изв*par[6] + Дол*par[9]
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

//***************************************************************************************************
model <Кп+Кп_кав+Изв+Дол+Анг>
equations
1/15 : АК = Кп*par[1] + Кп_кав*par[25] + Изв*par[4] + Дол*par[7] + Анг*par[10]
1/0.1 : ГГК = (Кп + Кп_кав)*par[2] + Изв*par[5] + Дол*par[8] + Анг*par[11]
1/0.05 : W = (Кп + Кп_кав)*par[3] + Изв*par[6] + Дол*par[9] + Анг*par[12]
500 : Кп + Кп_кав + Изв + Дол +Анг=1
in
АК, ГГК, W
out
Кп : [0.0, 0.3]
Кп_кав : [0.0, 0.3]
Изв : [0.0, 1.0]
Дол : [0.0, 1.0]
Анг : [0.0, 1.0]
end
program <Кп+Кп_кав+Изв+Дол+Анг> parameters <Константы>
Optimization('Кп+Кп_кав+Изв+Дол+Анг', 500, 1.0e-5)
АК_теор = Кп*par[1] + Кп_кав*par[25] + Изв*par[4] + Дол*par[7] + Анг*par[10]
ГГК_теор = (Кп + Кп_кав)*par[2] + Изв*par[5] + Дол*par[8] + Анг*par[11]
W_теор = (Кп + Кп_кав)*par[3] + Изв*par[6] + Дол*par[9] + Анг*par[12]
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

//***************************************************************************************************
model <Кп+Кг+Изв+Дол+Анг>
equations
1/15 : АК = Кп*par[1] + Кг*par[16] + Изв*par[4] + Дол*par[7] + Анг*par[10]
1/0.1 : ГГК = Кп*par[2] + Кг*par[17] + Изв*par[5] + Дол*par[8] + Анг*par[11]
1/0.05 : W = Кп*par[3] + Кг*par[18] + Изв*par[6] + Дол*par[9] + Анг*par[12]
500 : Кп + Кг + Изв + Дол + Анг = 1
in
АК, ГГК, W
out
Кп : [0.0, 0.3]
Кг : [0.0, 1.0]
Изв : [0.0, 1.0]
Дол : [0.0, 1.0]
Анг : [0.0, 1.0]
end
program <Кп+Кг+Изв+Дол+Анг> parameters <Константы>
Optimization('Кп+Кг+Изв+Дол+Анг', 500, 1.0e-5)
АК_теор = Кп*par[1] + Кг*par[16] + Изв*par[4] + Дол*par[7] + Анг*par[10]
ГГК_теор = Кп*par[2] + Кг*par[17] + Изв*par[5] + Дол*par[8] + Анг*par[11]
W_теор = Кп*par[3] + Кг*par[18] + Изв*par[6] + Дол*par[9] + Анг*par[12]
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

//***************************************************************************************************
model <Оценка Кг>
equations
1/15 : АК = Кп*par[1] + Кп_кав*par[25] + Кг*par[16] + Изв*par[4] + Дол*par[7] + Анг*par[10]
1/0.1 : ГГК = (Кп + Кп_кав)*par[2] + Кг*par[17] + Изв*par[5] + Дол*par[8] + Анг*par[11]
1/0.05 : W = (Кп + Кп_кав)*par[3] + Кг*par[18] + Изв*par[6] + Дол*par[9] + Анг*par[12]
300 : Кп + Кп_кав + Кг + Изв + Дол + Анг = 1
in
АК, ГГК, W, Изв, Дол, Анг
out
Кп : [0.0, 0.3]
Кп_кав : [0.0, 0.3]
Кг : [0.0, 0.3]
end
program <Кп+Кп_кав+Кг+Изв+Дол+Анг> parameters <Константы>
Optimization('Кп+Кп_кав+Изв+Дол+Анг', 500, 1.0e-5)
Optimization('Оценка Кг', 500, 1.0e-5)
АК_теор = Кп*par[1] + Кп_кав*par[25] + Кг*par[16] + Изв*par[4] + Дол*par[7] + Анг*par[10]
ГГК_теор = (Кп + Кп_кав)*par[2] + Кг*par[17] + Изв*par[5] + Дол*par[8] + Анг*par[11]
W_теор = (Кп + Кп_кав)*par[3] + Кг*par[18] + Изв*par[6] + Дол*par[9] + Анг*par[12]
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

//***************************************************************************************************
Program <Уравнение баланса>
SUM = Кп + Кп_кав + Кг + Изв + Дол + Анг + Арг + Кварц + Бит
End

//***************************************************************************************************
program <Интерпретация по комплексу моделей> parameters <Константы>
if Depth <= 1773 then
Optimization('Кп+Изв+Дол+Анг', 500, 1.0e-5)
АК_теор = Кп*par[1] + Изв*par[4] + Дол*par[7] + Анг*par[10]
ГГК_теор = Кп*par[2] + Изв*par[5] + Дол*par[8] + Анг*par[11]
W_теор = Кп*par[3] + Изв*par[6] + Дол*par[9] + Анг*par[12]
Кп = 0
endif
if Depth > 1773 then
Optimization('Кп+Кп_кав+Изв+Дол+Анг', 500, 1.0e-5)
Optimization('Оценка Кг', 500, 1.0e-5)
АК_теор = Кп*par[1] + Кп_кав*par[25] + Кг*par[16] + Изв*par[4] + Дол*par[7] + Анг*par[10]
ГГК_теор = (Кп + Кп_кав)*par[2] + Кг*par[17] + Изв*par[5] + Дол*par[8] + Анг*par[11]
W_теор = (Кп + Кп_кав)*par[3] + Кг*par[18] + Изв*par[6] + Дол*par[9] + Анг*par[12]
endif
SetDescription(АК_теор, 'DTP')
SetDescription(ГГК_теор, 'RHOB')
SetDescription(W_теор, 'NGP')
end

.....С уважением.........