JP6980242B1

JP6980242B1 - 汚染物質流出による環境影響の予測方法、及び汚染物質流出による環境影響の予測プログラム

Info

Publication number: JP6980242B1
Application number: JP2021047210A
Authority: JP
Inventors: 宏史肴倉; 裕之横山
Original assignee: National Institute for Environmental Studies
Current assignee: National Institute for Environmental Studies
Priority date: 2021-03-22
Filing date: 2021-03-22
Publication date: 2021-12-15
Anticipated expiration: 2041-03-22
Also published as: JP2022146314A

Abstract

【課題】所定の汚染源からの汚染物質の流出による環境影響を好適に予測する技術を提供する。【解決手段】本開示の汚染物質流出による環境影響の予測方法は、コンピュータが、解析格子１０、２０において、固相解析領域１１、２１と液相解析領域１２、２２とを設定することと、固相解析領域１１、２１について、汚染物質１の拡散を計算することと、解析格子１０、２０において、固相と液相との間の汚染物質の分配を、固液界面の平衡モデルを用いて計算することと、液相解析領域１２、２２について、単位格子間の汚染物質の移行を、移流分散解析モデルを用いて計算することと、を実行する。【選択図】図２

Description

本発明は、汚染物質流出による環境影響の予測方法、及び汚染物質流出による環境影響の予測プログラムに関する。

近年、工場跡地の再開発などにおいて実施される土壌調査により、土壌汚染、またはそれに伴う地下水汚染が顕在化する事例が増えている。そして、汚染土に対して、従来から各種技術を利用した浄化工事が実施されている。しかしながら、例えば、浄化工事として汚染土の掘削除去が実施される場合、掘削された汚染土の処分が新たな問題となり得る。

また、土木工事において、資源の有効活用の観点から、リサイクル材料が盛土材や埋土材として利用される場合がある。しかしながら、例えば、鉄鋼スラグ等のリサイクル材料にはアルカリ性物質が含まれていることがあり、雨水等によりリサイクル材料からアルカリが流出することで、周囲の環境を汚染してしまう事態が生じ得る。

ここで、特許文献１には、アルカリ性資材を陸上現場で使用する際に、該陸上現場でのアルカリ流出リスクを予め予測するための技術が開示されている。当該技術によれば、陸上現場の土壌のアルカリ吸着を考慮した移流分散解析手法において、陸上現場の土壌の最大アルカリ吸着能を測定するだけの簡便な方法で各種定数を予測し、陸上現場におけるアルカリ流出リスクを予測することができる。

特開２０１７−６６８２７号公報

地球規模での資源枯渇の問題から、鉄鋼スラグ等のリサイクル材料や、自然由来に存在する重金属を含んだ自然由来汚染土、または人為的な汚染土として掘削除去された後に洗浄された土などを、土木工事等に有効活用するニーズは大きい。しかしながら、これらの再利用にあたって、周囲の環境を汚染してしまうリスクが問題となり、結局のところ廃棄されてしまう場合がある。ここで、特許文献１に記載の技術によれば、陸上現場の土壌のアルカリ吸着を考慮した移流分散解析手法により、陸上現場におけるアルカリ流出リスクが予測される。これによれば、アルカリ性物質を含む鉄鋼スラグ等のリサイクル材料から流出するアルカリによる環境汚染のリスクを事前に把握できるようにも思われる。しかしながら、リサイクル材料からの汚染物質の溶出メカニズムまでは十分に説明されておらず、また、現場土壌の実測値を用いるなどの経験則によるところもあり、汚染源からの汚染物質の流出による環境影響を汎用的に評価できるとは言い難い。

そして、従来技術によれば、汚染源からの汚染物質の流出による環境影響について、該汚染物質によって汚染される被汚染側の固相表面での吸着平衡を考慮した移流分散解析モデルで評価される場合が多く、被汚染側の固相表面での汚染物質の挙動に加えて被汚染側の固相内部での汚染物質の挙動まで考慮したモデルは見られない。例えば、従来技術の解析モデルによれば、汚染源からの汚染物質の流入が停止した後においては、その汚染土壌の状態が、１００年前から汚染されている土壌と、１日前に汚染された土壌と、で同一に扱われてしまう問題がある。実際の現象としては、これら土壌を比較すると、その汚染の程度が異なることがわかっているため、従来技術の解析モデルでは、汚染源からの汚染物質の流出による環境影響を適切に評価できないことになる。そして、このような環境影響を事前に詳細に予測することができれば、上記のリサイクル材料や自然由来汚染土などを、土木工事等に適切に有効活用することができる。このように、所定の汚染源からの汚染物質の流出による環境影響を適切に評価したうえで資源の有効活用を図る技術については、未だ改善の余地を残すものである。

