Evaluation of Characteristics of TOPMODELParameters using SCE-UA Method

九州大学学術情報リポジトリ

Kyushu University Institutional Repository

Evaluation of Characteristics of TOPMODEL Parameters using SCE-UA Method

田齋, 秀章

九州大学大学院生物資源環境科学府

平松, 和昭

九州大学大学院農学研究院

森, 牧人

九州大学大学院農学研究院

原田, 昌佳

九州大学大学院農学研究院

https://doi.org/10.15017/8857

出版情報：九州大学大学院農学研究院学芸雑誌. 61 (2), pp.261-272, 2006-10-27. Faculty of Agriculture, Kyushu University

バージョン：

権利関係：

――

!"#$ %&

'()$* + %&,- )$* ./012 345 6789:;<=>?@AB CD:EF GHIJKL,-)$* IM0NO CDPQRJ067ST+U J

VWX YZ[PQR'()$*+, -)$*\] [^_`ab +c d,-)$*!" IM0 efJ0[VW ghiDjk lm

% n , o p q J 0r sYDghtuvwx yzqcy{|0}~0 Y€!" ‚ƒ'()$

*+c„JY…†%‡Y€ˆ‰ŠIy‹Œ

…†%‡Ž,,‘’“”•+ , -)$*+c„&+J&[–y—˜ +%

™Y€'()$* [^IM0š›34I… cœyZžŸ0YZ5 ¡¢:£¤Ÿ0

¥¦X˜J,-)$* %™…Ÿ0 +§¨+0c©c%0ªY«

XI I¬­®¢¯ °

J+J&±²³´rµ¶!""#xŠ

·¸y¹|0+Iº°»O

¼½¾±² ¿À©nrÁ µ!""$xÂÃĆJ >?@AB ¿ )ÅÆ:Ç È’Éwr % sÊ

€%xqJÁµ ±² ËÌ K˜ÍÎÏÐy~J0

ÑÃÄI %&±² ¿%ÒÓc ÔÏ ÕÖyיØÙÚÛ'ÜÝa ÕÖy_`abÞÍÎßàáIM0&'%

