JP3131661B2

JP3131661B2 - 水質汚染監視装置及びその方法

Info

Publication number: JP3131661B2
Application number: JP20311492A
Authority: JP
Inventors: 直樹原; 幹雄依田; 文智呉; 昭二渡辺; 捷夫矢萩
Original assignee: Hitachi Ltd
Current assignee: Hitachi Ltd
Priority date: 1992-07-07
Filing date: 1992-07-07
Publication date: 2001-02-05
Anticipated expiration: 2016-02-05
Also published as: JPH0627014A

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、湖沼水、海水、河川
水、ダム水などの水圏の汚染を監視する水質汚染監視装
置及びその方法に関する。

【０００２】

【従来の技術】湖沼や内湾などの閉鎖性水域に窒素、り
ん等の栄養塩類が流入すると、これらは水域に蓄積さ
れ、藻類やその他の水性生物が増殖し、富栄養化に至
る。富栄養化は、都市下水や食品工場排水などが流れ込
むとそれだけ早く進行する。湖沼、内湾などの閉鎖性水
域の状態は、水温、容存酸素濃度、ＰＨ、透明度、ＣＯ
Ｄ等の水質分析値やプランクトンをはじめとする生物の
生態等で表すことができ、特に、水域に生息するプラン
クトン等の微生物は水質の変化を最も正確に反映する。
たとえば、湖ではプランクトンの量は富栄養化の一つの
指標であり、植物性プランクトンの容積が３〜５ｃｍ/
ｍ³になるかどうかが貧栄養湖と富栄養湖とのおよその
境目とされている。従って、汚染監視対象となる水域で
は、微生物の種類やその出現量を定量的に計測し、水圏
監視に反映させる必要がある。閉鎖性水域の監視方法と
しては、水域内に設置した水質、気象計測器やプランク
トンの顕微鏡観察が実施されている。プランクトンの出
現種やその量の計測は、顕微鏡観察に依存しており、赤
潮などプランクトンの異常繁殖の監視対策は後手になら
ざるをえない。顕微鏡観察の場合、一般には定期的に船
上からサンプリングした特定個所の水を持ち帰り、顕微
鏡で目視観察しており、人手と労力を要する。また、種
類を見分ける専門家が必要なため、常時監視が困難であ
る。一方、微生物の常時監視のため、水中カメラと画像
処理装置を組み合わせた装置（特願平０２−１８２６３
０)で微生物を定量化する方式が考案されている。

【０００３】

【発明が解決しようとする課題】通常、湖沼や河川のプ
ランクトン数は高々数万個/ｍｌと希薄である。プラン
クトンを形状から分類可能な倍率（対物レンズ４倍以
上）で撮影すると、撮影画面内にプランクトンの存在す
る確率は非常に小さい。水中カメラと画像処理装置を組
み合わせた装置（特願平０２−１８２６３０)は、微生
物を出来るだけ自然な状態で撮影する方法で、微生物濃
度が高い液体中で使用可能である。即ち、培養槽や富栄
養価化が既に進行した水域など微生物濃度が高い状態で
は有効であるが、貧又は中栄養湖など現状汚染がさほど
進んでおらず、微生物濃度が希薄な水域の常時監視には
適さない。また、従来の装置では、湖沼や河川の現状汚
染を監視できるが、現状汚染が将来どのように変動する
かを予知し、この予知変動に対する浄化対策を推定する
ことに配慮がなされていない。本発明の目的は、微生物
濃度が希薄な水域内の微生物状態をオンライン計測し、
汚染状況を連続監視し、水質汚染状況を的確に把握する
と共に、水質汚染の変動を予知し、浄化対策を推定する
に好適な水質汚染監視装置及びその方法を提供すること
にある。

【０００４】

