海水導入工を利用した港内埋没防止技術の開発
独立行政法人水産工学研究所

水産総合研究センター
水産土木工学部

水理研究室

中山哲嚴

水産土木工学部

水理研究室

新井雅之

水産土木工学部

開発システム研究室 大村智宏

１．緒言
砂浜の漁港は、漂砂による埋没に悩まされている。沿岸漂砂系を遮断せずに漁港への漂砂侵入を抑制する事は、漁港機能を維持できるだけでなく、周辺の砂浜保全にも大きく貢献する。本研究は、港内埋没の主要なメカニズムとして、高濃度浮遊砂の長周期波等による港内への移流・拡散・沈降過程に着目し、漁港内への強制的な海水導入を行うことにより、浮遊砂の港内への侵入を抑制する工法を開発することを目的とする。
２．調査方法検討は以下の項目について行った。
２．１．海水導入工による港内埋没効果把握、適切配置（平成１６～１７年）
平面水理模型実験によって海水導入工による効果を検討した。波浪、サーフビート、海浜流による浮遊砂の挙動を明らかにするための数値計算手法を開発し、海水導入工による港内埋没防止効果について検討した。
２．２．海水導入工への砂侵入防止（平成１７～１８年）
大縮尺模型による海水導入工周辺の流体現象、浮遊砂の侵入状況を明らかにし、浮遊砂を最小限に抑えるための構造の検討を行った。同時に VOF 法による流体計算をベースに浮遊砂の挙動を計算する手法を開発し、その再現性について水理模型実験結果と比較した。数値計算の再現性を比較検討後、現地スケールによる海水導入工への浮遊砂侵抑制のための構造検討を数値計算により、行った。
２．３．港内埋没評価モデルの現地適用性（平成１７～１８年）
新たに開発した浮遊砂による港内埋没評価のための数値計算手法を実際に港内埋没に苦慮している漁港に適用し、その再現性について検討した。
３．調査結果３．１．平面水理模型実験
平面水槽に縮尺 1/40 の小規模な漁港を作成し、海水導入工の有無によって、漁港周辺の流動がどのように変化するかを明らかにした。同時にビーズを浮遊砂挙動に模し、波浪作用後のビーズの残留率によって、海水導入工の効果を把握した。

75

水理模型実験はピストン型不規則波造波機を有す
導水工
37.5cm×2 基る平面水槽（幅 10m×長さ 35m×深さ 0.8m）で行った．図-1 に示すように，1/20 勾配地形に導水工を有

主防波堤先端部
水深 h=-14.2cm

1/20

港口幅75cm

37.5
消波工

するモデル漁港を設置して固定床とした．また，導水外港港口部

外 港工には山本ら（1987）によって開発された潜堤方式を内 港用いた．潜堤の諸元は，法勾配 1/2，天端幅 20cm，潜堤長 37.5cm，遊水部幅 17.5cm，導水口の直径 3.5cm

内港口幅62.5cm

泊 池
水深-7.5cm

200cm

内港港口部

1/20

である．導水量の効果を比較できるようにあらかじめ2 基分設置してある．潜堤天端は静水面と一致させており，導水工の沖側から港内へビーズが侵入しないよ
うに 1.2mm×1.2mm の網を遊水部の上部全面に取り

汀線400cm
消波材（ヘチマロン）
岸壁（直立壁：天端高 S.W.L.+5cm）高さ 10×幅 5×長さ 225cm

図－1

付けた．実験波は，修正 Bretschneider・光易型スペクトルの 1 方向不規則波（有義波高 9.0cm，有義周期

水理模型実験で用いたモデル漁港採取範囲1.58s）とした．ここでの模型縮尺は 1/40 を想定してい

投入位置

AB

る．

25cm

漁港内外の流況や渦の挙動を可視化して調べるため，Cビーズを用いた．造波開始 2 分後にビーズを図－2 に示D1m

す外港港口部の港外側及び港内側と外港部中央に投入1m

