海岸工学 テキストhome.hiroshima-u.ac.jp/kiyosi/coastalEng.pdf海岸工学 テキスト Coastal Engineering 広島大学工学研究科 海岸工学研究室 Graduate School

海岸工学　テキスト

Coastal Engineering

広島大学工学研究科　海岸工学研究室Graduate School of Engineering, Hiroshima Univ.River, Estuarine & Coastal Engineering Laboratory

川西　澄Kiyoshi Kawanishi

Office: A2-423

E-mail: [email protected]

http://home.hiroshima-u.ac.jp/kiyosi

Tel & Fax: 082-424-7817

目　次

2.4 微小振幅波理論 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1

2.5 波長と波速（分散関係） . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2

2.7 水粒子の運動速度と軌跡 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

2.8 水中圧力 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

2.9 波のエネルギーとその伝達率 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

2.10 重複波 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10

2.11 有限振幅波. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

3.2 浅水変形 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

3.3 波の屈折 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

3.5 波の反射と透過 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

3.6 砕波 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

3.7 波の減衰 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18

4.3 高潮 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20

4.4 津波 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21

4.5 セイシュ・副振動と湾水振動 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

５ 波の統計的性質と波浪推算 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23

5.1 不規則波の解析方法 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23

5.2 波別解析法による不規則波の表示. . . . . . . . . . . . . . . . . . . . . . . . . . 23

5.3 スペクトル解析法による不規則波の表示 . . . . . . . . . . . . . . . . . . . . . . 24

5.4 風波の発生と発達 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

5.5 風波の発達過程 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

5.6 風波の推算法 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25

６ 沿岸海浜過程 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

6.2 海浜形状 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

6.4 波による底質の移動 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

6.5 漂砂量の算定 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29

6.6 海浜流とその発生メカニズム . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

— i —

2.4 微小振幅波理論

０．渦なし流れとみなせる理由［微小振幅波の仮定（7）の成立理由］

波による底面付近の流速は周期的に流向が反転するため，底面境界層で発生する渦

度も周期的に符号が反転し，正負の渦度が互いに打ち消し合うため，底面近傍を除い

て流れ場のほとんどは渦なし場（ポテンシャル流れ）とみなされる．

進行波の流れと底面渦度

１．基礎方程式と境界条件

（１）基礎方程式

教科書 16ページの微小振幅波の仮定（7）より，水平，鉛直速度は

u = −∂φ∂x, w = −∂φ

∂z(2.1)

となり，これらを連続式∂u

∂x+∂w

∂z= 0 (2.2)

に代入すると，次のラプラス方程式となる．

∂2φ

∂x2+∂2φ

∂z2= 0 (2.3)

また，非粘性の渦なし流れでは次式のベルヌーイの式（圧力方程式）が，流れ場のい

たるところで成立する．

−∂φ∂t+p

ρ+ gz +

1

2|∇φ|2 = f(t) (2.4)

ここで， f(t) は時間のみの関数であり， φ に含めて考えることも出来る．微小振幅波

では， |∇φ|2 を２次の微少量として無視することができる．(2.3)式から速度が，(2.4)

式から圧力が求まる．

（２）境界条件

底面での運動学的条件

　海底面通過する浸透流がなければ，底面での鉛直速度は −h のラグランジュ微分

で与えられる．すなわち

w = −DhDt.

底面が時間的に変化しなければ

w = −u∂h∂x

=∂φ

∂x

∂h

∂x. (2.5)

— 1 —

底面が水平なら

w = −∂φ∂z

¯z=−h

= 0 (2.6)

が底面での運動学的境界条件になる．

水面での運動学的条件

　水面の水粒子が水面から離脱しなければ，

w

µ≡ −∂φ

∂z

¶=Dη

Dt=∂η

∂t+ u

∂η

∂x≈ ∂η

∂t, z = η ≈ 0 (2.7)

水面での力学的条件

ベルヌーイの式（圧力方程式）から， |∇φ|2 を２次の微少量として無視し，大気圧

を 0とすると，水面 z = η では

η ≈ 1g

∂φ

∂t(2.8)

となる．

時間的に変化する水面境界を有する場を解析的に解くことは困難なので，微小振幅

波理論では，水面での境界条件が z = 0 で適用できるものと考える．

（３）速度ポテンシャルの解

波の観察から，速度ポテンシャルとして，次式の変数分離型の解が予想される．

φ = Z(z) sin(kx− ωt) (2.9)

これを (2.3)式に代入すると

d2Z

dz2− k2Z = 0 (2.10)

底面での運動学的条件と水面での力学的条件のもとで上式を解くと

φ(x, z, t) = −Hg2ω

cosh k(h+ z)

cosh khsin(kx− ωt) (2.11)

が得られる．

2.5 波長と波速（分散関係）

水面での運動学的条件 (2.7)式に (2.11)式の φ の解と η = H2 cos(kx− ωt) を代入す

るとωH

2sin(kx− ωt) = gkH

2ω

sinh kh

cosh khsin(kx− ωt) (2.12)

が得られる．これより

ω2 = gk tanh kh (2.13)

L =gT 2

2πtanh

µ2πh

L

¶(2.14)

繰り返し計算の必要がない近似式は

L ≈ gT2

2πtanh

(2π

sh

gT 2

Ã1 +

sh

gT 2

!)(2.15)

— 2 —

C = L/T = ω/k から

C =g

ωtanh kh =

gT

2πtanh

µ2πh

L

¶(2.16)

C =

rg

ktanh kh =

sgL

2πtanh

µ2πh

L

¶(2.17)

kh→∞ で tanh kh→ 1 なので (2.16)式より

C → gT

2π= C0 （深海波） (2.18)

L0 = C0T → gT 2

2π（深海波） (2.19)

従って，深海波の波速 C0 と波長 L0 は周期 T で決まり，水深の影響は受けない．

kh→ 0 で tanh kh→ kh なので (2.17)式より

C →pgh （長波） (2.20)

L = CT →pghT （長波） (2.21)

従って，長波の波速 C は水深 h のみで決まる（分散性のない波）．

【例題 2.1】の要点

分散関係を使って波の周期と水深から波長と波速を求める．