【課題を解決するための手段】上記目的は、湖沼水、海
水、河川水、ダム水等の水質を自動的にサンプルする手
段と、サンプル中の浮遊物質を濃縮する手段と、濃縮液
中の濁質を撮像する手段と、撮像画像を画像処理し、濁
質を抽出する手段と、該抽出濁質の特徴を計算する手段
と、該特徴から浮遊物質を複数の形状に分類する手段を
具備し、浮遊物質の形状分類毎に出現量を計測すること
によって、達成される。また、湖沼水、海水、河川水、
ダム水等の対象水域を複数の監視ブロックに分割し、監
視ブロック毎の流速、水温、水域への流入・流出情報に
基づいて水域全体の対流を計算し、この対流計算値か
ら、水質情報が監視ブロック間で拡散する状態をモデル
化し、次いで、現在を初期値として一定時間後の水質情
報とプランクトン発生量を計算し、これらの計算値を基
に栄養塩指標、有機物指標、プランクトン指標を計算
し、一定時間後の水域全体の状態をシミュレ−ション計
算し、水質変動を予知することによって、達成される。

【０００５】

【作用】湖沼や河川中のプランクトン数は、通常、高々
数万(細胞数／ｍｌ)と希薄である。そのため、サンプリ
ング手段、濃縮手段及び移送手段により、希薄な懸濁浮
遊物質を濃縮し、画像認識の効率を高める。また、拡大
手段と照明手段及び撮像手段により構成される撮像装置
は、移送流動する液中の浮遊物質を完全に保持し、静止
画面を得る。この撮像装置において、プランクトン抽出
手段は、形状認識画像を基に、画像間の座標や形状の特
徴量などを考慮してプランクトンの種類や個数、浮遊物
質の面積や粒径を計測し、さらに形状毎に分類してプラ
ンクトンの分布やその質量を計測する。これらの画像計
測情報と水質情報を連続計測し、情報の時系列変化、気
象情報及び水域の地形情報から水質の汚染状態とその範
囲、原因を判定し、表示する。さらに、シミュレーショ
ン手段により汚染状況を推定し、汚染水域の浄化対策を
推定する。このように、連続計測情報から水質汚染地域
の同定、及び、事前にその予測を行なうことができると
共に汚染水域の浄化対策を講ずることが可能になる。

【０００６】

【実施例】以下、本発明の実施例を図面により説明す
る、図１は、本発明の一実施例であり、水質監視装置の
全体構成のブロック図を示す。本実施例の水質監視装置
は、サンプリング装置２、濃縮装置３、撮像装置４、画
像処理装置５、水質計測装置６、気象情報計測装置７、
計算機８、表示装置９から構成する。

