土木学会論文集B1(水工学)水工学論文集,第60巻,2016年2月
Vol.72, No.4, I̲925-I̲930, 2016.

プレート型振動センサを用いた掃流砂量
計測手法に関する研究
STUDY ON BEDLOAD MEASUREMENT BY PLATE-TYPE VIBRATIONSENSOR
小柴孝太1・角 哲也2・堤 大三3
Takahiro KOSHIBA, Tetsuya SUMI and Daizo TSUTSUMI
1学生会員 学士（工）京都大学防災研究所 水資源環境研究センター（〒611-0011 宇治市五ヶ庄）
2正会員

博士（工）京都大学防災研究所 水資源環境研究センター（〒611-0011 宇治市五ヶ庄）
3正会員 博士(農)京都大学防災研究所流域災害研究センター (〒506-1422高山市奥飛騨温泉郷中尾436-13)

Measuring the bedload transport rate is utmost interest to understand and quantify all sediment transport
related phenomena. Three existing bedload measuring systems are 1) the hydrophone, 2) the geophone,and 3) the plate microphone as their modified system. All systems estimate transport rate based on
sediment impact, whereas each has respective disadvantages. In this paper, as the other improved system,
plate-type vibration sensor is discussed. Calibration experiment is conducted and the system is comparedwith the other plate hydrophone. In the experiments, impact data is perceived as both vibration pressureand vibration pulse, and analyzed the same as plate microphone. Although it has been clarified that finesediments lower than 10 mm are not detectable with plate hydrophone, plate-type vibration sensor
detected vibration pulses of them clearly. Moreover, detection rate can be linked to the sediment bulk
density and the gravel jump length.
Key Words : plate-type vibration sensor, bedload measurement , flume experiment, grain size

１．はじめに
日本には約3,000箇所のダムが所在するが，堆砂問題が顕在化しており，適切な貯水池土砂管理手法の開発が急務である．その中では，上流から流下する土砂をダム下流へとバイパスさせることにより貯水池内の堆砂を防止する排砂バイパストンネルは世界の中で日本とスイスがリードする最先端技術であり，日本においては，新宮川水系旭ダム，天竜川水系美和ダム等で実施されており，同じ天竜川水系の小渋ダム，松川ダムでも導入が進められている．しかし，トンネル内部の土砂動態をどのように把握するのかは未確立であり特に，長期的な維持管理に直結するトンネルインバートの摩耗対策を進めるためには，トンネル内部の水理特性と土砂輸送機構の解明が必須である．また，ダム下流に排出された土砂がその後の洪水でどのように輸送されるかは，河道の治水および環境管理の両面で重要事項である．土
砂の輸送機構のうち，ウォッシュロードや浮遊砂については，濁度計を用いてその動態が把握されてきている一方，掃流砂の現地観測は極めて難しいのが現状である．

これまで掃流砂に関しては，例えばハイドロフォンを用いた観測が京大穂高砂防観測所などで実施されてきている 1)．ハイドロフォンは鋼製パイプ内にマイクロフォンを導入し，砂礫の衝撃を発生する音圧とパルスで計測するものである．これに対して，日本と同様な山岳河川を有するスイスにおいては，鋼製プレートに振動センサを設置したジオフォンが開発され，スイスのSolisダムの排砂バイパストンネ
ルに導入され，平成25年より観測が開始された 2)．
ハイドロフォンは，2mm程度の粒径から感知が可能であり，小粒径に対する感度が高いといえる一方，大粒径の砂礫が衝突した際にパイプが変形してしまい，正確な計測が困難になるという脆弱性を持っている．これに対して，ジオフォンは巨石の衝突に対して変形しにくいという長所があるものの，検知可能な最小粒径が1cm程度とされ，感度の低さが欠点となっている．
そこで，これらの長所と欠点を鑑みて著者らはプレートマイクロフォンを提案した．プレートマイク
ロフォンはジオフォン同様の鋼製プレートマイクロフォンを設置した装置であり，計測原理はマイクロフォンに則っている．プレートマイクロフォンは，

I_925

高感度と巨石に対する強度を併せ持つことを想定して提案され，人工水路を用いたキャリブレーション実験を行った 3)．その結果，高流速による掃流砂量の計測可能性は明らかになってきているものの，現段階ではマイクロフォンに迫る細粒土砂までの感度の高さは確認されていない．そこで，新たに掃流砂量計測手法として提案された鋼製プレートにジオフォンとは異なる原理の振動センサを設置し，このプレート型振動センサに関するキャリブレーション実験を行った．本研究では当実験結果をもとに，プレート型振動センサの基本的な特性を明らかにし，従来の計測手法との比較を行う．

