この作品は、水質改善目標の指定された一連の農業保全慣行の最適な(最も低コストの)配置を解くために、進化的アルゴリズムを用いた最適化コンポーネントを使用した水質モデルの統合を示します。ソリューションは、トレードオフの明確な定量化を可能にする、マルチ客観的なアプローチを使用して生成されます。
Method Article
この作品は、水質改善目標の指定された一連の農業保全慣行の最適な(最も低コストの)配置を解くために、進化的アルゴリズムを用いた最適化コンポーネントを使用した水質モデルの統合を示します。ソリューションは、トレードオフの明確な定量化を可能にする、マルチ客観的なアプローチを使用して生成されます。
風景渡る特定の水質目標を達成するためのターゲティング保全実践投資のコスト効率の高い( すなわち 、最も低コストの)方法を見つけることは流域管理における最も重要である。流域コンテキスト( 例えば 、5,12,20)で最も低コストの解決策を見つけるの伝統的な経済学の方法は、外部への影響を正確に生成された敷地内の汚染の割合として記述することができることを想定しています。このようなアプローチは、汚染源の影響はしばしば複雑な生物物理学的プロセスによって決定される流域における実際の汚染プロセスを代表することはほとんどありません。現代の物理ベース、空間的に分散水文学的シミュレーションモデルの使用は、プロセス表現の面でリアリズムの大きい程度を可能にしますが、モデルは最適化の不可欠な一部となるシミュレーション最適化フレームワークの開発が必要である。
進化進アルゴリズムは流域シミュレーション·最適化問題のコンビナトリアル性質に対処するとフル水質モデルの使用を可能にすることが特に有用な最適化ツールであるように見えます。進化的アルゴリズムは、候補解として、流域の保全慣行の特定の空間的配分を扱い、最適化の目標に対して改善を見つけるために反復的に選択、交叉および突然変異の確率的演算子を適用する解候補の集合(集団)を利用する。この場合の最適化の目的は同時に保全慣行のコストを最小限に抑え、流域における非点、源汚染を最小にするためである。研究の最近の拡大とセットは同様の方法を使用しようとすると広義の進化的最適化手法3,4,9,10,13-15,17-19,22,23,25と水質モデルを統合しています。このアプリケーションでは、Rabotyagovらの適切に従うプログラムを実証ACHと保全の実践と保全の実践と、ユーザーが指定した水質のコストとの間の完全なトレードオフのフロンティアを検索するため、その費用の現代および一般SPEA2 26多目的進化アルゴリズムとSWATの水質モデル7を使用し、ユーザが指定したセットを統合目的。フロンティアは、様々な水質改善の目標に関連するコストの完全な範囲を提示することによって、流域管理者が直面してトレードオフを定量化する。プログラムは、指定された水質改善の目標と保全慣行の最適化された配置のマップの生産を達成する流域の構成を選択することができます。
1。流域モデルを準備して、最適化のための入力データを提供
候補解の合計コストは、流域単位( "水文応答ユニット"、またはHRUs)に適用保全慣行のコストの合計です。最適化プログラムは、流域内のすべての農地HRUの保全慣行の特定のセットから単一の保全の実践の最適割り当てを考慮します。のHRUに割り当てられている可能性保全の実践のETSは、対立遺伝子セットと呼ばれます。
2。最適化パラメータを選択する

