"土壌等に含まれる自然由来の重金属元素等の水系への溶出量を長期的に推定するための推定方法であって、前記土壌等から採取された試料の化学分析などから前記水系の少なくともpH情報を含む水質変化をシミュレートする水質変化シミュレーションと併せて、前記化学分析から得られた前記重金属元素等の含有情報から前記水質変化に対応させて前記重金属元素等の各々の理想溶存量及び水酸化鉄の沈殿量を求めるステップと、前記沈殿量に対応して吸着される前記重金属元素等の吸着量を求めるステップと、前記理想溶存量から前記吸着量を減じて前記重金属元素等の前記溶出量を推定するステップと、を含むことを特徴とする自然由来の重金属元素等の溶出量推定方法。前記土壌等に粘土鉱物が含まれる場合において、前記水質変化に対応させて前記粘土鉱物の沈殿量を求めこれに対応させて吸着される前記重金属元素等の吸着量を前記理想溶存量から更に減じるステップを含むことを特徴とする請求項1記載の溶出量推定方法。前記含有情報は前記重金属元素等の元素種であることを特徴とする請求項1又は2に記載の溶出量推定方法。前記含有情報は前記重金属元素等の元素種毎の含有量であることを特徴とする請求項3記載の溶出量推定方法。前記理想溶存量は、前記pH情報に対応させて求められることを特徴とする請求項4記載の溶出量推定方法。前記理想溶存量は、前記水質変化シミュレーションの中で行われることを特徴とする請求項3記載の溶出量推定方法。前記重金属元素等は、少なくともヒ素、鉛、セレンの1つ以上を含むことを特徴とする請求項1乃至6のうちの1つに記載の溶出量推定方法。非特許文献1に挙げられた地球化学用シミュレータ以外にもいくつかの水質変化シミュレーションが知られており、土壌等の化学分析などの実測値から水系の少なくともpH情報を含む水質変化を求めることができる。かかるシミュレーションでは、鉱物種の水系への溶出量を求めるのと同様に、重金属元素等の情報を入力することで、理想的な重金属元素等の水系への溶出量についても求め得る。本発明は、上記したような状況に鑑みてなされたものであって、本発明の目的とするところは、土壌及び/又は岩石からの自然由来の重金属元素等の水系への溶出量を推定するための推定方法において、水質変化シミュレータと併せて重金属元素等の水系への溶出量を長期的に推定するための推定方法を提供することである。本発明による推定方法は、土壌等に含まれる自然由来の重金属元素等の水系への溶出量を長期的に推定するための推定方法であって、前記土壌等から採取された試料の化学分析などから前記水系の少なくともpH情報を含む水質変化をシミュレートする水質変化シミュレーションと併せて、前記化学分析から得られた前記重金属元素等の含有情報から前記水質変化に対応させて前記重金属元素等の各々の理想溶存量及び水酸化鉄の沈殿量を求めるステップと、前記沈殿量に対応して吸着される前記重金属元素等の吸着量を求めるステップと、前記理想溶存量から前記吸着量を減じて前記重金属元素等の前記溶出量を推定するステップと、を含むことを特徴とする。かかる発明によれば、土壌等からの自然由来の重金属元素等の各々の理想溶存量を補正して該重金属元素等の水系への溶出量を長期的に推定できるのである。すなわち、本発明者は、重金属元素等が水中へ溶出する場合、水系に存在する水酸化鉄のような沈殿物によって重金属元素等が吸着され、残部だけが水中に溶出して存在することに注目し、本発明に至ったのである。上記した発明において、前記土壌等に粘土鉱物が含まれる場合において、前記水質変化に対応させて前記粘土鉱物の沈殿量を求めこれに対応させて吸着される前記重金属元素等の吸着量を前記理想溶存量から更に減じるステップを含むことを特徴としてもよい。かかる発明によれば、粘土鉱物を含む土壌等であっても、土壌等からの自然由来の重金属元素等の各々の理想溶存量を補正して該重金属元素等の水系への溶出量を長期的に高精度で推定できるのである。上記した発明において、前記含有情報は前記重金属元素等の元素種であることを特徴としてもよい。かかる発明によれば、簡便に重金属元素等の水系への溶出量を長期的に推定できるのである。上記した発明において、前記含有情報は前記重金属元素等の元素種毎の含有量であることを特徴としてもよい。また、前記理想溶存量は、前記pH情報に対応させて求められることを特徴としてもよい。かかる発明によれば、簡便に重金属元素等の水系への溶出量を長期的に推定できるのである。上記した発明において、前記理想溶存量は、前記水質変化シミュレーションの中で行われることを特徴としてもよい。かかる発明によれば、より正確に重金属元素等の水系への溶出量を長期的に推定できるのである。上記した発明において、前記重金属元素等は、少なくともヒ素、鉛、セレンの1つ以上を含むことを特徴としてもよい。かかる発明によれば、より正確に重金属元素等の水系への溶出量を長期的に推定できるのである。本発明による土壌等に含まれる自然由来の重金属元素等の水系への溶出量、特に、ヒ素、鉛、亜鉛の水系への溶出量について長期的に推定する推定方法について図1乃至3を用いて説明する。まず、該推定方法の概略を図1に沿って説明する。掘削工事の施工予定地域の土壌及び/又は岩石をサンプル採取しS1、これに含まれる元素についての化学分析を行うS2。かかる化学分析の結果から、公知の水質変化シミュレーション、例えば、背景技術の欄で述べたPHREEQCのような地球化学用シミュレータによって、水質のpH変化や、少なくともFeイオンを含む各種化学種の濃度を求めるS3。この水質のpH及びFeイオン濃度から水酸化鉄の沈殿量を求めS4、この沈殿量に対応して吸着され得る重金属元素等の吸着量を求めるS5。一方、化学分析により検出された重金属元素等、特に、ヒ素、鉛、亜鉛といった元素の量と、上記した水質変化シミュレーションによって得られたpHの値から、この水質で溶出し得る理想溶存量を求めるS7。若しくは、重金属元素等の元素量比から上記した水質変化シミュレーションにより理想溶存量を求める。そして、重金属元素等の理想溶存量から吸着量を減じることで溶出量を求め得るのである。ここで、土壌等に粘土鉱物が含まれる場合には、粘土鉱物の沈殿量を水酸化鉄の沈殿量とともに求め、これに吸着され得る重金属元素等の吸着量を求め、理想溶存量から更に減じて溶出量を求めるS8。次に、図1に沿って、各ステップの処理について、その詳細を説明する。掘削工事の施工予定地域から土壌及び/又は岩石をサンプル採取しS1、破砕粉砕して粉体試料にする。粉体試料は蛍光X線分析法や湿式分析法などで全岩化学分析され、土壌等の中に含まれる元素を特定するS2。つまり、後述する地殻を構成する主要元素や、As,PbSeなどの重金属元素等の各元素の含有量比を求め得る。ここで、塩酸のような所定の溶媒を用いることで全岩化学分析の結果から、土壌等の中に含まれる重金属元素等を含む各元素の含有量を求めることも出来得る。これに関しては後述する。次に、上記した公知のシミュレーションによれば、全岩化学分析S2の結果から、少なくとも水系のpH及び含まれ得る元素量などの水質情報を得ることができるS3。かかるシミュレーション方法の一例について以下に述べる。例えば、上記した全岩化学分析により、地殻を構成する主要元素としてのFe,Mn,Ti,Ca,K,P,Si,Al,Mg,Naの存在比率、微量元素としてAs,PbSeなどの存在比率が測定されたとする。かかる存在比率から、各主要元素に対応する酸化物Fe2O3,FeO,MnO,TiO2,CaO,K2O,P2O5,SiO2,Al2O3,MgO,Na2Oの比率に換算するS31。次に、主要元素の元素量比に基づいて複数のノルム鉱物に系統的に分配するノルム計算によって、土壌中のノルム鉱物の存在量比を決定するS32。土壌中のノルム鉱物の水への溶出については、図3に示すよう、不可逆的な溶解及び沈殿の平衡S41と、可逆的な溶解及び沈殿の平衡S42と、の2つを考慮され得る。ノルム鉱物の水への不可逆的な溶解及び沈殿の平衡S41では、ノルム鉱物A、B、C、D、...について、nAAnBB...=nCCnDD...式1の反応を生じるとする。各ノルム鉱物の含有量A,B,C,D,...について、平衡定数K’は、K’=CnCDnD/AnABnB式2と表される。つまり、式2で示すように、平衡定数K’は、ノルム鉱物の存在量比から計算できる。一方、ノルム鉱物の水への可逆的な溶解及び沈殿の平衡S42では、飽和度指数SIについて、平衡溶解度積Kspを用いて、式2の平衡定数K’から、SI=LogK’/Ksp式3と定義される。つまり、式3から飽和度指数SIが正ならば沈殿反応が生じ、飽和度指数SIが負ならば溶解反応が生じるのである。なお、飽和度指数SIが零ならば平衡状態である。各ノルム鉱物の水への溶解及び沈殿の平衡S43について、反応速度Rkmols1kgw1は、Rk=rkA0mkn/Vm0kn式4で表される。ここで、rkmolm2s1を定数、A0m2を各ノルム鉱物の初期の表面積、mkmolを特定時間におけるノルム鉱物のモル数、Vを溶液量kgw、m0kmolをノルム鉱物の初期モル数とする。更に、nは各ノルム鉱物の形状因子定数である。例えば、ノルム鉱物の形状が球又は立方体ならば、このnは2/3である。各ノルム鉱物の水への溶解及び沈殿の平衡における反応速度Rkから所定時間後の水中に存在する全元素イオンの含有量を求め得る。また、これらから水のpHなどの水質情報を求め得るS44。また、pHに対応して、重金属元素等を含むノルム鉱物の溶解度積などから重金属元素等の濃度Ccalcmについても求め得る。以上は、1つのシミュレーションの方法について述べたが、少なくとも水のpH及び地殻を構成する主要元素の1つであるFeイオンの水中での含有量を求めることができるシミュレーションであれば、本発明のために用い得る。ところで、上記したシミュレーションの反応速度Rkから得られる重金属元素等の水中での濃度Ccalcmは、実測される値とは差異を生じる。また、pHから簡易的に重金属元素等の水中での濃度Ccalcmを求め得るが、やはり差異を生じる。つまり、かかるシミュレーションによる水中での濃度Ccalcmは理想溶出量である。この原因の1つとして、上記した地殻を構成する主要元素の1つであるFeイオンや、粘土質鉱物が重金属元素等、特に、As,PbSeなどを吸着することにある。水に溶出したこれらの重金属元素等のうちAsなどは、オキソ酸陰イオンとなるが、これは水酸基と酸素を共有して結合し易い。このような理由により、水酸基を有する水酸化鉄や粘土質鉱物は重金属元素等を吸着すると考えられる。重金属元素等の吸着の補正は、図1に示すように、吸着物質である沈殿物の量を求めS4、かかる沈殿物に吸着され得る重金属元素等の量Csorpmを求めS5、理想溶出量Ccalcmからこれを減じるS8のである。まず、吸着物質の水酸化鉄については、シミュレーションでの水質情報のpH値からこの水中に存在し得る理想のFeイオン量Fe3idealについて、図4に示すような水酸化鉄のイオン反応における溶解度積から、Fe3ideal=10-37.2/10pH-14式5と表すことが出来る。つまり、水酸化鉄の沈殿量FeOH3Sは、シミュレーションでの水質情報のFe3calcとともに、FeOH3S=Fe3calc-Fe3ideal=Fe3calc10-37.2/10pH-14式6のように表される。重金属元素等の水酸化鉄への吸着量Csorpmは、水酸化鉄の沈殿量FeOH3SとpHとの関数であるから、Csorp,m=fFeOH3S、pH式7のように表される。かかる重金属元素等の水酸化鉄への吸着量Csorpmは、実験的に得ることが出来る。また、フロイドリッヒの式によれば、aを定数、nを1以上の定数とすると、Csorp,m=apH×FeOH3S1/n式8のように簡便にも表し得る。また、分配係数による方法や、ラングミュアの式による方法など、公知の経験式を用い得る。つまり、重金属元素等の水酸化鉄への吸着量Csorpmと、シミュレーションで得られる重金属元素等の濃度Ccalcmとから、重金属元素等の溶出量Caqmは、Caqm=Ccalcm-Csorpm式9と補正され得る。ここで、重金属元素等の濃度Ccalcmは、全岩化学分析により重金属元素等の含有量を求め、シミュレーションから得られたpHの値を用いて求めることも出来得る。なお、水酸化鉄以外の吸着物質についても、例えば、土壌等が粘土鉱物を含む堆積岩からなる場合などには、粘土鉱物が吸着物質となる。かかる場合にあっても、上記同様に粘土鉱物の沈殿量を求めて、重金属元素等の吸着量を得て補正を行っても良い。上記した実施例によれば、土壌等からの自然由来の重金属元素等の各々の理想溶存量を補正して該重金属元素等の水系への溶出量を長期的に正確に推定できる。以上、本発明による実施例及びこれに基づく変形例を説明したが、本発明は必ずしもこれに限定されるものではなく、当業者であれば、本発明の主旨又は添付した特許請求の範囲を逸脱することなく、様々な代替実施例及び改変例を見出すことができるであろう。"