本開示の目的は、所定の汚染源からの汚染物質の流出による環境影響を好適に予測する技術を提供することにある。

本開示は、所定の汚染源からの汚染物質の流出による環境影響を予測する方法であって、コンピュータが、解析格子において、固相を模擬した固相解析領域と、該固相と接触する液相を模擬した液相解析領域と、を設定することと、前記固相解析領域について、前記固相の内部での前記汚染物質の拡散を計算することと、前記解析格子において、前記固相と前記液相との間の前記汚染物質の分配を、固液界面の平衡モデルを用いて計算することと、前記液相解析領域について、単位格子間の前記汚染物質の移行を、移流分散解析モデルを用いて計算することと、を実行し、前記移流分散解析モデルには、前記固相の全体の体積に基づいて算出される前記固相からの溶出量が考慮される。

上記の汚染物質流出による環境影響の予測方法では、固相の内部での汚染物質の拡散（固相内拡散）と、固液界面での汚染物質の吸脱着（吸脱着平衡）と、を接続して、所定の汚染源からの汚染物質の流出による環境影響を予測することができる。なお、上記の汚染源とは、例えば、鉄鋼スラグ等のリサイクル材料や重金属を含んだ汚染土などである。ここで、固相内拡散は、例えばフィックの第１法則によって計算され、算出された固相濃度が、固液界面でヘンリー型、ラングミュアー型、またはフロイントリッヒ型の平衡モデルを用いて平衡分配される。そして、このように固相内拡散と吸脱着平衡とが接続された本開示の解析モデルによれば、汚染物質の経時的な濃度変化を合理的に説明することができ、以て、汚染源からの汚染物質の流出による環境影響を適切に評価することができる。また、本開示によれば、汚染源側と被汚染側とで同じ解析モデルを用いて汚染物質の挙動を解析することができるため、汚染源からの汚染物質の流出による環境影響を汎用的に評価できる。

そして、上記の汚染物質流出による環境影響の予測方法において、コンピュータが、前記液相解析領域について、前記固相との境界に位置する領域である液相境界層領域と、それ以外の領域である一般液相領域と、に分割することと、前記一般液相領域について、単位格子間の前記汚染物質の移行を、前記移流分散解析モデルを用いて計算することと、前記液相境界層領域について、前記汚染物質の拡散を計算することと、を更に実行してもよい。

また、上記の汚染物質流出による環境影響の予測方法において、コンピュータが、前記固相解析領域について、前記液相との境界に位置する領域に所定の微小な厚さの吸着層を設定することと、前記固相と前記液相との間の前記汚染物質の分配について、前記吸着層と前記液相解析領域との間で前記汚染物質の平衡分配を計算することと、を更に実行してもよい。

また、本開示は、汚染物質流出による環境影響の予測プログラムの側面から捉えることができる。すなわち、本開示は、所定の汚染源からの汚染物質の流出による環境影響を予測するプログラムであって、コンピュータに、解析格子において、固相を模擬した固相解析領域と、該固相と接触する液相を模擬した液相解析領域と、を設定することと、前記固相解析領域について、前記固相の内部での前記汚染物質の拡散を計算することと、前記解析格子において、前記固相と前記液相との間の前記汚染物質の分配を、固液界面の平衡モデルを用いて計算することと、前記液相解析領域について、単位格子間の前記汚染物質の移行を、移流分散解析モデルを用いて計算することと、を実行させ、前記移流分散解析モデルには、前記固相の全体の体積に基づいて算出される前記固相からの溶出量が考慮される。

本開示によれば、所定の汚染源からの汚染物質の流出による環境影響を好適に予測することができる。

第１実施形態における、汚染物質の移行メカニズムを模式的に説明するための第１の図である。第１実施形態における、汚染物質の移行メカニズムを模式的に説明するための第２の図である。第１実施形態に係る汚染物質流出による環境影響の予測方法の処理フローを示すフローチャートである。ユーザによってユーザ端末に入力される解析条件を説明するための図である。第２実施形態における、汚染物質の移行メカニズムを模式的に説明するための図である。第２実施形態に係る汚染物質流出による環境影響の予測方法の処理フローを示すフローチャートである。

