Evaluating the effectiveness of a geostatistical approach with groundwater flow modeling for three-dimensional estimation of a contaminant plume

3次元地下汚染分布評価に対する地下水流動を考慮した地球統計学的アプローチの有効性検証

高井 静霞 ; 島田 太郎 ; 武田 聖司 ; 小池 克明*

Takai, Shizuka; Shimada, Taro; Takeda, Seiji; Koike, Katsuaki*

放射性核種や化学物質により汚染された地下環境に対し、リスクを評価し除染対策を立てるためには、測定データから汚染物質の濃度分布とその不確かさを正確に推定する必要がある。汚染の放出時刻歴が明らかでない場合でも、これを時間方向の相関を考慮して測定濃度から逆解析することで、既知の汚染源から放出された汚染分布とその不確かさを地下水流動を考慮した地球統計学的手法によって推定することができる。しかし、従来の手法には3つの課題があった:(1)実際の汚染事例に対する3次元的な検証例がない、(2)推定値に制限がないため、負値の発生及び不確かさの増大につながりうる、(3)複数の汚染物質を対象とした検証例が少ない。これらを解決するために、本研究では地下水流動を考慮した地球統計学的手法に、ギブスサンプリングによる非負の制限を導入した。そして本手法を、Gloucester処分場(カナダ、オタワ)における地下汚染事例に適用した。評価対象は、水に可溶な3つの汚染物質(1,4-ジオキサン,テトラヒドロフラン,ジエチルエーテル)とした。その結果、1982年の測定濃度(66点)から推定した汚染分布は、全ての汚染物質に対し測定値と高い相関(相関係数R0.7)を示し、本手法の有効性が確認された。特に1,4-ジオキサンに対しては、1978年の実際の大規模投棄に応じた放出ピークが、最小エントロピー法による先行研究よりも正確に推定された。同様の放出ピークは他の汚染物質に対しても、有機炭素含有量からの遅延係数の推定範囲で概ね再現された。

When assessing the risk from an underground environment that is contaminated by radioactive nuclides and hazardous chemicals and planning for remediation, the contaminant plume distribution and the associated uncertainty from measured data should be estimated accurately. While the release history of the contaminant plume may be unknown, the extent of the plume caused by a known source and the associated uncertainty can be calculated inversely from the concentration data using a geostatistical method that accounts for the temporal correlation of its release history and groundwater flow modeling. However, the preceding geostatistical approaches have three drawbacks: (1) no applications of the three-dimensional plume estimation in real situations, (2) no constraints for the estimation of the plume distribution, which can yield negative concentration and large uncertainties, and (3) few applications to actual cases with multiple contaminants. To address these problems, the non-negativity constraint using Gibbs sampling was incorporated into the geostatistical method with groundwater flow modeling for contaminant plume estimation. This method was then tested on groundwater contamination in the Gloucester landfill in Ontario, Canada. The method was applied to three water soluble organic contaminants: 1,4-dioxane, tetrahydrofuran, and diethyl ether. The effectiveness of the proposed method was verified by the general agreement of the calculated plume distributions of the three contaminants with concentration data from 66 points in 1982 (linear correlation coefficient of about 0.7). In particular, the reproduced large spill of organic contaminants of 1,4-dioxane in 1978 was more accurate than the result of preceding minimum relative entropy-based studies. The same peak also appeared in the tetrahydrofuran and diethyl ether distributions approximately within the range of the retardation factor derived from the fraction of organic carbon.