２/３
修士論文発表会

伊勢湾における小型底びき網漁船を用いた環境計測及び資源分布予測
学籍番号 47166619
指導教員 多部田 茂 教授

岩場 公利

背景

1

伊勢湾における小型底引き網漁業の漁獲対象種の動態モデル
格子の生息適正の評価

各環境因子jの選好強度
Pj=0.00(不適)～1.00(最適)
各環境因子の選好強度の積を複合選好強度Pjとして算出
出典：本宮ら，伊勢湾におけるマアナゴの動態モデルの開発，2015水産工学
出典：鈴木ら，魚類動態モデルを利用したマアナゴ漁業シミュレータの開発，2013土木学会論文
集B2（海岸工学）

本研究室では，これまで伊勢湾におけるマアナゴ・シャコなどの小型底びき網漁業の漁獲対象種について魚類動態モデルをつくってきた

これらのモデルは水温・溶存酸素（DO）などの環境因子に基づいて動態をモデル化したものである

背景モデルの参考データ
環境因子
観測点：数十カ所程度
観測頻度：月に一度

2

動態モデル
水温やDOの選好性を決定する際，室内実験データなどを基にしている
出典：鈴木ら，魚類動態モデルを利用したマアナゴ漁業シミュレータの開発，2013
土木学会論文集B2（海岸工学）

30

しかし，小型底びき網漁船の漁場は日毎に変化しており，短期的な環境因子の変動が関わっていることが伺われる．
環境因子をより詳しく把握する為にも，実海域における頻繁かつ詳細な調査が必要である

4小型底びき網漁船に観測器をつけて直接測定すれば，より詳細な環境因子，さらに漁獲から実海域における動態のデータが得られるのでは？

π∏
選好強度

515102025π∏
選好強度0
出典：三重県水産研究所，貧酸素情報H2806

123

研究概要

3

漁船を用いて測定することで，これまでのデータよりも観測頻度や観測点を増やすことができ，貧酸素水塊の境目など環境因子の詳細を把握することが可能となる

１曳網毎に漁獲量を記録してもらうことで漁獲データを得ることができ，同時に測定した環境データと合わせて，魚種別に実海域における環境因子に対する選好性を推定できる

目的

4l 貧酸素水塊を含めた詳細な環境因子を把握する

l マアナゴ・シャコなど小型底びき網漁業の漁獲対象種について，実海域の環境因子に対する動態を把握する

l 漁獲対象種について資源分布予測を行う

手法
小型底びき網漁船の概要とデータ収集方法

5測定する船：
有滝地区の小型底びき網漁船
計３隻

手法

6

測定機器 2015年度試験的に数回行い，2016年6月以降現在も継続して測定
測定頻度：月に８回程度（3～7日に一度）
測定点：曳網中３分毎に測定（曳網は1時間程度を約8回，3隻で測定）
標本船操業野帳（手作業記録）
対象種：マアナゴ・シャコなど
記録内容：１曳網毎の魚種別漁獲量
船の曳網方向・曳網速度

２０１６年度から導入
塩分ロガー（伝導度計）
Onset社製：Hobo U-24
計測頻度：３分

水深ロガー
（圧力計）
Onset社製：Hobo U-20
計測頻度：３分

計５種類のロガー

測定場所

7有滝地区
伊勢湾における
小型底びき網漁業の主な漁港のひとつ

引用：
https://www.google.co.jp/maps/@34.7802929,136.8538846,10.25z?hl=ja&authuser=0

研究の流れ8漁船で計測して環境因子および漁獲の基データを得る
データ補正（センサの応答時間補正）
GPS・DO・水温・水深のデータをまとめて水質データセットを作成
漁獲データセットを作成し，水質データと1曳網毎に対応させる
各魚種毎にCPUE（Catch Per Unit Effort）を計算描画，グラフ化
考察データ整理9

水質データと漁獲データの対応作業
作成した水質データセットと漁獲データセットは別々に存在しており，CPUEの計算や，環境因子に対する動態把握の為にも，これらのデータを対応させる必要がある
水質データ
水質データの速度平均[kt] 水質データの平均方向

水質データLabel2

水質データの曳網時間 コメント

漁獲データ
曳網方向

kt

Tim e2.85117

2016/9/1 18:18

39130

2.8452.71324

2016/9/1 19:30

543102.8502.72140

2016/9/1 20:54

45

130

2.850

2.84307

2016/9/1 22:0954300

2.845

2.79136

2016/9/1 23:39

51

130

2.845

2.843142016/9/2 1:00

45

290

2.850

2.92992016/9/2 2:15

39130

2.850

2.88302

2016/9/2 3:2751

250

2.855

2.86304

2016/9/14 16:42

30

290

2.945