図－2

してビデオ撮影した．また，投入 30 分後の漁港内外のビーズ採取範囲及び投入位置

ビーズ量を図中に示す領域ごとに採取し，計量した．
その結果、図-３より、港奥部（図-２で C の領域）では、海水導入量が増加すると著しくビーズの量が減少し、海水導入工が有効であることがわかった。
実験において、漁港の防波堤の隅角部や先端に発生する渦がビーズの輸送に大きく影響していることがが海水導入工によってどのように変化するかを検討するために、電磁流速計を多点配置し計測した。その結果、海水導入工によって港内側にビーズの移動を促進する渦成分は小さくなることがわかった（写真１、図-４参照）。

採取重量/投入総重量

目視観測によって確認されたので、その周辺の流動場

0.5

導入量00.4

1000 (cm3/s)
減少0.30.20.10.0A

B

図－3：導入量 0cm3/s：導入量 1000 cm3/s

500

C

ビーズ残留量

O
Ｌ

回転方向
40cm

Ｍ
P
Ｋ

回転方向

図－４

渦の移動経路の違い

写真－1

76

渦によるビーズの輸送D

３．２．平面の数値計算
まず、サーフビービートが港内への浮遊砂輸送に大きな影響Exp.

を与えていると考え、サーフビートによる浮遊砂輸送を明らか

内港

内突堤

外港

内港

内突堤にすることとした。非線形長波方程式と非平衡性を考慮した準外港

３次元浮遊砂モデルを組み合わせて、漁港周辺の浮遊砂輸送を計算により検討した。その結果、防波堤先端や港内の波除堤周辺で顕著な堆積侵食が生じた。海水導入工がある場合には港奥Cal.

への浮遊砂輸送は激減した。しかし、漁港内は波浪、サーフビート、海浜流が共存しており、サーフビートのみでは定量的に評価できないことがわかった。そこで、波浪場、サーフビート、海浜流を同時に計算出来る非線形分散性波動モデルの一つである修正ブシネスク波動方程式と３次元浮遊砂モデルを組み合わせた新たな港内埋没評価モデルを開発した。本モデルの特徴は、
波、サーフビート、海浜流を一つ図－５短周期波高分布(H1/3/Ho’)

サーフビートのみの場合ほど効果がないのモデルで表現出来ること、浮遊

港周辺の波高分布、サーフビート

内港内突堤

-100B

A

Ｔ型突堤
導水工2基入工の効果を評価した。まず、漁-50
Ｔ型突堤導水工0基

験に用いた漁港に適用し、海水導0通常突堤導水工2基ている。このモデルを水理模型実

50

通常突堤
導水工1基

従来のモデルとは大きく異なっ

A.内港泊地B.外港航路

100

通常突堤
導水工0基

浪場に反映させている点であり、150

堆積量 [m3]

砂輸送による地形変化を逐一波

外港

振幅分布、海浜流流分布及び水位効果的である

変動のスペクトル分布について、
水理模型実験結果と比較し、流動非線形分散波
導水工なし・通常突堤場の再現性を検討した。
その結果、

単位：10 cm

地形変化量
導水工2基・通常突堤

単位：10 cm

航路

航路

数値計算結果は若干過小評価傾向があるものの十分な精度を有（スポンジ層）（スポンジ層）していることがわかった（図-５

海水導入

内港（スポンジ層）

外港

漁港周辺の地形変化を図-6 に示

地形変化量
導水工なし・Ｔ型突堤す。同図上段は、各ケースの港奥航路

単位：10 cm

内港

内突堤

同モデルによって、計算された内突堤参照）。

（スポンジ層）

外港

地形変化量
導水工2基・Ｔ型突堤

単位：10 cm

航路部（A）、航路部（B）の体積量を示している。
これから海水導入工（スポンジ層）（スポンジ層）２基＋T 型突堤（港奥部の波除堤

海水導入

内港

（スポンジ層）

外港

内港

