地下流體與多孔隙介質礦物顆粒間產生地球化學反應，伴隨著孔隙率與水力傳導係數變化，此現象在科學、工程與工業領域皆為重要之研究議題。本研究進行一系列批次表面反應實驗及管柱流動實驗，進而瞭解流體流經碳酸鈣產生溶解與沈澱反應對孔隙率與水力傳導係數之變化。利用鹽酸HCl(aq) 及硫酸H2SO4(aq) 兩種溶液作為地下流體以加速反應速率。以測壓計測量水位，記錄管柱兩邊壓力差以達西定律計算水力傳導係數；測量出流水中鈣離子與硫酸根離子之濃度，以便計算孔隙率變化。管柱流動實驗結束，將管柱之碳酸鈣礦石填充床取出，進一步進行碳酸鈣微觀型態與礦物組成分析。研究結果發現碳酸鈣之粒徑與溶解/沈澱彼此競爭以致孔隙率與水力傳導係數有震盪的情形；SEM觀測得知石膏沈澱皆於進流口，以上現象皆與理論相符。

Investigation of chemical reactions between a fluid and the porous medium through which it flows, and of the subsequent changes in the medium porosity and hydraulic conductivity, is an area of interest for many scientific, industrial and engineering. Dissolution and precipitation are two of the most important processes affecting fluid chemistry, and they can significantly modify the physical and chemical properties of porous media. In this study, we investigate the effect of particle size of carbonate rock on temporal evolution of porosity and hydraulic conductivity induced by dissolution and precipitation by performing a series of laboratory surface reaction experiment and column flow experiment by injecting HCl/H2SO4. The temporal evolution of overall hydraulic conductivity is calculated by measurement of change in hydraulic head gradient and application of Darcy’s law. Variation in porosity is obtained by analyzing the effluent acid for Ca2+ and and invoking the law of mass balance. Furthermore, after each experiment, the column sample is retrieved and sectioned in order to study the micro-morphology and mineral composition. Results show that significant fluctuations in porosity and hydraulic conductivity occur due to competition between dissolution and precipitation. The patterns of fluctuations in porosity and hydraulic conductivity are influenced by particle size of calcium carbonate. SEM analysis shows that gypsum precipitates in the inlet of the column. These findings are in qualitative agreement with conceptual understanding of such phenomena.

[1]Dijk, P. and B. Berkowitz, “Precipitation and dissolution of reactive solutes in fractures’’, Water Resources Research, 34, 457-470, 1998.
[2]Bolton, E. W., A. C. Lasaga, and D. M. Rye, “A model for the kinetic controlof quartz dissolution and precipitation in porous media flow with spatially variable permeability: Formulation and examples of thermal convection’’, Journal of Geophysical Research, 89, 4009-4025, 1984.
[3]Palmer, A. N., “Rates of limestone dissolution and calcite precipitation in cave streams of east-central New York State’’, Geological Society of America, 89, 227-264, 1996.
[4]Dijk, P. and B. Berkowitz, “Buoyancy-driven dissolution enhancement in rock fractures’’, Geology , 28 , 1051– 1054, 2000.
[5]Darmody, R. G., C. E. Thorn, R. L. Harder, J. P. L. Schlyter, and J. C. Dixon, “Weathering implications of water chemistry in an arctic-alpine environment, north Sweden’’, Geomorphology , 34, 89–100, 2000.
[6]陳文山，「台灣的岩石」，2002年岩盤工程研討會論文集，21-22， 中華大學，台灣新竹，2002年。
[7]何春蓀，「普通地質學」，三版，五南圖書出版公司，民國93年。
[8]陳錦芳，「恆春地區石灰岩地層之地化模擬」，國立臺灣大學，碩士論文，民國81年。
[9]林財富，「土壤與地下水污染：原理與應用」，中華民國環境工程學會，民國97年。
[10]張元耀，「孔隙率與化學反應之交互作用之溶質傳輸歷程」，國立臺灣大學，碩士論文，民國95年。
[11]Liu, C.W. and J.F. Chen, “The simulation of geochemical reactions in the Heng-Chung limestone formation’’, Applied Mathematical Modeling, 20, 549–558, 1996.
[12]De Graaff, J. W. M., and R. D. Schuiling, “Geochemical engineering in chalk: Neutralization of waste sulfuric acid’’, Annales Societatis Geologorum , 5, 345– 350 , 1997.
[13]Chadam, J., P. Ortoleva and A. Sen, “ Reactive-infiltration instabilities” , IMA Journal of Applied Mathematics, 36, 207–220, 1986.
[14]Sherwood, “ Stability of a plane reaction front in a porous medium” , Chemical Engineering Science , 42, 1823–1829, 1987.
[15]Hinch, E. J. and B. S. Bhatt, “ Stability of an acid front moving through rock” , Journal of Fluid Mechanics , 212, 279–288, 1990.
[16]Liu, X., A. Ormond, K. Bartko, Y. Li, and P. Ortoleva, “ A geochemical reaction-transport simulator for matrix aciding analysis andd esign” , Journal of Petroleum Science &Engineering, 17, 181–196, 1997.
[17]Hoefner, M. L. and H. S. Fogler r, “Influence of transport and reaction on wormhole formation in porous media”, AIChE journal, 44, 1933-1949, 1998.
[18]Fogler, H. S. and S. D. Rege, “ Porous dissolution reactors’’,
Chemical Engineering Communications, 42, 291–313, 1986.
[19]Singurindy, O. and B. Berkowitz, “Carbonate dissolution and precipitation in coastal environments: Laboratory analysis and theoretical consideration’’, Water Resources Research, 41, doi:10.1029/2005WR004433, 2005.
[20]Singurindy, O. and B. Berkowitz, “The role of fractures on coupled dissolution and precipitation patterns in carbonate rocks” , Advances in Water Resources, 28, 507–521, 2005.
[21]Singurindy, O. and B. Berkowitz, “Evolution of hydraulic conductivity by precipitation and dissolution in carbonate rock”, Water Resources Research, 39, doi: 10.1029/2001WR001055, 2003.
[22]行政院環境檢測所，水中鹼度檢測方法-滴定法(NIEA W449.0B)，民國92年。
[23]Krauskopf, K. B., “Introduction to Geochemistry’, McGraw-Hill, New York, 1967.
[24]行政院環境檢測所，水中總硬度檢測方法—EDTA滴定法（NIEA W208.51A），民國95年。
[25]Hoefner, M. L, H. S. Fogler, P. Stenius, and J. Sjoblom, “ Role of aciddiffusion in matrix acidizing of carbonates’’, Journal of petroleum technology, 32, 203– 208, 1987.
[26]Sudmalis, M. and R. Scheikholeslami, “Coprecipitation of CaCO3 and CaSO4”, Canadian Journal of Chemical Engineering, 78, 21–31, 2000.
[27]Hall, C. and D. C. Cullen, “Scanning force microscopy of gypsum dissolution and crystal growth”, AIChE Journal, 42, 232–238, 1996.
[28] Helfferich, F. G., “The theory of precipitation/dissolution waves”, AIChE Journal, 35, 75– 87, 1989.