以下、図面に基づいて、本開示の実施の形態を説明する。以下の実施形態の構成は例示であり、本開示は実施形態の構成に限定されない。

（ハードウェア構成）
本開示の汚染物質流出による環境影響の予測方法を実現するハードウェア資源の構成について説明する。本開示の汚染物質流出による環境影響の予測方法は、サーバ３００によって実行される。

サーバ３００は、データの取得、生成、更新等の演算処理及び加工処理のための処理能力のあるコンピュータ機器であればどの様な電子機器でもよく、例えば、パーソナルコンピュータ、サーバ、メインフレーム、その他電子機器であってもよい。すなわち、サーバ３００は、ＣＰＵやＧＰＵ等のプロセッサ、ＲＡＭやＲＯＭ等の主記憶装置、ＥＰＲＯＭ、ハードディスクドライブ、リムーバブルメディア等の補助記憶装置を有するコンピュータとして構成することができる。なお、リムーバブルメディアは、例えば、ＵＳＢメモリ、あるいは、ＣＤやＤＶＤのようなディスク記録媒体であってもよい。補助記憶装置には、オペレーティングシステム（ＯＳ）、各種プログラム、各種テーブル等が格納されている。

サーバ３００は、機能部として通信部３０１、記憶部３０２、制御部３０３を有しており、補助記憶装置に格納されたプログラムを主記憶装置の作業領域にロードして実行し、プログラムの実行を通じて各機能部等が制御されることによって、各機能部における所定の目的に合致した各機能を実現することができる。ただし、一部または全部の機能はＡＳＩＣやＦＰＧＡのようなハードウェア回路によって実現されてもよい。

ここで、通信部３０１は、サーバ３００をネットワークに接続するための通信インタフェースである。通信部３０１は、例えば、ネットワークインタフェースボードや、無線通信のための無線通信回路を含んで構成される。サーバ３００は、通信部３０１を介して、その他の外部装置と通信可能に接続される。

記憶部３０２は、主記憶装置と補助記憶装置を含んで構成される。主記憶装置は、制御部３０３によって実行されるプログラムや、当該制御プログラムが利用するデータが展開されるメモリである。補助記憶装置は、制御部３０３において実行されるプログラムや、当該制御プログラムが利用するデータが記憶される装置である。

制御部３０３は、サーバ３００が行う制御を司る機能部である。制御部３０３は、ＣＰＵなどの演算処理装置によって実現することができる。

なお、サーバ３００は、ネットワーク２００を介して所定のユーザ端末４００から操作されてもよい。ここで、ユーザ端末４００は、本開示の処理を実行するユーザが保有する電子機器であればよく、例えば、スマートフォン、タブレット端末、ウェアラブル端末、パーソナルコンピュータ等、その他端末機器であってもよい。

ユーザ端末４００は、機能部として通信部４０１、入出力部４０２、記憶部４０３を有している。通信部４０１は、ユーザ端末４００をネットワーク２００に接続するための通信インタフェースであり、例えば、ネットワークインタフェースボードや、無線通信のための無線通信回路を含んで構成される。入出力部４０２は、通信部４０１を介して外部から送信されてきた情報等を表示させたり、通信部４０１を介して外部に情報を送信する際に当該情報を入力したりするための機能部である。記憶部４０３は、サーバ３００の記憶部３０２と同様に主記憶装置と補助記憶装置を含んで構成される。

入出力部４０２は、更に、表示部４０２１、操作入力部４０２２を有している。表示部４０２１は、各種情報を表示する機能を有し、例えば、ＬＣＤ（Liquid Crystal Display）ディスプレイ、ＬＥＤ（Light Emitting Diode）ディスプレイ、ＯＬＥＤ（Organic Light Emitting Diode）ディスプレイ等により実現される。操作入力部４０２２は、ユーザからの操作入力を受け付ける機能を有し、具体的には、タッチパネル等のソフトキーあるいはハードキーにより実現される。

また、ネットワーク２００は、例えば、ＩＰネットワークである。ネットワーク２００は、ＩＰネットワークであれば、無線であっても有線であっても無線と有線の組み合わせであってもよく、例えば、無線による通信であれば、ユーザ端末４００は、無線ＬＡＮアクセスポイント（不図示）にアクセスし、ＬＡＮやＷＡＮを介してサーバ３００と通信してもよい。また、ネットワーク２００は、これらの例に限られず、例えば、公衆交換電話網や、光回線、ＡＤＳＬ回線、衛星通信網などであってもよい。