図1は、最適化の目的とパラメータの設定。
最適化パラメータを選択する:
プリセット:最適化する流域を選択します。 "Apply"をクリックすると、この画面上のコントロールの値を埋めるためにプリセットファイル "流域presets.csv"からエントリを選択します。
出力変数:最適化のための環境目標を選択します。 (Nアウトレット、アウトレットP)を選択したように、これは3次元の目的関数を定義しています:窒素(有機N + NO3 + NH4 + NO2)はリン(オルガ、出口での5年間の平均NICのP +ミネラルP)は出口で5年間の平均化され、保全実務の総費用。これは3次元のトレードオフのフロンティアを作成することに注意してください。多目的プログラム({出力変数}、トータルコスト)を最小化するためにある別の出力変数は、選択することができます。
人口規模:設定初期人口規模。これは、解候補数の初期値を決定します。オプションの "各アレルによるシード"が選択された場合、流域内のすべての耕作地HRUsに設定アレルで指定された各保全実践の統一的適用を表す解候補が最初に作成されています。残りの候補解は耕作地HRUsに設定アレルから保全慣行の無作為割り当てによって作成されます。オプションの "各アレルによるシード"を選択した場合は、1は、初期集団のサイズが対立遺伝子セット内の対立遺伝子の数(このデモでは23)と少なくとも同じ大きさであることを確認する必要があります。
世代の > 数:最適化のラン(ランが再開できることに注意してください)の世代の希望の番号を(反復)を設定します。
交叉確率:2つの候補解が新たな解候補を作成するために選択された場合は、交叉確率は異なる新たな解が(このデモのために1に設定)作成されている確率を指定します。
一時的な人口の大きさ:これは、作成した新しい解候補の数を決定します。この値はプロセッサスレッドの数の整数倍(このデモで16)である場合にプロセッサのリソースを最も効率的に使用されています。
突然変異確率:対立遺伝子セットから別の省エネ実践にHRU割り当てのランダムな変化の確率を指定します。 (それはこのデモのために0.03に設定されています)。
数スレッド:使用プロセッサまたはスレッドの数を選択します。 16は、このデモで使用されています。
カーブ番号校正係数:これは、SWATのモデルのキャリブレーションから提供される。
テキストファイル内の人口を保存:これは、指定した数だけ繰り返された後に実行する最適化を再起動するように一つ願いが完了したかどうか選択することが重要です。このオプションをチェックすると、すべての存続候補解(個体)のすべてのHRUのアレルの値を指定したテキストファイルを生成します。これは実行を再開し、続行します後ろに読み取ることができます。
二次最適化パラメータ
最初の年は:過去の気象情報の開始年後に設定されており、遅くとも7年以上このデータの終わりの前にしなければなりません。
トウモロコシの価格:肥料削減のコストを見積もるために歩留まり損失の式と組み合わせて使用します。
採点方法:SPEA2アーカイブ。採点は存続個人がクロスオーバーするために選択されることになっている方法おそらく決定します。
パージ方法:支配した。すべて3次元で悪化している個人が支配し、パージされます。
HUCソース: "指定した場所"に設定し、次のフィールドから"7100006"の値を意味する"流域HUC"はHUCゾーンテーブルの行を検索するために使用されます。値は、 "07100006"はアライグマ流域のための8桁のHUCコードです。
コストソース:CRP以外の費用は上記の練習のコストテーブルに郡のFIPSコードによって決定されるであろうことを示すために、 "郡(HRUロケーション·コード)"に設定されています。
コストソースCRP:CRPのコストは上記の練習のコストテーブルに郡のFIPSコードによって決定されるであろうことを示すために、 "1場所"に設定してください。
3。代表的な結果
GeneticiSWAT.exeはすべての解候補(個人)だけでなく、最適化の実行を再開するために使用することができます最終的なアルゴリズムの反復との結果をコードする "保存"ファイルの設定と結果を示すログファイルを生成します。
この時点で、1次のステップに従うことで、パレート効率的なソリューション(トレードオフのフロンティア)のセット全体を視覚化することができます:

図2:Scが3次元可視化のためのフロンティア "スナップショット"を作成するためのreenshot。
出力は、POV-Rayのプログラムを使用して、次に "レンダリング"、 "ファイルキュー"を選択して画像ファイルに一度にレンダリングされたすべてのことができるファイルのシリーズです。画像は自分で使用したり、アルゴリズムの進行を示すムービーにまとめることができます。

図3はトレードオフのフロンティアの静的可視化。
必要であれば、アルゴリズムの進行を示すムービー "はFramescanner.exe"を実行し、次の手順を実行することによって作成できます。
フロンティア内の各ポイントは、流域構成(景観上の保全慣行の特定の割り当て)を表します。これらの構成のマップは、次の手順を実行して、全体のフロンティアとして見ることができる:

図4:最後のフロンティア内の個々のマップを作成するのスクリーンショット。
関心のある特定の流域構成(個人)をエクスポートします。
多くの場合、目的の質問は特定の水質目標を達成する具体的な流域の構成(個人)を選択することです。例えば、人は個人ラージヒルを見つけたがっているかもしれませんベースラインの負荷に対して20%の30%、リンによって窒素を低減フロンティアでidual。 MapSWATは1つが指定された目標に対する最小ユークリッド距離を使用して個々のためにフロンティアを検索することができます。これは、次の操作を実行して行うことができます。

図5:水質目標に基づくフロンティアで特定の個人を検索のスクリーンショット。