図-1 実験水路に設置されたプレート型振動センサ
表-1 実験材料諸元

２．検討方針本研究ではプレート型振動センサによって計測される振動値の基礎的な性質を検証するため，プレートマイクロフォンと同条件，同解析でのキャリブレーション実験を行った．振動センサによって計測されるデータは，プレートマイクロフォンと同じく以下の3項目であり，これらの計測結果を用いて比較を行う．
(1) 生 波 形 ： 土 砂 の 衝 突 に よ り 生 じ た 振 動 （ プ
レートマイクロフォン，ハイドロフォンの場
合は音）を電圧に変換し，波形として表したもの．
(2) パルス数：電圧値に関して6段階の閾値を設定し，その値を越える値の数をカウントした数．
(3) 電圧平均値：電圧値の一定時間毎の時間平均値．
パルス数は生波形に現れたピークの数であり，衝突した礫の数に比例すると考えられる．しかし，プレートマイクロフォンを用いた実験では，衝突礫数が多くなるに伴いピーク同士が干渉しあうことによってピーク数が過小評価されることがわかった 4)．
そこで，礫が衝突した際のエネルギーをより正確に反映する値として波形の時間積分値が提案された．
これは一定間隔でサンプリングされている電圧値を単位時間（10秒）間隔に区切って平均値を求めたものである．しかし，電圧平均値においても，衝突礫が重複するにつれ衝突礫数との比例関係よりも減少することがプレートマイクロフォンに関して明らかになっており，このような特性についてプレート型振動センサについても検証する．また，比例関係からの減少を説明するパラメータとしてパルス数を用いた検出率 5)を導入し一定の関連性を確認したが，本研究においても同様にプレート型振動センサにおける検出率の有効性についても確認を行う．

図-2 実験に使用した礫（Ds=100mm）する手法であり，従来の計測手法であるジオフォンと異なる周波数特性をもつ振動センサを用いている．
本実験で使用した振動センサはKEYENCE社のGA313A（応答周波数40～8kHz）また，変換器として同社GA-223（応答周波数100～8kHz）を用いている．
よって双方を考慮した応答周波数は100～8kHzとなる．振動センサは厚さ2cm，長辺49cm，短辺36cmの鋼製プレートに内蔵されている(図-1参照)．一般
的に，短辺に平行な方向に流水があるように設置され，本実験でも流路と短辺が平行となるよう実験用水路に設置した．
(2) 実験材料
実験に用いる石材には，2mm，5mm，10mm，50mm，100mm，5種類の粒径の礫を用いた．それぞれの粒径と平均重量は表-1に示す通りである．平
均重量より計算した50gあたりの粒子個数は粒径2mmは3125個，粒径5mmは350個である．また，礫の密度は2.7g/cm3である．実験材料の礫の形は様々であり，図-2の例にあるように，球形に近いものから扁平なものまで多様である．

３．実験方法(1) プレート型振動センサ
プレート型振動センサは間接的に掃流砂量を計測(3) 実験水路
本実験はスイス連邦工科大学チューリヒ校(ETH)
水理・水文・氷河学研究所(VAW)の人工水路を用

I_926

り礫を投入した．投入時よりデータロガーを作動させ，プレート型振動センサに衝突した礫のエネルギーが電圧値としてロガーに記録した．
(4) 記録方法
全ての実験において，計測された衝突エネルギー
は生波形，パルス数，電圧平均値の３通りのパラメータである．それぞれの計測原理を以下に示す．
a) 生波形
礫の衝突による電圧は20μs毎に計測され，波形データとして記録された．図-4は記録波形データの例であり，縦軸が電圧値，横軸が時間である．

図-3 実験水路諸元

