第四章 討論
4.6 建議的對比流程
由理論與本研究的經驗對比流程包含三個部分：檢查資料、井測資料對比海
水面變遷曲線和濾波。檢查井測資料的目的是要檢查資料的完整性與採樣頻率，
以確認該資料是否有測量上的錯誤及判斷資料是否可以進行濾波(即採樣頻率為 定值)。此外，井測資料必須要反映穩定的沉積環境才有對比的價值。由對比成 果決定進行濾波的材料(井測資料或海水面變遷曲線)及濾除的頻帶。反覆進行對 比，最後選擇出最佳的對比成果來進行地質學解釋。
對比流程(圖 4.1)如下:
第一步 檢查井測資料是否完整。
第二步 確認井測資料反映穩定沉積環境中的沉積相變化。
第三步 井測資料對比。
第四步 決定需要濾波的材料與濾波頻帶；若對比成果已經很好，則 不需要濾波，完成對比工作。
第五步 濾波與對比。
反覆進行第四步與第五步。
第六步 獲得最佳的對比成果，完成對比。
圖 4.1 對比流程
4.7
未來值得進一步探討未來值得進一步探討未來值得進一步探討的項目未來值得進一步探討的項目的項目的項目五、井測資料對比海水面邊遷曲線整合研究計畫
本研究採用國際海洋鑽探計畫與中油公司提供的井測資料與年代模式，雖然 已是相當好的觀測資料與可信的年代模式，仍有不足之處：其一是沒有證據 顯示珈瑪井測資料忠實反映沉積物粒度變化，其二是沒有高解析度的年代模 式。若未來有機會進行深入的探討，期望能在穩定的沉積環境進行地球物理 井測、採標本進行粒度分析並建立高解析度的年代模式。如果能對年代模式 及岩層粒度變化有充份的掌控，就可以確切地證明井測資料對比全球海水面 變遷曲線的方法是否可靠。除了海相地層的研究之外，也希望能測試本研究 的對比技術是否能應用在陸相的沉積地層。
六、尋找其他適合與全球海水面變遷曲線對比以建立年代模式地質資料。
第五章 第五章 第五章
由本研究井測資料對比全球海水面變遷曲線的成果，基本上可以發現：在大 陸邊緣具有穩定沉積速率的地層中，大多數探測井的珈瑪射線紀錄，皆會依循著
大尺度的海水面升降週期而改變。至於對比性較差的探測井，推測是受到所處位
置沉積物受地形阻隔、深海海流造成的沉積流沖刷、重力崩落、受侵蝕作用等等 狀況之影響，使得原有的沉積層序受到破壞，影響其對比性。總合來說，珈瑪井 測紀錄應該可以反映全球海水面的升降變遷。即便理論基礎或牽涉到的變因都不 是很清楚；就經濟效益觀點，此技術可以快速且有效率地決定地層的沉積年代，
應該可以廣泛地運用到全球的探勘工作上。
參考文獻 參考文獻 參考文獻 參考文獻
Bleuer, N. K., 2004, Slow-logging subtle sequences : the gamma-ray log character of glacigenic and other unconsolidated sedimentary sequences, Bloomington, Indiana University : Indiana Geological Survey.
Boggs, S., 1995, Principles of sedimentology and stratigraphy 2nd edition, Englewood Cliffs, N.J., Prentice Hall. 774 p.
Chen, C.-S., and Jeng, Y., 2013, Natural logarithm transformed EEMD instantaneous attributes of reflection data: Journal of Applied Geophysics, v. 95, no. Aug., p.
53-65.
Ellis, D. V., and Singer, J. M., 2007, Well logging for earth scientists, Dordrecht, The Netherlands, Springer. 692 p.
Fleming, A. H., 1994, The hydrogeology of Allen County, Indiana: a geologic and ground-water atlas, Allen County Dept. of Planning Services, Indiana Geological Survey Special Report, 111 p.
Haq, B. U., Hardenbol, J., and Vail, P. R., 1987, Chronology of Fluctuating Sea Levels Since the Triassic: Science, v. 235, no. 4793, p. 1156-1167.
Haykin, S. S., and Van Veen, B., 2005, Signals and systems 2nd edition, Hoboken, NJ, John Wiley & Sons. 820 p.
Hays, J. D., Imbrie, J., and Shackleton, N. J., 1976, Variations in the Earth's Orbit:
Pacemaker of the Ice Ages: Science, v. 194, no. 4270, p. 1121-1132.
Huang, N. E., Shen, Z., Long, S. R., Wu, M. C., Shih, H. H., Zheng, Q., Yen, N.-C., Tung, C. C., and Liu, H. H., 1998, The empirical mode decomposition and the
Jeng, Y., Lin, M.-J., Chen, C.-S., and Wang, Y.-H., 2007, Noise reduction and data recovery for a VLF-EM survey using a nonlinear decomposition method:
Geophysics, v. 72, no. 5, p. F223-F235.
Jervey, M. T., 1988, Quantitative geological modeling of siliciclastic rock sequences and their seismic expression: Sea-Level Changes: An Integrated Approach:
SEPM, Special Publication, v. 42, p. 47-69.
Keigwin, L. D., Rio, D., Acton, G. D., Bianchi, G. G., Borowski, W., Cagatay, N., Chaisson, W. P., Clement, B. M., Cortijo, E., Dunbar, G. B., Flood, R. D., Franz, S.-O., Giosan, L., Gruetzner, J., Hagen, S., Haskell, B., Horowitz, M. J., Laine, E. P., Lund, S. P., Okada, M., Poli, M.-S., Raffi, I., Reuer, M. K.,
Ternois, Y. G., Williams, T., Winter, D. M., Yokokawa, M. E., and Swanson, S.
E., 1998, Proceedings of the Ocean Drilling Program; Initial reports;
Northwest Atlantic sediment drifts, covering Leg 172 of the cruises of the drilling vessel JOIDES Resolution, Charleston, South Carolina, to Lisbon, Portugal, Sites 1054-1064, 14 February-15 April, 1997: Ocean Drilling Program.
Kendall, C. G. S. C., and Lerche, I., 1988, The rise and fall of eustasy: Sea-Level Changes: An Integrated Approach: SEPM, Special Publication, v. 42, p. 3-17.
Laskar, J., Fienga, A., Gastineau, M., and Manche, H., 2011, La2010: a new orbital solution for the long-term motion of the Earth: Astronomy & Astrophysics, v.
532, A89.
Laskar, J., Joutel, F., and Boudin, F., 1993, Orbital, precessional, and insolation quantities for the Earth from -20 Myr to+ 10 Myr: Astronomy and Astrophysics, v. 270, p. 522-533.
Laskar, J., Robutel, P., Joutel, F., Gastineau, M., Correia, A., and Levrard, B., 2004, A long-term numerical solution for the insolation quantities of the Earth:
Astronomy & Astrophysics, v. 428, no. 1, p. 261-285.
Lisiecki, L. E., and Lisiecki, P. A., 2002, Application of dynamic programming to the correlation of paleoclimate records: Paleoceanography, v. 17, no. 4, 1049.
Lisiecki, L. E., and Raymo, M. E., 2005, A Pliocene-Pleistocene stack of 57 globally distributed benthic δ18O records: Paleoceanography, v. 20, no. 1, PA1003.
Mallat, S. G., 2009, A wavelet tour of signal processing : the sparse way 3rd edition, Boston, Elsevier/Academic Press. 832 p.
Milankovitch, M., 1969, Canon of insolation and the ice-age problem (Kanon der Erdbestrahlung und seine Anwendung auf das Eiszeitenproblem) Belgrade, 1941, Jerusalem, Israel Program for Scientific Translations; [available from U.S. Dept. of Commerce, Clearinghouse for Federal Scientific and Technical Information, Springfield, Va.].
Miller, K., Mountain, G., Wright, J., and Browning, J., 2011, A 180-million-year record of sea level and ice volume variations from continental margin and deep-sea isotopic records: Oceanography, v. 24, no. 2, p. 40-53.
Miller, K. G., Kominz, M. A., Browning, J. V., Wright, J. D., Mountain, G. S., Katz, M. E., Sugarman, P. J., Cramer, B. S., Christie-Blick, N., and Pekar, S. F., 2005, The Phanerozoic record of global sea-level change: science, v. 310, no.
5752, p. 1293-1298.
Mitra, S. K., 2001, Digital Signal Processing: A Computer-based Approach 2nd edition, Boston, McGraw-Hill Higher Education. 896 p.
Raffi, I., 2002, Revision of the early-middle pleistocene calcareous nannofossil biochronology (1.75–0.85 Ma): Marine Micropaleontology, v. 45, no. 1, p.
25-55.
Rider, H., 2002, The Geological Interpretation of Well Logs 2nd Edition, Houston, Rider-French Consulting Limited, 280 p.
Tian, J., Zhao, Q., Wang, P., Li, Q., and Cheng, X., 2008, Astronomically modulated Neogene sediment records from the South China Sea: Paleoceanography, v. 23, no. 3, p. PA3210.
Vail, P. R., Mitchum, R. M., Jr., and Thompson, S., III, 1977a, Seismic Stratigraphy and Global Changes of Sea Level: Part 3. Relative Changes of Sea Level from Coastal Onlap: Section 2. Application of Seismic Reflection Configuration to Stratigrapic Interpretation: Seismic Stratigraphy--Applications to Hydrocarbon Exploration Memoir 26, p. 63-81.
Vail, P. R., Todd, R. G., and Sangree, J. B., 1977b, Seismic Stratigraphy and Global Changes of Sea Level: Part 5. Chronostratigraphic Significance of Seismic Reflections: Section 2. Application of Seismic Reflection Configuration to Stratigraphic Interpretation: Seismic Stratigraphy--Applications to
Hydrocarbon Exploration Memoir 26, p. 99-116.
Wu, Z., and Huang, N. E., 2005, Ensemble empirical mode decomposition: a noise assisted data analysis method, Calverton, MD, Center for
Ocean-Land-Atmosphere Studies, COLA Technical Reports, v. 193. 49p.
附錄 附錄
附錄 附錄一 一 一 一 使用 使用 使用 MATLAB 繪製頻譜 使用 繪製頻譜 繪製頻譜 繪製頻譜
% 以繪製 SL05 曲線的頻譜(本論文的圖 3.1)為例：miller2005_imf.xls 中 miller imf 分頁的第一列為年代、第二列是海水面高度。
% 輸入 xls 格式(excel)資料。紅色為選取的 xls 檔，藍色為選取的分頁。
A = xlsread('miller2005_imf.xls','miller imf');
age=A(:,1); %將載入資料的第一列定為 age
% 準備繪製頻譜 Y = fft(SL,NFFT)/L;
% 進行快速傅立葉轉換；紅色為輸入資料，藍色為 FFT 的點數
% 將 FTT 處理成果除以資料點數目(L)，將成果存為 Y。
f = Fs/2*linspace(0,1,NFFT/2+1); % 算出頻譜資料的 x 軸座標 subplot(2,1,2); % 繪製圖型的下半部
plot(f,2*abs(Y(1:NFFT/2+1))); % 繪製頻譜。
grid; % 自動繪製網格
title('Spectrum'); % 次標題命名為 Spectrum xlabel('Frequency (times per ma)');
ylabel('|Y(f)|');
% x 軸座標標記為 Frequency (times per ma)，y 軸座標標記為|Y(f)|
% 儲存圖型
saveas(gcf,'miller 2005.fig'); %將圖型存為 miller 2005.fig (MATLAB 圖型格式) saveas(gcf,'miller 2005.jpg'); %將圖型存為 miller 2005.jpg
附錄二 附錄二
spectrogram(m1,64,60,64,1/0.005,'yaxis'); % 繪製時頻圖
% spectrogram(輸入資料,視窗的點數,重疊點數,FFT 的點數,採樣頻率,頻率在座 標軸)。
ylabel('frequency (times per Ma)');xlabel('AGE (Ma)');
% y 座標標為 frequency (times per Ma)；y 座標標為 AGE (Ma)
colorbar; % 自動產生圖例
colormap('hsv'); % 採用 Hue-saturation-value color map title('SPECTROGRAM of SL05'); % 標題為'SPECTROGRAM of SL05
alpha(0.8); % 調整透明度
grid minor; set(gca,'GridLineStyle', '-'); % 設定網格形式
附錄三 附錄三 附錄三
save SL05_IMF.txt -ASCII M1
經驗模組拆解所獲得成分資料，可以透過各種軟體繪圖。
為了處理大量資料方便，本研究採用師大地球物理實驗室開發的 MATLB 批次 經驗模組拆解與繪圖指令，進行資料處理與繪圖，如果讀者有需要請洽台灣師 範大學地球科學系鄭懌教授。