図6:検索出力のスクリーンショット
files/ftp_upload/4009/4009fig7.jpgの "alt ="図7 "/>
図7。フロンティアで選択された個々のサンプルを記述したマップのスクリーンショット。 拡大図を表示するには、ここをクリックしてください 。
さらに分析するために、地図データをエクスポートすると、以下の手順に従うことにより可能である:
| プログラムの名称 | ソース | 説明 |
| 回転子 | カード | |
| Swat2005GA.exe | 米国農務省草地、土壌·水研究所 | 流域シミュレーションモデル |
| i_SWAT.exe | カード | SWATのデータベースインターフェイス |
| GeneticISWAT.exe | カード | 進化アルゴリズムのSWATのコントローラ。からGALib組み込んhttp://lancet.mit.edu/ga/ 。 |
| MapSWAT.exe | カード | i_SWATデータベースとシェープファイルを読み込んで、世代や個人の画像を生成します。 |
| POV-Rayは | Povray.org | ビジョンレイトレーサーの永続性。 |
| Framescanner.exe | トッドキャンベル | AVIファイルコンバータへのPNG画像 |
| Windows Liveムービーメーカー | マイクロソフト | AVを圧縮するために使用私はWMVへ |
表1に必要なプログラムの表。
| ファイルの名前 | タイプ | 説明 |
| アライグマGA.mdb | Accessデータベース | アライグマ流域の構造と管理の説明。 GeneticiSWATとMapSWATによる読み込み。 |
| 流域presets.csv | テキスト | GeneticiSWAT.exeとMapSWATの設定プリセット |
| Alleles.csv | テキスト | 進化的アルゴリズムのための対立遺伝子セットの一覧。 |
| アライグマアレルHRU.txt | テキスト | GeneticISWATが農地であると判断、それらの対立を記載することによって作成されたファイル。 MapSWATによる読み込み。 |
| subbasin Josh.mdbの訓練費 | Accessデータベース | 管理の実践や郡によってコストを削減できます。 |
| テラスZones.mdb | Accessデータベース | 表[HUCデータ]流域用のテラスと歩留まりゾーン番号を保持しています。 |
| NRI Budgets.mdb | Accessデータベース | この実行で使用されていない作物&マシンのテーブルのGeneticISWAT.exeによる読み込み。 |
| phucrp 2008-12-15.dat | テキスト | この実行で使用されていない植物の熱ユニットルックアップテーブル、。 |
| Management.mdb | Accessデータベース | この実行で使用されていない回転ルックアップテーブル、。 |
| アライグマジョージア2011年9月28日1313.log、アライグマジョージア2011年9月29日0732.log、アライグマジョージア2011年10月7日0644.log | テキスト | GeneticISWAT実行のログファイル。 |
| アライグマGA.wmv | アニメーション | 世代別の個人の3Dディスプレイ |
| Subbasin.shp > | ESRIシェイプファイル | 流域におけるsubbasinsの概要。 |
| アライグマMap.wmv | アニメーション | フロンティア上の個々の各subbasinための支配的な対立の表示。 |
表2に必要なサンプルファイルの表。
Access restricted. Please log in or start a trial to view this content.
私たちは、流域レベルの栄養素削減目標の範囲を達成するために農業保全慣行の最低コストミックスと場所を含む流域構成のパレート効率的なセットを検索するための統合されたシミュレーション·最適化フレームワークを構築します。シミュレーション·最適化システムの概念図を図8に示されている。農業保全慣行の水質への影響をシミュレートすることを含む、流域シミュレーションは、Windowsベースのデータベース管理システム、i_SWAT 6,8と相まって水文モデル、SWAT2005によって処理されます。最適化コンポーネントは、SWATの水文応答単位(HRUs)上で動作し、同時に栄養負荷(N、P、または両方)と保全慣行のコストを最小限に抑える保全慣行の配分を見つけるために進化アルゴリズム26のロジックを採用しています。アルゴリズムの反復が、個々を存続のセットを終えられた後にualsは、おおよそのトレードオフのフロンティアを表しています。 2つの栄養素を同時に(硝酸塩-N及び全リン)が検討されているので、我々は3次元のトレードオフのフロンティアを入手してく...
Access restricted. Please log in or start a trial to view this content.
特別な利害関係は宣言されません。
この研究は、米国環境保護庁のターゲット流域グラント·プログラム(プロジェクト#WS97704801)、結合自然と人間系の国立科学財団のダイナミクス(プロジェクト#DEB1010259-CARD-クリン)、及び米国司法省から受け取ったサポートから一部で賄われていたFoodand農業協定の農業プロジェクトの農研究所(プロジェクト#20116800230190-CARD-)の。
Access restricted. Please log in or start a trial to view this content.
Access restricted. Please log in or start a trial to view this content.
Request permission to reuse the text or figures of this JoVE article
Request Permission