図-4 生波形データの例b) パルス数
生波形データの絶対値を用いてパルス数を求めた．
パルス数は電圧値に一定の閾値を設定し，その値を超過した電圧値の数を勘定した値である．具体的には，生波形の絶対値を2倍，4倍，16倍，64倍，256倍，1024倍の6段階に増幅しそれぞれ2Vの閾値を超える電圧値の数をパルス数として勘定した．パルス数は，礫の衝突の重複がない場合に衝突エネルギーに比例し，またそれぞれの増幅率に応じて衝突礫サイズが推定できると考えられる．すなわち，小さな礫は衝突エネルギーも小さいので，高増幅しなければパルスとして数えられないが，大粒径の礫が衝突した場合は低増幅率でもパルスとして現れる．
c) 電圧平均値
電圧平均値は，図-5に示された手順で計算された．
まず，計測された波形データは絶対値変換され，バンドパスフィルタを用いて特定周波数帯のデータのみが抽出される．その波形の振幅データのみ（包絡線）を抽出し，電圧値とした．この電圧値の時間平均値（本実験では10秒間）を電圧平均値とした．

図-5 電圧平均値の計算過程
表-2 実験条件
FlowCaseNo.
12345678910Q[m3/s]
0.100.100.100.100.100.180.180.180.180.18V[m/s]
2.5 (Low)
2.5 (Low)
2.5 (Low)
2.5 (Low)
2.5 (Low)
4.5 (High)
4.5 (High)
4.5 (High)
4.5 (High)
4.5 (High)

Sediment
Ds
W
[mm]
[g]
250
550
10505020stones
100 20stones
250
550
1050
50
20stones
100 20stones

いて行った．実験水路の諸元は図-3に示す通りである．全長8.00m，幅0.50m，高さ0.60m，底面はコンクリート，また側面はガラス製の水路でありプレー
ト型振動センサは上流から7.56mの地点に設置されている．上流のポンプから水が供給され，下流プレート型振動センサ直下には排水口及び礫の集石ネットが設置されている．勾配はなく，流速は上流側のジェットボックスにより水路断面積を変化させることにより調節する．
実験では，まず，ジェットボックスを用いて実験
条件の流速の流水を流し，ジェットボックス直下よ(5) 実験条件
本 実 験 の 実 験条 件 を 表-2 に 示 す ． 水 深は 一 定（0.08m）に保たれ，流速Ｖ及び流量Ｑは高流速（4.5m/s）と低流速（2.5m/s）2種類が用いられた．
ただ，一般的に排砂バイパストンネル内の土砂流下
流速は10m/sを越え6)，十分であるとは言えないが，流速の違いによる特性の比較を行うことができた．
Wは投入土砂量であり，粒径2mm，5mm，10mmの礫 に つ い て は 50g を 量 り 投 入 し た ． ま た ， 粒 径50mm，100mmについては各20個を，それぞれの実験において重量を計測して投入した．これらを踏まえ，それぞれのケースについて50回ずつ，合計500回の実験を行った．

４．実験結果及び考察(1)

実験結果a)

生波形データ

図-6，図-7の下段の波形は計測した衝突データの一例である．図-6はCase10（粒径100mm，高流速4.5m/s）の結果，図-7はCase6（粒径2mm，高流速速4.5m/s）の実験結果である．

I_927

6

図-6

Case10生波形（粒径100mm, 高流速4.5m/s）

10510410310210110010

Case6生波形（粒径2mm, 高流速4.5m/s）

10510410310210110010(a)
2mm
5mm
10mm
50mm
100mm

2 41664 256 1024
Amplitude [-]6図-7

また，参考としてそれぞれ上段にはプレートマイク
ロフォンの同条件での実験結果を示している．これ
より，振動センサの結果はマイクロフォンによる波形と比べてピークが鋭く明確である．また，Case6の粒径2mmの実験結果を比べると，プレートマイクロフォンでは波形が流速のノイズに隠れているのに対し，振動センサでは明確にピークを示している．

b)

図-8

2mm
5mm
10mm
50mm
100mm

2 41664 256 1024
Amplitude [-]パルス数と増幅率の関係．(a)低流速実験(2.5m/s)，(b)高流速実験(4.5m/s)
1000

(a)

100

パルス数

図-8に全実験における増幅率別のパルス数の結果を示す．(a)が低流速実験，(b)が高流速実験である．
これらより，増幅率が増大するにつれパルス数も増大することが明らかである．増幅率16倍以上では，平均パルス数の大きさの順が10mm<5mm<2mm<50mm<100mmとなっている．これより，10mm以下
では礫の投入個数がパルス数に大きく反映され，50mm以上では，礫1つあたりの衝突エネルギーの影響が卓越して反映されていると考えられる．

c)

(b)

1010101000