内突堤抑制効果があることがわかる。同

内突堤を T 型にする）が最も港内埋没図下段に示した地形変化からも同様な結果を得た。

図-6 評価モデルによる漁港周辺の地形変化

77（スポンジ層）

外港

３．３．導入水に含まれる砂を最小限にするための仕切壁設置

大縮尺水理模型実験

吸水口改良
160

が平面水理模型実験及び数値計算により、明らかに1405.8された。この評価の前提条件は、導入水に多量の砂が混入していないことである。このためには、通常

1：2
潜堤部排砂口（穴あけ）の海水導入工では十分ではなく、導入水に含まれる浮遊砂濃度を極力少なくすることが不可欠である。

図-７

導水管

海水導入工の断面図

このために大縮尺の２次元断面実験を行った。想定[ H=50cm，T=4.0s ]海水導入工より沖側に中央粒径 0.17mm の砂を敷き詰め、波高 50cm、周期４秒程度（現地換算 2.5m、9 秒）の規則波を作用させ、海水導入周辺の波高、各ケース/改良なしれた潜堤方式とし、その断面形状を図-７に示す。

平均導水流速2.02.0

改良なし
導水管上向き仕切壁
潜堤穴あけ

1.5

各ケース/改良なし

縮尺は 1/5 で、対象波は現地スケールで 2.5m の波を作用させた。海水導入工の形式は水工研で開発さ

14.68

直立壁

海水導入工による港内埋没防止効果は高いこと

1.00.50.0

流動、浮遊砂濃度を計測した。また、導水管背後に

平均浮遊砂濃度
改良なし
導水管上向き仕切壁
潜堤穴あけ

1.51.00.50.0

図-８

捕砂器を設置し、波浪作用時間に導水により輸送さ模型実験結果(流速・浮遊砂濃度)れる全浮遊砂量を計測した。
図-８に各ケースでの平均導水流速及び導水に含まれる浮遊砂濃度（改良無しの場合で無次元化している）を示す。①改良しない場合（直線の導水管）には、導入水の浮遊砂濃度はそれほど大きくないもの

水部の渦が導水口に到達するために導水中の浮遊砂濃度が高くなるとみられたが、仕切壁を設けても遊水部が短いために効果は小さい結果となった。
④潜堤に排砂口を設けることにより遊水部での堆入射波

模型実験
数値計算

直立壁
導水管
潜堤

-100

0

砂が著しく低減する効果が見られた。
３．４．導入水に含まれる砂を最小限にするため

H=50cm，T=4.0s

図-９の数値計算結果

100 200 300 400
法先からの距離x(cm)ることは適当ではない。このような理由から、数値計算を用いてさらなる検討を加えた。用いた数値計算法は、２次元ではあるが流体現象を比較的精度良く表現出来るとされる
VOF 法をベースに浮遊砂の移流拡散方程式を加え、浮遊砂

150

導水流速
4.E-06120

3.E-06902.E-06

浮遊砂濃度

60

1.E-06300.E+000
模型実験

3く依存するので、水理模型実験結果をそのまま現地に換算す

H=50cm，T=4.0s

3

浮遊砂濃度[m /m ]1/5 とスケールが小さい。浮遊砂の現象は、スケールに大き600

模型実験と数値計算の比較(波高)5.E-06

水理模型実験は、縮尺が大きいとはいえ、現地に比べると500

導入流速[m/s]

導入水中の浮遊砂濃度の低減に効果がある。③遊80706050403020100

数値計算

上向きに改良すると導入流速が 3 割程度減少し、波高 Have(cm)の港内への浮遊砂の侵入が見られた。②導水口を濃度変化を計算出来るようにした。まず、水理模型実験と同

図-10 模型実験と数値計算の比較

じスケールで、流れ及び浮遊砂濃度の再現性を検討した。図(導水流速と浮遊砂濃度)-９に海水導入工周辺での計算と実験結果の波高変化の例を78示す。図-10 に導水流速及び浮遊砂濃度対数表示