（ソフトウェア構成）
次に、本開示の汚染物質流出による環境影響の予測方法を実現するソフトウェア構成について、第１実施形態、第２実施形態、および第３実施形態に基づいて説明する。

＜第１実施形態＞
本実施形態における汚染物質流出による環境影響の予測方法の概要について、図１および図２に基づいて説明する。

図１は、本実施形態における、汚染物質の移行メカニズムを模式的に説明するための第１の図である。本実施形態は、所定の汚染源からの汚染物質の移行の態様を説明するものである。ここで、上記の汚染源とは、例えば、鉄鋼スラグ等のリサイクル材料や、自然由来に存在する重金属を含んだ自然由来汚染土、または人為的な汚染土として掘削除去された後に洗浄された土などである。そして、このような汚染源に含まれる重金属等の汚染物質が、雨水等により溶出し隣接する被汚染側の地盤に流入する。

ここで、図１（ａ）に示す汚染源および被汚染地盤の部分拡大図を表す図１（ｂ）に示すように、汚染源側および被汚染側の夫々には、固相と液相とが含まれることになる。そして、本実施形態では、これら固相と液相とを含んだ微小体積に解析格子が設定される。

図２は、本実施形態における、汚染物質の移行メカニズムを模式的に説明するための第２の図である。図２では、２つの単位格子を接続した例として、汚染源側から被汚染側への汚染物質の移行メカニズムを表している。図２に示すように、所定の汚染源からの汚染物質１が流出する汚染源側の解析格子１０においては、所定の拡散係数が設定されることで固相（汚染源）の内部の構造を模擬した固相解析領域１１と、該固相と接触する液相（雨水等）を模擬した液相解析領域１２と、が設定される。また、汚染物質１が流入する被汚染側の解析格子２０においては、所定の拡散係数が設定されることで固相（被汚染物質である土壌）の内部の構造を模擬した固相解析領域２１と、該固相と接触する液相（雨水等）を模擬した液相解析領域２２と、が設定される。このように、本実施形態では、汚染源側と被汚染側において、同様にして解析領域が設定される。

ここで、固相解析領域１１および固相解析領域２１では、汚染物質１が拡散するとしてモデル化される。また、固相解析領域１１と液相解析領域１２との界面である固液界面１３、および固相解析領域２１と液相解析領域２２との界面である固液界面２３では、固相と液相との間の汚染物質１の吸脱着がモデル化される。

また、液相解析領域１２および液相解析領域２２について、単位格子間の汚染物質１の移行がモデル化される。そうすると、図２に示すメカニズムでは、汚染源側から被汚染側へ汚染物質１が移行することになる。なお、図２における線Ｌ１０および線Ｌ２０は、汚染物質１の量を模式的に表していて、汚染源側の固相解析領域１１に含まれる汚染物質１の一部が、汚染源側の液相解析領域１２に溶出している。また、汚染源側の液相解析領域１２に含まれる汚染物質１の一部が、被汚染側の液相解析領域２２に移行し、該液相解析領域２２に含まれる汚染物質１の一部が、被汚染側の固相解析領域２１に滲入している。

次に、本実施形態における汚染物質流出による環境影響の予測方法の処理フローについて、図３に基づいて説明する。図３は、本実施形態に係る汚染物質流出による環境影響の予測方法の処理フローを示すフローチャートである。本フローは、サーバ３００によって実行される。