図-9に実験より得られた流砂量と電圧平均値の関連性を示す．(a)は低流速実験，(b)が高流速実験である．これらより粒径が大きくなるに連れ1mVから1000mVまで幅広く電圧平均値が増大している. ただし，粒径2mm及び5mmでは，ほとんど0mVに近く，検出できていない．50mm以上では，大きくそれぞれの粒径のグループにわかれた．10mmとその他の区別は付いているが，50mmと100mmでは明確な差がでておらず，波形が干渉し飽和していると予測される．これらの結果は，プレートマイクロフォンを用いた実験結果と大きな差異は見られなかった．
実験結果の分析a)

検出率

全投入礫数に対しどの程度のパルスが検知されたかを明らかにするために，検出率Rdを導入した．検出率は検知されたパルス数を投入土砂数で除した値である．図-10は検出率と増幅率の関係を表しており，(a)は低流速実験，(b)が高流速実験である．

110210310410510

Bedload tranport rate [g/sec]

流砂量と電圧平均値の関係性(2)2mm5mm10mm50mm100mm(b)

1001011002mm
5mm
10mm
50mm
100mm
101102103

104105

Bedload tranport rate [g/sec]
図-9 流砂量と電圧平均値の関係性(a)低流速(b)高流速

増幅率を大きくするに伴い検出率も単調に増加している．小粒径（2mm，5mm）については増幅率2，4倍では極めて低い検出率であるが，中高増幅率では大粒径と類似性のある明確な単調増加がみられ，衝突が検知できていると考えられる．これらの関係
性は高流速，低流速共に見られ，検知可能な最低増幅率は16倍程度である．

I_928

高流速，低流速を比べると低流速の結果の方が全体的に大きな値となっている．これは，低流速であるほど1つの礫が複数回プレートに衝突する可能性が高くなるからと考えられる．実際，単一礫を流下，衝突させると複数回の衝突が起きることが生波形から確認されている．そこで，衝突回数はプレート上に石が存在する時間に比例，つまり流速の逆数に10510410310210110010(a)10-20105104103

飽和度

検出率を定量的に算出することが可能であれば，実験結果と算出された検出率を用いてプレート上を流下した土砂量を再現することが可能である．そこ
で検出率の定量化を目的とし，これに関係するパラメータとして飽和度を設定した．飽和度は一度にプレート上を流れる土砂の密度を表し，この値と検出率は負の相関を持つと仮定できる．そこで，その密度を表す指標として飽和度RSを導入した．飽和度は以下の式で求められる．
∆t × Pn(1)T
ここに，∆t：単一礫の衝突によって生じるパルスの時間幅，Pn ：一度の実験での投入礫数，T：一度の実験での，1つ目の礫の衝突から最後の衝突までの間の時間である．図-12に各パラメータを図示する．
図-13に粒径50mm，低流速における飽和度の計算結果を示す．この結果より，増幅率1024倍では明確な検出率と飽和度の反比例性がみられるが，その他の増幅率では不明瞭であった．これは，生波形データからもわかるように，プレートマイクロフォンによる波形はピークの減衰が緩やであるのに対し，プレート型振動センサの波形はピークが鋭く，減衰も早いことが原因と考えられる．
Rs =2mm5mm
10mm
50mm
100mm-110b)

12101010Amplitude [-]3410(b)10210110

図-102mm5mm
10mm
50mm
100mm

010-110-2100123410101010
Amplitude [-]検出率と増幅率の結果(a)低流速(b)高流速

低流速（変換後）

図-12

飽和度を構成するパラメータ

高流速Rd [-]

図-11 低流速実験結果を低流速に対する高流速の比で除した値と高流速実験結果との比較比例すると仮定し，低流速実験結果を低流速に対する高流速の比1.8で除した値と高流速実験結果とを比較した（図-11参照）．これにより，増幅率16倍以上では概ね双方の値が一致していることがわかる．
しかし，2種類の流速だけでは十分な結果とは言えず，今後より多様な流速で実験を行い確認する必要がある．

Rs [-]
図-13 飽和度と検出率の関係（粒径50mm，低流速）増幅率は(a)1024倍(b)256倍(c)64倍(d)16倍

I_929

c)

衝突率

５．結語

