宜蘭清水地區是台灣在1981年曾經發展過地熱發電的地區，但由於後期出水孔隙阻塞嚴重，熱水與蒸汽量的減少，導致清水地熱區發電終止運作。為了解深部是否有較佳的地熱潛能條件，本研究將使用能探測深部地層（大於1公里）的方法－大地電磁法(Magnetotellurics, MT)，對清水地區做的深部地層的調查。
MT使用上多有限制，特別易受人為電磁雜訊影響，本研究除了依循傳統的MT方法，並研究改善資料品質，包括延長資料收集時間、遠端參考站的使用以及GB張量分解等。
MT結果並參考前人研究與地質的分析，可以看出地熱區與低電阻構造呈現相關性，由此推估測區的地熱，極可能是來自清水地區南方深部8公里，但溫泉水的來源，應該只是天水深入地層受地熱加溫，再沿岩層裂隙出露地表，並非直接來自深層地熱水。测區有明顯熱源的集中，極具發電潛能，此外本研究還看到更深層的熱源在20公里以下，推測此深部的熱源可能來至台灣地區的板塊運動推擠有關，有待後續調查。

The Chingshui geothermal area in Yilan has ever been developed as a geothermal power plant in Taiwan, 1981. The power station was stopped operation because the water pipe had been blocked seriously and the hot water and steam reduced greatly. We did MT survey in order to understand whether there is potential geothermal energy in depth. For deep investigation, we used the magnetotelluric (MT) method that could analyze the deep structure more than 1 km.
To increase the signal to noise ratio, besides following the traditional MT data collection, we have tested to record MT data by prolonging time as long as three days and we also use the remote-reference for improving the quality of data. The result shows that three days record and remote-reference are needed for high quality MT data acquisition in Taiwan area. Then we use the Groom and Bailey decomposition analysis (GB decomposition) for the geoelectric strike determinations.
The result of MT study in Chingshui area shows the relation between geothermal source and low resistivity structure. We can suggest that the heat source comes from the south of the region and its depth is about 8 km. The hot spring could be the meteoric water which permeated into deep strata and was heated by geothermal source and then appeared on the earth’s surface along the fracture of strata. Thus, the Chingshui geothermal area is high potential for power plant. Lastly, we also find a deeper heat source that down to the depth of 20 km which is perhaps due to the collision of plates in Taiwan area.

Bai, D., Meju, M.A., and Liao, Z., Magnetotelluric images of deep crustal structure of the Rehai geothermal field near Tengchong, southern China, Geophys. J. Int., 147, 677-687, 2001.
Cagniard, L., Basic theory of the magneto-telluric method of geophysical prospecting, Geophysics, 18, 605-635, 1953.
Chang, Chung-Pai, Strain and stress field in Taiwan oblique convergent system: constraints from GPS observation and tectonic data, Earth and Planetary Science Letters, 214, 115-217, 2003.
Chang, C.R.Y., Ramey, H.J., Well interference test in the Chingshui geothermal field. Paper presented at Fifth Geothermal Reservoir Engineering Workshop, Stanford University, Stanford, CA, 64–76, 1979.
Gamble, T.D., Goubau, W.M., and Clark, J., Magnetotelluric with a remote magnetic reference, Geophysics, 44, 53-68, 1979.
Gianni, V., Adele, M., Adolfo F., Investigation of geothermal structures by magnetotellurics (MT): an example from the Mt. Amiata area, Italy, Geothermics, 32, 131-145, 2003.
Groom R.W., and Bailey, R.C., Decomposition of magnetotelluric impedance tensor in the presence of local three-dimensional galvanic distortion, J. Geophy. Res., 94, 1913-1925, 1989.
Groom, R.W., and Bahr, K., Correstion of near surface effects: Decomposition of magnetotelluric impedance tensor and scaling correction for regional resistivities, Geophysics, 13, 341-380, 1992.
Harinarayana, T., Abdul Azeez, K.K., Naganjaneyulu, K., Manoj,C., Veeraswamy, K., Murthy, D.N., and Rao, S.P.E., Exploration of geothermal structure in Puga geothermal field, Ladakh Himalayas, India by magnetotelluric studies, Journal of Applied Geophysics, 58, 280-295, 2006.
Jones, A. G., Chave, A.D., Egbert, G., Auld, D, Bahr, K., A comparison of techniques for magnetotelluric response function estimation. J. Geophy. Res., 94, 14201-14213, 1989.
Jones, A. G. and Jodicke, H., Magnetotelluric transfer function estimation improvement by a coherence-based rejection technique, paper presented at 54th Annual International Meeting, Society of Exploration Geophysicists, Atianta, Dec. 2-6, 1984.
Price, A. T., The theory of magnetotelluric methods when the source field is considered, J. Geophy. Res., 67, 5, 1907-1918, 1962.
Rodi, W., Mackie, R. L., Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion. Geophysics, 66, 174-187, 2001.
Vozoff, K., Hasegawa, H., and Ellis, R.M., Results and limitation of magnetotelluric surveys in simple geologic situation, Geophysics, 28, 778, 1963.
Vozoff, K., and Ellis, R., Magnetotelluric measurements in southern Alberta, Geophysics, 31, 1153-1157, 1966.
Vozoff, K., The Magnetotellruic Method, Chapter 8, Electromagnetic method in applied geophysics-Applications part A and part B, edit by Corbett, J.D., published by Scociety of Exploration Geophysicists, 641-711, 1991.
Wen Xiao, Magnetotelluric Exploration in the Rocky Mountain Foothills, Alberta, master thesis of UNIVERSITY OF ALBERTA, 2004.
Yen, Tsang-Po, etc., A stratigraphical study on the Tananao Schist in northern Taiwan, Bulletin of the Geological Survey of Taiwan, 12, 53-66, 1960.
林啟文、林偉雄，五萬分之一台灣地質圖：南澳圖幅，經濟部中央地質調查所，1993。
林啟文、林偉雄，五萬分之一台灣地質圖：三星圖幅，經濟部中央地質調查所，1995。
季松青，大地電磁資料多站多頻分析於台灣中部及金門地區地殼電性構造，國立中央大學博士論文，145頁，2006。
姜智文，應用大地電磁法研究台灣中部地區之電性構造，國立中正大學碩士論文，台灣嘉義，85頁，2004。
張麗旭，臺灣變質區第三系基於小型有孔蟲之生物地層學研究（四：中央山脈最北部），中國地質學會會刊，17，85-93，1974。
曾長生，宜蘭縣清水及土場區地質及地熱產狀，台灣石油地質，15，11-23，1978。
何春蓀，臺灣地質概論，中華民國經濟部，118頁，1975。