本フローでは、先ず、Ｓ１０１において、解析条件が取得される。本実施形態では、サーバ３００は、ユーザによってユーザ端末４００に入力された解析条件を通信により取得する。ここで、図４は、ユーザによってユーザ端末４００に入力される解析条件を説明するための図である。本開示の処理を実行するユーザは、ユーザ端末４００を用いて、図４に例示する画面ＳＣ１を表示部４０２１に表示させることができる。図４に例示する画面ＳＣ１には、解析領域ＳＣ１１、入力項目ＳＣ１２、入力フィールドＳＣ１３、および備考欄ＳＣ１４が示される。そして、ユーザは、ユーザ端末４００の操作入力部４０２２を用いて、入力フィールドＳＣ１３に解析条件を入力することができる。例えば、解析領域ＳＣ１１が固相で入力項目ＳＣ１２が質量の入力フィールドＳＣ１３には、固相の質量が入力される。また、本実施形態では、固相が複数の微小な球体または板であるとして解析される。そのため、解析領域ＳＣ１１が固相で入力項目ＳＣ１２が形状の入力フィールドＳＣ１３には、固相を複数の微小な球体として解析する場合には０が、固相を複数の微小な板として解析する場合には１が入力される。また、本実施形態では、固相と液相との間の汚染物質の吸脱着を計算する固液界面の平衡モデルとして、ヘンリー型、ラングミュアー型、またはフロイントリッヒ型の平衡モデルを利用できる。そのため、解析領域ＳＣ１１が固液界面で入力項目ＳＣ１２が平衡モードの入力フィールドＳＣ１３には、ヘンリー型を利用する場合には０が、ラングミュアー型を利用する場合には１が、フロイントリッヒ型を利用する場合には２が入力される。なお、このような解析条件は、解析格子夫々について、設定することができる。

次に、Ｓ１０２において、解析格子が設定される。本実施形態では、上記の図１および図２の説明で述べたように、全体の構成のうち固相と液相とを含んだ微小体積に解析格子が設定される。例えば、上記の図２に示したメカニズムでは、汚染源側の解析格子１０において、固相解析領域１１と液相解析領域１２とが設定され、被汚染側の解析格子２０において、固相解析領域２１と液相解析領域２２とが設定される。なお、格子分割数は適宜設定することができる。また、このような解析格子が設定される１つの微小固体と、固相全体と、の関係を下記式１で表すことができる。

φ：空隙率
Ｖ：全体の固相体積（空隙含む）
ｖ：１つの微小固体の体積
Ｎ：微小固体の個数
なお、全体の固相体積Ｖは、Ｓ１０１の処理で取得される固相の質量と密度とに基づいて、算出することができる。また、１つの微小固体の体積ｖは、Ｓ１０１の処理で取得される微小固体の形状と半径（又は板厚）とに基づいて、定義することができる。

次に、Ｓ１０３において、固相解析領域について、固相の内部での汚染物質の拡散（固相内拡散）が計算される。これは、フィックの第１法則によって計算され、固相解析領域における拡散方程式を下記式２に表す。

Ｆ：固相濃度
ρｓ：固相密度
Ｄｓｍ：固相内分子拡散係数
ここで、固相内分子拡散係数Ｄｓｍは、例えば、拡散溶出試験で得ることができる。

次に、Ｓ１０４において、固液界面での汚染物質の吸脱着（吸脱着平衡）が計算される。本実施形態では、上述したように、固液界面の平衡モデルとして、ヘンリー型、ラングミュアー型、またはフロイントリッヒ型の平衡モデルが利用され得る。そして、上記式２により算出される固相濃度Ｆが、固液界面でヘンリー型のモデルにより平衡分配される場合を下記式３に、ラングミュアー型のモデルにより平衡分配される場合を下記式４に、フロイントリッヒ型のモデルにより平衡分配される場合を下記式５に表す。

Ｃ：液相濃度
ｋｄ：ヘンリー型の分配係数

Ｃ：液相濃度
ｋＬ：ラングミュアー型の平衡定数
ｑｓａｔ：ラングミュアー型の飽和吸着量

Ｃ：液相濃度
ｋＦ：フロイントリッヒ型の平衡定数
ｎＦ：フロイントリッヒ型のべき数

次に、Ｓ１０５において、隣接する単位格子間の汚染物質の移行が計算される。本実施形態では、Ｓ１０３の処理により算出される液相濃度Ｃを用いて、下記式６に示す移流分散方程式によって上記の移行が計算される。

ｕ：ダルシー流速
Ｄｌ：液相の分散テンソル
Ｑ：固相からの溶出量
なお、ダルシー流速ｕは、液相を流れる流体の流量と、その通水面積と、に基づいて、算出することができる。また、固相からの溶出量Ｑは、微小体積における固相内拡散と、微小固体の個数Ｎと、に基づいて、算出される。また、液相の分散テンソルＤｌは、下記式７によって表される。

αＬ：縦分散長
αＴ：横分散長
Ｄｌｍ：液相内分子拡散係数
τ：屈曲率