r()**!sµ(â)**$sµ(â!"""x

%0>?@AB ÍÎÏy{J

>?@AB ãO,äåqc

>?@AB [P†%‡º°P˜Jæçc

>?@ABÍÎÏ è89 Œš›%%

0àáä{n é.êëìír x

î+)ï î!ð !+)!,!r!""+x

&'%yqJ >?@AB [Pñò

µ ¶ ó ô

õö ÷ ø ù

û ü ý õË µ þ

é.ë.ëêëÃÄëYXë ëÃÄ r!""+W$-!""+W,!#Dx

( &

( &'% . / .

01(/

. /

3 .4 2

5.( 6( %( ( 0( (' 6((-)!-$-)7

)é.ë.ëCëëYXë ëÃÄ

)3 . 4 2

5.( .( 5 (.5 0( ('

ú& .( r8 (93 ( ((:x

１解析対象流域

!"#$%&

'()*+,-./0 123045,-./

./6789: ;<=./0 .>? @ABCDEB;F G6:HIJ KLMHINOPQRST./ UVWX YZ ./ [D\]^_#` abcWdefg#0Z'Sh./67 ijklm'n'opqro04s Z ./ tu HvwxyuWzy{|

}uq~y~tu €‚ƒ„…†‡

ˆ‰#„…q~yŠ‹Œ

２解析対象データ

.2Ž XY&'<‘%&' Z `+’S\+“~.W”•

–—2Žj˜™v™š™W›œl“

žŸ6 hŒS'‹ŒsXY qO89 ¡¢£h%&'‹Œ

¤¥ WŒZ ¥¦¢§¨C

©hª—«¬­®’SM¯>+s0°±B­®

0 h²S'‹ŒZ’S³‹Œ

´µj04Bp.>¦¢¶

H·¸0¹º&'¶H·¸ WŒ

３流出解析手法

PT UA

Sq~»¼&'P

½ª¾S T¿HIŽÀÁ.>ÂÃ04 ¿HIŽÀÁ<‚.>pqr<Ä+ÅÆs 0WÇj‚ƒpqr„…HHk&'ÈÉÊË Ì HIŽÀÁ‹£Íh<Ä+ÎÏ PÐÑ/T +ÒÓÔdŽÀÁ‹ÕOZWÖ

׋ŒZ Âà Ø:./WÙÚÛÈÉÊ Ë0HkÜÈÉÊË Ý`3’S²S' PT’S£Í&'<

ÞßÝP<ÞØà3 !"#$TW¶

./‚ƒ„… áâ·ãB«ä­®W£Í<‚.

CDWâ·ã͎0å‘4Bpæçè 0éê LÞÈÉÊËWë

P T tìí

0<Ä+ .'Úí0‚?0 å¶ ‹Œ

PT

¾ î ï ð S

図 1 ()*+,-./

!"#$%&'(!"#$

)*+,-./012345 678&9%&:;<6 :)78 9*+,-:=

> >?@A:?BCDEF6

GHIJK#LIMNOPCQRS9TUV6:W X Y&Z?[I OP\]^_`ab?:Vcd

:^6^e8fg?hiDjklV6 abPmnV67:7opq79r oc6

s;\*+,-8t7?suv!"#

$;\V6:wxpc6

8su!"#$y 6t7 8zRS{TUV6:BCcEF6 8su[IOP\|:BCD 7?}FE?~V6:q7?€6q7TUV 6:9 ‚ƒ„6

8 f … † ‡ = ˆ

8‰s;\Š

6)s;\Š8

‹Xpc6t7?}F6:78ŒŽ

Pc6

8‡=ˆ

:‘dc8s;\‡=ˆ6t 789‡=ˆ’“:”•O^ ’

“–—c6:?˜DEF69su[I

OP\|^pd"7™š›œ8žŸ 0 c ¡¢9£ocdŠop

!"#$Pe¤69¥Cc6 9¦§Š¨8©ªOP¦]P«¬V6 8©ªO­Nb$

%$:{‘dcEF6

®¯°±

8²ŽP³´s

&,s''&s'

&A*+µ¡?W¶V6cp

*+µ¡u·¸®¯8ŒŽP<Ž

¹º³´s·¸®¯

»¼½¾·¸8³´s¿WX6 iÀP^Á¡8,sPÂIOPÃÄ c,s¡PÅXpc6³

´s8YC Æ:DEÇ£c·¸eÈ»¼½

¾8q7®¯c6

®¯HIuJK#LIÉ

³´s*+,-²%opÊ9Ë ³´s¡Ì}ͬ

»¼½¾¼½¾ 6H¼½¾8 ):V6³´sPÎÏÀ9½ ÐVcd,sP(¡V6A) H? , s P ( ¡ D :V6

¹ºs·¸®¯

s¡*+8s Y:DEÇ£c ()*+Ñ?}F 12345ÒÓÔÕ

図 2 ÖT

ޓroot zone (㓸ਛቯᢙဳ㧕

unsaturated zone (ಽᏓቯᢙဳ)

saturated zone (㓸ਛቯᢙဳ) EA R

EXi

⴫ጀ࿯ფ

࿾ਅ᳓

⾂⇐ Qb

Si SRMAX

SRZ

SUZi

UZi

!

"#$%&'()*

$%&+,%&-./0123 4563745$89:7;

<=>? $%&%&@A B9:CDEFGHIJKLM /NOPQR)*ST)*UV1 WXY45$8Z[\]

!T

^ _`abcCDEF >)*ST>deMTfgCDEF

%&h3ST+,ijkl m I9:3>?6n

+,ijkop q

KT

>rsSTtuTv5 1wxyz{|}~9:t%&€

P9:t‚}DƒCz{|}~9:t

„…23 †d‡S9:M w ˆ€>M

‰Š‹ ŒR>de3 flm IŽR

R>?3f

R3 ‘’“M O”•–>M

—^)*˜CDEF™š3

 ) * 1 w › œ ”  ž 3 Ÿ   ¡RŸ ¢ d”)*™š£¤bM>‰Š‹¥¦

§¨©ª«¬ŽR—^$%&)*

cCDEF¢™o

"#lm I)*

­®–3¥¦§¯°E±0

Rlm I>²³Tlm I´µ )* ™š

t I>¶·RCDEF¸

¹dº»¼1w9:t½¾¿dT CDEF‘’ÀÁ£<@ÂÃ9: t IÄ9:tCDEF­

‘’ÀÁ£<tCDEF‘’„…ÅÆ tÇdCDEFPM

" #XYÈÉ$8)*

—^ʚ+,%&$%&)*1 ˙%&'( 1ː

­XYÈÉ$8 )*

3)*UV1̔d)

*UV1R >h3Íʚ3Î1 Ï.•–>M

—^ʐ­ÐÑz{|}

~Z[“  ÊÒ —^Ó3

_`a

_`a)* !Ô

§ÕDÖ×Øٍ—gÚ ÛÜÝ

Ô§ÕDÖ×c?£_CDEF ÞߐÊÒàá£<£<£<âãg

šãžäCDEFÊaå>ÎæçCDEF bŽKdèCDEFw aåé?£_CDEFêëR CDEF£<@„…ÅÆ

ìK

íÞߣ<

íàá£<

£<@îï†C DEF­ï† Oð

t£<ñJKLt£<

@®òM£<óåôõöžR Ù ÷ ø ù

!"#

$%&'(

)*+,"#-'./01$234(

5678

9:;:<= *5678

*23>?@ ABCDEFG

))HIJ KLM5NO*PQR SHTU2V SHWXYZ [ HT\] ^

_S`S

)) Ha !NI SH !b.NO H !cdNO H.*e]f gJ HhN*Oij C HT\hN*NkV CH1l ^mWXYZ[$Nnop=q r4(=m)-'*stu@

!

! ))Hvw*xyz {H|}~

:€y;:V {‚ C Ha !M5NO Hƒ„V H…†-'*aV H‡J Hˆ

‡ H‰ HŠ‹‰Œ H a‹‰Œ @T\]*23$

@

))HlJ *Žg*a ! H a !NIKLM5NO Ha

!UJ ‘^

４モデルパラメータの感度分析

’“”$•–—˜™š*P›œ23]%&

žŸJ ¡@¢£J¤¥¦)) §*–—˜™š""q¨©ªy- '§*23]""""œ

«'(=¬­*–—˜™š*PQ £J" ® '

"

""

""

""

m¯°"$–—˜™š*PQR23

]*PQR*i"$23]PQœ±

