[httpresearch.ncku.edu.twrearticlesc200812261.html.pdfVIP

成大研發快訊 - 文摘 成大研發快訊 第七卷 第二期 - 2008年十二月二十六日 [ .tw/re/articles/c1.html ] 碎形非均勻隨機介質場之溶質傳輸研究 * 徐國錦 、陳冠志 國立成功大學工學院資源工程學系 kchsu@.tw Chen, K.-C. and K.-C. Hsu, A General Fractal Model of Flow and Transport in Randomly Heterogeneous Porous Media, Water Resources Research, 43, W12501, doi:10.1029/2007WR005934, 2007.(SCI) 水力傳導係數為影響地下水流及溶質傳輸最為重要的物 理參數之一。由於水力傳導係數的高度空間變異性及實務 上資料取得的困難，傳統的定率式分析方法無法對異質場 進行分析，取而代之的是序率式分析方法。過去諸多研究 指出，現地水力特性參數在空間上之變異分佈適合以高斯 碎形模式來加以描述，其中碎形高斯噪音(fGn)為最普遍被 使用的一種碎形模式，fGn之累加即成為碎形布朗運動 (fBm)。近年來已有越來越多地球物理井測資料及現地岩心 證據顯示，現地水力特性參數之空間變異性往往呈現非高 斯分佈型態，而較適用以非高斯碎形模式例如碎形Levy運動(fLm)及碎形Levy噪音(fLn)來加以描述。FLn 指的是變數的空間分佈呈現Levy-stable機率分佈型態，而fLm所指的是變數的差值呈現Levy-stable機率分 佈型態。Levy-stable機率分佈類似於高斯分佈但具有更大的使用彈性，其機率分佈之尾部具有冪指數之效 應，可解釋為在空間尺度上統計特徵具有向上與向下延伸的特性，稱之為重尾效應。此外，相較於高斯分 佈，Levy-stable機率分佈函數在均值處具有較高之機率值。重尾效應使得Levy模式更適用於描述變異劇烈 之水力特性參數，但也使得Levy-stable機率分佈函數在使用上遭遇到具有無窮大變異數的困擾。 由於碎形模式在數學上具有變異數無限大的特性，使得碎形模式在實際使用上會產生困難，因此，可採用 下列兩種方式來加以解決：(1)給予Levy-stable機率密度函數一個門檻值，依此產生出Levy碎形模式之地質 隨機場，以進行水流及溶質傳輸之數值模擬；(2)給予fBm之變差函數一個門檻值，再藉由最小擾動法來進 行水流及溶質傳輸之數值模擬，門檻值之上限相當於研究場址的範圍大小，門檻值的下限相當於採樣的精 度。 上述方法遭遇的困難包括Levy-stable機率密度函數門檻值的定義並無有一定準則，加上傳統高斯碎形模式 無法描述高度非均勻性地質場的限制，有必要統合上述方法發展出通用碎形模式，使得序率理論能夠更靈 活地運用於地下水流及溶質傳輸行為之描述。 碎形模式主要用以描述一個具有自我仿射特性的序率過程，fLm為fBm之通式，同時也具有統計上之自我 仿射特性及尺度不變性，統計特性可表示如下式 or (1) 式中α為Levy指數，C 及C 分別代表分隔距離為sh及h時之Levy-stable機率分佈函數之寬度。當H =1/α sh h 時，fLm為一個單純隨機過程；當1/αH 1時，為一個持續性過程；當0H 1/α時，則為一個反持續性過 程。在此，為了滿足0H 1，α值需限制在1～2之間。此外，fLm之冪指數型變差函數(power variogram) 可表示為 1 of 5 成大研發快訊 - 文摘 (s)=a sH (2) 1 式中 為常數，H為H