図３に示すフローでは、Ｓ１０３〜Ｓ１０５の処理が実行されることで、固相の内部での汚染物質の拡散と固液界面での汚染物質の吸脱着、および隣接する単位格子間の汚染物質の移行が計算される。なお、Ｓ１０３〜Ｓ１０５の処理は、周知の技術に基づいて、所定の演算周期で繰り返し実行される。そして、本実施形態に係る汚染物質流出による環境影響の予測方法の処理フローでは、固相の初期濃度の設定により、汚染源側と被汚染側の解析格子を模擬することができ（例えば、上記の図２に示したメカニズムでは、固相の初期濃度が高く設定された解析格子が汚染源側の解析格子１０として、固相の初期濃度が低くまたは０に設定された解析格子が被汚染側の解析格子２０として模擬される。）、汚染源側と被汚染側とで同じ解析モデルを用いて汚染物質の挙動を解析することができる。そのため、汚染源からの汚染物質の流出による環境影響を汎用的に評価できる。また、本実施形態では、上記式２および上記式３〜５に表されるように、固相内拡散と吸脱着平衡とが接続されるため、汚染物質の経時的な濃度変化を合理的に説明することができる。

そして、本実施形態に係る汚染物質流出による環境影響の予測方法の処理フローでは、上記のＳ１０３〜Ｓ１０５に例示される処理による、固相の内部での汚染物質の拡散と固液界面での汚染物質の吸脱着、および隣接する単位格子間の汚染物質の移行の計算が、予め定められた回数繰り返し実行される。これにより、汚染物質の経時的な濃度変化を算出することができる。そうすると、例えば、上記の図２に示したメカニズムにおいて、液相を流れる流体の流量が通水開始から所定時間経過後には０になる場合を仮定すると、汚染源から被汚染側への汚染物質の移行が停止した後においては、従来技術の解析モデルによれば、固相と液相とが瞬時平衡となるため、汚染源側および被汚染側での汚染物質の濃度変化が生じないのに対して、本実施形態の固相内拡散と吸脱着平衡とが接続された解析モデルによれば、汚染源側および被汚染側での汚染物質の経時的な濃度変化が生じることになる。ここで、汚染源からの汚染物質の流入が停止した後の土壌について、１００年前から汚染されている土壌と、１日前に汚染された土壌と、では、その汚染の程度が異なることがわかっている。そして、本実施形態では、上記の如く汚染物質の経時的な濃度変化を算出することができるため、実際に汚染土から固相内拡散によって染み込んだ汚染物質が徐々に放出されるような事象を再現することができ、以て、汚染源からの汚染物質の流出による環境影響を適切に評価することができる。

以上に述べた処理によれば、環境影響を事前に詳細に予測できることで、循環資源の環境影響や長期的な地盤汚染の解析において、有用性が高い。これにより、リサイクル材料や自然由来汚染土などを、土木工事等に適切に有効活用することができる。

そして、本実施形態では、サーバ３００によって、汚染源となり得るリサイクル材料や汚染土などの利用可能性が判定されてもよい。従来技術によれば、汚染源からの汚染物質の流出による環境影響を適切に評価し難いため、これらの再利用にあたって、周囲の環境を汚染してしまうリスクが問題となり、結局のところ廃棄されてしまう場合がある。一方、本実施形態の固相内拡散と吸脱着平衡とが接続された解析モデルによれば、汚染物質の経時的な濃度変化を予測できるため、例えば、或る期間における汚染物質の濃度の変化を算出し、その結果と（土壌汚染対策法等により規定される）所定の閾値とを比較することで、上記の利用可能性を判定することができる。

詳しくは、サーバ３００は、例えば、上記の図２に示したメカニズムの汚染源側の解析格子１０および被汚染側の解析格子２０において、所定期間（例えば、１０年間）における汚染物質の濃度の経時的変化を算出する。このとき、流体が液相を流れるタイミングを任意に複数回設定してもよいし、その流量を任意に設定してもよい。そして、サーバ３００は、算出された汚染物質の濃度が所定の閾値よりも小さくなるタイミングを、所定の汚染源およびそれが埋土等された土壌の利用可能開始時期として判定する。なお、これは、汚染源から汚染物質が徐々に放出されるような事象や、汚染源からの汚染物質が染み込んだ土壌から該汚染物質が徐々に放出されるような事象を想定したものであって、このような場合、長い期間をかけて固相と液相とで汚染物質の濃度が平衡にいたることや、微量の汚染物質が地下水に混入し所定の基準値以下の汚染水として流出していくことで、全体としての汚染物質の濃度が低下する。したがって、全体としての汚染物質の濃度が所定の閾値よりも低下したタイミングを、上記の利用可能開始時期として判定することができる。