【０００７】次に、これら装置の詳細を説明する。サン
プリング装置２は、湖沼等の閉鎖水域１の貯溜水をサン
プリングする。水域１内の懸濁浮遊物質、特に植物プラ
ンクトンは、位置、時間、季節により時々刻々変化する
ため、この変化に対処したサンプリングを実施する必要
がある。図２において、貯溜水をサンプリングするに
は、位置調節装置２０２により導通管２０３を上下方向
に操作し、導入口２０４の水深位置を変化させ、サンプ
リングポンプ２０１がサンプル液を導入口２０４から吸
引、取水して、濃縮装置３に送る。ここで、水深位置操
作は、時間帯、季節等を考慮して自動的あるいは遠隔地
より調節する。導入口２０４を水面付近にすると、アオ
コ、赤潮の原因とされる浮上性の植物性プランクトンの
サンプリングが可能になる。サンプリング装置２からの
サンプル液は濃縮装置３に導かれ、懸濁浮遊物質が濃縮
される。湖沼やダムの水域１は、雨天時には数千ｐｐｍ
の濁質となる場合もあるが、通常は数十ｐｐｍ以下と清
澄である。また、植物プランクトンは、大量発生でも数
万（細胞数／ｍｌ）であり、平均数百（細胞数／ｍｌ）
とされている。これを対物レンズ４倍の顕微鏡によって
観察しても、１画面（観察容積約０．２ｍｍ³）に１細
胞あるか否かであり、効率的な監視とは云えない。そこ
で、効率的な監視を行うにはサンプル液を濃縮すること
が必要となる。図３に、濃縮装置３の一例を示す。サン
プル液は、導管３０２により散水され、網状の分離膜３
０３において懸濁物質が濃縮分離され、濃縮槽３０４に
貯溜する。サンプル液量と濃縮液量の比で濃縮倍率が決
まる。懸濁物質が分離された濾過液は排水される。ま
た、濃縮液は濃縮槽３０４から引抜き管３０６を介して
ポンプ３０５により採取する。なお、懸濁物質を濃縮分
離する方式として、遠心力を利用して分離する遠心濾過
方式、密度差や重量差を利用した遠心分離方式あるいは
遠心沈降方式を採用してもよい。濃縮装置３からの濃縮
液は撮像装置４に導かれる。図４に、撮像装置４の一例
を示す。撮像装置４は、濃縮液の導入管４０９に対して
第１筐体４０１と第２筐体４０２から構成される。各々
の筐体が対向する導入管４０９の一部に透明ガラス製の
接液窓４０５と４０６を配置する。接液窓４０５、４０
６はその一方を平型に、他方を凹形にし、両者が接した
ときに形成される空間部をサンプル室とし、濃縮液の一
部を保持する。第１筐体４０１にはサンプル室に焦点を
合わせた拡大光学レンズ４０４とＩＴＶ（テレビカメ
ラ）４０３を内蔵し、一方、第２筐体４０２には集光レ
ンズを有する照明装置４０２と、接液窓４０６の位置を
調節する駆動装置４０８を設置する。このような構成に
おいて、駆動装置４０８は後述の画像処理装置５からの
指令あるいは外部からのタイマ−により制御され、接液
窓４０６を上下動操作する。図４は、サンプル室開放時
に接液窓４０６を下降させた状態であり、サンプル室の
液は濃縮液の流れにより入れ替えらる。そして接液窓４
０６を再び上昇させることにより、新たな濃縮液がサン
プル室に保持される。すなわち、サンプル室は接液窓４
０６上昇動作時に形成される。サンプル室内の濃縮液
は、透過光方式で照明され、拡大光学レンズ４０４を介
してＩＴＶ４０３が受光し、電気信号に変換される。こ
の時、サンプル室内の保持液は、導入管４０９内を流通
する濃縮液の影響を受けずに静止状態にある。これらの
操作は連続的に行われ、自動的に画面が更新される。な
お、接液窓４０６を下降させた状態で、導入管４０９内
にワイパ−（図示せず）を設置し、両接液窓面の洗浄
と、サンプル室保持液の強制的入れ替えを実行すること
もできる。撮像装置４において電気信号（映像信号）に
変換された画像は、画像処理装置５に入力され、懸濁浮
遊物質と植物プランクトンの形状を画像処理抽出する。
また、画像処理装置５からの情報に基づいて画像の取り
込み指令を行い、撮像装置４の動作制御を実施する。

【０００８】図５に、画像処理装置５の一実施例を示
す。画像処理装置５は、ＣＰＵ（中央処理装置）５０
１、主メモリ５０２、記憶装置５０３、通信インタフェ
ース５０４、システムバス５０５、画像処理部５０６か
ら構成する。ＣＰＵ５０１は、主メモリ５０２に格納さ
れているプログラムを実行し、画像処理装置５全体を制
御する。記憶装置５０３には、画像処理装置５を制御す
るプログラムと画像処理結果が格納される。プログラム
は、画像処理装置のイニシャルスタート時にＣＰＵ５０
１のマイクロプログラムによって主メモリ５０２に転送
される。通信インタフェース５０４は他計算機、計測装
置とのデータ送受信機能を持ち、ＣＰＵ５０１から通信
インタフェース５０４を介して、計算機８への画像計測
情報の送信及び撮像装置４の駆動装置４０８への制御信
号の送信を実施する。画像処理部５０６は、画像メモリ
５０７、画像処理プロセッサ５０８、映像インタフェー
ス５０９からなる。画像メモリ５０８には、縦方向２５
６画素、横方向２５６画素で輝度階調が２５６（８ビッ
ト）の濃淡画像メモリと、縦方向２５６画素、横方向２
５６画素で輝度階調が２（１ビット）の２値メモリがそ
れぞれ複数個配されている。画像処理プロセッサ５０８
は、画像メモリ５０８の濃淡画像並びに２値画像に対し
て、濃淡画像処理演算、２値化処理、形状特徴量抽出な
どを高速に実施する複数のＬＳＩから構成されている。
映像インタフェース５０９は、複数の入力と出力を有
し、入力にはＡ／Ｄ変換器、出力にはＤ／Ａ変換器を配
しており、アナログの映像信号をデジタル信号に変換す
る。画像処理装置５内部ではデジタル信号を用いて処理
する。撮像装置４からの映像信号は、映像インタフェー
ス５０９においてデジタル信号に変換され、１つの濃淡
画像として画像メモリ５０７にリアルタイムで格納さ
れ、画像処理プロセッサ５０８は画像メモリ５０７に格
納される該濃淡画像に対して各種画像処理を実行する。
なお、画像処理のタイミングはＣＰＵ５０１が制御す
る。映像インタフェース５０９に接続したモニタテレビ
５１０では、画像メモリ５０７に格納されている画像や
撮像装置４の映像信号を表示する。画像処理プロセッサ
５０８は駆動装置４０８を駆動させて接液窓４０９を下
降、上昇動作させる。この下降、上昇動作によりサンプ
ル室の濃縮液が交換されて、新たな画像が撮像装置４に
取り込まれる。