2.7 水粒子の運動速度と軌跡

速度ポテンシャル φ が (2.11)式のように求まっているので

dx

dt= u

µ≡ −∂φ

∂x

¶=H

2

gk

ω

cosh{k(h+ z)}cosh(kh)

cos(kx− ωt)dz

dt= w

µ≡ −∂φ

∂z

¶=H

2

gk

ω

sinh{k(h+ z)}cosh(kh)

sin(kx− ωt) (2.22)

分散関係 (2.13)式を使うと

u =Hω

2

cosh{k(h+ z)}sinh(kh)

cos(kx− ωt)

w =Hω

2

sinh{k(h+ z)}sinh(kh)

sin(kx− ωt) (2.23)

進行波の水面変動と水平流速

— 3 —

水粒子の運動軌跡

(2.23)式の連立常微分方程式を解けば水粒子の運動が求まるが，(2.23)式は右辺に未

知関数 x, z を含む非線形方程式であり解析解を求めることは困難である．

水粒子の平均位置を (x0, z0) とすると，微小振幅波では x − x0, z − z0 は微小で無

視できる．従って

dx

dt=Hω

2

cosh{k(h+ z0)}sinh(kh)

cos(kx0 − ωt)dz

dt=Hω

2

sinh{k(h+ z0)}sinh(kh)

sin(kx0 − ωt) (2.24)

これはすぐに時間 t で積分できて

x = x0 − H2