平均浮遊砂濃度の計算と実遊水部で渦の発生

験結果を示す。導水流速及び浮遊砂濃度については、数値
天端幅 8m

計算結果が実験結果に比較して若干小さめであるが、ほ遊水部幅 7m

ぼ妥当な結果を与えていることがわかる。
図-11 に数値計算結果の一図-11 浮遊砂濃度計算の一例(改良無し)

例を示す。この図から遊水部での大きな渦が浮遊砂濃度分布に大きな影

H=4.0m，T=10.0s

浮遊砂濃度

高 4m、周期 10 秒の規則波を用いた。
模型スケールでの数値計算及び水理模型

0.0E+00

遊水部幅2.0倍

遊水部幅1.5倍

天端幅2.0倍

天端高+0.5m

2.0
天端幅1.5倍来襲するので、現地スケールでの検討では波5.0E-05

潜堤排砂口行った。なお、現地では 2.5m 以上の波浪が1.0E-04

仕切壁

いて、現地スケールでの海水導入工の検討を導水口上向き以上から、ここで開発した数値計算法を用

2.5
1.5E-04

改良無し

ることが確認されている。

3.02.0E-04
平均浮遊砂濃度 [ m 3 /m 3 ]型実験でも遊水部内に強い渦が発生してい以降、導水口上向き

1.5
1.0

導水流速 [ m/s ]響を与えていることがわかる。なお、水理模

導水流速

0.50.0

図-12 各改良の効果

実験結果から、排砂口、遊水部の長さ、潜堤天端幅が、主に導入水に含まれる浮遊砂濃度に影響を与えていると考えられたので、海水導入工の天端幅の拡導水口上向き

潜堤部排砂構造を変化させて、数値計算を行った。その結果を図-12 に示す。この図から、模型実験で検討した改良は、現地スケールの計算でも類似した傾向

導水口上向き潜堤部排砂を示し、潜堤に排砂口を設けた場合は遊水部での堆砂が低減された。遊水部幅を広くすると渦が沖遊水部幅の拡側に移動し導水口付近の浮遊砂は減少する。導入
図-13 改良案のイメージ

水中の浮遊砂を低減させるには、導水口を上向きにする、潜堤に排砂口を設ける、遊水部幅を広くする、天端幅を広くする等が挙げられる(図-13)。構造により高波浪時にも 0.17mm 程度の粒径の浮遊砂濃度は 10-5(m3/m3)以下となり、例えば 100 万 m3 の海水導入に対して侵入する砂は数 m3 程度と小さくすることができる。
３．５．港内埋没評価モデルの現地適用性開発した浮遊砂による港内埋没評価のための数値計算を実際に港内埋没に苦慮している漁港に適用し、その再現性について検討した。
本モデルが良い再現性を有することが確認された。さらに対策に関して、検討を行った。図-14（左）に示すように、実際の深浅測量の変化と、計算結果の傾向が良く合うことがわかる。さらに同漁港への対策として防波堤を検討し、効果のあることがわかった（図-14 右）。しかしながら、依然として港口部周辺には堆積領域が存在するために港口部周辺に潜堤、海水導入工などを配置したところ一定の改善が79見られた。今後、さらに検討していく予定である。
①現況

深浅図比較結果
2005.10～2005.12
赤：堆積
青：侵食

堆積量が減少航路での堆積が見られない
堆積

侵食
この領域の侵食が減少

堆積

侵食

