JPH0619476B2 - 潮流転流時期予測方法 - Google Patents
潮流転流時期予測方法Info
- Publication number
- JPH0619476B2 JPH0619476B2 JP24501885A JP24501885A JPH0619476B2 JP H0619476 B2 JPH0619476 B2 JP H0619476B2 JP 24501885 A JP24501885 A JP 24501885A JP 24501885 A JP24501885 A JP 24501885A JP H0619476 B2 JPH0619476 B2 JP H0619476B2
- Authority
- JP
- Japan
- Prior art keywords
- tidal current
- time
- commutation
- predicted
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 238000000034 method Methods 0.000 title claims description 19
- 238000012937 correction Methods 0.000 claims description 29
- 238000011156 evaluation Methods 0.000 claims description 15
- 238000012545 processing Methods 0.000 description 23
- 238000010586 diagram Methods 0.000 description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 description 4
- 238000005259 measurement Methods 0.000 description 2
- 244000017020 Ipomoea batatas Species 0.000 description 1
- 235000002678 Ipomoea batatas Nutrition 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
【発明の詳細な説明】 〔産業上の利用分野〕 本発明は、海峡又は狭水道における潮流の転流時期を予
測する潮流転流時期予測方法に関するものである。
測する潮流転流時期予測方法に関するものである。
海峡又は狭水道を通航する船舶は潮流変化、特に流れの
方向が換わる転流の前後で、自由な操船が不可能となる
ことがあり、他船への衝突、岩礁への乗場等海難事故の
危険性がある。これに対する対策として、従来は海上保
安庁から発表される潮汐表をもとに潮流の変化を予測し
たり、潮流信号として海流の方向により機械的に信号板
が反転するもので潮流の方向のみを表示していた。
方向が換わる転流の前後で、自由な操船が不可能となる
ことがあり、他船への衝突、岩礁への乗場等海難事故の
危険性がある。これに対する対策として、従来は海上保
安庁から発表される潮汐表をもとに潮流の変化を予測し
たり、潮流信号として海流の方向により機械的に信号板
が反転するもので潮流の方向のみを表示していた。
しかしながら、上記潮汐表により潮流の変化を予測する
方法では海峡又は狭水道の局所の潮流変化、特に転流時
期の予測はできないという問題があった。また、海流の
方向により機械的に信号板が反転するものでは潮流の方
向のみを表示する方法では潮流変化の予測、特に転流時
期の予測が極めて困難であという欠点があった。
方法では海峡又は狭水道の局所の潮流変化、特に転流時
期の予測はできないという問題があった。また、海流の
方向により機械的に信号板が反転するものでは潮流の方
向のみを表示する方法では潮流変化の予測、特に転流時
期の予測が極めて困難であという欠点があった。
これに対する対策として本出願人は、先に潮流予測方法
を開発し特許出願している｛特願昭６０−１９７２９９
号（特開昭６２−５８１８６号公報）｝。この潮流予測
方法は過去の潮流の観測データ（潮流データ）より、そ
れらのデータが示す近似式を作成し該近似式による予測
値と観測データからの値とから誤差を算出し、一定時間
にその誤差の量に対応したグレード評価を行って補正角
を求め、転流判定角度（適用海域により決まる）に至る
時刻（転流時期）の予測を行ない、転流時直前の予測値
予報の不安定性（後に詳述する）を除去するため、転流
予測時期の予め定められた範囲に入った時その予測を固
定するものである。
を開発し特許出願している｛特願昭６０−１９７２９９
号（特開昭６２−５８１８６号公報）｝。この潮流予測
方法は過去の潮流の観測データ（潮流データ）より、そ
れらのデータが示す近似式を作成し該近似式による予測
値と観測データからの値とから誤差を算出し、一定時間
にその誤差の量に対応したグレード評価を行って補正角
を求め、転流判定角度（適用海域により決まる）に至る
時刻（転流時期）の予測を行ない、転流時直前の予測値
予報の不安定性（後に詳述する）を除去するため、転流
予測時期の予め定められた範囲に入った時その予測を固
定するものである。
しかしながら上記潮流予測方法は、転流時直前の予測値
予報の不安定性を除去するという作用効果があるもの
の、後に詳述するように潮流の特異状況により転流の予
測時期と実際の転流時期に大きな誤差が発生することが
あり、航行船舶の適正な運航に障害を与える可能性があ
るという問題があった。
予報の不安定性を除去するという作用効果があるもの
の、後に詳述するように潮流の特異状況により転流の予
測時期と実際の転流時期に大きな誤差が発生することが
あり、航行船舶の適正な運航に障害を与える可能性があ
るという問題があった。
本発明は上述の点に鑑みてなされたもので、転流時直前
の予測値予報の不安定性を除去すると共に、潮流の特異
状況により予測転流時期と実際の転流時期との間に大き
な誤差の発生することのない潮流転流時期予測方法を提
供することにある。
の予測値予報の不安定性を除去すると共に、潮流の特異
状況により予測転流時期と実際の転流時期との間に大き
な誤差の発生することのない潮流転流時期予測方法を提
供することにある。
上記問題点を解決するため本発明は、潮流を観測しその
観測結果より潮流データを作成し、一定時間該潮流デー
タを記憶し、該記憶された過去の潮流デーにより流速と
角度換算した流向を示す近似式を作成し、該近似式によ
り予測される観測データ取得時の流向値の予測値と観測
データによる流向角の地との誤差を算出し、一定時間毎
にグレード評価を行い、該グレードにより補正角を算定
し、該補正角により転流時期を算出予測すると共に、潮
流の流速が零となる時期付近では、流速の予測値と観測
値の比較を行い、両者の差が所定の範囲外の場合流速の
観測値から流速と時刻の近似式を新たに作成し、両者の
差が所定の範囲内になるまで該近似式の更新を行い、前
記グイレード評価による補正角の継続更新を、前記更新
された近似式による転流時期の算出予測に変更するよう
に構成した。
観測結果より潮流データを作成し、一定時間該潮流デー
タを記憶し、該記憶された過去の潮流デーにより流速と
角度換算した流向を示す近似式を作成し、該近似式によ
り予測される観測データ取得時の流向値の予測値と観測
データによる流向角の地との誤差を算出し、一定時間毎
にグレード評価を行い、該グレードにより補正角を算定
し、該補正角により転流時期を算出予測すると共に、潮
流の流速が零となる時期付近では、流速の予測値と観測
値の比較を行い、両者の差が所定の範囲外の場合流速の
観測値から流速と時刻の近似式を新たに作成し、両者の
差が所定の範囲内になるまで該近似式の更新を行い、前
記グイレード評価による補正角の継続更新を、前記更新
された近似式による転流時期の算出予測に変更するよう
に構成した。
本発明によれば、近似式により流向角の予測される観測
データ取得時の予測値と該取得時に観測データより得ら
れた流向角との誤差を算出し、該誤差を一定時間毎にグ
レード評価を行い補正角を算定し、該補正角により転流
時期を算出予測すると共に、潮流の流速が零となる転流
時期付近では流速の予測値と観測値の比較を行い、両者
の差が所定の範囲外の場合には所定の範囲になるまで流
速の観測値から流速と時刻の近似式を新たに作成し、グ
レード評価による補正角の継続更新を、更新された近似
式による転流時期の算出予測に変更するので、転流時期
の予測値を実際の転流時期に略一致させることが可能と
なる。
データ取得時の予測値と該取得時に観測データより得ら
れた流向角との誤差を算出し、該誤差を一定時間毎にグ
レード評価を行い補正角を算定し、該補正角により転流
時期を算出予測すると共に、潮流の流速が零となる転流
時期付近では流速の予測値と観測値の比較を行い、両者
の差が所定の範囲外の場合には所定の範囲になるまで流
速の観測値から流速と時刻の近似式を新たに作成し、グ
レード評価による補正角の継続更新を、更新された近似
式による転流時期の算出予測に変更するので、転流時期
の予測値を実際の転流時期に略一致させることが可能と
なる。
以下、本発明の一実施例を図面に基づいて説明する。
第２図は本発明の潮流転流時期予測方法に用いるデータ
処理装置のシステム構成を示すブロック図である。同図
に示すように、本データ処理装置は、海峡又は狭水道１
の所定箇所の改定に沈設した潮流センサ２からの潮流の
流速及び流向の観測信号をデータ処理装置３で処理し、
信号装置４に送り、該信号処理装置４において、前記海
峡又は狭水道を航行中の船舶から観測できる表示盤５に
潮流の流向、流速及び転流時刻等の潮流予測データを表
示したり、或いは音声再生装置６及びスピーカ７を通し
て音声で航行中の船舶に潮流予測データを知らせる。
処理装置のシステム構成を示すブロック図である。同図
に示すように、本データ処理装置は、海峡又は狭水道１
の所定箇所の改定に沈設した潮流センサ２からの潮流の
流速及び流向の観測信号をデータ処理装置３で処理し、
信号装置４に送り、該信号処理装置４において、前記海
峡又は狭水道を航行中の船舶から観測できる表示盤５に
潮流の流向、流速及び転流時刻等の潮流予測データを表
示したり、或いは音声再生装置６及びスピーカ７を通し
て音声で航行中の船舶に潮流予測データを知らせる。
第３図はデータ処理装置３の構成を示すブロック図であ
り、図示するように、データ処理装置３は前記潮流セン
サ２からの潮流の流速及び流向きの観測信号を中央処理
装置（ＣＰＵ）２２に入力したり、該中央処理装置２２
からの信号を信号処理装置４に出力する入出力（Ｉ／
Ｏ）回路２１、記憶装置２３、印字装置２４、操作盤２
５、ＣＲＴ２６等を具備し、前記中央処理装置で潮流セ
ンサ２からの流速及び流向の観測信号を収集し、前記記
憶装置２３に該観測データを格納し、これら観測データ
に基づき転流時刻等の潮流予測値を算出し、表示盤５及
びスピーカ７を介して航行中の船舶に予報を発する。海
峡又は狭水道１の潮の流れは、第２図に示すように、矢
印Ａ或いはＢに示すように西向きになったり東向きにな
ったり、更には矢印Ｃ或いはＤに示すように渦になった
りする。
り、図示するように、データ処理装置３は前記潮流セン
サ２からの潮流の流速及び流向きの観測信号を中央処理
装置（ＣＰＵ）２２に入力したり、該中央処理装置２２
からの信号を信号処理装置４に出力する入出力（Ｉ／
Ｏ）回路２１、記憶装置２３、印字装置２４、操作盤２
５、ＣＲＴ２６等を具備し、前記中央処理装置で潮流セ
ンサ２からの流速及び流向の観測信号を収集し、前記記
憶装置２３に該観測データを格納し、これら観測データ
に基づき転流時刻等の潮流予測値を算出し、表示盤５及
びスピーカ７を介して航行中の船舶に予報を発する。海
峡又は狭水道１の潮の流れは、第２図に示すように、矢
印Ａ或いはＢに示すように西向きになったり東向きにな
ったり、更には矢印Ｃ或いはＤに示すように渦になった
りする。
第４図は前記海峡又は狭水道１の潮流の流向及び流速を
示す図である。図示するように、潮の流れは時間ととも
に変化し流速が零となる時刻ｔ1,t2,t3, ……の点で潮
流の方向が変わる、即ち転流が起こる。この潮流の流速
や流向は、下式で近似できる。
示す図である。図示するように、潮の流れは時間ととも
に変化し流速が零となる時刻ｔ1,t2,t3, ……の点で潮
流の方向が変わる、即ち転流が起こる。この潮流の流速
や流向は、下式で近似できる。
Y(t)=Σ{aisin(cit+k)bicos(di+k)}・・・・(1) ここで、流速はＹ（ｔ）の絶対値、流向は第４図の位相
角度である。この位相角度をθとすると、 θ=sin-1{Y(t)／MaxY(t)}±2nπ =sin-1{Σ{aisin(cit+k)+bicos(di+k)}／ Max{Σaisin(cit+k)+Σbicos(di+k)}〕 ±2nπ ここでΣｂｉｃｏｓ（ｄｉ＋ｋ）は時刻と無関係で、周
期波のバイアス値（交流に加えた直流分）に相当する。
多くの場合この値は零に近く、この場合θは、 θ=sin-1〔Σaisin(cit+k)／MaxΣaisin(cit+k)〕±2n
π また、Y(t)は大略同様の周期で繰り返されることが多い
ので、ａｉ，ｃｉ共にそれぞれａ，ｃで代表されるとす
ると θ=sin-1{sin(ct+k)／1} ±2nπ ＝ct+k±2nπ となる。
角度である。この位相角度をθとすると、 θ=sin-1{Y(t)／MaxY(t)}±2nπ =sin-1{Σ{aisin(cit+k)+bicos(di+k)}／ Max{Σaisin(cit+k)+Σbicos(di+k)}〕 ±2nπ ここでΣｂｉｃｏｓ（ｄｉ＋ｋ）は時刻と無関係で、周
期波のバイアス値（交流に加えた直流分）に相当する。
多くの場合この値は零に近く、この場合θは、 θ=sin-1〔Σaisin(cit+k)／MaxΣaisin(cit+k)〕±2n
π また、Y(t)は大略同様の周期で繰り返されることが多い
ので、ａｉ，ｃｉ共にそれぞれａ，ｃで代表されるとす
ると θ=sin-1{sin(ct+k)／1} ±2nπ ＝ct+k±2nπ となる。
但し、ｋ，ｃｉ，ｄｉは天体の運行で定まる固定定数、
ａｉ，ｂｉは過去の観測データより決定される係数、ｔ
は時間である。ま、ｋ，ｃｉ，ｄｉ，ａｉ，ｂｉは時間
と共に変化するので、ｉ＝１，２，３，……n回目（周
期）とし、データと区別している。潮流の流速や流向が
上記近似式で表すことができるとすれば、過去の一定期
間(i＝１，２，３，……）の観測データから上記係数ａ
ｉ，ｂｉを求め、潮流センサ２からの観測データにより
潮流の流れが変わる時間、所謂転流時刻が予測できる。
ａｉ，ｂｉは過去の観測データより決定される係数、ｔ
は時間である。ま、ｋ，ｃｉ，ｄｉ，ａｉ，ｂｉは時間
と共に変化するので、ｉ＝１，２，３，……n回目（周
期）とし、データと区別している。潮流の流速や流向が
上記近似式で表すことができるとすれば、過去の一定期
間(i＝１，２，３，……）の観測データから上記係数ａ
ｉ，ｂｉを求め、潮流センサ２からの観測データにより
潮流の流れが変わる時間、所謂転流時刻が予測できる。
中央処理装置２２は過去の観測データから上記係数
ａｉ，ｂｉを算出し、所定時間毎の観測データをもとに
観測時を基準に潮流の流速、流向を予測する。その際、
流向については前記所定の観測時間毎に後述するグレー
ド評価方法により、予測値の補正を行なっている。
ａｉ，ｂｉを算出し、所定時間毎の観測データをもとに
観測時を基準に潮流の流速、流向を予測する。その際、
流向については前記所定の観測時間毎に後述するグレー
ド評価方法により、予測値の補正を行なっている。
第６図は転流予測のための処理の流れを示す図である。
潮流センサ２からの観測データより流向角度θを求める
（ステップ２０１）。次に上記（１）式に基づいて、予
測値の流向角θ′を求め（ステップ２０２）。これら予
測値の流向角θ′と観測データよりの流向角θとを比較
し両者の角度差、即ち角度誤差θｅ（θ′−θ＝θｅ）
を求める（ステップ２０３）。該角度誤差θｅを後述す
るグレード評価により、補正角度Δθを求め（ステップ
２０４）、前記予測値の流向角θ′を補正した補正角
θ″＝θ′＋Δθを求め（ステップ２０５）、前に予測
した転流予測時刻の変更を行なう（ステップ２０６）。
所定時間（例えば１５分）が経過したら（ステップ２０
７）、前記ステップ２０１に戻り上記処理を繰り返す、
即ち予測のための処理を所定時間間隔で定期的に行ない
転流時刻の予測精度を上げる。
潮流センサ２からの観測データより流向角度θを求める
（ステップ２０１）。次に上記（１）式に基づいて、予
測値の流向角θ′を求め（ステップ２０２）。これら予
測値の流向角θ′と観測データよりの流向角θとを比較
し両者の角度差、即ち角度誤差θｅ（θ′−θ＝θｅ）
を求める（ステップ２０３）。該角度誤差θｅを後述す
るグレード評価により、補正角度Δθを求め（ステップ
２０４）、前記予測値の流向角θ′を補正した補正角
θ″＝θ′＋Δθを求め（ステップ２０５）、前に予測
した転流予測時刻の変更を行なう（ステップ２０６）。
所定時間（例えば１５分）が経過したら（ステップ２０
７）、前記ステップ２０１に戻り上記処理を繰り返す、
即ち予測のための処理を所定時間間隔で定期的に行ない
転流時刻の予測精度を上げる。
第７図はグレード評価と補正角の関係を示す図である。
図示するようにグレードは角度誤差θｅの大きさにより
Ａ，Ｂ，Ｃ，Ｄ，Ｅの５段階に分け、それぞれグレード
Ａ〜Ｅに応じて補正角Δθを与えている。即ち、第６図
の処理フローのステップ２０４において、記憶装置４に
格納されている第７図に示すテーブルを参照して補正角
Δθを求め、ステップ２０５において、予測補正、即ち
該補正角Δθにより、前記予測値の流向角度θ′を補正
した補正角θ″＝θ′＋Δθを求める。
図示するようにグレードは角度誤差θｅの大きさにより
Ａ，Ｂ，Ｃ，Ｄ，Ｅの５段階に分け、それぞれグレード
Ａ〜Ｅに応じて補正角Δθを与えている。即ち、第６図
の処理フローのステップ２０４において、記憶装置４に
格納されている第７図に示すテーブルを参照して補正角
Δθを求め、ステップ２０５において、予測補正、即ち
該補正角Δθにより、前記予測値の流向角度θ′を補正
した補正角θ″＝θ′＋Δθを求める。
上記の如く、過去の一定期間の観測データにより、これ
ら過去の観測データの示す近似式を求め、観測により得
られた観測データからの流向角と、前記近似式で求めら
れた予測値の流向角の角度誤差を求め、該角度誤差をグ
レード評価し転流時刻を予測する方法は、潮流の流向が
第４図に示すように時間と共に、正確に変化すれば予測
も正確であるが、実際は第５図（ａ）及び（ｂ）に示す
ように潮流の特異状況により予測転流時期が実際の観測
転流時期より速くなる場合と遅くなる場合がある。第５
図（ａ）及び（ｂ）は第４図のＡ部分の拡大図であり、
同図（ａ）は予測転流時期が実際より速い場合を示し、
Ｔａは従来のグレード評価時で転流時刻を予測する時
期、Ｔｂはそれにより予測される予測転流時期、Ｔｃは
実際に観測した観測転流時期を示す。同図（ｂ）は予測
転流時期が実際より遅い場合を示し、Ｔｄは従来のグレ
ード評価法で転流時刻を予測する時期、Ｔｅはそれによ
り予測される予測転流時期、Ｔｆは実際に観測した観測
転流時期を示す。
ら過去の観測データの示す近似式を求め、観測により得
られた観測データからの流向角と、前記近似式で求めら
れた予測値の流向角の角度誤差を求め、該角度誤差をグ
レード評価し転流時刻を予測する方法は、潮流の流向が
第４図に示すように時間と共に、正確に変化すれば予測
も正確であるが、実際は第５図（ａ）及び（ｂ）に示す
ように潮流の特異状況により予測転流時期が実際の観測
転流時期より速くなる場合と遅くなる場合がある。第５
図（ａ）及び（ｂ）は第４図のＡ部分の拡大図であり、
同図（ａ）は予測転流時期が実際より速い場合を示し、
Ｔａは従来のグレード評価時で転流時刻を予測する時
期、Ｔｂはそれにより予測される予測転流時期、Ｔｃは
実際に観測した観測転流時期を示す。同図（ｂ）は予測
転流時期が実際より遅い場合を示し、Ｔｄは従来のグレ
ード評価法で転流時刻を予測する時期、Ｔｅはそれによ
り予測される予測転流時期、Ｔｆは実際に観測した観測
転流時期を示す。
そこで本実施例では更に第５図（ａ）における予測電流
時期Ｔｂを実際に観測した観測転流時期Ｔｃに間は同図
（ｂ）における予測転流時期Ｔｅを実際に観測した観測
転流時期Ｔｆに近付けるめ、処理を行なう。
時期Ｔｂを実際に観測した観測転流時期Ｔｃに間は同図
（ｂ）における予測転流時期Ｔｅを実際に観測した観測
転流時期Ｔｆに近付けるめ、処理を行なう。
第１図は本発明に係る潮流転流時期予測方法の処理の流
れを示すフローチャートである。同図において、ステッ
プ１０１のグレード評価補正角は第６図のステップ２０
４のグレード評価により補正角Δθを求めるのに相当
し、ステップ１０２の予測補正は第６図のステップ２０
５の予測補正、即ち補正角Δθにより予測値の流向角
θ′を補正した補正角θ″＝θ′＋Δθを求める予測補
正に相当する。
れを示すフローチャートである。同図において、ステッ
プ１０１のグレード評価補正角は第６図のステップ２０
４のグレード評価により補正角Δθを求めるのに相当
し、ステップ１０２の予測補正は第６図のステップ２０
５の予測補正、即ち補正角Δθにより予測値の流向角
θ′を補正した補正角θ″＝θ′＋Δθを求める予測補
正に相当する。
前記ステップ１０１において、グレード評価により補正
角Δθが得られたら転流予測時期の補正を実施する（ス
テップ１０２）。即ち上記補正角θ″＝θ′＋Δθを求
める。この場合、転流時期、即ち潮流の流速が零となる
時期付近までの各時刻と潮流速度の対応を予測する。次
にステップ１０３で潮流流速の前記各時刻で予測した予
測値ど実際にその時刻で観測した観測値とを比較する予
測・観測比較処理を行なう。次にステップ１０４におい
て、この予測値と観測値の一致性を確認し、一致してい
ればステップ１０５において予測値を表示する状態を継
続する。ステップ１０４で予測値と観測値が一致しない
場合はステップ１０６において、潮流速度の観測値のサ
ンプルを収集し、その時点以前の観測データの規定量を
含めて近似式の作成処理を実施する（ステップ１０
７）。この場合代表的な近似式の例としては、例えば ｙ＝ｋｘn である。ここで ｙ＝潮流速度変化 ｘ＝時刻 ｋ＝定数 ｎ＝定数 である。次にステップ１０７で得られた近似式を用いて
得られた各時刻における潮流の流速の予測値とその時刻
で実際に観測した潮流の流速の観測値との比較を行ない
（ステップ１０８）、その一致性を判定し（ステップ１
０９）、一致していれば前記ステップ１０５の表示に移
行し、不一致であれば前記ステップ１０６の潮流速度の
観測値のサンプル収集に戻り、前記ステップ１０６から
ステップ１０９までの処理を実行する。
角Δθが得られたら転流予測時期の補正を実施する（ス
テップ１０２）。即ち上記補正角θ″＝θ′＋Δθを求
める。この場合、転流時期、即ち潮流の流速が零となる
時期付近までの各時刻と潮流速度の対応を予測する。次
にステップ１０３で潮流流速の前記各時刻で予測した予
測値ど実際にその時刻で観測した観測値とを比較する予
測・観測比較処理を行なう。次にステップ１０４におい
て、この予測値と観測値の一致性を確認し、一致してい
ればステップ１０５において予測値を表示する状態を継
続する。ステップ１０４で予測値と観測値が一致しない
場合はステップ１０６において、潮流速度の観測値のサ
ンプルを収集し、その時点以前の観測データの規定量を
含めて近似式の作成処理を実施する（ステップ１０
７）。この場合代表的な近似式の例としては、例えば ｙ＝ｋｘn である。ここで ｙ＝潮流速度変化 ｘ＝時刻 ｋ＝定数 ｎ＝定数 である。次にステップ１０７で得られた近似式を用いて
得られた各時刻における潮流の流速の予測値とその時刻
で実際に観測した潮流の流速の観測値との比較を行ない
（ステップ１０８）、その一致性を判定し（ステップ１
０９）、一致していれば前記ステップ１０５の表示に移
行し、不一致であれば前記ステップ１０６の潮流速度の
観測値のサンプル収集に戻り、前記ステップ１０６から
ステップ１０９までの処理を実行する。
上記処理の具体的例を以下に説明する。
第５図のＴａにおける流速をｙａ、その次の観測におけ
る流速をｙａ１とすると、この場合の流速の差Δｙａ１
は、 Δｙａ１＝ｙａ−ｙａ１ 流速観測の時間間隔を単位時間と見做すと、ｘ＝１で、 Δｙａ１＝ｋ を得る。次の観測値として流速ｙａ２を得た場合、ｙａ
との差Δｙａ２は、 Δｙａ２＝Δｙａ１２2 （∵ｘ＝２） となる。つまり logΔｙａ２＝logΔｙａ１＋ｎlog２ となり、ここでΔｙａ２，Δｙａ１は既知であることか
ら、 ｎ＝（logΔｙａ２−logΔｙａ１）／log２ としてｎが求められる。つまい上記ｙａ，ｙａ１，ｙａ
２の観測により流速変化の近似式が得られる。これで第
３回目の実測値ｙａ３に対する予測が行えることになる
が、予測の結果、実測値ｙａ３との誤差が運用上問題と
ならない範囲として認められる場合は、この近似式を第
４回目の予測に用いることとし、若しこの誤差が大きい
と判断される場合は、第２回目と第３回目の実測を用い
上記と同様の手順で近似式を更新する。この操作を繰返
し予測値を実際の転流時期に近づける。
る流速をｙａ１とすると、この場合の流速の差Δｙａ１
は、 Δｙａ１＝ｙａ−ｙａ１ 流速観測の時間間隔を単位時間と見做すと、ｘ＝１で、 Δｙａ１＝ｋ を得る。次の観測値として流速ｙａ２を得た場合、ｙａ
との差Δｙａ２は、 Δｙａ２＝Δｙａ１２2 （∵ｘ＝２） となる。つまり logΔｙａ２＝logΔｙａ１＋ｎlog２ となり、ここでΔｙａ２，Δｙａ１は既知であることか
ら、 ｎ＝（logΔｙａ２−logΔｙａ１）／log２ としてｎが求められる。つまい上記ｙａ，ｙａ１，ｙａ
２の観測により流速変化の近似式が得られる。これで第
３回目の実測値ｙａ３に対する予測が行えることになる
が、予測の結果、実測値ｙａ３との誤差が運用上問題と
ならない範囲として認められる場合は、この近似式を第
４回目の予測に用いることとし、若しこの誤差が大きい
と判断される場合は、第２回目と第３回目の実測を用い
上記と同様の手順で近似式を更新する。この操作を繰返
し予測値を実際の転流時期に近づける。
以上説明したように本発明に係る潮流転流時期予測方法
は、最終転流時期付近までは近似式により予測される観
測データ取得時の流向角の予測値と観測データからの値
との誤差を算出し、該誤差の一定時間毎にグレード評価
を行ない、該グレードにより補正角を算定し、該補正角
により転流時期の算出予測を行ない、潮流の流速が零と
なる転流時期付近では、流速の予測値と観測値の比較を
行ない、両者の差が所定の範囲内になるまで近似式の更
新を継続し、該更新された近似式で転流時期の予測を行
なうから、グレート評価による補正角の継続更新のみに
よる転流時期の算出予測を行なう場合より、潮流転流時
期が正確に予測でき海峡又は狭水道を航行する船舶の安
全な航行に極めて優れた効果を発揮する。
は、最終転流時期付近までは近似式により予測される観
測データ取得時の流向角の予測値と観測データからの値
との誤差を算出し、該誤差の一定時間毎にグレード評価
を行ない、該グレードにより補正角を算定し、該補正角
により転流時期の算出予測を行ない、潮流の流速が零と
なる転流時期付近では、流速の予測値と観測値の比較を
行ない、両者の差が所定の範囲内になるまで近似式の更
新を継続し、該更新された近似式で転流時期の予測を行
なうから、グレート評価による補正角の継続更新のみに
よる転流時期の算出予測を行なう場合より、潮流転流時
期が正確に予測でき海峡又は狭水道を航行する船舶の安
全な航行に極めて優れた効果を発揮する。
第１図は本発明に係る潮流転流時期予測方法の処理の流
れを示すフローチャート、第２図は潮流予測装置のシス
テム構成を示すブロック図、第３図は第２図のデータ処
理装置のシステム構成を示すブロック図、第４図は海峡
又は狭水道の潮の流方向及び流速を示す図、第５図
（ａ），（ｂ）は第４図のＡ部分（転流時期近傍）の拡
大図、第６図は転流予測のための処理の流れを示す図、
第７図はグレード評価と補正角の関係を示す図である。 図中、１……海峡又は狭水道、２……潮流センサ、３…
…データ処理装置、４……信号装置、５……表示盤、６
……音声再生装置、７……スピーカ、２１……入出力回
路、２２……中央処理装置、２３……記憶装置、２４…
…印字装置、２５……操作盤、２６……ＣＲＴ。
れを示すフローチャート、第２図は潮流予測装置のシス
テム構成を示すブロック図、第３図は第２図のデータ処
理装置のシステム構成を示すブロック図、第４図は海峡
又は狭水道の潮の流方向及び流速を示す図、第５図
（ａ），（ｂ）は第４図のＡ部分（転流時期近傍）の拡
大図、第６図は転流予測のための処理の流れを示す図、
第７図はグレード評価と補正角の関係を示す図である。 図中、１……海峡又は狭水道、２……潮流センサ、３…
…データ処理装置、４……信号装置、５……表示盤、６
……音声再生装置、７……スピーカ、２１……入出力回
路、２２……中央処理装置、２３……記憶装置、２４…
…印字装置、２５……操作盤、２６……ＣＲＴ。
───────────────────────────────────────────────────── フロントページの続き (72)発明者 茨田 和生 東京都港区虎ノ門１丁目７番12号 沖電気 工業株式会社内 (56)参考文献 特公 平４−31358（ＪＰ，Ｂ２)
Claims (1)
- 【請求項１】潮流を観測しその観測結果より潮流データ
を作成し、一定時間該潮流データを記憶し、該記憶され
た過去の潮流データにより流速と角度換算した流向を示
す近似式を作成し、該近似式により予測される観測デー
タ取得時の流向角の予測値と観測データによる流向角の
値との誤差を算出し、一定時間毎にグレード評価を行
い、該グレードにより補正角を算定し、該補正角により
転流時期を算出予測すると共に、潮流の流速が零となる
時期付近では、流速の予測値と観測値の比較を行い、両
者の差が所定の範囲外の場合流速の観測値から流速と時
刻の近似式を新たに作成し、両者の差が所定の範囲内に
なるまで該近似式の更新を行い、前記グレード評価によ
る補正角の継続更新を、前記更新された近似式による転
流時期の算出予測に変更することを特徴とする潮流転流
時期予測方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP24501885A JPH0619476B2 (ja) | 1985-10-31 | 1985-10-31 | 潮流転流時期予測方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP24501885A JPH0619476B2 (ja) | 1985-10-31 | 1985-10-31 | 潮流転流時期予測方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS62105076A JPS62105076A (ja) | 1987-05-15 |
| JPH0619476B2 true JPH0619476B2 (ja) | 1994-03-16 |
Family
ID=17127349
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP24501885A Expired - Lifetime JPH0619476B2 (ja) | 1985-10-31 | 1985-10-31 | 潮流転流時期予測方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0619476B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2700871B1 (fr) * | 1993-01-28 | 1995-04-07 | Roger Jean Michel | Feux de signalisation pour chaussée maritime submersible. |
| JP2016206067A (ja) * | 2015-04-24 | 2016-12-08 | 株式会社東芝 | センサー監視システム、災害監視システム、及び災害監視方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60197299A (ja) * | 1984-03-16 | 1985-10-05 | Sanyo Chem Ind Ltd | 脱水用助剤および脱水法 |
-
1985
- 1985-10-31 JP JP24501885A patent/JPH0619476B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPS62105076A (ja) | 1987-05-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6558760B2 (ja) | 海象推定装置及び海象推定方法 | |
| CA1093662A (en) | Shipboard reference for an aircraft navigation system | |
| KR20200024094A (ko) | Gm 계산 시스템, 방법 및 프로그램, 및 횡파 주기 예측 시스템, 방법 및 프로그램 | |
| JPH0619476B2 (ja) | 潮流転流時期予測方法 | |
| JP3935828B2 (ja) | 航海支援装置及び航海支援システム | |
| JP2786494B2 (ja) | 電波航法装置 | |
| CN120949265B (zh) | 一种船舶动态高精度定位与异常行为识别方法 | |
| Lu et al. | Ocean vehicle inertial navigation method based on dynamic constraints | |
| Trump | Estimating absolute current velocities by merging shipboard Doppler current profiler data with LORAN-C data | |
| JPS5956108A (ja) | 航法装置 | |
| JPS59125014A (ja) | 複合船首方位測定装置 | |
| Watson et al. | Ultra-short baseline acoustic tracking system | |
| Severance Jr | Optimum filtering and smoothing of buoy wave data | |
| JPS5653428A (en) | Water temperature measuring system onboard ship | |
| JPH0431358B2 (ja) | ||
| Touya et al. | Automatic on‐line measurement of ship's attitude by use of servo‐type inclinometers | |
| JPS5956107A (ja) | 航法装置 | |
| JPS60133316A (ja) | 航法装置 | |
| JP2000088585A (ja) | 距離計算装置、位置計算装置 | |
| JPS61247917A (ja) | 自動航海装置 | |
| RU2057680C1 (ru) | Измеритель фактической остойчивости судна | |
| JP2736325B2 (ja) | 微小定数積算追尾測定方法 | |
| JPS5937415A (ja) | 船舶の横すべり速度解析器 | |
| Franke | Sensitivity of the error in multivariate statistical interpolation to parameter values | |
| Wilson et al. | The acoustic Doppler current profiling system at AOML |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EXPY | Cancellation because of completion of term |