【０００９】ここに、画像処理装置５の処理手順の詳細
を図６、図７、図８を用いて説明する。本発明者らの観
察によれば、透過照明方式の撮像装置４の拡大画像で
は、植物プランクトンは一部光が透過して輝度が高いも
のの、液相部よりは暗い。また、ゴミ状の物質はプラン
クトンより輝度が低い。図６の濃淡画像(a)は、濃縮液
の植物プランクトン画像の一例である。濃淡画像は、同
図(a)に示すように、明るい液相部Ｚの領域と、液相部
より暗い植物性プランクトンＡ、Ｂ、Ｃと、ゴミ状の物
質Ｄを含んでいる。このように濃淡画像は、輝度(明る
さ)に応じた濃淡情報を持っており、該濃淡画像を２値
化処理して、植物性プランクトンのみを分離抽出する。
すなわち、特定の輝度レベルｈ₁を基準にして、輝度レ
ベルｈ₁より低輝度の領域を１、ｈ₁より高輝度の領域を
０にする２値化処理を実行する。この０、１の２値情報
に変換された２値画像を図６（ｂ）に示す。まず、液相
部輝度と懸濁物質（プランクトンとゴミ状の物質）の輝
度間に２値化の輝度レベルｈ₁を設定し、懸濁物質を１
に、液相部を０にした懸濁物質の２値画像を得る。次
に、ゴミ状物質部輝度とプランクトンの輝度間に２値化
の輝度レベルｈ₂を設定し、ゴミ状物質を１に、他を０
にしたゴミ状物質の２値画像を作成し、前記懸濁物質の
２値画像と差分することで図６（ｂ）の２値画像を得
る。図７に、この２値化処理の手順を示す。まず、内部
のシステム時刻を読み出し、予め設定した起動時刻ある
いは起動周期と一致した場合画像処理を実行する（５４
０）。画像処理プロセッサ５０８は、映像インタフェー
ス５０９を介して得られたプランクトン映像を画像メモ
リ５０７内の濃淡画像メモリに格納後、濃淡画像処理演
算、２値化処理を施し、２値メモリに格納する（５４
１）。２値化処理の詳細については、図８において説明
する。さらに、この２値メモリに対し、ラべリング処理
を実行し、画像全体の物体総個数、全面積を計算し（５
４２）、続いて、粒径分布、平均粒径等（５４３）を計
算する。次に、プランクトン個々について、面積、周囲
長、形状係数、２次モーメント、端点数、交点数、穴の
数などの形状特徴量を計算し、主メモリ５０２に格納す
る（５４４）。ＣＰＵ５０１は、これら各種画像処理計
測値を通信インタフェース５０４を介して計算機８へ送
信する（５４５）。さらに図８に、２値化処理の詳細な
手順を示す。撮像装置４のプランクトン映像信号は、濃
淡画像メモリＧ１に格納（５５０）後、輝度情報を輝度
−画素のヒストグラムとして抽出する（５５１）。輝度
−画素ヒストグラムは、横軸に輝度（０〜２５５）、縦
軸に画素数をとると画素数の曲線として表現できる。こ
の曲線から、例えば、曲線の最大値からある固定定数を
減算する方式を用いて植物性プランクトンのみを分離抽
出する特定の輝度レベルｈ₁を計算し（５５２）、輝度
レベルｈ₁を基準にして、輝度レベルｈ₁より低輝度の領
域を１、ｈ₁より高輝度の領域を０にする２値化処理を
実行、２値メモリＢ１に格納する（５５３）。２値抽出
されたプランクトンのみの輝度情報を計算するため、濃
淡画像メモリＧ１と２値メモリＢ１とＡＮＤ演算を実行
（マスク処理）し、背景液相部の輝度を０にした画像を
濃淡画像メモリＧ２に格納する（５５４）。濃淡画像メ
モリＧ２に対し、先のステップ５５１、５５２と同様の
方法で、輝度−画素ヒストグラム抽出し（５５５）、ゴ
ミ状物質部とプランクトンを分離する２値化の輝度レベ
ルｈ₂の計算（５５６）を行い、輝度レベルｈ₂による２
値化を実行し、２値メモリＢ２に格納する（５５７）。
なお、上記２値化処理手順において、濃淡画像（ａ）を
直接処理したが、各画像の背景画像（液相部のみの画
像）を予め撮像し、その画像と差分処理した画像を対象
に画像処理を実行しても良い。この前処理を実行するこ
とにより、画像全体の輝度ムラ、すなわち、照明光によ
る明るさの影響をなくし、良好な抽出画像を得ることが
できる。また、２値化処理のための輝度レベルの設定は
固定値並びに各画像の輝度分布（ヒストグラム）を考慮
して変化させる自動２値化法を用いることができる。