漁獲データは１曳網につき１データしか存在していない
一方，水質データは曳網中の３分毎にデータがある為，1曳網毎に区切り，平均値で比較した漁獲データに記載された・曳網方向・船速・曳網時間と比較して対応させた

結果①（三重水研の貧酸素情報と比較）

10

三重県の水産研究所が公表している
貧酸素情報と比較
l 月に一度公開されている貧酸素情報データ（2016年6月6日）と近い日に,漁船で測定を行っていた2016年6月5日のデータを水平分布化して可視化し，比較を行う
l 漁船測定データは3隻全てのデータを統合したものである（これ以降の結果も同様）
l 比較を行うことで，・測定データの妥当性（貧酸素情報と大きなずれが生じていないか）
・測定データがより詳細にデータを得ているかどうかを確認する

結果①（三重水研の貧酸素情報と比較）
６/６
三重水研の貧酸素情報（黒点の場所が観測点）

６/５
漁船を用いた実測データ（色の付いた点1つ1つが観測点）
[km]

11

結果②（環境・漁獲データの水平分布）12

マアナゴ・シャコCPUEおよびDO・水温の水平分布図
l マアナゴ・シャコのCPUEおよび，環境因子であるDO（溶存酸素量）・水温について，観測した日付毎に水平分布化を行い，地図上のどの場所で曳網したかを判別できるようにした
l データは2016年7月30日～8月11日までのものを示すl これまでの月に一度の測定では捉えにくかった，短期的な物理場変動の詳細を確認する
l 環境因子と漁獲対象種のデータを同時に可視化することで，環境因子に対する動態やエサ生物との関係性や漁業者の漁場選択について考察する

結果②（環境・漁獲データの水平分布）
７月３０日のマアナゴ・シャコ
CPUEとDO・水温の水平分布図
貧酸素の縁辺部であると考えられる低DOの場所でマア
ナゴ・シャコ共に漁獲が多い
低DOの場所は低水温となっている

13

結果②（環境・漁獲データの水平分布）
８月６日のマアナゴ・シャコ
CPUEとDO・水温の水平分布図

漁業者は日毎に曳網場所を変えている貧酸素が少し解消されている。

14

結果②（環境・漁獲データの水平分布）
８月１１日のマアナゴ・シャコ
CPUEとDO・水温の水平分布図

貧酸素が再び生じている。
このことから，水塊構造は数日スケールで大きく変動している15

結果③（水質とCPUEの散布図）

16

DO・水温と漁獲対象種CPUEの関係
l DO（溶存酸素量）を横軸，水温を縦軸にとり，魚種別にCPUEを色の濃淡で示した散布図を作成した．例としてマアナゴ・シャコ・マダイを示す
l 水温・DOの値は１曳網の間の平均値をとったものである
l 各漁獲対象種がどのような水質条件の時に多く漁獲されているかをグラフから確認
し，既往研究で得られている知見と比較を行った
l 漁獲記録が全曳網の１割以上で確認された月のデータをまとめたグラフを，漁獲対象種別に作成した

結果③（水質とCPUEの散布図）

17

室内実験で仮死状態とされている2.01mg/Lを下回る場所でも多く漁獲がみられる
既往研究で適正水温とされている23℃以下で主に漁獲されている丸で囲んだ低DO・低水温で漁獲が多い

結果③（水質とCPUEの散布図）

18

室内実験で生存閾値とされている1.62mg/Lを下回る場所でも多く漁獲がみられる
既往研究で適正水温とされている20℃以下で主に漁獲されている

丸で囲んだ
低DO・低水温で漁獲が多い

結果③（水質とCPUEの散布図）

19
室内実験で行動異常値とされている2.89mg/Lを下回る場所での曳網は確認されなかった
既往研究で適正水温とされている8～28℃の間で漁獲されていた丸で囲んだ水温が17~19℃
DOが6~8mg/L
で漁獲が多い

結果④（機械学習を用いた分布予測）

20

マアナゴの資源分布予測
l ニューラルネットワーク（NN）を用いて，主な漁獲対象種であるマアナゴのCPUEの予測を行った
l NNはRの単一中間層ニューラルネットワークパッケージ「nnet」を使用した
l NNの学習には実測データを使用し，マアナゴの漁獲がみられた2016年6月～9月のデータを利用した
l 各値は０～１となるよう正規化して用いた

結果④（機械学習を用いた分布予測）

マアナゴCPUEの予測方法
入力・DO・水温・水深・シャコCPUE（エサ環境）の４つ
出力・マアナゴCPUE
なお，入力および出力の値は０～１となるよう正規化した

月毎に学習データを与えて予測を行っている．
例えば，6月の曳網は3隻合計で158回あった．
これを機械に学習させる為のデータと，検証用のデータに振り分ける必要がある
そこで，26回を検証用データとして利用し，その他を学習用の教師データとした21