飽和度の他に検出率を左右するパラメータとして，流下する礫の飛躍長に関係する衝突率が挙げられる．
ここでは，これらの飛躍長を考慮した衝突率を以下の式7)に則り計算した．
Lp= 251θ
Ds(2)ここに，Lp ：飛躍長，Ds ：礫の粒径，θ：シールズ数であり，シールズ数は以下の式で求めた．
θ=τbρG𝑠 gDs(3)ここに，G𝑠 ：礫の比重，g：重力加速度，τb は以下の式で計算した．
τb = ρu∗ 2(4)ここに，ρ：水の密度，u∗ ：摩擦速度でありここで，u∗ = √ghib(5)

ib ：水路勾配（本実験では水路勾配がないためエネルギー勾配を用いている）である．これらをもとに，衝突率P(Lp )をB
P(Lp ) =Lp(6)で計算した．ここで，Bはプレート長である．
図-14に増幅率256倍における検出率に対する衝突率の計算結果を示す．図よりどの粒径においても流速が大きくなると衝突率が減少していることがわかる．これは，流速が速いほど飛躍長が大きくなり，プレートに衝突する確率も低くなるので，それに伴い検出率が減少することを示している．

Rd [-]

P(Lp) [-]
図-14

増幅率256倍における検出率と衝突率の関係

本研究では新たな流砂量計測装置として提案されたプレート型振動センサの高流速下でのキャリブレーション実験を行った．電圧平均値については既存のプレートマイクロフォンによる実験結果と大きな差異は現れなかったが，生波形及びパルス数の結果より，粒径2mm程度の小粒径への適用の可能性を示すことができた．また，飽和度，衝突率と検出率との関係について一定の方向性を明らかにすることができた．今後，実験結果から掃流砂量を再現することを目指すに当り，さらなる実験を通して飽和度や衝突率から定量的に検出率を算出する手法を提案することが必要である．プレート型振動センサ及びプレートマイクロフォンは平成28年度より運用が開始される長野県小渋ダムの排砂バイパストンネル出口に設置される予定となっており，今後は，現地での土砂流下実験も併せて行うことにより，実機でのキャリブレーションを進める予定である．
謝辞：本研究では，スイス連邦工科大学チューリヒ校(ETH)水理・水文・氷河学研究所(VAW)より実験水路の提供を受けた．また，本研究はJSPS科研費26257304および国土交通省河川砂防技術開発地域課
題研究の助成を受けた．ここに記して謝意を表する．
参考文献
1) 堤大三，野中理伸，水山高久，藤田正治，宮田秀介，
市田児太朗: 掃流砂観測におけるプレート型ジオフォ
ンとパイプ型ハイドロフォンの比較，京都大学防災研究所年報, (57), pp. 385-390, 2013.
2) Hagmann, M., Albayrak, I. and Boes, R.M.: Field research:
Invert material resistance and sediment transportmeasurements, Proc. Int. Workshop on Sediment Bypass
Tunnels, VAW-Mitteilung 232 (R. Boes, ed.), ETH Zurich,Switzerland, pp. 123-135, 2015.
3) 小 柴 孝 太 ， 角 哲 也 ， 竹 門 康 弘 ， 堤 大 三 : Flume
Experiment on Bedload Measurement with a Plate
Microphone, 京 都 大 学 防 災 研 究 所 年 報 B= Disaster
Prevention Research Institute Annuals, B 58, 2016.
4) 鈴木拓郎, 水野秀明, 小山内信智, 平澤良輔, 長谷川祐治: 音圧データを用いたハイドロフォンによる掃流砂量計測手法に関する基礎的研究, 砂防学会誌, 62(5),pp.18-26, 2010.
5) 堤 大 三 ， 平 澤 良 輔 ， 水 山 高 久 ， 志 田 正 雄 ， 藤 田 正
治: 山地流域における音響法を用いた流砂量観測, 京都
大 学 防 災 研 究 所 年 報 B= Disaster Prevention ResearchInstitute Annuals, B 53(B), pp. 537-544, 2010.
6) Auel, C. and Boes, R.M.: Sediment bypass tunnel design –review and outlook, Proc. ICOLD Symposium Dams underchanging challenges (A.J. Schleiss & R.M. Boes, eds.), 79th
Annual Meeting, Lucerne. Taylor & Francis, London, UK,
pp. 403-412, 2011.

7) Auel, C.: Flow characteristics, particle motion andinvert abrasion in sediment bypass tunnels.Doctoral
dissertation, ETH Zürich, No. 22008, 2014.
（2015.9.30受付）

I_930