【００１０】計算機６は、画像処理装置５からの画像情
報と、気象情報計測装置７からの気象情報と、水質計測
装置６からの水質情報を連続的に受信し、水域の汚染状
態を監視する。ここで、水質情報とは、水温、濁度、ｐ
Ｈ、ＤＯ（容存酸素濃度）、電気伝導度、化学的酸素要
求量ＣＯＤ、流速、水位などであり、気象情報とは、気
温、風向、風力、日射量、雨量などである。画像処理装
置５、気象情報計測装置７、水質計測装置６は、監視対
象水域の地理的条件を考慮し、監視に最適な位置に設置
する。水域を複数の監視ブロックに分割し、各ブロック
の代表地点に各計測装置を配置することが望ましい。こ
こで、図９に、計算機６の構成例を示す。ＣＰＵ（中央
処理装置）８０５は、プログラム記憶装置８０３に格納
されている各種プログラムを主メモリ８０６に転送して
実行する。実行タイミングは、映像インタフェース８０
８に接続した表示装置９からの起動指令信号、タイマ
ー、入出力インタフェース８０７からの情報などによ
る。計測値データーベース８０１には入出力インタフェ
ース８０７を介して受信した画像情報、水質情報、気象
情報及びシミュレーション結果が格納され、地形データ
ーベース８０２には水域の形状、流入水量、流入箇所、
流出水量、流出箇所等のデータが格納される。プログラ
ム記憶装置８０３には、異常診断プログラム、水質シミ
ュレーションプログラム、データ送受信プログラム、表
示装置操作プログラム等が格納されている。知識ベース
８０９には、過去の現象に基づくルール群など知識工
学、ファジー、ニューラルネットに必要な情報が格納さ
れる。異常診断プログラムは、各種計測情報並びに知識
ベ−ス８０９の格納ル−ルに基づいて水域の汚染状況と
その要因及び浄化対策を診断する。まず、プランクトン
総体積、個数などを用いて植物プランクトン増減傾向を
把握し、植物プランクトンが影響要因であれば、その種
類を同定し、知識ベ−ス８０９のル−ルから増殖要因を
診断する。種類の同定は、形状の特徴から求める。例え
ば、アオコや赤潮の原因となるミクロキスチスやウロゲ
レナは円形状係数や粒径、さらに細胞と中空（穴）部の
面積比等の計測情報を用いて同定する。また、糸状や矩
形状のプランクトンは長軸・短軸長比、面積と周囲長
比、端点や交点数等から同定する。さらに、星形状や連
環状のプランクトンは端点や交点数、細胞と中空（穴）
部の面積比、穴数、細胞と中空（穴）部の面積比、面積
と周囲長比等から判定する。これらの判定には、クリス
プル−ルやファジ−ル−ルによる知識工学的手法、並び
に、各種計測情報を入力層に、計測情報に対応する特徴
値を教師デ−タとするニュ−ラルネットワ−ク手法を用
いることもできる。各種ル−ルや特徴値及びそれに対応
するプランクトン名等は知識ベ−ス８０９に格納し、必
要に応じて呼出して推論を実行する。このように、水域
の汚染状態をプランクトンの出現量やその出現種を連続
計測して監視する。

