설명
국제석유의 비상계획슬롯 나라(CPM) 생산, 운송, 정제가 공동으로 개발되었습니다. 스탠포드 대학교, Sobotka and Company 및 미국학과 에너지. 이 슬롯 나라은 다양한 효과를 조사하기 위해 설계되었습니다. 만일의 사태(원유 공급 중단) 및 정책(예: 금지) 알래스카 수출). 다음 표현은 단일 3장에서 논의된 지역 버전.
소형 슬롯 나라 유형 :NLP lp
카테고리 : 슬롯 나라 모델 라이브러리
메인 파일 : srcpm.gms
$title 단일 지역 비상 계획 슬롯 나라(SRCPM,SEQ=52)
$onText
국제석유의 비상계획슬롯 나라(CPM)
생산, 운송, 정제가 공동으로 개발되었습니다.
스탠포드 대학교, Sobotka and Company 및 미국학과
에너지. 이 슬롯 나라은 다양한 효과를 조사하기 위해 설계되었습니다.
만일의 사태(원유 공급 중단) 및 정책(예: 금지)
알래스카 수출). 다음 표현은 단일
지역 버전은 3장에서 논의됩니다.
Manne, A S, Nelson, C R, So, K C 및 Weyant, J P, CPM: 우발 상황
국제석유시장 기획슬롯 나라, 국제에너지
프로그램 보고서. 기술. 대표, 스탠포드 대학교, 1982.
키워드: 비선형 계획법, 선형 계획법, 석유 생산,
생산계획, 에너지경제학
$offText
세트
c '모든 상품' / ls-light '저유황 경질유'
ls-medium '저황 중간 원유'
ls-heavy '저유황 중질유'
hs-light '고유황 경질유'
hs-medium '고유황 중질유'
hs-heavy '고유황 중질유'
ngl '천연가스 액체'
ls-feed '저유황 피드'
hs-feed '고유황 피드'
가솔린 '모터 가솔린'
증류액 '난방유+경유+등유'
나프타 '1차 나프타 LPG + 기타'
ls-resid '저유황 연료유 + 윤활유 및 왁스'
hs-resid '고유황 연료유+아스팔트' /
cr(c) '원유 유형' / ls-light, ls-medium, ls-heavy, hs-light
hs-중형, hs-중형, ngl /
ci(c) '중간제품' / ls-feed, hs-feed /
cf(c) '최종 제품' / 가솔린, 증류액, 나프타, ls-잔유, hs-잔유 /
p '공정' / d-lsl '저유황 경질유 증류'
d-lsm '저황 중질 원유 증류'
d-lsh '저유황 중질유 증류'
d-hsl '고유황 경질유 증류'
d-hsm '고유황 중질 원유 증류'
d-hsh '고유황 중질유 증류'
ngl '천연가스 액체 처리'
개혁 '개혁'
ls-hc '저유황 중말단 변환'
hs-hc '고유황 중말단 변환'
탈황 '탈황' /
pt '이송 공정' / t-gas '가솔린을 증류액으로 이송'
t-nap '나프타를 증류수로 옮김'
t-lsr '저유황 잔류물의 이송'
t-hsr '고유황 잔류물의 이송' /
m '생산 단위' / p-still '1차 증류'
ngl 'ngl 시설'
개혁가
he-conv '무거운 끝 변환기'
탈황 '탈황제' /
l '시설의 비용 수준' / b '기본 또는 저가',
m '한계 또는 높은 비용' /;
테이블 a(c,l,p) '입력 출력 계수'
d-lsl d-lsm d-lsh d-hsl d-hsm d-hsh ngl 개혁 ls-hc hs-hc desulf
ls-light.(b,m) -1.0
ls-medium.(b,m) -1.0
ls-heavy.(b,m) -1.0
hs-light.(b,m) -1.0
hs-medium.(b,m) -1.0
hs-heavy.(b,m) -1.0
ngl.(b,m) -1.0
ls-feed.b .35 .50 .68 -1.0 .99
ls-feed.m .34 .49 .67 -1.0 .99
hs-feed.b .38 .48 .66 -1.0 -1.0
hs-피드.m .38 .47 .65 -1.0 -1.0
가솔린.b .20 .96 .67 .61
가솔린.m .20 .96 .36 .25
증류기.b .28 .25 .20 .26 .23 .18 .07 .02
증류기.m .28 .25 .20 .26 .23 .18 .18 .10
나프타.b .35 .23 .10 .33 .27 .14 .80 -1.0 .04 .06
나프타.m .35 .23 .10 .33 .27 .14 .80 -1.0 .03 .04
ls-resid.b .23
ls-resid.m .42
hs-resid.b .30
hs-resid.m .60 ;
테이블 tr(c,pt) '전송 프로세스'
t-가스 t-nap t-lsr t-hsr
휘발유 -1.0
증류액 1.0 .93
나프타 -1.0
ls-피드 -1.0
hs-피드 -1.0
ls-잔류 1.0
hs-잔류 1.0;
매개변수 loss(p,l) '프로세스 손실';
손실(p,l) = 합계(c, a(c,l,p));
디스플레이 손실;
테이블 b(m,p) '용량 활용 매트릭스'
d-lsl d-lsm d-lsh d-hsl d-hsm d-hsh ngl 개혁 ls-hc hs-hc desulf
p-아직 1 1 1 1 1 1
영어 1
개혁자 1
그 변환 1 1
탈황 1;
테이블 k(m,l) '초기 용량(mbd)'
비엠
p-스틸 50.5 7.5
영어 3.4
리포머 7.1 .8
그 변환 7.3 2.9
탈황 2.7.3;
매개변수
hecr(l,p) '무거운 끝 변환 제한' / b.ls-hc -.83, (b,m).hs-hc 1.0 /
hecc(l) '헤비 엔드 변환 용량(mbd)' / b 3.9, m 2.5 /
trb(pt) '전송 범위(mbd)' / t-gas 1.7, t-nap 1.9, (t-lsr,t-hsr) inf /
trc(pt) '이송 비용(b당 $)' / t-gas -.3 /
mcpe(c) 'cpe + 드로우다운(mbd)에서 가져오기' / hs-light .7, distillat .5, hs-resid .1 /;
테이블 crd(cr,*) '원유 데이터'
가격 제한
* (mbd) ($/b)
ls-라이트 13.6 32.0
ls-매체 1.1 32.0
ls-무거움 1.0 32.0
hs-라이트 16.2 32.0
hs-중간 8.9 32.0
hs-heavy 4.4 32.0
NGL 3.1 32.0;
테이블 pc(l,p) '공정 비용(배럴당 미국 달러)'
d-lsl d-lsm d-lsh d-hsl d-hsm d-hsh ngl 개혁 ls-hc hs-hc desulf
b .45 .45 .45 .50 .50 .50 .41 .27 .32 .9
m .50 .50 .50 .55 .55 .55 .50 .45 .28 1.0;
테이블 ddat(cf,*) '수요 데이터'
ref-q 참조-p 참조-t elas
* (mbd) ($/b) ($/b)
휘발유 12.8 45.0 0.0 -.20
증류액 13.8 41.0 0.0 -.20
나프타 8.3 43.0 0.0 -.20
ls-잔류 4.2 35.3 0.0 -.30
hs-잔류 8.6 31.0 0.0 -.30;
매개변수
dema(cf) '수요 매개변수 a'
demb(cf) '수요 매개변수 b';
demb(cf) = 1/ddat(cf,"elas") + 1;
dema(cf) = ddat(cf,"ref-p")/demb(cf)/ddat(cf,"ref-q")**(demb(cf) - 1);
dema, demb, ddat를 표시합니다.
변수
z(p,l) '프로세스 수준(mbd)'
u(c) '원유 공급(mbd)'
w(pt) '이동 활동(mbd)'
xf(c) '최종제품 출하(mbd)'
pcost '생산 비용(일당 밀$)'
scps '소비자 및 생산자 잉여(일당 밀 $)';
양수 변수 z, u, w, xf;
방정식
mb(c) '물질수지(mbd)'
cc(m,l) '용량 제약 조건(mbd)'
ccr(l) '용량 제한 - 무거운 끝(mbd)'
pcdef '생산 비용 정의(일당 밀 $)'
dcps '잉여 정의(일일 밀 $)';
mb(c).. mcpe(c) + sum((p,l), a(c,l,p)*z(p,l)) + u(c)$cr(c) + sum(pt, tr(c,pt)*w(pt)) =g= xf(c)$cf(c);
cc(m,l).. sum(p, b(m,p)*z(p,l)) =l= k(m,l);
ccr(l).. sum(p, hecr(l,p)*z(p,l)) =l= hecc(l);
pcdef.. pcost =e= sum((p,l), pc(l,p)*z(p,l)) + sum(pt, trc(pt)*w(pt)) + sum(cr, crd(cr,"price")*u(cr));
dcps.. scps =e= sum(cf, dema(cf)*xf(cf)**demb(cf) - ddat(cf,"ref-t")*xf(cf)) - pcost;
u.up(cr) = crd(cr,"한계");
w.up(pt) = trb(pt);
슬롯 나라
refnlp '세계 참조 슬롯 나라 nlp' / 모두 /
reflp '세계 참조 슬롯 나라 lp' / all - dcps /;
xf.fx(cf) = ddat(cf,"ref-q");
lp를 사용하여 pcost를 최소화하는 reflp를 해결합니다.
xf.lo(cf) = 2;
xf.up(cf) = +inf;
nlp를 사용하여 scps를 최대화하는 refnlp를 해결합니다.