설명
이 슬롯 커뮤니티은 최적의 산림 관리 계획을 찾아 기존 산림을 전환합니다. 숲을 고도로 관리되는 숲으로 만듭니다. 다양한 형태의 지속수익률 조건이 부과될 수 있습니다.
대형 슬롯 커뮤니티 유형 :LP
카테고리 : 슬롯 커뮤니티 모델 라이브러리
메인 파일 : 슬롯 커뮤니티gms
$title 안탈리아 임업 슬롯 커뮤니티 - 동적(TFORDY,SEQ=62)
$onText
이 슬롯 커뮤니티은 기존 숲을 전환하여 최적의 산림 관리 계획을 찾습니다.
숲을 고도로 관리되는 숲으로 만듭니다. 다양한 형태의 지속수익률
조건을 부과할 수 있습니다.
Bergendorff, H, Glenshaw, P 및 Meeraus, A, 투자 계획
제지 산업 프로그램. 기술. 대표자, 세계은행, 1980.
키워드: 선형 계획법, 임업, 투자 계획, 산림 관리 계획
$offText
세트
c '상품' / 펄프로그, 제재목, 잔류물, 펄프, 제재목 /
cf(c) '최종 제품' / 펄프, 제재목 /
cl(c) '로그 유형' / 펄프로그, 톱로그 /
s '종' / nigra, brutia /
k '사이트 클래스' / 좋음, 중간, 나쁨 /
'나무 나이'에서 / a-10, a-20, a-30, a-40, a-50, a-60, a-70, a-80, a-90
a-100, a-110, a-120, a-130, a-140, a-150, a-160, a-170 /
u(at) '초기 연령' / a-10, a-20, a-30, a-40, a-50, a-60, a-70, a-80, a-90 /
p '프로세스' / 펄프-pl, 펄프-sl, 펄프-rs, 톱질 /
m '생산 단위' / 펄프 공장, 제재소 /
te '확장된 범위' / 기간-1*기간-12 /
t(te) '슬롯 커뮤니티 범위' / 기간-1*기간-9 /;
별칭(t,tp);
매개변수
scd(k) '사이트 클래스 분포' / 좋음 .25, 중간 .50, 나쁨 .25 /
토지 '사용 가능한 토지(1000ha)' / nigra 143.679, brutia 227.58 /;
테이블은(at,s) '기존 산림의 초기 연령 분포(비율)'입니다.
니그라 브루티아
(a-10,a-20,a-30) .2 .2
(a-40,a-50,a-60) 8.5 12.7
(a-70,a-80,a-90) 24.6 20.6;
테이블 yef(at,s,cl) '기존 숲의 수확량(ha당 m3)'
nigra.pulplogs nigra.sawlogs brutia.pulplogs brutia.sawlogs
a-10 38.8 1.2 17.8 3.2
a-20 48.4 8.6 16.8 19.1
a-30 43.4 26.6 15.4 32.6
a-40 34.4 45.6 16.0 41.0
a-50 27.8 59.2 18.2 46.8
a-60 28.5 66.5 16.8 53.2
a-70 27.3 72.7 17.5 55.5
a-80 25.2 79.8 17.3 54.7
a-90 26.4 83.8 17.0 54.0
a-100 27.1 85.9 16.3 51.7
a-110 27.8 88.2 14.9 47.1
a-120 28.3 89.7 12.0 38.0
a-130 28.8 91.2 9.1 28.9
a-140 28.8 91.2 6.2 19.8
a-150 28.3 89.7 3.1 9.9
a-160 27.8 82.2 1.2 3.8
가-170 27.1 85.9 .5 1.5;
테이블 ymf(at,k,s,cl) '관리림의 수확량(ha당 m3)'
nigra.pulplogs nigra.sawlogs brutia.pulplogs brutia.sawlogs
a-10.좋음 17.5
a-10.medium
a-10.불량
a-20.좋음 120.0 66.8
a-20.중간 95.0 51.1
a-20.poor 80.0 37.8
a-30.좋음 132.6 37.4 91.3 25.7
a-30.중간 120.2 14.8 81.4 10.0
a-30.poor 115.0 71.3
a-40.좋음 121.0 99.0 91.3 74.7
a-40.중간 115.5 59.5 86.5 44.5
a-40.poor 119.0 21.0 90.1 15.9
a-50.좋음 108.0 162.0 76.0 114.0
a-50.중간 112.0 108.0 77.0 74.0
a-50.poor 112.2 57.8 92.0 47.6
a-60.좋음 104.0 221.0 76.0 116.0
a-60.중간 106.0 159.0 76.0 113.0
a-60.poor 110.0 90.0 95.2 77.8
a-70.좋음 105.0 270.0 78.0 200.0
a-70.중간 98.0 207.0 72.0 153.0
a-70.나쁨 97.0 128.0 88.0 116.0
a-80.좋음 102.0 323.0 76.0 240.0
a-80.중간 105.0 235.0 80.0 177.0
a-80.나쁨 92.0 163.0 84.0 148.0;
매개변수
yw(te,at,s,cl) '기존 산림의 수확량(ha당 m3)'
yv(te,te,s,cl,k) '관리 산림의 수확량(ha당 m3)'
iad(at,s) '초기 연령 분포(비율)';
iad(u,s) = is(u,s)/sum(at, is(at,s))/100;
yw(t,at,s,cl)$u(at) = yef(at + (ord(t) - 1),s,cl);
loop(at, yv(t,t + ord(at),s,cl,k) = ymf(at,k,s,cl));
yv(te,te+3,s,cl,k) = ymf("a-30",k,s,cl);
yw, yv, iad를 표시합니다.
세트
wpos(u,te) '가능성'
vpos(t,te) 'v 가능성';
wpos(u,t) = yes$sum((s,cl), yw(t,u,s,cl));
vpos(t,te) = yes$sum((s,cl,k), yv(t,te,s,cl,k));
테이블 a(c,p) '입력 출력 행렬'
펄프-pl 펄프-sl 펄프-rs 톱질
펄로그 -1.0
톱로그 -1.0 -1.0
잔차 -1.0 0.4
펄프 .207 .207 .207
제재목 0.6;
테이블 b(m,p) '용량 활용도'
펄프-pl 펄프-sl 펄프-rs 톱질
펄프 공장 1 1 1
제재소 1;
매개변수
pc(p) '프로세스 비용(m3 입력당 미국 달러)' / (pulp-pl,pulp-sl,pulp-rs) 5.96
톱질 6.00 /
pd(cf) '판매가격(개당 US$)' / 펄프 147.0, 제재목 70.0 /
nu(m) '투자 비용(입력 m3당 US$)' / 펄프 공장 37.8, 제재소 61.5 /
age(at) '나무의 나이(년)'
avl(t,te) '식물은 생리 기간에 산다'
delt(t) '할인 요인';
스칼라
mup '식재 비용(ha당 US$)' / 150.0 /
muc '절감 비용(m3당 US$)' / 7.0 /
생명 '식물수명(년)' / 30 /
rho '할인율' / .1 /
sgm '자본회수율';
나이(at) = 10*ord(at);
avl(t,t) = 1;
avl(t,t-1) = 1;
avl(t,t-2) = 1;
sgm = rho/(1 - (1 + rho)**(-수명));
delt(t) = (1 + rho)**(-10*(ord(t) - 1));
연령, avl, sgm, delt를 표시합니다.
$sTitle 슬롯 커뮤니티 정의
방정식
efs(s,k,u) '기존 산림자원(1000ha)'
pfs(s,k,t) '식재된 산림자원(1000ha)'
lbal(cl,te) '로그 잔액'
bal(c,t) '목재 가공의 물질 수지'
cap(m,t) '목재 가공 용량'
sy1(te) '지속적인 수율 제약 - 모든 로그'
sy2(cl,te) '지속적인 수율 제약 - 로그 유형'
sy3(cl,te) '지속적인 수율 제약 - 펄프 통나무'
sy4(c,t) '지속적인 수율 제약 - 펄프'
wbnd '절단 제한'
ainvc(t) '투자비용'
aproc(t) '프로세스 비용'
asales(t) '판매 수익'
acutc(t) '비용 절감'
aplnt(t) '식재 비용'
이익;
변수
w(s,k,u,te) '기존 산림 벌채(연간 1000ha)'
v(s,k,t,te) '신규림 관리(연간 1000ha)'
r(c,te) '산업에 원목 공급(연간 1000m3)'
z(p,t) '프로세스 수준(연간 1000m3 입력)'
h(m,t) '생산능력 확장(연간 1000m3 투입)'
x(c,t) '최종 출하량(연간 1000개)'
phik(t) '투자비용(연간 1000us$)'
phir(t) '처리 비용(연간 1000us$)'
phix(t) '판매 수익(연간 1000us$)'
phil(t) '비용 절감(연간 1000us$)'
phip(t) '식재 비용(연간 1000us$)'
phi '총 이익(할인된 비용)';
양수 변수 w, v, r, z, h, x;
efs(s,k,u).. 10*sum(t$wpos(u,t), w(s,k,u,t)) =l= iad(u,s)*scd(k)*land(s);
pfs(s,k,t).. sum(te$vpos(t,te), v(s,k,t,te)) =l= sum(u$wpos(u,t), w(s,k,u,t))
+ 합계(tp$vpos(tp,t), v(s,k,tp,t));
lbal(cl,te).. r(cl,te) =e= sum((k,s,t), yv(t,te,s,cl,k)*v(s,k,t,te))
+ 합계((k,s,u), yw(te,u,s,cl)*w(s,k,u,te));
sy1(te-1).. sum(cl, r(cl,te)) =g= sum(cl, r(cl,te-1));
sy2(cl,te-1)..r(cl,te) =g= r(cl,te-1);
sy3("pulplogs",te-1).. r("pulplogs",te) =g= r("pulplogs",te-1);
sy4("펄프",t-1).. x("펄프",t) =g= x("펄프",t-1);
wbnd.. sum((s,k,u,te)$wpos(u,te), w(s,k,u,te)$((ord(u) + ord(te)) <= 5)) =e= 0;
bal(c,t).. sum(p, a(c,p)*z(p,t)) + r(c,t)$cl(c) =g= x(c,t)$cf(c);
cap(m,t).. sum(p, b(m,p)*z(p,t)) =l= sum(tp$avl(t,tp), h(m,tp));
ainvc(t).. phik(t) =e= sgm*sum(tp$avl(t,tp), sum(m,nu(m)*h(m,tp)));
aproc(t).. phir(t) =e= sum(p, pc(p)*z(p,t));
asales(t)..phix(t) =e= sum(cf, pd(cf)*x(cf,t));
acutc(t)..phil(t) =e= muc*sum(cl, r(cl,t));
aplnt(t).. phip(t) =e= mup*sum((s,k,te)$vpos(t,te), v(s,k,t,te));
이점.. phi =e= sum(t, delt(t)*( phix(t) - phik(t) - phir(t) - phil(t) - phip(t)));
슬롯 커뮤니티
antala '기본 사례' / all - sy1 - sy3 - sy4 /
antalb 'case b' / all - sy1 - sy3 - sy4 - wbnd /
antalc 'case c' / all - sy2 - sy3 - sy4 - wbnd /
antald 'case d' / all - sy1 - sy2 - sy4 - wbnd /
antale 'case e' / all - sy1 - sy2 - sy3 - sy4 - wbnd /
antalf 'case f' / all - sy1 - sy2 - sy3 - wbnd /;
$sTitle 대소문자 정의 및 보고서
매개변수
rp(*,* ,s,k) '새 포리스트의 순환 기간(%)'
tc(s,k) '총 신규 벌채량(1000ha)'
csum(*,t,s) '오래된 숲의 요약 요약';
phi를 최대화하는 lp를 사용하여 antala를 푼다.
tc(s,k) = sum((t,tp), v.l(s,k,t,tp));
rp("case-a",at,s,k) = (100*sum(t, v.l(s,k,t,t + ord(at)))/tc(s,k))$tc(s,k);
rp("case-a","total-cut",s,k) = tc(s,k);
csum("case-a",t,s) = sum((k,u), w.l(s,k,u,t));
phi를 최대화하는 lp를 사용하여 antalb를 해결합니다.
tc(s,k) = sum((t,tp), v.l(s,k,t,tp));
rp("case-b",at,s,k) = (100*sum(t, v.l(s,k,t,t + ord(at)))/tc(s,k))$tc(s,k);
rp("case-b","total-cut",s,k) = tc(s,k);
csum("case-b",t,s) = sum((k,u), w.l(s,k,u,t));
rp, csum 표시;