【００１１】図１０に、この異常診断プログラムの処理
の一例を示す。ＣＰＵ８０５は、計算機８の起動と同時
に異常診断プログラムを主メモリ８０６にロードし、実
行する。まず、異常診断プログラムは、システム時刻を
読みだし、予め登録した時刻または周期と比較し、起動
時刻と判定した場合、ステップ８５１に移行する（８５
０）。続いて、水域の監視ブロック毎に、栄養塩指標Ｉ
ｋ、有機物指標Ｉｃ、プランクトンＩｐ指標を計算する
（８５１）。計算手順は後述する図１１、図１２におい
て説明する。次に、水質情報、画像情報及び該３つの指
標を初期値として、今後１カ月の水質情報、画像情報及
び該３つの指標をシミュレーションする（８５２）。シ
ミュレーションの手順は後述する図１３において説明す
る。現状値またはシミュレーション値が基準範囲を超え
た場合、水質情報や出現プランクトンの種類及び量が年
間変動と一致しないときに異常発生と判定し（８５
３）、異常発生場所、異常要因となる流入水をシミュレ
ーション結果もしくは過去の経験則から推定する（８５
４）。異常要因に対する対策は、知識ベースに格納して
おき、発生場所、発生要因、季節、気象などから推論す
る（８５５）。図１１に、栄養塩指標Ｉｋ、有機物指標
Ｉｃの計算方法の一例を示す。水質情報からＤＯ、濁度
Ｔｕ、リンＰ、窒素ＴＮ、化学的酸素要求量ＣＯＤ、ｐ
Ｈを読み出し（８６１）、栄養塩指標Ｉｋ、有機物指標
Ｉｃをそれぞれ（１）式、（２）式により計算する（８
６２）。Ｉｋ＝ｗ₁×Ｐ＋ｗ₂×ＴＮ（１）Ｉｃ＝ｗ₃×ＤＯ＋ｗ₄×ＣＯＤ＋ｗ₅×ｐＨ（２）ただし、ｗ₁，ｗ₂，₃ｗ₄，ｗ₅は重み係数。降雨、台風などの要因では、水質の変動が予想されるの
で、降雨量、水域への流入量、濁度Ｔｕを変数に持つ係
数ｗ₆，ｗ₇を用いて（３）、（４）式により補正する
（８６３）。Ｉｋ＝ｗ₆× Ｉｋ（３）Ｉｃ＝ｗ₇× Ｉｃ（４）以上のように計算した指標は、主メモリ８０６と計測値
データベース８０１に格納する（８６４）。図１２に、
プランクトン指標Ｉｐの計算手順の一例を示す。まず、
画像情報を基にプランクトンの分類を形状の特徴から数
個（Ｎ）に分類した（８７０）後、分類毎にプランクト
ンの体積ＰＶ（ｉ）、個数ＰＮ（ｉ）を計算（８７１）
する。糸状や矩形状のプランクトンは、長軸・短軸長
比、面積と周囲長比、端点や交点数等から同定できる。
さらに、星形状や連環状のプランクトンは、端点や交点
数、細胞と中空（穴）部の面積比、穴数、細胞と中空
（穴）部の面積比、面積と周囲長比等から判定できる。
次に、アオコや赤潮の原因となる特定種のプランクトン
の存在有無は形状特徴から判定する（８７２）。例え
ば、アオコや赤潮の原因となるミクロキスチスやウロゲ
レナは円形状係数や粒径、さらに細胞と中空（穴）部の
面積比等の情報を用いて同定する。存在を確認した場
合、異常指標のプランクトン発生を主メモリ８０６及び
計測値データベースに格納する（８７４）。プランクト
ン発生量ＶＶとプランクトン指標Ｉｐを（５）、（６）
式により計算し、主メモリ８０６及び計測値データベー
ス８０１に格納する（８７４、８７６）。ただし、ｗ₈（ｉ）：体積計算時のプランクトン分類ｉ
の重み係数Ｃ：係数