結果④（機械学習を用いた分布予測）22

中間層の数別のマアナゴCPUE予測精度変化
RMSE（平均二乗誤差）は低いほど予測精度が高い
よって，中間層は4層が最も精度が良いと判断できるため，4層を採用する

結果④（機械学習を用いた分布予測）
NNによる月別のマアナゴCPUE予測結果

予測値

予測精度が低い例：8月

23
予測値

予測精度が高い例：9月

横軸実測値から求めた
マアナゴCPUE
縦軸
NNによって予測した
マアナゴCPUERMSEの月別結果
6月：0.0939
7月：0.0866
8月：0.101
9月：0.0337

実測値

8月のNNによる予測は，CPUEが高い時に精度が悪い
理由として，学習データに高CPUEのデータが多く存在していないことが考えられる

実測値

9月は高CPUEのデータも揃っていた為，予測精度が高かったと考えられる

結果④（機械学習を用いた分布予測）

24

各月のマアナゴCPUEの予測に最も寄与している因子の検討
入力
l DO
l 水温
l 水深
l シャコCPUE（エサ環境）
４つの内，1つのデータを一定値としてNNの学習に反映させないようにする

出力l マアナゴCPUE
月毎に学習データを与えて，予測精度の変化をみることで各入力の出力に対する寄与を確認する→予測精度が悪くなった項目が最も予測に寄与している

結果④（機械学習を用いた分布予測）
月別のマアナゴの動態に寄与している因子の特定結果

マアナゴの動態に最も寄与している
因子は月によって異なる結果となった

25

6月

7月

8月

9月

結果④（機械学習を用いた分布予測）

各月のマアナゴCPUEの予測に最も寄与している因子の整理
結果
６月：DO
７月：シャコ（エサ環境）
８月：シャコ（エサ環境）
９月：水温

既往研究
田島ら(2016)：
マアナゴは7月・8月のエサ環境が悪い
黒木(2016)：
マアナゴは水温が低いほど漁獲が多かった（本研究でも低水温・低DOで漁獲が多い）
7月・8月は少ないエサを優先し，6月は低水温の場所が低DOである傾向が見られていたことから，低水温を好みながらも貧酸素に移動を阻まれ貧酸素水塊の縁辺部に蝟集していたと考えられる

26

結果④（機械学習を用いた分布予測）

27

マアナゴCPUEの水平分布予測
入力
l DO
l 水温
l 水深
学習データは6月～9月のデータを全て使用した
水平分布予測を行う日付のデータは学習データから取り除いた

出力l マアナゴCPUE
NNによって，環境因子から予測した結果が水平分布的にどの程度整合しているかを確認する
現実的には，事前にシャコの分布を把握できない為，入力からシャコCPUEを取り除いた結果④（機械学習を用いた分布予測）
NNによる
７月３０日のマアナゴCPUEの予測結果と実測データの比較
実測データでは貧酸素の縁辺部と考えられる場所で漁獲が多かった
予測結果でもCPUEの高い場所を予測できている28

結果④（機械学習を用いた分布予測）
NNによる
８月１１日のマアナゴCPUEの予測結果と実測データの比較
実測データでは全体的に漁獲量が少ない予測結果では全体的にCPUEが高いと予測しており，あまり整合していない

29

結論
三重県の貧酸素情報との比較•

実測データの妥当性を示すと共に，計測点が大幅に増えたことを示した

DO・水温とCPUEの水平分布図•

これまで捉えるのが難しかった短期的な水塊構造の変化を捉えており，環境因子は数日スケールで大きく変化していることがわかった

DO・水温とCPUEの散布図•

マアナゴ・シャコにおいて，実海域における行動特性は，既往研究の室内実験で得られている知見と異なる可能性が示唆された

NNを用いた因子特定• マアナゴの動態に最も寄与している因子は月によって異なっていることを示した
NNを用いた水平分布予測•

一部の条件以外ではマアナゴCPUEが高い場所を予測できており，NNを用いた資源分布予測の有用性を示した

30

今後の課題31

実測範囲の拡大•

有滝地区の漁船３隻に計測を依頼した，今後計測する船を増やすことができれば，より広範囲のデータを収集することができる

データ精度の向上•

漁獲データの記録は手作業で記録している．記録の自動化等を行うことで，データの記録漏れ防止や精度向上が期待できる•

漁獲量は水質データの1曳網間の平均値と対応させている．水中カメラ等により曳網中の網に入ったタイミングを判断できれば，実海域の動態についてより詳細なデータを得ることができる

分布予測精度の向上•

１曳網毎のデータを与えて予測したが，局所的なデータだけでなく，曳網場所全体が貧酸素など，全面的な環境因子の情報も予測に反映することで，より高精度な予測ができると考えられる