cosh{k(h+ z0)sinh(kh)

sin(kx0 − ωt)

z = z0 +H

2

sinh{k(h+ z0)sinh(kh)

cos(kx0 − ωt) (2.25)

ξ = x− x0, ζ = z − z0 とおくと µξ

A

¶2+

µζ

B

¶2= 1 (2.26)

ここで,

A =H

2

cosh{k(h+ z0)}sinh(kh)

, B =H

2

sinh{k(h+ z0)}sinh(kh)

以上の結果から，以下のことがわかる．

・波にともなう水粒子の運動軌道は楕円である．

・深海波 (kh > π) の変動振幅は水深方向に指数関数的に減少する．

・長波 (kh < π/10) の変動振幅は水深方向に変化しない．鉛直方向成分は底面から

の高さに比例する．

【例題 2.2】の要点

まず分散関係式から波長を求める．この波長を速度ポテンシャルから求まる速度と加

速度式（教科書の 2.30式）に代入して水粒子の速度，加速度を求める．

2.8 水中圧力

(2.7)式の 2次の微小項を無視した圧力方程式に，速度ポテンシャルの解 (2.11)式を

代入すれば，水中圧力が求まる．

p = ρ∂φ

∂t− ρgz

=1

2γ0H

cosh{k(h+ z)}cosh(kh)

cos(kx− ωt)− γ0z (2.27)

ここに， γ0 は海水の単位体積重量 ρg である．(2.27)式から圧力は，右辺第 1項の波

の運動によって発生する動圧 pd と第 2項の平均水位に対する静水圧 ps の重ね合わ

せで表されることが分かる．水面変動 η = 12H cos(kx− ωt) を用いて書き直すと，

pd = p− ps = Kpγ0η, Kp =cosh{k(h+ z)}cosh(kh)

(2.28)

— 4 —

ここで， Kp は圧力係数と呼ばれる．上式から，水中圧力を測定して， η を算定する

ことができる．

η =N(p− ps)γ0Kp

=Npdγ0Kp

(2.29)

ここで， N は補正係数で約 1.35.

【例題 2.3】の要点

周波数から周期を求め，分散関係式を使って波長（すなわち波数）を求める．波数が

わかると圧力係数が求まる．

— 5 —

2.9 波のエネルギーとその伝達率

波のエネルギー =位置エネルギー+運動エネルギー

2.9.1 位置エネルギー

単位底面積当たりの水柱の持つ位置エネルギーは

dEp1 =

Z η

−hρgz dz =

1

2ρg(η2 − h2) (2.30)

波のない静水時の位置エネルギーを dEp2 とすると

dEp2 =

Z 0

−hρgz dz = −1

2ρgh2 (2.31)

となる．波の位置エネルギー dEp は水柱の持つ位置エネルギーから静水時の位置エネ

ルギーを引くことによって得られる．

dEp = dEp1 − dEp2 =1

2ρgη2 (2.32)

規則波の水位変化 η = (H/2) cos(kx − ωt) を代入して一波長にわたって積分し，波

長平均値を求めると

Ep =1

2ρg1

L

Z x+L

x

H2

4cos2(kx− ωt) dx = 1

16ρgH2 (2.33)

となる．これが，単位幅当たりの波の（平均）位置エネルギーである．

2.9.2 運動エネルギー

単位底面積当たりの水柱の持つ運動エネルギーは

dEk =1

2ρ

Z η

−h(u2 + w2) dz (2.34)

単位幅当たりの波長平均値は

Ek =1

2

ρ

L

Z x+L

x

Z η

−h(u2 + w2) dxdz =

1

16ρgH2 (2.35)

となり，位置エネルギーと同じ大きさになる．

2.9.3 波の全エネルギー

前述の位置エネルギーと運動エネルギーを加えた単位幅当たりの波の全エネルギー

は

E = Ep +Ek =1

8ρgH2 = ρgη2 (2.36)

ここで， η2 は水位変動の２乗平均である．一波長間の総エネルギーは

EL =1

8ρgH2L (2.37)

で与えられる．

— 6 —

2.9.4 エネルギー伝達率（単位時間当たりのエネルギーの輸送量）

エネルギー伝達率は仕事率で表される．

仕事率 =力×速度

波の場合

力 = Pd（動圧）, 速度 = u（水粒子の水平速度）

であるから，一周期平均したエネルギー伝達率 W は

W =1

T

Z t+T

t

Z 0

−hPd u dzdt

=ρgH2

8

ω

k

1

2

½1 +

2kh

sinh(2kh)

¾(2.38)

となり，エネルギー伝達率はエネルギーと群速度の積になっている．すなわち，

W = ECg = ECn (2.39)

ここで， n は群速度と波速の比で

n =CgC=1

2

½1 +

2kh

sinh(2kh)

¾(2.40)

深海波（沖波）で n = 0.5，極浅海波（長波）で n = 1.0 となり，群速度は波速以下で

ある．

(2.39)式より 波のエネルギー E は群速度 Cg で運ばれる ことがわかる．

群速度とは

— 7 —

— 8 —

群速度の実験

間隔の少し異なる縦縞を描いた図AとBを重ねると第３の縞模様（スダレ，モアレ）

が見える．モアレの黒い部分は２つの波の山と谷が重なって波が消えた所，白い部分は

山と山，谷と谷が重なって，波の振幅が大きくなった所と解釈できる．振幅最大のとこ

ろが動く速度が群速度である．

さて，Aと Bの縦縞を平行に重ねたまま，縞に直交する方向に紙を動かして以下の

性質を確認する．

(1) AとBを同じ速度（位相速度）で動かせば，当然のことながら，モアレも同じ速度で

動く．このことは，波に分散性がなければ，群速度は位相速度に等しいことを示して

いる．

(2) Aを Bより速く動かせば，モアレは Bよりも遅く動く．このことは，波数の小さい

波の位相速度の方が，波数の大きな波の位相速度より速ければ，群速度は，両者の位

相速度より遅いことを示している．

(3) Aを Bより遅く動かせば，モアレは Bよりも速く動く．このことは，波数の小さい

波の位相速度の方が，波数の大きな波の位相速度より遅ければ，群速度は，両者の位

相速度より速いことを示している．

水の波に関する群速度の性質

・ 群速度は波速（位相速度）以下である．

・ 群速度と波速の比は水深の減少とともに増加し，1.0に漸近する（群速度は波速に近

づいて行く）．深海波では Cg/C = 1/2 である．

・ 波速が水深とともに一様に減少するのに対して，群速度は一時的に増加する．

水深による群速度の変化

【例題 2.4】の要点

波長と周期から波速を計算し，(2.40)式を使って群速度を求める．

— 9 —

2.10 重複波

重複波

2.10.1 完全重複波（定在波）

振幅が等しく，進行方向が反対の波が重なってできる．

η = η1 + η2 =H

2cos(kx− ωt)| {z }

x の正の方向へ進む波

+H

2cos(kx+ ωt)| {z }

x の負の方向へ進む波

= H cos kx| {z }x によって波の振幅が変化

cosωt = H cos2π

Lx cos

2π

Tt (2.41)

水位変動と水粒子速度の関係は進行波と異なり，

波の腹： kx = 0,π, 2π, · · ·x = 0, L/2, L, · · ·

水粒子速度 =水平速度の振幅は 0，鉛直速度の振幅は最大

波の節： kx = π/2, 3π/2, 5π/2, · · ·x = L/4, 3L/4, 5L/4, · · ·

水粒子速度 =水平速度の振幅は最大，鉛直速度の振幅は 0

完全重複波の水面変動と流速

— 10 —

2.10.2 部分重複波

η = η1 + η2 =H12cos(kx− ωt)| {z }

x の正の方向へ進む波

+H22cos(kx+ ωt)| {z }

x の負の方向へ進む波

=1

2(H1 −H2) cos(kx− ωt)| {z }

進行波

+ H2 cos kx cosωt| {z }定在波

=1

2(H1 +H2) cos kx cosωt| {z }

定在波

+1

2(H1 −H2) sin kx sinωt| {z }

定在波

(2.42)

最大波高 Hmax = H1 +H2, x =mL

2

最小波高 Hmin = H1 −H2, x =2m+ 1

2L

ここで， m = 0,±1,±2, · · ·

H1 =Hmax +Hmin

2, H2 =

Hmax −Hmin2

KR =H2H1

=Hmax −HminHmax +Hmin

Hmax と Hmin を測定すれば，反射率 KR が求まる（Healyの方法）．

— 11 —

2.11 有限振幅波

波の非線形効果を評価する指標

長波では，Ursell数 Ur =HL2

h3深海波では，波形勾配 H/L

— 12 —

2.11.3 波による質量輸送

微小振幅波では水粒子の運動は楕円軌道を描くが，実際の波では水粒子の運動は閉

じた楕円軌道からはずれ，水粒子は徐々に波の進行方向へ進んでいく．

【問題】

(1) 完全重複波（定在波）の速度ポテンシャルを求め，速度場が図 2.12に示したようにな

ることを確かめよ．

(2) 完全重複波（定在波）の位置エネルギーと運動エネルギーはどのように表されるか．

— 13 —

3.2 浅水変形

狭い意味での浅水変形とは，水深変化による波形の変化を意味する．

エネルギー収支は

(ECg)I − (ECg)II =Wloss (3.1)

エネルギー損失 Wloss を無視すると

(ECg)I = (ECg)II (3.2)

E = ρgH2/8 の関係を用いると，

HIIHI

=

s(Cg)I(Cg)II

(3.3)

断面 Iを深海域にとり Cg = nC の関係を用いると，深海波の n は 1/2だから

H

H0=

rC02nC

= Ks（浅水係数） (3.4)

の関係がある．分散関係 C = C0 tanh kh より浅水係数は

Ks =1√

2n tanh kh

=1q

tanh kh¡1 + 2kh

sinh 2kh

¢ (3.5)

となる．

【例題 3.1】の要点

分散関係⇒波長⇒波数⇒浅水係数⇒波高

3.3 波の屈折

波の進行方向（波向）と波速の変化する方向とが異なる場合，スネルの法則にしたが

って波向が変化する．sin θ1C1

=sin θ2C2

(3.6)

b1cos θ1

=b2

cos θ2(3.7)

(EbCg)0 = EbCg (3.8)

H

H0=

µCg0Cg

¶1/2µb0b

¶1/2= KsKr (3.9)

Kr =

µb0b

¶1/2=

µcos θ0cos θ

¶1/2[屈折係数] (3.10)

— 14 —

波高変化を浅水変形によるものと屈折によるものに分けると

H

H0=H

H 00

H 00

H0= KsKr (3.11)

H = H 00Ks, H 0

0 = H0Kr (3.12)

H 00 を相当（換算）沖波波高と呼ぶ．一般に波がどのような屈折経路をとって観測地

点に到達したかわからないので， H 00 の沖波が浅水変形のみを受けて波高 H の波に

なったと考える．(3.6)式から，屈折角は次式のように波速の比を使って求められる．

sin θ

sin θ0=C

C0(3.13)

等深線が平行な直線状の海岸に波が斜めに入射する場合，(3.6) と (3.7) 式から屈折係

数は

Kr =

µcos θ0cos θ

¶1/2=

µ1− sin2 θcos2 θ0

¶−1/4

=

∙1− (C/C0)2 sin2 θ0

cos2 θ0

¸−1/4=

"1 +

(1−

µC

C0

¶2)tan2 θ0

#−1/4(3.14)

上式を使えば，入射角 θ0 から屈折係数が求まる（屈折角を求める必要はない）．

【例題 3.2】の要点

分散関係から波長と波速を求め，屈折係数と浅水係数を求める．

— 15 —

3.5 波の反射と透過

3.5.1 規則波の反射率と透過率

透過性の構造物に波が入射する場合のエネルギー収支は

(ECg)I = (ECg)R + (ECg)T +Wloss (3.15)

ここで，添字 I,R,T はそれぞれ入射波，反射波および透過波に関する量を示す．

水深の変化がなければ，入射波，反射波，透過波の群速度は等しいから，

H2I = H

2R +H

2T +

µWloss18ρgCg

¶(3.16)

となる．入射波高に対する比で表すと

1 =

µHRHI

¶2+

µHTHI

¶2+

µWloss

18ρgH

2ICg

¶= K2

R +K2T +Kloss (3.17)

である．ここで， KR =反射率，KT =透過率，Kloss =エネルギー損失率である ．

3.5.2 砕波をともなう斜面（海浜）からの反射率

KR = χ1 χ2 (3.18)

ここで， χ1 :斜面の粗度や浸透性に関係する係数 ， χ2 :砕波の有無（斜面勾配と波形

勾配）に関係する係数 である．

係数 χ1 は，不透過滑面で χ1 = 1，不透過粗面で χ1 = 0.7～0.9，砂浜で χ1 ≈ 0.8，捨石で χ1 = 0.3～0.6 の値をとる．

係数 χ2 は，砕波によるエネルギー損失に関係し，Michieにより以下のように与え

られている．

砕波が起こる限界波形勾配 (H0/L0)max は次式で与えられる．µH0L0

¶max

=

r2β

π

sin2 β

π(3.19)

ここで，β は水平面と斜面とのなす角度（ラジアン）である．波形勾配が (H0/L0)max

より大きいとき，

χ2 =(H0/L0)maxH0/L0

(3.20)

と与えられる．

— 16 —

3.6 砕波

3.6.1 砕波の形式

〔崩れ砕波〕

勾配がゆるい海岸に，波形勾配の大きな波が入射する場合に見られる．波形の非対称

性は小さい．

〔巻き砕波〕

崩れ砕波の場合より急勾配の海岸に，波形勾配の小さい波が入射するときに発生す

る．波形の非対称性は大きく，エネルギー消散は急激に起こる．

〔砕け寄せ波砕波〕

さらに急勾配の海岸に，波形勾配の非常に小さい波が入射するときに発生する．波形

の非対称性は大きい．

3.6.2 砕波波高と砕波水深－砕波指標

沖波の条件 (H 00と T ) が与えられたとき，砕波がどのような条件 (Hbと hb) で起こ

るかを推定する．

浅海波：HbLb

= 0.142 tanh2πhbLb

(3.21)

深海波：

µH0L0

¶max

= 0.142 (3.22)

極浅海波：γ =Hbhb≈ 0.83 (3.23)

参考

深海波では kh→∞, tanh kh→ 1 から (3.22)式が得られる．

極浅海波では kh→ 0, tanh kh→ kh から

HbLb

= 0.142× 2πhbLb→ Hb

hb= 0.142× 2π = 0.89

砕波指標（教科書の図 3.20～図 3.23）

図 3.20： 沖波の周期（波長）→設計水深で砕波が発生するときの波高を求める．

図 3.21： 沖波の周期（波長）→設計水深で砕波が発生するときの水底から波頂まで

の高さを求める．

図 3.22： 沖波の波高と周期（波長）→砕波波高を求める．

図 3.23： 沖波の波高と周期（波長）→砕波水深を求める．

3.6.3 浅水変形と砕波条件

浅水変形計算によって求められる砕波点での波高と，砕波指標から求められる波高と

は一致しない．これは，浅水変形計算式は微小振幅波理論から求められたもので，有限

振幅の効果や砕波点付近での波形変化の効果が考慮されていないためである．

【例題 3.5】の要点

微小振幅波理論から求めた砕波水深と砕波波高を，砕波指標から求められるものと

比較する．

— 17 —

3.7 波の減衰

エネルギー損失により，波は伝播するにつれ波高が減少する．エネルギー損失の原因

としては，（１）内部粘性，（２）海底摩擦，（３）海底面での粘性浸透，（４）砕波，（５）

逆風などがあげられる．

3.7.1 海底摩擦による波高減衰

（１）波の乱流境界層

底面剪断応力 τb は

τb = fρub|ub| (3.24)

ここで， ub =底面での水粒子速度 ， f =底面摩擦係数 である．

摩擦によって失われるエネルギー Wf は底面剪断応力のなす仕事率として

Wf = τbub = fρ|ub|3 (3.25)

と表されるから，半周期平均のエネルギー損失は

Ef =2

T

Z T/4

−T/4τbub dt =

2

Tfρ

Z T/4

−T/4u3b dt

=4

3π2

fρ

sinh3 kh

µH

T

¶3=4

3πfρub

3 (3.26)

ここで， ub は海底面での水粒子の最大平均速度である．

エネルギーフラックスの保存則は

∂

∂x(ECg) = −Ef (3.27)

となり，この微分方程式を解けば，海底摩擦による波高の変化が求められる．波のエネ

ルギーと群速度を代入すると

∂

∂x(ECg) =

ρgH

8

k

ω

µ1 +

2kh

sinh 2kh

¶dH

dx= −Ef

µ≡ − 4

3πfρub

3

¶(3.28)

これに分散関係 ω/k = g/ω tanh kh = gT/(2π) tanh kh を代入すると

ρgH

8

gT tanh kh

2π

µ1 +

2kh

sinh 2kh

¶dH

dx= − 4

3πfρub

3 (3.29)

これを整理して，浅水係数の２乗

Ks2 =

∙tanh kh

µ1 +

2kh

sinh 2kh

¶¸−1(3.30)

を用いると

−dHH2

=64π3fKs

2

3g2 T 4 sinh3 khdx (3.31)

∆x の区間で水深が一定とみなし， x = 0 で H = H1 ， x = ∆x で H = H の条

件で積分すると

H

H1=

∙1 +

64π3fKs2H1∆x

3g2 T 4 sinh3 kh

¸−1(3.32)

（２）海底摩擦係数

Jonnsonは底面での摩擦応力 τb の振幅 τb と水平速度の振幅 ub の関係を

τb =ρ

2fwub

2 (3.33)

とおいた．ここで， fw は摩擦係数で Reynolds数 ubAb/ν と相対粗度で決まる．代表

長さ Ab は水粒子の軌道振幅， ds は相当粗度である．

— 18 —

3.7.2 浸透による波高減少

　水底が透水性の物質で構成されている場合，波がその上を進行すると水底の水圧

変化によって水の浸透または湧き出しが起こる．透水層内の流れは粘性の作用によって

エネルギーの一部を失う．単位面積，単位時間当たりのエネルギー損失 Ep は，

Ep =π

4ρg

KH2

L cosh2 2kh

1− exp(−8πD/L){1 + exp(−4πD/L)}2 (3.34)

ここで， D は透水層の厚さ， K は透水係数である．

透水層の厚さが D > 0.3L なら D の影響はほとんどなく，

Ep =π

4ρg

KH2

L cosh2 2kh(3.35)

と表される．

【例題 3.6】の要点

分散関係から水深 25 mでの波長（波数）を求め，(3.33)式から浅水係数の 2乗を計

算する．海底での摩擦損失を考慮した (3.35)式を用いて，波が 5 km進んだ時の波高を

求める．

— 19 —

4.3 高潮

１．発生原因

（１）気圧低下，（２）強風による海水の吹き寄せ

２．発生条件

（１）湾口が台風の襲来方向に開いている．

（２）台風の進行方向東側危険半円域に入る．

（３）湾奥に向かって強風が吹く．

（４）湾内の水深が浅い．

３．高潮の計算

（１）気圧低下による海面上昇

1

ρ

∂p

∂r+ g

∂η

∂r= 0 (4.1)

1

ρg(p∞ − p0)| {z }

∆p

+ (η∞ − η0)| {z }−ηps

= 0 (4.2)

ηps =1

ρg∆p (4.3)

mb = 103 dyn/cm2 だから，cm-mb単位系では

1

ρg=

1000

1.03× 980 = 0.991

（２）風の吹き寄せによる海面上昇

水面勾配による圧力勾配と剪断応力が釣り合っている場合，運動方程式は

−g ∂η∂x

+∂

∂z

µKz

∂u

∂z

¶= 0 (4.4)

境界条件は

ρKz∂u

∂z= τs　 for z = 0 (4.5)

ρKz∂u

∂z= −τb, u = 0　 for z = −h (4.6)

これらの条件のもとで運動方程式を z で積分すると，海面勾配は

I =∂η

∂x=τs + τbρgh

(4.7)

［海面に働く風応力と海底摩擦応力の合力に釣り合うように海面勾配が決まる］

海面に働く風応力 τs は

τs = γs2ρaU

2; γs2 ≈ 1.2× 10−3 (4.8)

海底摩擦応力 τb は

τb = γb2ρub

2; γb2 ≈ 2.6× 10−3 (4.9)

— 20 —

から計算される．ただし，海底付近の流速 ub が分からないので， τb は不明である．

τb が τs に比例すると考え， τb = λτs とおくと

∂η

∂x=

1

ρgh(1 + λ)γs

2ρaU2 (4.10)

0 ∼ F まで x で積分すると

η = γs2 ρaρ(1 + λ)

U2

hgF = k

F

hU2 (4.11)

【例題 4.1】の要点

水深が変化する場合は式 (4.10)を積分して海面上昇量を求める必要がある．

（３）砕波による海面上昇

Wave setupにより，

ηB = cH1/3; c ≈ 0.1 (4.12)

の海面上昇が生じる．ここで H1/3 は有義波高である．

（４）高潮の経験式

高潮による最大海面上昇量は半経験的に次式で推算される．

ηM = a∆p+ b U2 cos θ + cH1/3 (4.13)

ここで， ∆p = 最大気圧降下量 [mb] ， U = 最大風速 [m/s] ， θ = 海面上昇が最大とな

る主風向と最大風速時の風向とのなす角 [◦]， a, b, c =海域に固有な経験定数 である．

4.4 津波

津波は微小振幅理論の長波として扱うことが出来る．

(a)津波の変形

H2H1

=

µCg1Cg2

¶1/2| {z }

浅水変形

µb1b2

¶1/2| {z }

屈折

(4.14)

長波の群速度は Cg = C =√gh だから，

H2H1

=

µh1h2

¶1/4µb1b2

¶1/2(4.15)

上式の関係は Greenの法則と呼ばれる．

(b)津波の打ち上げ高と進行速度

一定水深と一定斜面からなる海岸に長波が入射するときの静水面からの打ち上げ高

R は教科書の理論式 (4.9)または実験式 (4.10)で求められる．

陸上に打ち上げた段波上の津波先端の進行速度 U は長波の波速から

U = kpgh (4.16)

— 21 —

で推定される．ここで， k は陸上地表面の条件によって決まる係数で，0.7～2.0程度の

値をとる．

【例題 4.2】の要点

沿岸等深線が凹状であるところにV字形の湾がある．この湾に津波が進入する．Green

の法則を使って，湾口の波高と湾内の波高を求める．

【例題 4.3】の要点

長波と考え，津波の波長を求める．平均勾配から斜面長 l を計算し，教科書の理論

式 (4.9)または実験式 (4.10)を使って打ち上げ高を求める．

4.5 セイシュ・副振動と湾水振動

閉鎖水域や一部開放水域には固有振動周期が存在し，これが外部擾乱の周期と一致

すると共振を起こし，波高が非常に大きくなる．この固有振動をセイシュまたは副振動

と呼ぶ．

鉛直壁では波が完全反射され，完全重複波が形成される．鉛直壁の位置は重複波の腹

の位置と一致する．この境界条件を満足するのは，

両端が鉛直壁の場合：水域の長さ l が 1/2波長の倍数となる

１端のみが鉛直壁の場合：水域の長さ l が 1/4波長の倍数となる

場合である．

一般に，セイシュ・副振動は長波と考えられるので，その周期は

T = L/C = L/pgh (4.17)

で与えられる．波長 L は上述したように，

両端が鉛直壁の場合：L = 2l/n n = 1, 2, 3, · · · (4.18)

１端のみが鉛直壁の場合：L = 4l/(2n− 1) n = 1, 2, 3, · · · (4.19)

である． k = 1 の場合が最も起こりやすい．(4.18), (4.19)式を (4.18)式に代入すれば，

固有振動の周期が求められる．ただし，１端のみが鉛直壁の場合の L は，回折効果に

より，(4.19)式で与えられるものより若干大きい．

エネルギー損失を無視した場合，湾口の波高 H0 に対する湾奥での波高 H1 の増幅

率 R は

R =H1H0

=1

cos kl(4.20)

である．共振状態では kl = 2πL l =

2π4l l =

π2 となるから， R →∞ ．現実には粘性など

の影響で有限にとどまる．

— 22 —

５ 波の統計的性質と波浪推算

不規則波では確率論的な取り扱いが必要となる．すなわち，平均値，分散（標準偏差）

などの統計量を使って不規則波を取り扱う．

5¯.1 不規則波の解析方法

(a) 波別解析法

１波ごとの波高や周期などの発生確率分布を算定する．非線形の強い波に適用され

る．

(b) スペクトル解析法

不規則波を様々な周期と振幅をもった成分波の合成と考え，各成分のエネルギー分布

によって不規則波を記述する．線形波に適用される．

5¯.2 波別解析法による不規則波の表示

5.2.1 波の定義法

(a) ゼロアップクロス法：水位が上昇しながら平均水位を切る時刻を区切りとする．

(b) ゼロダウンクロス法：水位が下降しながら平均水位を切る時刻を区切りとする．

5.2.2 代表波とその定義

N 個の波群中から，波高の大きい順に N/n 波を選び，その波高と周期を平均した

ものを 1/n 最大波高： H1/n ， 1/n 最大周期： T1/n と呼ぶ．また， H1/n と T1/n

を持つ波を 1/n 最大波と呼ぶ．

(i) 最大波 (Hmax, Tmax)

(ii) 1/10最大波 (H1/10, T1/10)

(iii) 1/3最大波 (H1/3, T1/3) 有義波

(iv) 平均波 (H, T )

5.2.3 波の統計的性質

(1) 水面変動 η の確率密度分布

水面変動 η は

p

µη

ση

¶=

1√2πexp

(−12

µη

ση

¶2)(5.1)

のガウス分布に従う．

(2) 波高 H の発生確率密度分布

波高 H は

p(H) =H

4ση2exp

(−18

µH

ση

¶2)(5.2)

の Reyleight分布に従う．

波群中で波高 H より大きい波高が出現する確率（超過発生確率 P ）は

P (H) =

Z ∞H

p(H) dH = exp

(−18

µH

ση

¶2)(5.3)

となる．

— 23 —

(3) 代表波間の関係

平均波高 H は

H =

Z ∞0

Hp(H) dHZ ∞0

p(H) dH

=√2π ση (5.4)

2乗平均波高 Hrms は

Hrms =

qH2=

⎧⎪⎪⎨⎪⎪⎩Z ∞0

H2p(H) dHZ ∞0

p(H) dH

⎫⎪⎪⎬⎪⎪⎭1/2

=√8ση (5.5)

有義波高との関係は

H1/3 = 1.597H = 1.416Hrms = 1.597√2πση = 1.416

√8ση = 4.004ση (5.6)

(4) 最大波高

波群中の最大波高は Hmax/H1/3 の度数分布の最多値と平均値によって予測される．

HmaxH1/3

の最多値 ≈ 0.706√lnN (5.7)

HmaxH1/3

の平均値 ≈ 0.706½√

lnN +γ

2√lnN

¾(5.8)

ここで， N =波群中の波の数，γ = Eulerの定数 (= 0.5772)．

一般に，構造物の設計などにあたっては，設計条件や構造物の重要度などを考慮して

Hmax = (1.6 ∼ 2.0)H1/3 (5.9)

が用いられている．

【例題 5.1】の要点

波高の超過確率を， Hrms =√8ση の関係を使って Hrms で表し，波高が 2Hrms 以

上となる確率を求める．

5.3 スペクトル解析法による不規則波の表示

不規則波は様々な振幅と周波数を持つ波が重なり合ったものと考えられる．不規則波

を周波数のふるいにかけて不規則波を構成している成分波を取り出し，各成分波の振

幅の 2乗値（エネルギー）を示したものをエネルギースペクトル，周波数 f を中心と

した単位周波数当たりのエネルギー（振幅の 2乗）を周波数エネルギースペクトル密度

関数と呼ぶ．

5.3.1 周波数スペクトルと代表波との関係

周波数の関数としてエネルギースペクトル密度を表したものが周波数エネルギース

ペクトル密度関数 S(f) である．従って，スペクトルと振幅の 2乗平均値の間には

η2 = ση2 =

Z ∞0

S(f) df (5.10)

— 24 —

の関係がある．上式でスペクトルから振幅の 2乗平均値が求まれば，代表波高は (5.4)

～(5.6),(5.9)式の関係を用いて計算される．

5.3.2 風波の周波数スペクトル

風波の周波数スペクトルに大きく関係する物理量としては，風速（水面での摩擦速

度），吹送時間と吹送距離がある．

5.3.3 代表的な風波のスペクトル

〔１〕Pierson-Moskowitzのスペクトル

外洋で完全に発達した風波に対するもので，吹送時間と吹送距離には関係しない．

〔２〕光易のスペクトル

吹送距離が有限の場合，いいかえれば発達段階の風波のスペクトルを表現したもの．

5.4 風波の発生と発達

5.4.1 波の発生と発達の理論

静かな水面に微風が吹くと水面に局所的な乱れが生じ，風速が増加するとさざ波へ

と発達する．風波の発生，発達機構を説明する理論として，フィリップスの共鳴理論，

マイルズのせん断流理論がある．

5.5 風波の発達過程

風波の発達は，風速 U ，吹送時間 tw と吹送距離 F に支配されている．

5.5.1 代表波の変化

風域の風上端 A で発生した波が，風域の風下端に到達するのに要する時間 tA よ

り長く風が吹き続けても，風波の波高と周期は一定となり風波の発達は止まる．すなわ

ち，風波の波高と周期は風域の大きさ（吹送距離）によって制限される． tA を最小吹

送時間という．最小吹送時間は風上端 A で発生した波群が，風域の風下端に達する時

間であるから，

tA =

Z F

0

dx

Cg(5.11)

で計算される．

t < tA の場合には，風波の発達は吹送時間 t によって制限される． tC(< tA) の時

間風が吹き続けたとすると，波群の伝播距離は

— 25 —

Fmin =

Z t=tC

0

Cg dt (5.12)

で与えられる．これを最小吹送距離と呼ぶ．実際の風域の大きさ（吹送距離）がこの最

小吹送距離より大きい場合，風波の発達に関係するのは (5.12) 式で計算される最小吹

送距離である．

5.6 風波の推算法

(1) 風域の設定

(2) 波の発達と減衰の計算

5.6.1 風域の設定

(1) 地上天気図・気象データ

(2) 計算期間の決定

(3) 天気図から傾度風の計算

(4) 傾度風から海上風の推定

〔１〕地衡風と傾度風

地衡風と傾度風は等圧線に沿って吹く．

地衡風：気圧傾度と Coriolis力が釣り合っている．

傾度風：気圧傾度，Coriolis力と遠心力が釣り合っている．

〔２〕海上風

地衡風と傾度風は，海面近くでは海面での摩擦力の影響を受け，海上風となる．

海上風：気圧傾度，Coriolis力，遠心力，摩擦力が釣り合っている．

5.6.2 SMB法（有義波法）

波の発達は風速 U と吹送距離 F または吹送時間 t のどちらかによって支配され

るので， U と F ， U と t の組み合わせに対して H1/3 と T1/3 を観測資料を整理

して作成された図表から求め，小さい方の値の組を推算値として採用する．

【例題 5.3】の要点

SMB法による風波の予測曲線を用いて，風速と吹送距離，風速と吹送時間の２つの

組み合わせから得られる有義波高のうち小さい方を選ぶ．（吹送距離と吹送時間のうち

のどちらが風波の波高と周期を規定しているか考える）

【例題 5.5】の要点

SMB法との波高差 0.3mは， T1/3 から L1/3 を求め，波数が得られるので，式 (3.32)

を用いて計算できる．

— 26 —

６ 沿岸海浜過程

海浜での底質移動現象と移動している底質自体をともに漂砂 (littoral drift) という．

波・流れと漂砂・海浜地形の相互作用，それらの変化過程を総合して，海浜過程 (coastal

processes)と呼んでいる．

6.2 海浜形状

沖波波形勾配 H 00/L0 と底質の粒径と比重が海浜縦断形状を決める主要因である．す

なわち，

・粒径 d の大きい底質ほど安定力が大きく，海浜勾配は急になる．

・同一の底質粒径に対しては波の大きい海岸の方が前浜勾配がゆるい．また，波が小

さく粒径が大きいほど前浜勾配は急となる．

〔静穏な波浪条件〕

汀線が前進し，前浜勾配が急になる．

〔高波浪条件〕

前浜が削られ，砕波点付近に堆積（沿岸砂州の形成）する．前浜勾配は緩やかになる．

6.4 波による底質の移動

〔粒子の抵抗力〕

FR = (W − FL) tanφ ≈W tanφ (6.1)

ここで，W = π(ρs − ρ)gd3/6（粒子の水中重量）， FL =粒子に働く揚力 ， φ =静止

摩擦角，d =粒径 である．

〔粒子に働く流体力〕

τbm = fρubm2 (6.2)

ここで， f は摩擦係数， ubm は底面付近の最大流速である．粒子に働く水平力 FT

は

FT ≈ τbmπ

µd

2

¶2=π

4fρd2ubm

2 (6.3)

(6.3)と (6.1)式の比をとると

FTFR∝ fρubm

2

{(ρs − ρ)}gd tanψ =(u∗bm)

2

sgd

1

tanψ(6.4)

ここで， u∗bm =pτbm/ρ =

√fubm （最大底面摩擦速度）， s = (ρs − ρ)/ρ （粒子の水

中比重）である．

ψm =(u∗bm)

2

sgd

を Shields数と呼ぶ．この無次元数は，海底に作用する底面剪断力と底質の剪断抵抗の

比で，Shields数が 0.05～0.06の大きさの時，底質が移動を開始する．これを限界Shields

数と呼ぶ．

摩擦係数 f は２つの Reynolds数の関数と考えられる．すなわち，

f = a

µubmd

ν

¶p µubmδ

ν

¶q(6.5)

— 27 —

境界層の厚さ δ は式より δ ≈³2νω

´1/2=³2νTπ

´1/2だから，

f = a

µubmd

ν

¶p "ubm

µ2

νω

¶1/2#q= a

µubmd

ν

¶p "ubm

µT

νπ

¶1/2#q(6.6)

粗面乱流では f は ν に関係しないから， p = −q/2 となる．従って，

f = a

µubmT

πd

¶q/2(6.7)

以上のように，基本的には，底質の移動限界は Shields数とReynolds数によって決定

されるが，さらに加速度の効果等を考慮する必要がある．また，工学的には，海岸と底

質の特性と入射してくる波の特性がわかったときに，底質の有意な移動はどの程度の水

深から始まるのかを推定する必要がある．このような観点から，以下の予測式が提案さ

れている．H 00

L0= α

µd

L0

¶nµsinh

2π hiL

¶µH 00

H

¶(6.8)

ここで， hi = 移動限界水深 ， α, n = 経験定数 である．なお， LとH は hi から決

まる．

【例題 6.1】の要点

水深 hi を仮定し，分散関係式から波長 L を求め，浅水係数から波高 H を計算す

る．なお，沖波波長 L0 は与えられた周期から求める．得られた諸量が移動限界予測

式を満足しているかどうかチェックし，満足していない場合は水深 hi を仮定し直して

再計算する．

— 28 —

6.5 漂砂量の算定

6.5.1 掃流漂砂量

qsw0d

= 12.5

µτbmsρgd

¶3= 12.5ψm

3 (6.9)

ここで， qs =半周期間の単位幅，単位時間当たりの平均掃流砂量，d =底質粒径，w0 =

底質粒子の沈降速度，s =底質の水中比重，g =重力加速度 ．最大底面剪断応力 τbm は

τbm = fwρubm2/2 で計算される．

6.5.2 浮遊漂砂量

浮遊砂濃度 c は，通常の移流拡散方程式に砂粒子の沈降フラックスを加えた次式を

解いて求められる．

∂c

∂t+

∂(uc)

∂x+∂(vc)

∂y+∂(wc)

∂z| {z }移流フラックスの発散

=∂

∂x

µεx∂c

∂x

¶+

∂

∂y

µεy∂c

∂y

¶+

∂

∂z

µεz∂c

∂z

¶| {z }

拡散フラックスの発散

+ w0∂c

∂z| {z }沈降フラックスの発散

(6.10)

周期平均した流れ場と濃度場が定常で水平方向に一様とみなせれば，

d

dz

µεzdc

dz

¶+ w0

dc

dz= 0 (6.11)

参考：保存方程式

簡単のため，１次元の場合を考える．

c(x, t) を濃度などの単位体積当たりの物質量，f(x, t) を物理量のフラックスとする．

区間 [a, b] における物理量の変化率は

∂

∂t

Z b

a

c(x, t) dx = f(a, t)− f(b, t)

となる．

f(a, t)− f(b, t) = −Z b

a

∂f(x, t)

∂xdx

から∂

∂t

Z b

a

c(x, t) dx = −Z b

a

∂f(x, t)

∂xdx

区間 [a, b] は任意であるから，保存則は

∂

∂tc(x, t) +

∂

∂xf(x, t) = 0

— 29 —

の様に表される．3次元の場合も同様に考えて，

∂c

∂t+∇ · f = 0

時間変化率＋フラックスの発散＝０

上式が述べていることは， c の時間変化はフラックスの収束，発散によって引き起こ

されると言うことである．

化学，生物過程を含み， c が保存されない場合は

∂

∂tc(x, t) +

∂

∂xf(x, t)− P +B = 0 (6.12)

ここで， P は物質の生成率， B は消滅率である．

フラックスの種類：

・移流フラックス

・拡散（分散）フラックス

・沈降フラックス

6.6 海浜流とその発生メカニズム

6.6.1 海浜流

海流や潮流に加えて沿岸域では一般に密度流や吹送流が発生する．さらに，砕波帯付

近では波が直接の原因となって生ずる海浜流と呼ばれる流れが存在する．沿岸部に存在

する流れは，海岸流（海岸線に平行：海流，潮流）と海浜流に大別される．海浜流には，

沿岸流（汀線に平行）と離岸流（砕波帯を貫通して沖方向）がある．

6.6.2 Radiation応力

波の存在によってもたらされる過剰な運動量流束は

d(Me)x = (p+ ρu2)dz − ρgzdz (6.13)

全水深に対する波の１周期平均を計算すると

(Me)x =1

T

Z T

0

Z η

−h(p+ ρu2)dzdt− 1

2ρgh2 (6.14)

となる．運動量流束は応力であるから， (Me)x は，波の進行方向に垂直な単位幅の鉛

直面に働く応力［過剰運動量流束（比力×単位重量）］である．波の進行方向に平行な

面に働く応力は

(Me)y =1

T

Z T

0

Z η

−hpdzdt− 1

2ρgh2 (6.15)

となる．こうした波の運動にともなって発生する応力を radiation応力と呼び，テンソ

ル表示によって Sij と表す． (Me)x = Sxx ， (Me)y = Syy である．

微小振幅波の水粒子速度と水中圧力を代入して積分すると

Sxx = E(2n− 1/2), Syy = E(n− 1/2) (6.16)

— 30 —

が得られる．極浅海域では n = C/Cg = 1，深海域では n = 1/2 だから，

極浅海域では Sxx =3

16ρgH2 (6.17)

深海域では Sxx =1

16ρgH2 (6.18)

となる．Radiation応力は波高の２乗に比例して大きくなる．

6.6.3 平均水位の変化

等深線が汀線に平行で一定勾配の直線海岸に，波が直角に入射する場合を考える．入

射方向を x とすると，運動量保存則は

Fx+dx − Fx − Fb = 0 (6.19)

となる．ここで

Fx = Sxx +1

2ρg(h+ η)2 (6.20)

Fx+dx = Sxx +1

2ρg(h+ η)2 + dx

d

dx

½Sxx +

1

2ρg(h+ η)2

¾(6.21)

Fb =pdx

cosβsinβ = ρg(h+ η)

dh

dxdx (6.22)

これらを式 (6.19)に代入すると

dSxxdx

+ ρg(h+ η)dη

dx= 0 (6.23)

η ¿ h ならdη

dx∼= − 1

ρgh

dSxxdx

(6.24)

砕波帯外では，浅水変形によって波高が増大するので，radiation 応力 Sxx も波の

進行方向に向かって増加し， dη/dx < 0 となって平均水位は岸方向に低下する wave

setdownが生ずる．砕波帯内では，砕波によって波の進行とともに波高が減少するので，

dSxx/dx < 0， dη/dx > 0 となって平均水位が上昇する wave setupが生ずる．

— 31 —