【００１２】図１３に、シミュレ−ション手順の詳細を
示す。まず、水域への流入・流出情報、監視ブロック毎
の水質情報、気象情報、地形データベース８０２から水
域の地形情報を読み出し（８８０）、監視ブロック毎
の、流速、水温、水域への流入・流出情報に基づき水域
全体の対流を計算する（８８１）。この対流計算値と水
温から、水質情報が監視ブロック間で拡散する状態をモ
デル化する。次に、現在の水質情報を初期値として一定
時間後のＤＯ、濁度Ｔｕ、リンＰ、窒素ＴＮ、化学的酸
素要求量ＣＯＤ、ｐＨを計算し、続いて、画像情報と栄
養塩、水温、ＤＯ、ＣＯＤを初期値とし、プランクトン
繁殖モデルを用いてプランクトン発生量を計算する（８
８２）。これらの計算値を基に栄養塩指標Ｉｋ、有機物
指標Ｉｃ、プランクトン指標Ｉｐを計算する（８８
３）。ステップ８８２、８８３を監視ブロック数分繰返
す（８８４）と一定時間後の水域全体の状態が計算でき
る。さらに、時間を更新し、１カ月間程度の計算を繰り
返し（８８５）、結果を主メモリ８０６及び計測値デー
タベース８０１に格納する（８８６）。このように、水
質シミュレ−ションプログラムは、当日の画像情報、水
質情報、気象情報及び水域に流入する河川等の流入量、
水質並びに水圏固有の地形、大きさなどの地理情報を基
にシミュレ−ション計算を行い、汚染発生の原因となる
場所の推定や翌日の水域内対流、拡散、反応などの水質
変動を解析する。すなわち、翌日の水質変動を予知する
ことにより、その対策を早期に行うことができる。診
断、シミュレーションにより異常を検知した場合、空気
又は酸素吹き込みによる曝気操作、撹拌操作、紫外線照
射などの浄化手段を起動する。浄化手段は、水域内複数
箇所に固定設置する方式と汚染箇所に移動する方式があ
る。表示装置９は、計算機８内の情報を表示し、計測デ
ータ、汚染状態、シミュレーション結果などを水域の水
平方向と垂直方向について表示する。また、マウス、ラ
イトペン、キーボードなどの信号を計算機８に送信す
る。

【００１３】

【発明の効果】本発明によれば、湖沼や河川中のプラン
クトン数が希薄であっても、懸濁浮遊物質を濃縮するこ
とにより、画像認識の効率を高めことができる。また、
閉鎖性水域の汚染状態をプランクトン出現量と水質情報
から連続的かつ定量的に計測できるので、水域の水質監
視及び浄化を効率的に実行することが可能となる。ま
た、シミュレーションにより、汚染状況を事前に予測す
ることができると共に汚染水域の浄化対策を講ずること
が可能になる。

【図面の簡単な説明】

【図１】本発明の一実施例を示す図。

【図２】サンプル装置の構成を示す図。

【図３】濃縮装置の構成を示す図。

【図４】撮像装置の構成を示す図。