以上に述べたように、本実施形態によれば、所定の汚染源からの汚染物質の流出による環境影響を好適に予測することができる。

＜第２実施形態＞
第２実施形態における汚染物質流出による環境影響の予測方法について、図５および図６に基づいて説明する。

図５は、本実施形態における、汚染物質の移行メカニズムを模式的に説明するための図である。図５では、２つの単位格子を接続した例として、汚染源側から被汚染側への汚染物質の移行メカニズムを表している。本実施形態では、図５に示すように、液相解析領域１２および液相解析領域２２について、固相との境界に位置する領域である液相境界層領域１２ａ、２２ａと、それ以外の領域である一般液相領域１２ｂ、２２ｂと、に分割される。そして、固相解析領域１１と液相境界層領域１２ａとの界面である固液界面１３、および固相解析領域２１と液相境界層領域２２ａとの界面である固液界面２３において、固相と液相との間の汚染物質１の吸脱着がモデル化される。また、一般液相領域１２ｂおよび一般液相領域２２ｂについて、隣接する単位格子間の汚染物質１の移行がモデル化される。

そして、液相解析領域では、液相境界層領域１２ａ、２２ａにおいて、汚染物質１が拡散するとしてモデル化される。そうすると、汚染物質１の量を模式的に表す線Ｌ１０および線Ｌ２０に示されるように、液相境界層領域１２ａ、２２ａにおいて、拡散による汚染物質１の濃度変化が生じることになる。

次に、本実施形態における汚染物質流出による環境影響の予測方法の処理フローについて、図６に基づいて説明する。図６は、本実施形態に係る汚染物質流出による環境影響の予測方法の処理フローを示すフローチャートである。本フローは、サーバ３００によって実行される。なお、本実施形態において、上述した第１実施形態と実質的に同一の処理については、その詳細な説明を省略する。

図６に示すフローでは、Ｓ１０４の処理の後、Ｓ２０１において、液相境界層領域について、液相境界層での汚染物質の拡散が計算される。これは、フィックの第１法則によって計算され、液相境界層領域における拡散方程式を下記式８に表す。

Ｃｂ：液相境界層内の液相濃度
Ｄｌｍ：液相内分子拡散係数
そして、Ｓ２０１の処理の後、サーバ３００はＳ１０５の処理へ進む。なお、本フローにおいても、Ｓ１０３〜Ｓ１０５の処理が、周知の技術に基づいて、所定の演算周期で繰り返し実行される。

以上に述べた処理によっても、所定の汚染源からの汚染物質の流出による環境影響を好適に予測することができる。

＜第３実施形態＞
第３実施形態における汚染物質流出による環境影響の予測方法について説明する。上述した第１実施形態では、固相内拡散によって、固相内部での汚染物質の濃度の経時的変化を算出することができるものの、従来の固液界面での平衡分配においては、固相の全体の体積についての濃度に対する分配が計算され得る。そして、このような場合には、実現象と差異が生じ得ることが新たに見出された。

そこで、本実施形態では、固相解析領域について、液相との境界に位置する領域に所定の微小な厚さの吸着層が設定される。ここで、微小な厚さとは、例えば、０．００１ｍｍである。そして、本実施形態では、固相と液相との間の汚染物質の吸脱着について、上記の吸着層と液相解析領域との間で汚染物質の平衡分配が計算される。

詳しくは、サーバ３００は、式３、式４、または式５における固相濃度Ｆについて、固相全体に対する値を上記の吸着層の体積に基づく値に換算して計算する。つまり、固相内部での汚染物質の濃度の経時的変化が算出された固相濃度Ｆにおいて、微小体積の領域である吸着層に存在する汚染物質が平衡分配に関与するものとして、固相と液相との間の汚染物質の吸脱着が計算される。これにより、実現象を再現し易くなる。

そして、以上に述べた処理によっても、所定の汚染源からの汚染物質の流出による環境影響を好適に予測することができる。

＜その他の変形例＞
上記の実施形態はあくまでも一例であって、本開示はその要旨を逸脱しない範囲内で適宜変更して実施しうる。例えば、本開示において説明した処理や手段は、技術的な矛盾が生じない限りにおいて、自由に組み合わせて実施することができる。

