(57)【特許請求の範囲】
【請求項１】
予測対象海域におけるフィールド調査水質データとして、水中の懸濁物質（ＳＳ）、強熱減量（ＩＬ）、懸濁態有機炭素（ＰＯＣ）及びクロロフィル濃度を求めるステップＡ、
ステップＡで求めた水質データを基に、懸濁物質（ＳＳ）に含まれる有機物（ＳＳorg）と無機物（ＳＳin）の量を下記式により算出するステップＢ、
[数１]
ＳＳorg＝ＳＳ×ＩＬ
ＳＳin＝ＳＳ−ＳＳorg
波浪や流れにより底泥から巻き上がった懸濁物質（ＳＳbm）に含まれる無機成分（ＳＳin-bm）を下記式により抽出するステップＣ、
[数２]
ＳＳin-bm＝ＳＳin−ＳＳin-off
（ここでＳＳin-offは巻上げの影響の小さい沖合表層の無機物（ＳＳin）の平均値をバックグランドとしたもの）
巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を下記式により抽出するステップD、
[数３]
ＳＳorg＝ａ・Ｃｈｌａ＋ｂ
（ここでａは係数、ｂは定数で、これらは実験値から求める。Ｃｈｌａはクロロフィル量）
ステップＤで抽出した一次式を用いて有機物（ＳＳorg）をバックグラウンド（ＳＳorg-off）として、巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を下記式により算出するステップE、
[数４]
ＳＳorg-bm＝ＳＳorg−ＳＳorg-off
上記で求めた値に基づき、底泥からの巻上げ懸濁物質濃度（ＳＳbm）を下記式により算出するステップF、
[数５]
ＳＳbm＝ＳＳin-bm＋ＳＳorg-bm
上記で求めたＳＳbmの値に基づき下記式により巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）の割合を算出するステップG、
[数６］
ＳＳorg-bm＝ｃ・ｌｏｇ（ＳＳbm）＋ｄ
又はｌｏｇ（ＳＳorg-bm）＝ｃ・ｌｏｇ（ＳＳbm）＋ｄ
（ここで、ｃは係数、ｄは定数で、これらは実験値から求める）
上記で求めたＳＳorg-bmの値に基づき下記式により巻上げ懸濁物質（ＳＳbm）の有機成分（ＳＳorg-bm）の割合の有機炭素量（ＰＯＣbm）への変換を行うステップＨ
[数７]
ＰＯＣbm／ＳＳorg-bm＝ｐ（ＳＳorg-bm）ｑ
ＰＯＣbm＝ＰＯＣ×（ＳＳorg-bm／ＳＳorg）
（ここで、ｐは係数、ｑはべき数でこれらは実験値より求める）
を有することを特徴とする巻上げ有機物量予測計算方法。
【請求項２】
予測対象海域におけるフィールド調査水質データとして、水中の懸濁物質（ＳＳ）、強熱減量（ＩＬ）、懸濁態有機炭素（ＰＯＣ）及びクロロフィル濃度を求めるステップＡ、
ステップＡで求めた水質データを基に、懸濁物質（ＳＳ）に含まれる有機物（ＳＳorg）と無機物（ＳＳin）の量を下記式により算出するステップＢ、
[数１]
ＳＳorg＝ＳＳ×ＩＬ
ＳＳin＝ＳＳ−ＳＳorg
波浪や流れにより底泥から巻き上がった懸濁物質（ＳＳbm）に含まれる無機成分（ＳＳin-bm）を下記式により抽出するステップＣ、
[数２]
ＳＳin-bm＝ＳＳin−ＳＳin-off
（ここでＳＳin-offは巻上げの影響の小さい沖合表層の無機物（ＳＳin）の平均値をバックグランドとしたもの）
巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を下記式により抽出するステップD、
[数３]
ＳＳorg＝ａ・Ｃｈｌａ＋ｂ
（ここでａは係数、ｂは定数で、これらは実験値から求める。Ｃｈｌａはクロロフィル量）
ステップＤで抽出した一次式を用いて有機物（ＳＳorg）をバックグラウンド（ＳＳorg-off）として、巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を下記式により算出するステップE、
[数４]
ＳＳorg-bm＝ＳＳorg−ＳＳorg-off
上記で求めた値に基づき、底泥からの巻上げ懸濁物質濃度（ＳＳbm）を下記式により算出するステップF、
[数５]
ＳＳbm＝ＳＳin-bm＋ＳＳorg-bm
上記で求めたＳＳbmの値に基づき下記式により巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）の割合を算出するステップG、
[数６]
ＳＳorg-bm＝ｃ・ｌｏｇ（ＳＳbm）＋ｄ
又はｌｏｇ（ＳＳorg-bm）＝ｃ・ｌｏｇ（ＳＳbm）＋ｄ
（ここで、ｃは係数、ｄは定数で、これらは実験値から求める）
上記で求めたＳＳorg-bmの値に基づき下記式により巻上げ懸濁物質（ＳＳbm）の有機成分（ＳＳorg-bm）の割合の有機炭素量（ＰＯＣbm）への変換を行うステップHと、
[数７]
ＰＯＣbm／ＳＳorg-bm＝ｐ（ＳＳorg-bm）ｑ
ＰＯＣbm＝ＰＯＣ×（ＳＳorg-bm／ＳＳorg）
（ここで、ｐは係数、ｑはべき数でこれらは実験値より求める）
ステップHで求めた有機炭素量（ＰＯＣbm）の値を用いて水質を予測計算するステップＩを有することを特徴とする巻上げ水質予測計算方法。
【請求項３】
水質を予測計算するステップＩにおいて、有機物生成・分解過程をモデル化した水質予測計算モデルを用いる請求項２の水質予測計算方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この出願の発明は、巻上げ有機物量予測計算方法および水質予測計算方法に関するものである。さらに詳しくは、この出願の発明は、沿岸海域等における水質の変動機構の解明や水質・底質環境の保全等に有用な、新しい巻上げ有機物量予測計算方法および水質予測計算方法に関するものである。
【０００２】
【従来の技術】
従来、波と流れにより水中へ巻き上げられる底泥の量（以下、巻上げ底泥量と呼ぶ）の予測計算は、海岸保存対策を目的とした漂砂による海浜の地形変化の予測計算や、航路の埋没対策（シルテーション対策とも呼ばれる）を目的とした粒径の小さな底泥の巻上げおよび輸送による航路埋没の予測計算などにおいて、盛んに研究・開発され、その実用化が進んでいる。
【０００３】
【発明が解決しようとする課題】
しかしながら、従来の巻上げ底泥量の予測計算は、予測対象物質が底泥中の懸濁物質（浮遊砂とも呼ばれる）であり、懸濁物質中の有機物のみを対象としてはいないといった問題があった。
【０００４】
すなわち、巻上げ底泥量の予測計算は具体的には波と流れにより水中へ巻き上げられる底泥からの懸濁物質の量（以下、巻上げ懸濁物質量と呼ぶ）を予測計算しており、他方、懸濁物質中には当然無機物も有機物も存在しているため、両者を含めた巻上げ懸濁物質量を予測計算しても、たとえば有機物生成・分解過程をモデル化した水質予測計算モデルを用いた水質予測計算では精度が十分に高くならないのである。
【０００５】
この出願の発明は、以上のとおりの事情に鑑みてなされたものであり、従来技術の問題点を解消し、底泥から巻き上げられる有機物の量（以下、巻上げ有機物量と呼ぶ）を予測計算することのできる、新しい巻上げ有機物量予測計算方法、ならびに、有機物生成・分解過程をモデル化した水質予測計算モデルによる水質予測計算をより精度良く行うことのできる、新しい水質予測計算方法を提供することを課題としている。
【０００６】
【課題を解決する手段】
この出願の発明は、上記の課題を解決するものとして、
予測対象海域におけるフィールド調査水質データとして、水中の懸濁物質（ＳＳ）、強熱減量（ＩＬ）、懸濁態有機炭素（ＰＯＣ）及びクロロフィル濃度を求めるステップＡ、
ステップＡで求めた水質データを基に、懸濁物質（ＳＳ）に含まれる有機物（ＳＳorg）と無機物（ＳＳin-bm）の量を下記式により算出するステップＢ、
[数１]
ＳＳorg＝ＳＳ×ＩＬ
ＳＳin＝ＳＳ−ＳＳorg
波浪や流れにより底泥から巻き上がった懸濁物質（ＳＳbm）に含まれる無機成分（ＳＳin）を下記式により抽出するステップＣ、
[数２]
ＳＳin-bm＝ＳＳin−ＳＳin-off
（ここでＳＳin-offは巻上げの影響の小さい沖合表層の無機物（ＳＳin）の平均値をバックグランドとしたもの）
巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を下記式により抽出するステップD、
[数３]
ＳＳorg＝ａ・Ｃｈｌａ＋ｂ
（ここでａは係数、ｂは定数で、これらは実験値から求める。Ｃｈｌａはクロロフィル量）
ステップＤで抽出した一次式を用いて有機物（ＳＳorg）をバックグラウンド（ＳＳorg-off）として、巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を下記式により算出するステップE、
[数４]
ＳＳorg-bm＝ＳＳorg−ＳＳorg-off
上記で求めた値に基づき、底泥からの巻上げ懸濁物質濃度（ＳＳbm）を下記式により算出するステップF、
[数５]
ＳＳbm＝ＳＳin-bm＋ＳＳorg-bm
上記で求めたＳＳbmの値に基づき下記式により巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）の割合を算出するステップG、
[数６]
ＳＳorg-bm＝ｃ・ｌｏｇ（ＳＳbm）＋ｄ
又はｌｏｇ（ＳＳorg-bm）＝ｃ・ｌｏｇ（ＳＳbm）＋ｄ
（ここで、ｃは係数、ｄは定数で、これらは実験値から求める）
上記で求めたＳＳorg-bmの値に基づき下記式により巻上げ懸濁物質（ＳＳbm）の有機成分（ＳＳorg-bm）の割合の有機炭素量（ＰＯＣbm）への変換を行うステップＨ
[数７]
ＰＯＣbm／ＳＳorg-bm＝ｐ（ＳＳorg-bm）ｑ
ＰＯＣbm＝ＰＯＣ×（ＳＳorg-bm／ＳＳorg）
（ここで、ｐは係数、ｑはべき数でこれらは実験値より求める）
を有することを特徴とする巻上げ有機物量予測計算方法（請求項１）を提供する。
【０００７】
また、この出願の発明は、予測対象海域におけるフィールド調査水質データとして、水中の懸濁物質（ＳＳ）、強熱減量（ＩＬ）、懸濁態有機炭素（ＰＯＣ）及びクロロフィル濃度を求めるステップＡ、
ステップＡで求めた水質データを基に、懸濁物質（ＳＳ）に含まれる有機物（ＳＳorg）と無機物（ＳＳin）の量を下記式により算出するステップＢ、
[数１]
ＳＳorg＝ＳＳ×ＩＬ
ＳＳin＝ＳＳ−ＳＳorg
波浪や流れにより底泥から巻き上がった懸濁物質（ＳＳbm）に含まれる無機成分（ＳＳin-bm）を下記式により抽出するステップＣ、
[数２]
ＳＳin-bm＝ＳＳin−ＳＳin-off
（ここでＳＳin-offは巻上げの影響の小さい沖合表層の無機物（ＳＳin）の平均値をバックグランドとしたもの）
巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を下記式により抽出するステップD、
[数３]
ＳＳorg＝ａ・Ｃｈｌａ＋ｂ
（ここでａは係数、ｂは定数で、これらは実験値から求める。Ｃｈｌａはクロロフィル量）
ステップＤで抽出した一次式を用いて有機物（ＳＳorg）をバックグラウンド（ＳＳorg-off）として、巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を下記式により算出するステップE、
[数４]
ＳＳorg-bm＝ＳＳorg−ＳＳorg-off
上記で求めた値に基づき、底泥からの巻上げ懸濁物質濃度（ＳＳbm）を下記式により算出するステップF、
[数５]
ＳＳbm＝ＳＳin-bm＋ＳＳorg-bm
上記で求めたＳＳbmの値に基づき下記式により巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）の割合を算出するステップG、
[数６]
ＳＳorg-bm＝ｃ・ｌｏｇ（ＳＳbm）＋ｄ
又はｌｏｇ（ＳＳorg-bm）＝ｃ・ｌｏｇ（ＳＳbm）＋ｄ
（ここで、ｃは係数、ｄは定数で、これらは実験値から求める）
上記で求めたＳＳorg-bmの値に基づき下記式により巻上げ懸濁物質（ＳＳbm）の有機成分（ＳＳorg-bm）の割合の有機炭素量（ＰＯＣbm）への変換を行うステップHと、
[数７]
ＰＯＣbm／ＳＳorg-bm＝ｐ（ＳＳorg-bm）ｑ
ＰＯＣbm＝ＰＯＣ×（ＳＳorg-bm／ＳＳorg）
（ここで、ｐは係数、ｑはべき数でこれらは実験値より求める）
ステップHで求めた有機炭素量（ＰＯＣbm）の値を用いて水質を予測計算するステップＩを有することを特徴とする巻上げ水質予測計算方法（請求項２）をも提供し、この予測計算方法における水質を予測計算するステップＩにおいて、有機物生成・分解過程をモデル化した水質予測計算モデルを用いること（請求項３）を提供する。
【０００８】
【発明の実施の形態】
この出願の発明は、巻上げ懸濁物質濃度と巻上げ懸濁物質中の有機物濃度との関係式が巻上げ懸濁物質濃度の関数として表すことができるという、この出願の発明の発明者等による全く新しい知見に基づいてなされたものであり、その関係式を用いることで巻上げ有機物量の予測計算を実現している。
【０００９】
上記関係式は、巻上げ量や水質予測の対象となる実際の海域毎にフィールド調査を行い、海域内における波や流れによる巻上げの大きな場所（たとえば水深の浅い場所）および小さな場所の水質データ（たとえば懸濁物質、懸濁物質の強熱減量、懸濁態有機炭素、クロロフィル濃度など）を解析することにより得ることができ、そしてそれを用いて、巻上げ有機物の巻上げ過程のモデル化を行い、予測対象海域毎の巻上げ有機物量を予測計算する。
【００１０】
より具体的には、まず、波と流れによる底泥からの懸濁物質の巻上げ・沈降過程を説明すると、たとえば図１に例示したように、巻上げ過程（図中Ｕｐ）では、底泥によりなる泥層から底面直上に懸濁物質が高濃度で含まれる底面直上層が形成され、その底面直上層から水中へ懸濁物質が巻き上げられる。
【００１１】
沈降過程（図中Ｗ１，Ｗ２）では、水中の巻上げ懸濁物質（ＳＳbm）に含まれる有機物は、土粒子等の無機物（ＳＳin-bm）と結合した有機物（ＳＳorg-bm1）と、無機物から遊離した有機物（ＳＳorg-bm2）に分かれて沈降し、流れにより輸送される。無機物と結合した有機物（ＳＳorg-bm1）の沈降速度は、無機物の沈降速度に左右され、無機物から遊離した有機物（ＳＳorg-bm2）の沈降速度よりも大きくなる。これらの有機物を含んだ巻上げ懸濁物質（ＳＳbm）では、輸送過程において時間の経過とともに、その有機物組成が変化し、無機物と有機物の混合粒径よりなる懸濁物質の平均的な沈降速度も変化していくと考えられる。
【００１２】
巻上げ有機物の巻上げ過程のモデル化は、この懸濁物質の巻上げ・沈降過程に基づき、上記関係式、および既存の巻上げ懸濁物質量予測計算モデル（たとえば底泥浮上量予測計算モデルや浮遊砂量予測計算モデルなど）と水質予測計算モデルとの整合性を踏まえて、たとえば次のように行う。
【００１３】
最初に、上記底面直上層は極薄い層であり、実際の海域では測定が難しく不確定な要素を多く含むため、仮想の層とする。次いで、フィールド調査水質データの解析により求められる巻上げ懸濁物質濃度と有機物濃度との関係式により、巻上げ懸濁物質を巻上げ有機物へ変換する。この変換には、変換のタイミングに従って、▲１▼水中への巻き上がり後に行う手法、▲２▼底面直上層で変換して水中への巻上げを計算する手法がある。これにより、図２に例示したような巻上げ過程モデルが得られる。
【００１４】
すなわち、この巻上げ有機物の巻上げ過程モデルが巻上げ有機物量の予測計算モデルであり、これに基づいて、上記既存の巻上げ懸濁物質量予測計算モデルにより巻上げ懸濁物質の濃度を予測計算し、それを上記関係式により巻上げ有機物量へ換算することで、底泥から巻き上げられる有機物量を正確に予測計算できるのである。
【００１５】
そしてさらに、この出願の発明によれば、上述のように予測計算された巻上げ有機物量を用いることで、有機物の生物化学的な生成・分解過程をモデル化した水質予測計算モデルによる水質予測計算を、より精度良く行うことができる。
【００１６】
なおこの場合、巻上げ有機物量の予測計算を、巻上げ有機物が水中起源の有機物と同様に分解・消費（被食）されるものとして、水質の予測計算と同時に行えるようにすべく、巻上げ有機物量の予測計算モデルは水質予測計算モデルに組み込んでもよい。
【００１７】
この出願の発明は、以上のとおりの特徴を有するものであるが、以下に、実施例を示し、さらに詳しくその実施の形態について説明する。
【００１８】
【実施例】
［実施例１］
ここでは、一実施例として、実際の予測対象海域におけるフィールド調査水質データに基づいた巻上げ懸濁物質濃度と有機物濃度との関係式の作成について説明する。
【００１９】
＜１．フィールド調査水質データ＞
関係式作成に用いたフィールド調査水質データは、下記のとおりである。
・調査場所：稲毛の浜および検見川の浜の砕波帯（測定点：Ｓ１〜Ｓ６，Ｎ１〜Ｎ５）および沖合水域（測定点：沖合１〜５）（図３参照）
・調査時期：平成１１年７月２８日、９月２６日、１０月２９日の計３回
・調査層：海底面上＋０．３ｍ或いは＋０．５ｍから海面までの２層或いは３層
・調査項目：水中の懸濁物質（ＳＳ）、強熱減量（ＩＬ）、
懸濁態有機炭素（ＰＯＣ）、クロロフィル濃度
【００２０】
＜２．懸濁物質中の有機物と無機物の算出＞
この水質データを基に、まず、懸濁物質（ＳＳ）に含まれる有機物（ＳＳorg）と無機物（ＳＳin）を算出する。これは、水中の懸濁物質（ＳＳ）を有機物（ＳＳorg）と無機物（ＳＳin）に区分し、各々次式により行う。
ＳＳorg＝ＳＳ×ＩＬ
ＳＳin＝ＳＳ−ＳＳorg
【００２１】
＜３．巻上げ懸濁物質の算出＞
続いて、巻上げ懸濁物質（ＳＳbm：波浪や流れにより底泥から巻き上がったＳＳ成分）を、巻上げの大きい場所と小さい場所の水質データを用いて、以下のように算出する。
【００２２】
＜３−▲１▼．巻上げ懸濁物質中の無機成分の算出＞
巻上げ懸濁物質（ＳＳbm）に含まれる無機成分（ＳＳin-bm）を抽出する。これは、巻上げの影響の小さい沖合表層の無機物（ＳＳin）の平均値をバックグラウンド（ＳＳin-off）とし、次式により行う。
ＳＳin-bm＝ＳＳin−ＳＳin-off
【００２３】
＜３−▲２▼．巻上げ懸濁物質中の有機成分の算出＞
巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を抽出する。有機物（ＳＳorg）には植物プランクトン由来の内部生産成分が含まれるが、クロロフィル量（Ｃｈｌａ）が植物プランクトン量の指標になるため、巻上げの影響の小さい沖合でのクロロフィル量（Ｃｈｌａ）と有機物（ＳＳorg）との関係を調べたところ、下記の一次式［ＳＳorg＝ａ・Ｃｈｌａ＋ｂ］で表せることがわかった（図４参照）。
【００２４】
７月    Ｙ＝０．０３８２Ｘ＋２．１６４  （Ｒ＝０．７７７）
９月    Ｙ＝０．０３０８Ｘ＋０．７１２１  （Ｒ＝０．９１４）
１０月  Ｙ＝０．１１７９Ｘ＋１．５９７１  （Ｒ＝０．８７６）
ここで、
Ｘ：Ｃｈｌａ濃度［μｇ／ｌ］
Ｙ：ＳＳorg濃度［ｍｇ／ｌ］
Ｒ：相関係数
【００２５】
この一次式を用いて計算した有機物（ＳＳorg）をバックグラウンド（ＳＳorg-off）として、巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）を次式により算出する。
ＳＳorg-bm＝ＳＳorg−ＳＳorg-off
【００２６】
＜４．巻上げ懸濁物質濃度の算出＞
以上から、底泥からの巻上げ懸濁物質濃度（ＳＳbm）は、
ＳＳbm＝ＳＳin-bm＋ＳＳorg-bm
として算出する。
【００２７】
＜５．巻上げ懸濁物質の濃度分布＞
このようにして得られた巻上げ懸濁物質濃度（ＳＳbm）の分布について解析すると、波による巻上げ現象が明瞭であった７月期における濃度分布の特徴は、以下のとおりであった。
【００２８】
すなわち、巻上げがないと仮定した沖合表層のＳＳbmは０〜１．６（ｍｇ／ｌ）程度の範囲であり、ＳＳbmの算出に際して、この程度の誤差は含まれるものと考えられる。突堤等の養浜構造の違いから波による巻上げが大きいと予想された稲毛の浜においてＳＳbmは高い値を示し、岸寄りの海域（汀線〜沖合２００ｍ付近）の下層が２０〜２００（ｍｇ／ｌ）程度の範囲にあった。一方、巻上げが小さいと予想された検見川の浜のＳＳbmは、稲毛の浜の１／１０程度であった。
【００２９】
＜６．巻上げ懸濁物質と有機成分の関係＞
次いで、巻上げ懸濁物質（ＳＳbm）に含まれる有機成分（ＳＳorg-bm）の割合を調べたところ、ＳＳbmの濃度が大きいほどＳＳorg-bmの割合が小さくなる傾向がみられ、このＳＳbmとＳＳorg-bmとの関係は、次式［ＳＳorg-bm＝ａ・ｌｏｇ(ＳＳbm)＋ｂまたはＬｏｇ（ＳＳorg-bm）＝ａ・ｌｏｇ(ＳＳbm)＋ｂ］で表せることがわかった（図５（ａ）（ｂ）（ｃ）参照）。
【００３０】
７月  Ｙorg-bm＝０．５１６３ｅｘｐ（１．４８１９Ｌｏｇ10（Ｘbm））
（Ｒ2＝０．６２９４）
９月   Ｙorg-bm＝１．５２７２Ｌｏｇ10（Ｘbm）＋０．１４７８
（Ｒ2＝０．４２６３）
１０月  Ｙorg-bm＝２．０５１４Ｌｏｇ10（Ｘbm）＋０．０５５７
（Ｒ2＝０．７２９６）
ここで、
Ｘbm：巻上げ懸濁物質濃度［ｍｇ／ｌ］
Ｙorg-bm：巻上げ懸濁物質中の有機物濃度［ｍｇ／ｌ］
Ｒ2：決定係数
【００３１】
＜７．巻上げ懸濁物質の有機成分と有機炭素量の関係＞
ここでは水質予測モデルとして有機物を炭素量で表現するモデルを用いているため、上記方法で算出した巻上げ懸濁物質（ＳＳbm）の有機成分（ＳＳorg-bm）を有機炭素量（ＰＯＣbm）に変換する必要がある。
【００３２】
予測対象海域の水質予測モデルの構築は巻上げ現象が明瞭であった７月を対象に行ったので、巻上げ懸濁物質（ＳＳbm）の有機成分（ＳＳorg-bm）の有機炭素量（ＰＯＣbm）への変換は、７月の水質データの整理・解析により得られた下記の関係式を用いた。
【００３３】
ＰＯＣbm／ＳＳorg-bm＝０．１５７６（ＳＳorg-bm）-0.2976
（Ｒ2＝０．４４２３）
ＰＯＣbm＝ＰＯＣ×（ＳＳorg-bm／ＳＳorg）
ここで、
ＳＳorg-bm：巻上げ懸濁物質の有機成分
ＰＯＣbm：巻上げ懸濁物質の有機成分中の有機炭素量
【００３４】
＜８．関係式＞
以上の＜６＞および＜７＞における関係式が、フィールド調査水質データに基づいて作成された巻上げ懸濁物質濃度（ＳＳbm）と巻上げ懸濁物質中の有機物濃度（ＳＳorg-bm）との関係式であり、巻上げ懸濁物質濃度（ＳＳbm）の関数として表されていることがわかる。したがって、以上と同様にして、予測対象海域毎に関係式を求めることで、巻上げ有機物量の正確な予測計算を実現することができる。
【００３５】
［実施例２］
ここでは、実施例１の関係式により得られた巻上げ有機物量に基づき、有機物生成・分解過程をモデル化した水質予測計算モデルを用いて水質予測を計算する場合において行う、巻上げ有機物量の予測計算モデルと水質予測計算モデルの組込みの一例について説明する。
【００３６】
＜１．水質予測モデルにおける有機物の生成・分解過程への組込み＞
＜１−▲１▼．懸濁態有機物の区分＞
まず、次式のように、水中の懸濁態有機物（ＰＯＭ）を水中起源の懸濁態有機物（ＰＯＭwa）と巻上起源の懸濁態有機物（ＰＯＭbm）に区分する。
ＰＯＭ＝ＰＯＭwa＋ＰＯＭbm
【００３７】
これら各懸濁態有機物（ＰＯＭwa，ＰＯＭbm）の生物化学的な生成・分解過程は次のとおりにモデル化することができる。
【００３８】
0004115109.tif 000002
【００３９】
0004115109.tif 000003
【００４０】
＜２．有機物の巻上げ過程の組込み＞
予測対象海域のフィールド調査水質データより得た巻上げ懸濁物質の巻上げ有機物量への換算式（実施例１参照）が水中の巻上げ懸濁物質を対象としたものであることから、予測対象海域の水質予測モデルの構築においては、巻上げ懸濁物質の巻上げ有機物への変換は水中の巻き上がり後に行う方法を用いる。
【００４１】
巻上げ有機物量の計算は、波浪と流れの場の条件を基にした、底面直上層の巻上げ懸濁物質量の計算結果を用いて次のとおりに行う。
【００４２】
＜２−▲１▼．水中の巻上げ懸濁物質量の計算＞
まず、底面直上層から水中への輸送は鉛直拡散によるものとし、次式により底面直上層から水中への懸濁物質の鉛直拡散量Ｄ、つまり巻上げ量Ｄの算出を行う。
【００４３】
Ｄ＝Ｋｚｂ（ｄＳＳbm／ｄｚ）
ここで、
Ｋｚｂ：鉛直拡散係数
ｄＳＳbm／ｄｚ：底面直上層と上層の水中間の鉛直方向の巻上げ懸濁物質（ＳＳbm）の濃度勾配
【００４４】
＜２−▲２▼．巻上げ有機物量の計算＞
そして、水中への巻上げ懸濁態有機物量（ＰＯＭbm）は、予測対象海域のフィールド調査水質データの解析により得られた関係式（実施例１参照）を用いて、上記懸濁物質の巻上げ量Ｄの換算により算出する。
【００４５】
以上のようにして、巻上げ有機物量の予測計算モデルを水質予測計算モデルに組み込むことで、巻上げ有機物量の計算を、巻上げ有機物が水中起源の有機物と同様に分解・消費（被食）されるものとして、水質の予測計算と同時に行うことができる。
【００４６】
［実施例３］
図６は、一実施例としてのこの出願の発明による巻上げ有機物量の予測計算から水質の予測計算までの全体の流れを例示したものである。本実施例では、突堤等の建造物があり、且つ干潟部分がある沿岸海域を予測対象海域としている。
【００４７】
＜１．巻上げ有機物量の予測計算＞
＜１−▲１▼＞まず、懸濁物質の巻上げを発生させる外力は波と流れであるので、本実施例では、その波と流れを、下記の既存の波浪の予測計算モデルと流れの場の予測計算モデルを用いて計算する。
【００４８】
Ａ．波浪の予測計算モデル
予測対象海域には突堤等の構造物があり、波の変形計算には構造物の反射、回折現象等を精度良く再現できる必要がある。そこで、本実施例では、非定常緩勾配方程式を基本として不規則波扱いに拡張した既存の予測計算モデルを用いる。
【００４９】
Ｂ．流れの場の予測計算モデル
潮流計算には干潟に対応した既存の２次元単層モデルを用いる。また、波浪の計算結果から既存の計算式を用いて海浜流を計算し、潮流との合成流を求める。
【００５０】
＜１−▲２▼＞次いで、波浪の計算結果と流れの場の計算結果を入力条件として、既存の底泥浮上量予測計算モデル（W.Bijker, Sedimentation is Channels and Treaanches., conf.coastal Eng. pp1708-1718, 1980を参照）により巻上げ懸濁物質量を予測計算する。
【００５１】
＜１−▲３▼＞そして、この巻上げ懸濁物質量の予測計算値を、予測対象海域におけるフィールド調査水質データを基に作成した巻上げ懸濁物質濃度と有機物濃度との関係式により、巻上げ有機物量に換算する。
【００５２】
以上から、巻上げ有機物量の正確な予測計算が実現される。
【００５３】
＜２．水質の予測計算＞
水質予測計算モデルとしては、有機物の生成・分解過程をモデル化した既存の生態系（低次）モデルタイプを用いる。干潟部等砂浜域において重要な泥系についても既存の生態系（低次）モデルタイプを採用する。なおここでの水質予測計算には底質の予測計算も含まれるものとする。
【００５４】
そして、この水質予測計算モデルに、底泥からの有機物の巻上げ過程モデル、つまり巻上げ有機物量の予測計算モデルを組み込み、この結合モデルに基づいて水質を予測計算する。
【００５５】
以上から、有機物の巻上げ量のみを考慮した高精度な水質予測計算が実現される。
【００５６】
もちろん、この発明は以上の例に限定されるものではなく、細部については様々な態様が可能である。
【００５７】
【発明の効果】
以上詳しく説明したとおり、この出願の発明によって、底泥から巻き上げられる有機物の量を予測計算することのできる、新しい巻上げ有機物量予測計算方法、および、有機物生成・分解過程の予測計算モデルによる水質予測計算をより精度良く行うことのできる、新しい水質予測計算方法が提供される。
【図面の簡単な説明】
【図１】底泥からの懸濁物質の巻上げ・沈降過程を例示した概念図である。
【図２】巻上げ有機物量の予測計算モデルを例示した概念図である。
【図３】一実施例としての予測対象海域におけるフィールド調査地点を示した模式図である。
【図４】一実施例としてのクロロフィル量と懸濁物質中の有機物との関係を示した図である。
【図５】一実施例としての巻上げ懸濁物質と有機成分との関係を示した図であり、（ａ）は７月調査データ、（ｂ）は９月調査データ、（ｃ）は１０月調査データのものである。
【図６】一実施例としてのこの出願の発明による巻上げ有機物量および水質の予測計算を説明するためのフローである。
【図１】
0004115109.tif 000004
【図２】
0004115109.tif 000005
【図３】
0004115109.tif 000006
【図４】
0004115109.tif 000007
【図５】
0004115109.tif 000008
【図６】
0004115109.tif 000009