【図５】画像処理装置のを示す図。

【図６】プランクトン画像の一例を示す図。

【図７】２値化処理の手順を示す図。

【図８】２値化処理の詳細な手順を示す図。

【図９】計算機の構成を示す図。

【図１０】異常診断プログラムの処理を示す図。

【図１１】計算方法の一例を示す図。

【図１２】計算方法の一例を示す図。

【図１３】シミュレ−ション手順の詳細を示す図。

【符号の説明】

１水域２サンプリング装置３濃縮装置４撮像装置５画像処理装置６水質計測装置７気象情報計測装置８計算機９表示装置

───────────────────────────────────────────────────── フロントページの続き (72)発明者呉文智茨城県日立市大みか町五丁目２番１号株式会社日立製作所大みか工場内 (72)発明者渡辺昭二茨城県日立市久慈町4026番地株式会社日立製作所日立研究所内 (72)発明者矢萩捷夫茨城県日立市久慈町4026番地株式会社日立製作所日立研究所内 (56)参考文献特開平３−61836（ＪＰ，Ａ) 実開昭50−114360（ＪＰ，Ｕ) (58)調査した分野(Int.Cl.⁷，ＤＢ名) G01N 15/06 C02F 1/00 C02F 3/00 G01N 33/18

Claims

(57)【特許請求の範囲】

【請求項１】湖沼水、海水、河川水、ダム水等の水質
を自動的にサンプルする手段と、サンプル中の浮遊物質
を濃縮する手段と、濃縮液中の濁質を撮像する手段と、
撮像画像を画像処理し、濁質を抽出する手段と、該抽出
濁質の特徴を計算する手段と、該特徴から浮遊物質を複
数の形状に分類する手段を具備し、浮遊物質の形状分類
毎に出現量を計測することを特徴とする水質汚染監視装
置。
【請求項２】請求項１において、分類毎の出現量分布
と、水質、季節、気象データから汚染状態を判定するこ
とを特徴とする水質汚染監視装置。
【請求項３】請求項１または請求項２において、分類
毎の出現量分布計測は、水域の水平方向及び水深方向に
ついて実施することを特徴とする水質汚染監視装置。
【請求項４】請求項１、請求項２または請求項３にお
いて、対象水域の地形、面積等の地理情報を記憶する手
段と、汚染状況をシミュレーションする手段を具備し、
汚染箇所を推定することを特徴とする水質汚染監視装
置。
【請求項５】請求項４において、シミュレーション結
果に基づいて汚染領域を判定し、該領域を浄化する手段
を具備することを特徴とする水質汚染監視装置。
【請求項６】請求項１から請求項５のいずれかにおい
て、それぞれの結果を表示する手段を具備することを特
徴とする水質汚染監視装置。
【請求項７】湖沼水、海水、河川水、ダム水等の対象
水域を複数の監視ブロックに分割し、該監視ブロック毎
に栄養塩指標、有機物指標、プランクトン指標を計算
し、これらの指標と水質情報及び浮遊物質の濁質を画像
処理した画像情報を基にシミュレーションし、該シミュ
レーション値が基準範囲を超えた場合、出現プランクト
ンの種類及び量が異常発生したと判定し、異常診断する
ことを特徴とする水質汚染監視方法。
【請求項８】請求項７において、異常診断の結果から
異常発生場所、異常要因となる流入水を推定し、該異常
要因に対する対策を推論することを特徴とする水質汚染
監視方法。
【請求項９】湖沼水、海水、河川水、ダム水等の対象
水域を複数の監視ブロックに分割し、監視ブロック毎の
流速、水温、水域への流入・流出情報に基づいて水域全
体の対流を計算し、該対流計算値から水質情報が監視ブ
ロック間で拡散する状態をモデル化し、次いで、現在を
初期値として一定時間後の水質情報とプランクトン発生
量を計算し、これらの計算値を基に栄養塩指標、有機物
指標、プランクトン指標を計算し、一定時間後の水域全
体の状態をシミュレ−ション計算し、水質変動を予知す
ることを特徴とする水質汚染監視方法。