また、上記の実施形態では、汚染源となり得るリサイクル材料や汚染土などの再利用を例にして説明したが、本開示の汚染物質流出による環境影響の予測方法が用いられる態様はこれに限定されない。例えば、本開示の汚染物質流出による環境影響の予測方法を用いて、工場等の人為的汚染源による環境影響を事前に評価することで、その工場跡地の再開発などが計画された場合に、浄化工事の必要性の有無や浄化工事の工法等を事前に検討することができる。

また、１つの装置が行うものとして説明した処理が、複数の演算処理装置によって分担して実行されてもよい。このとき、複数の演算処理装置同士が好適に協働可能に構成される。また、異なる装置が行うものとして説明した処理が、１つの装置によって実行されても構わない。コンピュータシステムにおいて、各機能をどのようなハードウェア構成（サーバ構成）によって実現するかは柔軟に変更可能である。

本開示は、上記の実施形態で説明した機能を実装したコンピュータプログラムをコンピュータに供給し、当該コンピュータが有する１つ以上のプロセッサがプログラムを読み出して実行することによっても実現可能である。このようなコンピュータプログラムは、コンピュータのシステムバスに接続可能な非一時的なコンピュータ可読記憶媒体によってコンピュータに提供されてもよいし、ネットワークを介してコンピュータに提供されてもよい。非一時的なコンピュータ可読記憶媒体は、例えば、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクドライブ（ＨＤＤ）等）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤディスク・ブルーレイディスク等）など任意のタイプのディスク、読み込み専用メモリ（ＲＯＭ）、ランダムアクセスメモリ（ＲＡＭ）、ＥＰＲＯＭ、ＥＥＰＲＯＭ、磁気カード、フラッシュメモリ、光学式カード、電子的命令を格納するために適した任意のタイプの媒体を含む。

１・・・・・・汚染物質
１０・・・・・汚染源側の解析格子
２０・・・・・被汚染側の解析格子
１１、２１・・固相解析領域
１２、２２・・液相解析領域
１３、２３・・固液界面

Claims

所定の汚染源からの汚染物質の流出による環境影響を予測する方法であって、
コンピュータが、
解析格子において、固相を模擬した固相解析領域と、該固相と接触する液相を模擬した液相解析領域と、を設定することと、
前記固相解析領域について、前記固相の内部での前記汚染物質の拡散を計算することと、
前記解析格子において、前記固相と前記液相との間の前記汚染物質の分配を、固液界面の平衡モデルを用いて計算することと、
前記液相解析領域について、単位格子間の前記汚染物質の移行を、移流分散解析モデルを用いて計算することと、
を実行し、
前記移流分散解析モデルには、前記固相の全体の体積に基づいて算出される前記固相からの溶出量が考慮される、
汚染物質流出による環境影響の予測方法。
請求項１に記載の汚染物質流出による環境影響の予測方法において、
コンピュータが、
前記液相解析領域について、前記固相との境界に位置する領域である液相境界層領域と、それ以外の領域である一般液相領域と、に分割することと、
前記一般液相領域について、単位格子間の前記汚染物質の移行を、前記移流分散解析モデルを用いて計算することと、
前記液相境界層領域について、前記汚染物質の拡散を計算することと、
を更に実行する。
請求項１又は請求項２に記載の汚染物質流出による環境影響の予測方法において、
コンピュータが、
前記固相解析領域について、前記液相との境界に位置する領域に所定の微小な厚さの吸着層を設定することと、
前記固相と前記液相との間の前記汚染物質の分配について、前記吸着層と前記液相解析領域との間で前記汚染物質の平衡分配を計算することと、
を更に実行する。
所定の汚染源からの汚染物質の流出による環境影響を予測するプログラムであって、
コンピュータに、
解析格子において、固相を模擬した固相解析領域と、該固相と接触する液相を模擬した液相解析領域と、を設定することと、
前記固相解析領域について、前記固相の内部での前記汚染物質の拡散を計算することと、
前記解析格子において、前記固相と前記液相との間の前記汚染物質の分配を、固液界面の平衡モデルを用いて計算することと、
前記液相解析領域について、単位格子間の前記汚染物質の移行を、移流分散解析モデルを用いて計算することと、
を実行させ、
前記移流分散解析モデルには、前記固相の全体の体積に基づいて算出される前記固相からの溶出量が考慮される、
汚染物質流出による環境影響の予測プログラム。