²m)"$–—˜™šœ³f23

]q³f)´Xu=

５最適化手法

µ¶Q·<$¸¹º¨¸»4(¨­$¼ :½¾¿|<À–ÁÂy<?@+,ÃÄÅ<

¨­$À <?@¸ÆÃÄÅ<

^ÇÈ$ÄÅÉFÊ*ËVÌÍd*"#mÎ ÏÐmœ'ËVœÑ4ºmÄÅÉÒ Ó›4Ô=ºq*^)*<$ÄÅɜ

ÌÍd*ÕÖ*×ØÙ(ÚÛ(œÜ¢Ð¹µ <@!"#"$%–—˜™š*ÝÞ ¡

!"#$

%&'()"*+, -./01234&56*789$

2&:5;#)<

&:=>?789@AB CDEFG& @2H(*$

IJ2K6

DG @LMNOPQ

@RS"*TUVWXYZ [2\]Y^_789ABC5

;`ab2cdefgh 4:ijklmaC&n`ab2o pijklm2qrJ&"*stuv./w

&]uv./wxyz2{|}

_~€4o+-z2VW‚

ƒ" D G L MNOPQ56*„…}_†‡ˆ‰2K;~

J}(Š‹$IŒ(C&]"

+Ž‘’“”†•‘’“”T5–—˜™š

‘’“”›œlm2K;žŸ (¡~

J7:(5;‘’“”2K;žŸ

¢ £K&$IŒ¤¥*K 5

;¦§2¨U©"K;&ª«£

HŒC($IŒ2%

DG @<¬­€

a®¯Mijklm2$I…<¬­€

u5;CT<¬­€°2±²

³VWa`VWT¥a2 T5–2´5;µ¶

·)¸&./$IŒ¹º‹t(

*»¼½¾(¿ÀÁY./

aDG2jÂÃQ ÄÅ"`T¥~J2ÆÇ a2~J$2È&"*

}(#€4

É2aVW‚ƒ`V W a 2ÊË#&

Ì`VW2+Ž LMNOPQ56*„…

}_

Í4*VWÊÎ2XÏY^_a 2~J$2È&"*}(#

€4†‡ˆ‰2†‡ˆ‰ÐÑ´

}nÒÓ"o;¥n­€É„Ë

­€ÌT(*`VW2„…}_^V W ÊÎa2o J} 5;ÔÕ LMNOPQ u5;CT<¬­€2±² Ö×ØÙÚ2µ¶

C

Ö×ÊÎ`uÛ܂ݪ«

µ¶žŸ2ƒUÞ*

E F ß à

図 3 @<¬­€

㐿ᆎ

ࡄ࡜ࡔ࡯࠲ᢙn,㓸࿅ᢙp,㓸࿅ౝὐᢙmࠍ⸳ቯ

ೋᦼὐ⟲ࠍ↢ᚑߒ,㑐ᢙ୯㗅ߦਗߴࠆ

ోὐ⟲ࠍp୘ߩ㓸࿅ߦಽഀߔࠆ CCEࠕ࡞ࠧ࡝࠭ࡓߦߡㅴൻߐߖࠆ

෼᧤್ቯ

⚳ੌ

YES

NO

図 4 LMNOPQ<¬­€

SCEࠕ࡞ࠧ࡝࠭ࡓ߆ࠄ

SCEࠕ࡞ࠧ࡝࠭ࡓ߳

ㅴൻ↪ὐ⟲qὐࠍㆬᛯ ᦨᖡὐUએᄖߩὐ⟲ߩ㊀ᔃGࠍ᳞߼ࠆ

ኻ⒓ὐRࠍ᳞߼ࠆ.ត⚝▸࿐ᄖߦ಴ߚ႐ว

࡜ࡦ࠳ࡓߦ1ὐࠍㆬ߮Rߣߔࠆ

fR㧨 fU

ਛὐCࠍ᳞߼ࠆ

fC㧨 fU

࡜ࡦ࠳ࡓߦ1ὐࠍㆬ߮Uߣ⟎߈឵߃ࠆ ㅴൻᓟߩὐߣ⟎឵ߒ㑐ᢙ୯㗅ߦਗߴࠆ

ǩ࿁➅ࠅ㄰ߒߚ߆ UࠍRߦ⟎឵

UࠍCߦ⟎឵

YES

YES NO

NO

NO YES

! "#$%

&' &("#$)*+, -./0 ')*12 3+, 456789:;<

= >?@"#$ABC7 D7EF%

EF"< GH4IJ%

KLMNOP

''JAQRST12 7

="#$ UV,56WX

YJZ =<[1X?\]

^_UV,

+< `abMNOP

=JcCde56

WXYJZ =+EF>

%>"#$ UV< fgMNO

P

=JcCde56

WXYJZ =<[1X

?\]^_UV,

+< `abMNOP

56hijkl<

m,3+,nop

89cCde12 "

#$\/7@6q[

56 562Jijkl<

rstuvwxvyz{

J?|$QRG}+ rst

uvw~ €$

&( J ‚t uvw+,~

Aƒ„/ &k …†‡&7@' ˆ‰,2Š‹#

Œ7BCZ =QR#ABCŽk7

‘QRJC A’“i#@]Ž '

”•–—Š‹#˜,Z =7A' 'J~ˆ‰,&™…†‡?

|$QRš‰,rstuvw›œ$+

œ$ž<2’“"#UV~

‚‰7Ÿtuvw|‰,

%UV ¡+¢£ ¡¤G¥¦‰7

|‰§¨+©ª

１模擬発生流量を用いた適用結果

|«$¬y­tuvw®$

+¯,' +°±¡&(²³´µ‰7 JUV, ¡¢£ ¡+¶<

4¬y­tuvwrs+,' ?|$Q R‚&(J·‚J

~Š‹#A J¸Š‹12 #AJk _#¤«< QR~+kjkl,i

#iJ~_#A¹º»º#i+

,¼#½be i¾¿,À··

Á‚‰7Âғ"#$

+ÄŒ“i#"ƞ<‚J~’“

"#$A’“i#ÇÅÈ7ÉÊ,&kQRË

~ÌÍJ G,J~Î7ÏÐJ

’“"#$AÑÒÓ$+kQRm³ÌÔ ÕÖ,&kQRËAÌ7+~×7A7

_#&(QRŒ7

J~_#AJ G,

‚J~QRAJk‚ØA (k@Ù]7'J&7@_#

+,‚+ÚÛ,Ü+'ÝBCÞß

~à Á_#¹º»º#

i+,’“"#$+’“"#ÄŒ“i#

"ƞ<

UVáTŒ7

¹º»º#i‚jkl, i#i+,ÂÂUVáT~

âÙ+ 'J¹º»º#

i,¼#½be i¾¿,À··

§¨'ã@Ñ+äåÞßAà Á_#¹º»º#i+,’“"

#$+’“"#ÄŒ“i#"ƞ<

2Á‚‰7ÂÃ

¬y­tuvw?|$QR§¨4IJ¥¦æ ,ž<&Á>çèA®$™<

''饦æ,¬y­tuvwé¬y

­tuvwé¬y­tuvw?\$

é ¬ y ­ t u v w ? B $ J  ‚‰7tuvwê뒓

!"#$%&

'()*+,-./

0 123 !"45 6 7897 :;<;9= 0 >?@AB 7897 :;<

;9=0 >?CDE0 '()*+,-

.FGH?IJ4KLM

NN OPQ 0 >?

R*+,-.STUVWXYZKLM STUV[T\]^4_ N A`2STUV[T_2'()*+,-.

Habcde045 010 STUV[T'()*+,-.

fYZ_gh

i j k l

10^–3 10^–2 10^–1 10⁰

10^–3 10^–2 10^–1 10⁰

10^–3 10^–2 10^–1 10⁰

10^–3 10^–2 10^–1 10⁰

0 15000 30000

Nunber of function evaluation

Objective function value

0 2500 5000

Nunber of function evaluation

Objective function value

SCE–UA GA

0 15000 30000

Nunber of function evaluation

Objective function value

GA

0 22500 45000

Nunber of function evaluation

Objective function value

GA

ត⚝ὐ㧦300ὐ

➅ࠅ㄰ߒ࿁ᢙ㧦100࿁

⸘▚ᤨ㑆㧦160ᤨ㑆

୘૕ᢙ㧦50୘

਎ઍ੤ઍᢙ㧦100࿁

⸘▚ᤨ㑆㧦18ᤨ㑆

୘૕ᢙ㧦50୘

਎ઍ੤ઍᢙ㧦900࿁

⸘▚ᤨ㑆㧦160ᤨ㑆

୘૕ᢙ㧦300୘

਎ઍ੤ઍᢙ㧦100࿁

⸘▚ᤨ㑆㧦110ᤨ㑆

図 5 mno9pqmn=9

表 1 rs*+,-.t_tu

OP OP OP OP OP

図 6

0 1

0 1

0 1

0 1

m T

t

SR

RV SRZ

UZ

Normalized parameter value

m T

t

SR

RV SRZ

UZ

m T

t

SR

RV SRZ

UZ

m T

t

SR

RV SRZ

UZ

SCE–UA ត⚝ὐ㧦300ὐ

➅ࠅ㄰ߒ࿁ᢙ㧦100࿁

⸘▚ᤨ㑆㧦160ᤨ㑆

Normalized parameter value Normalized parameter value Normalized parameter value

GA

୘૕ᢙ㧦50୘

਎ઍ੤ઍᢙ㧦100࿁

⸘▚ᤨ㑆㧦18ᤨ㑆

GA

୘૕ᢙ㧦300୘

਎ઍ੤ઍᢙ㧦100࿁

⸘▚ᤨ㑆㧦110ᤨ㑆

GA

୘૕ᢙ㧦50୘

਎ઍ੤ઍᢙ㧦900࿁

⸘▚ᤨ㑆㧦160ᤨ㑆

２実測流量を用いた適用結果

! "

#$%&'$%(( )*+,-.- /0./123 2.-4,5.678'-9:,';&

&'< '=>

/?@A,'BCDEFGHI$%JKLM NOPQ'RSJK4TU VWX$%JK4X 4YZ[\]Z^_XBCDEFGHI `X

３による再現性

"#A_a'BCDEF GHIbc,da'ef 4ghL,'ijklmFn !ij opmFn<LM q4 <rsL_atuvwxyz{|

}~_aX4€‚ƒ„X…'

q†‡4,ˆ‰Š‹ŒŽ‘’“”•

l–—X˜™8a

š ›

œ‘žŸ ™¡_a¢Q.£4|

¤Q.£¥a¦aŸ§¨©>ª_a'BCD4,

@|¤Q.£«¬BCD

­®_aŸœ¯°±²³´µ¶

U · X | ¸ ¹ _ a º

!"#$%& #'(&%»‘²¯°¼—½¾¿

À,y‚Á&8(ÂÃaq ‘Äz«¬ÅzÆÇ¢Q.

£4,ȑ’y !‘ÄzÉÊËÌ '’Í !ÎÏÅ||Ð_a'mÑÒkÓ4|

¤Q.£4,È4ŸÔÕ qq4

*‘Äz¢Q.£Ÿ§Ö×XØْ

Ú Û Ü Ý

図 7 BCDEFGHIYZ[\

–0.4 –0.2 0 0.2 0.4

m t0 t_d SR_MAX RV SRZ0UZ_i,MAX

Relative Sensitivity

図 8 ?@A,'BCDEFGHI7Þ 0

1

0 1

m T0 t_d SR_MAX RV SRZ0UZ_i,MAX

m T0 t_d SR_MAX RV SRZ0UZ_i,MAX Normalized parameter value

SCE–UA

GA Normalized parameter value

図 9 ijklmFn4ijopmFnºœß¿

0 5 10

0 1000 2000 3000 4000

10^–2 10⁰ Rainfall (mm per 10min) Discharge (mm per 10min)

Time (10min–unit)

ታ᷹୯

⸘▚୯㧔GA㧕

⸘▚୯㧔SCE–UAᴺ㧕

!"#$%& ' ()*+(, - # ./

01230456718)9&:;

,(<=>)?@:;ABCD0=>1E F1GHIJ; KLMNOPQ ARM1STU$V;WAXY:;AZ[\

]1ST^_]`aAKLbcdO0e Of ghijk1KLMl

<((m 1n9&op qA

pq1rstuvw1xub)yzghijk O(m{|,rs}

'()0zA~;utuvw1xub) /:ghijk0O0{|,rs )}J;A)f (cdrsJz) yz-€A !"ghijk1KLM O )  ‚ 0 e ( ƒ - , ) J 19&o„(…†`J‡KLM NOEF)e 10,ˆ‰Š‹Œ,†1Ž0 e '()-€A7‘ ’“

” J • O L f A – ] 1—˜™špq '(0n9

&›w1œ)6ž:;

Ÿ  

! !"

#!!$%%&'(

#)*+,!))- #%"). '' //+")-//+"01)%"23)"/, ++%"#) ,)/) 2-! ’4•56

¡¢£¤¥”¦§¡¥¨©ª¨«¬¥­ ® '7(

.¯°±²on1³´µ¶²o·y¸¹º»

’¥•564

¡¢£¤¥¨©ª¨«¬¥­®¥¼½¾ ¿

ÀÁ ÃÄÅM]Æ1ÇÈÉ1bÊË

·y¸¹º»’•56&67

ÌÍΑ ''' 8ÇÈ1ÏйÑ%±18ÒÓ¥

ÔÒÓÑÕÖ×ØÙÚ(&Û

Ü.ÝÞ¥ßà á ''' ÇÈ8Ÿ¹ÕÖ×

ØÙÚ6(6&Û

“”â㥡¢£¤¥¨©ª¨«¬¥­äå # ./Á ‹Œ

1 !æçè·y¸¹º»’•5(6 锿ê¥ë¡ 쥔¦Üíî¥ï ðñ .

ò/1óô623õ1 1 LfÑóô6öf !"(1÷øÑ8Ÿ¥

8ùú¹û»’4•566'4

”¦Üíî '' kçü !"1ycd

¸ýþ¹ûŸZ’4•56

”¦Üíî ‰þ¹Á =

>Vþ¹ûÙÚ&Û

”¥¡¢£¤¥¨©ª¨«¬¥­äå 4 # ./1óô62 3·y¸¹º»’•5&676

OfA ghijk1²o &

! "

!"

# ! $ %&

# # &! ' &!

"# % ! (# )*+$ % # ),

$ ! )*+$

! %! !

% )*+$ % !