srcpm.gms : 단일 지역 비상 계획 슬롯 나라

설명

국제석유의 비상계획슬롯 나라(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를 해결합니다.