6-1 結論
前人研究之邊界源解析解模式已被廣泛研究與討論，但現地場址通常為任意 位置多物種污染源注入情形，故在探討具多內部源溶質傳輸污染問題時容易受到 限制且不利廣泛應用。而過去發展的單物種傳輸模式無法模擬序列降解下多物種 的衰減與傳輸情形。本研究發展具多內部源二維多物種溶質傳輸解析解模式，並考 慮移流、延散傳輸、線性平衡吸附及一階衰減反應，以描述二維溶質的傳輸行為。
本研究推導解析解過程中，因為使用 finite Fourier cosine 轉換與廣義型積分轉 換消去空間微分項，故產生兩個無窮級數 l、n 累加，Pe 較大的情況下所需較多的 累加項數 l、n 才能達成收斂，反之 Pe 較小時達成收斂所需的累加項數 l、n 較少。
發展之解析解與有限差分數值解在不同 Pe 值下均有良好的驗證結果，並確認此模 式解析解推導的正確性。
本研究接著針對實際的含氯有機污染場址進行現地模擬，並模擬多物種污染 團序列降解反應下的濃度分布情形，其中由於污染源為 DNAPL 的殘留項，較不適 合做為邊界條件，故本研究考慮多個內部源注入較符合實際情況，模式可快速進行 長時間情境的模擬，且對於場址內污染物的擴散與整治情況能夠有初步的了解。由 模擬結果可知，發展的模式可表現多污染源釋出情形，並表現污染團濃度集中及隨 地下水流遷移的現象，但與 Dover site 現地數據比較發現無法準確預測污染團濃度 分布情形，亦顯示模式在模擬大尺度流場時，無法準確預測地下水流況，且 Pe 大 於 25 時解析解模式不容易達到收斂，模擬結果在 x ≥300 m 容易產生誤差。
6-2 建議
具多內部源溶質傳輸模式可以幫助了解現地場址污染源的釋出與傳輸情形，
並可作為土壤與地下水污染整治參考與依據，而現地情況多為三維流場，溶質傳 輸除了受到縱向與側向延散效應影響外，地下水流與污染物傳輸亦與流場厚度與 重力影響有關，故本研究發展的具多內部源二維多物種溶質傳輸解析解方法可擴 展至求解三維溶質傳輸解析解模式中，可更完整的描述污染物的溶質傳輸情形。
本研究發展的解析解模式考慮線性平衡吸附作用並假設延散係數為常數，可 進一步擴展至考慮非平衡吸附或尺度延散作用下的溶質傳輸模式，以描述各種不 同情況溶質傳輸情形。
具多內部源多物種溶質傳輸模式可嘗試應用在其他中小尺度或是 Pe≤25 之流 場，並與現地數據比較，討論發展模式的適用性。
參考文獻
Atteia, O., Franceschi, M., Dupuy, A., (2008). Validation of Reactive Model Assumptions with Isotope Data: Application to the Dover Case. Environmental Science &
Technology, 42, 3289-3295.
Aziz, C.E., Newell, C.J., Gonzales, J.R., Haas, P., Clement, T.P., Sun, Y., (2000).
BIOCHLOR: Natural attenuation decision support system, User’s manual version 1.0, EPA/600/R-00/008. USEPA Office of Research and Development, Washington, DC.
Batu, V., (1989). A generalized two-dimensional analytical solution for hydrodynamic dispersion in bounded media with the first-type boundary condition at the source.
Water Resources Research, 25 (6), 1125–1132.
Batu, V., (1993). A generalized two-dimensional analytical solute transport model in bounded media for flux-type finite multiple sources. Water Resources Research, 29(8), 2881-2892.
Bauer, P., Attinger, S., Kinzelbach, W., (2001). Transport of a decay chain in homogeneous porous media: analytical solutions. Journal of Contaminant Hydrology, 49, 217-239.
Chen, J.S., Lai, K.H., Liu, C.W., Ni, C.F., (2011). A novel method for analytically solving multi-species advective–dispersive transport equations sequentially coupled with first-order decay reactions. Journal of Hydrology, 420-421, 191-204.
Chen, J.S., Liang, C.P., Liu, C.W., Li, L.Y., (2016). An analytical model for simulating two-dimensional multispecies plume migration. Hydrology and Earth System Sciences, 20, 733-753.
Chiang, C.Y., Salanitro, J.P., Chai, E.Y., Colthart, J.D., Klein, C.L., (1989). Aerobic biodegradation of benzene, toluene, and xylene in a sandy aquifer-data analysis
and computer modeling. Ground Water, 27,823–834.
Cho, C.M., (1971). Convective transport of ammonium with nitrification in soil. Canadian Journal of Soil Science, 51 (3), 339-350.
Clement, T.P., (2001). Generalized solution to multispecies transport equations coupled with a first-order reaction network. Water Resources Research, 37:157-163.
Clement, T.P., (1997). RT3D: A modular computer code for simulating reactive multi-species transport in 3-dimensional groundwater systems, Pacific Northwest National Laboratory Report, PNNL-SA 11720.
Clement, T.P., Sun, Y., Hooker, B.S., Petersen, J.N., (1998). Modeling multi-species reactive transport in groundwater aquifers. Groundwater Monitoring Remediation, 18(2), 79-92.
Clement, T.P., Johnson, C.D., Sun, Y., Klecka, G.M., Bartlett, C., (2000). Natural attenuation of chlorinated ethane compounds: model development and field-scale application at the Dover site. Journal of Contaminant Hydrology, 42,113 -140.
Cotta, R.M., (1993). Integral transforms in computational heat and fluid flow. CRC Press, Boca Raton, FL.
Essaid, H.I., Bekins, B.A., Godsy, E.M., Warren, E., Baedecker, M.J., Cozzarelli, I.M., (1995). Simulation of aerobic and anaerobic biodegradation processes at a crude oil spill site. Water Resources Research, 31(12),3309-3327.
Lu, G., Clement, T.P., Zheng, C., Wiedemeier, T.H., (1999). Natural attenuation of BTEX compounds: model development and field-scale application. Ground Water, 37(5), 707-717.
Lunn, M., Lunn, R.J., Mackay, R., (1996). Determining analytic solution of multiple species contaminant transport, with sorption and decay. Journal of Hydrology, 180,
Newell, C.J., McLeod, R.K., Gonzales, J., (1996). BIOSCREEN: Natural Attenuation Decision Support System, User’s manual version 1.3, EPA/600/R-96/087. USEPA Office of Research and Development,Washington, DC.
Parlange, J.Y., Barry, D.A., Starr. J.L., (1985). Comments on “boundary conditions for displacement experiments through short laboratory soil columns.” Soil Science Society of America Journal, 49, 1325.
Parlange, J.Y., Starr, J.L., (1978). Dispersion in soil column: effect of boundary conditions and irreversible reactions. Soil Science Society of America Journal, 42, 15-18.
Parlange, J.Y., Starr, J.L., van Genuchten, M.Th., Barry, D.A., Parker, J.C., (1992). Exit condition for miscible displacement experiments in finite columns. Soil Science, 153, 165-171.
Park, E., Zhan, H., (2001). Analytical solutions of contaminant transport from finite one-, two-, three-dimensional sources in a finite-thickness aquifer. Journal of Contaminant Hydrology, 53 (1-2), 41–61.
Parker, J.C., van Genuchten, M.Th., (1984) Fluxaveraged and volume-averaged concentrations in continuum approaches to solute transport. Water Resources Research, 20, 866-872.
Pérez Guerrero, J.S., Pimentel, L.C.G., Skaggs, T.H., van Genuchten, M.Th., (2009).
Analytical solution of the advection-diffusion transport equation using a change-of-variable and integral transform technique. Int. J. Heat Mass Transf., 52, 3297- 3304.
Pérez Guerrero, J.S., Skaggs, T.H., van Genuchten, M.Th., (2010). Analytical Solution for Multi-Species Contaminant Transport in Finite Media with Time-Varying Boundary Conditions. Transp Porous Med., 85:171–188.
Quezada, C.R., Clement, T.P., Lee, K.K., (2004). Generalized solution to multi- dimensional multi-species transport equations coupled with a first-order reaction network involving distinct retardation factors. Adv. Water Resources, 27, 507-520.
Rifai, H.S., Bedient, P.B., Borden, R.C., Haasbeek, J.F., (1987).BIOPLUME II:
Computer model of two-dimensional contaminant transport under the influence of oxygen limited biodegradation in ground water. National Center for Ground Water Research, Rice University.
Rifai, H.S., Bedient, P.B., Wilson, J.T., Miller, K.M., Armstrong, J.M., (1988).
Biodegradation modeling at aviation fuel spill site. Journal of Environmental Engineering, 114(5), 1007-1029.
Semprini, L., Kitanidis, P.K., Kampbell, D.H., Wilson, J.T., (1995).Anaerobic transformation of chlorinated aliphatic hydrocarbons in a sand aquifer based on spatial chemical distribution. Water Resources Research, 31, 1051-1062.
Srinivina, V., Clement, T.P., (2007). Analytical solution for sequentially coupled one- dimensional reactive transport problems – Part I:Mathematical derivations.
Advances in Water Resources, 31, 203-218.
Srinivina, V., Clement, T.P., (2007). Analytical solution for sequentially coupled one- dimensional reactive transport problems – Part II:Special cases, implementation and testing. Advances in Water Resources, 31, 219-232.
Sun, Y., Clement, T.P., (1999). A decomposition method for solving coupled multi- species reactive transport problems. Transport in Porous Media, 37, 327-346.
Sun, Y., Petersen, J.N., Clement, T.P., Hooker, B.S., (1998). Effects of reaction kinetics on predicted concentration profiles during subsurface bioremediation. Journal of Contaminant Hydrology, 31, 359-372.
solutions for multispecies transport with serial and parallel reactions. Water Resources Research, 35, 185-190.
Van Genuchten, M.Th., Alves, W.J., (1982). Analytical solutions of the one-dimensional conv e ct i v e – di s pe rs i ve s ol ut e t r ans por t equ at i on . US D epa rt m ent o f Agriculture.Technical Bulletin Number 1661.
Van Genuchten, M.Th., (1985). Convective–dispersive transport of solutes involved in sequential first-order decay reactions. Computer & Geosciences, 11(2), 129-147.
Van Genuchten, M.Th., Parker, J.C., (1984) Boundary conditions for displacement experiments through short laboratory soil columns. Soil Science Society of America Journal, 48, 703-708.
田傑仁(2013)，「層狀地質介質溶質傳輸之解析解」，碩士論文，國立中央大學應用 地質研究所。
凃佑霖(2015)，「二維有限域溶質傳輸解析解模式發展」，碩士論文，國立中央大學 應用地質研究所。
附錄
此附錄提出推導正交化特徵函數K(l,X)之過程。根據方程式(3-33)與邊界條 件(3-34)、(3-35)定義特徵值問題如下：
0
其中
l為特徵值，則正交化特徵函數K (
l, X )
可表示為