図-14 実際の漁港へのモデルの適用及び対策案の検討
４．考察
本調査により、海水導入工を利用した港内埋没対策技術は十分効果のあることを確認した。導入水に含まれる砂を最小限に抑えるためには、遊水部の長さを通常の海水導入工の２倍程度にすること、遊水部に排砂口をもうけること、導水管を上向きにすることが必要である。現在、開発した数値計算法により、さらなる検討を行っているところである。
本工法は複雑な構造を持たず、波力など外力に対しても十分安定性を確保することができる。正確な整備コストは未検討であるが、比較的安価であると考えられる。既存の防波堤を改良しても設置可能であるので、非常に汎用性が高い。
また、本課題においては港内埋没を精度良く評価することが、海水導入工の開発に当たって不可欠であった。そのために開発した港内埋没評価モデルは再現性があることが確認された。本モデルは漁港周辺の地形変化のみならず、従来手法では評価できなかった港内の砂移動を評価できる。今後、港内埋没防止対策を検討するための基本ツールになると考えられる。
港内埋没対策の基本的な考え方は、過去の研究成果を踏まえ、今年度の全国漁港漁場整備技術研究発表会に発表する予定である。
摘要
以下のような課題が残されており、さらに検討していく予定である。
（１）海水導入工の現場への適用と効果確認及び改良（２）港内埋没評価モデルの積極的利用、問題点の抽出、改良のための現地調査（３）同モデルへの掃流砂、粒径分布の変化等の組み込み、高速化、汎用性
参考文献
例えば・中山哲嚴・牧野弘幸・新井雅之・大村智宏・小林学・田村仁・灘岡和夫・佐藤勝弘：港内埋没対策技
80

術と地形変化予測モデルの開発，海岸工学論文集，第 53 巻，pp.526-530．
・鈴木崇之・岡安章夫・柴山知也 (2004)：砕波による時空間的浮遊砂濃度変動の数値シミュレーション，
海岸工学論文集，第 51 巻，pp.381-385．
・中山哲嚴・新井雅之・大村智宏・小林学・田村仁・灘岡和夫・半沢稔 (2005)：海水導入工による港内埋没対策に関する実験的研究，海岸工学論文集，第 52 巻，pp.481-485．
・中山哲嚴・山本潤・松岡道男・池田正徳・竹内聖一・灘岡和夫・八木宏 (1994)：中小港湾港口部に発生する水平大規模渦についての現地観測，海岸工学論文集，第 41 巻，pp.321-325．
・灘岡和夫・田村仁・飯塚広泰 (2004)：中小港湾港口部に発生する水平大規模渦による浮遊砂流入特性に関する数値解析－浮遊砂流入から見たハーバーパラドックス－，海岸工学論文集，第 51 巻，pp.451-455．
・平山克也・上原功・永松宏一・平石哲也(1998)：珊瑚礁リーフにおける波と流れの計算法の適用性，

海岸工学論文集，第 45 巻，pp161-165
・平山克也・平石哲也 (2004)：ブシネスクモデルによる砕波・遡上計算法とその適用性，海岸工学論文集，第 51 巻，pp11-15
・松本朗・松岡道男・中山哲嚴・山本潤・灘岡和夫・八木宏(1993)：中小港湾埋没機構解明に向けての長期連続現地観測の試み，海岸工学論文集，第 40 巻，pp.491-495．
山本正昭・中泉昌光・間辺本文 (1987)：潜堤付防波堤による海水交流工法の開発，第 34 回海岸工学講演会論文集，pp.675-679．
・Nielsen, P. (1992):Coastal bottom boundary layers and sediment transport, World Scientific, RiverEdge, N.J.,p.324.
本研究で投稿した論文
・中山哲嚴・牧野弘幸・新井雅之・小林学・佐藤勝弘：導入水に含まれる浮遊砂量を低減するための海水導入工の構造に関する研究，海岸工学論文集，第 54 巻（投稿中）・中山哲嚴・牧野弘幸・新井雅之・大村智宏・小林学・田村仁・灘岡和夫・佐藤勝弘：港内埋没対策技
術と地形変化予測モデルの開発，海岸工学論文集，第 53 巻，pp.526-530．
・中山哲嚴・新井雅之・大村智宏・小林学・田村仁・灘岡和夫・半沢

稔 (2005)：海水導入工による港内埋没対策に関する実験的研究，海岸工学論文集，第 52 巻，pp.481-485．

81