时频分析之应用-用水行为分析.PDF

Visual Signal 技術文章專欄 ■ .tw ______________________________________________________2012 _Jan. 18 NO.8 【時頻分析之應用-用水行為分析】 我們要如何從地下水水位變化中，瞭解當地主要的人為抽水產業類型呢？本篇介紹如何利用Visual Signal 算出年平均水位變化及日平均水位變化，搭配頻譜與時頻分析來找出當地抽用水特徵。以下我們 以地下水位站-芳草為例，示範怎麼從原始水位訊號分析人為抽水行為。 首先我們匯入原始地下水訊號，如圖一 ，可以看到每年的五月為水位最低時，之後水位逐漸上升， 在十月左右達到一年中最高水位，之後又逐漸下降，此水位變化為季節性雨量豐枯所造成。 芳草(1) 16 10 2004 2006 2008 2010 Date 圖一 芳草(1)原始地下水位訊號 資料在紀錄的過程中常會因為機器故障、干擾而使數據產生異常值或是有發生缺漏的現象，此時可 以使用 FillNull 將缺值補遺，須注意使用缺值補遺時，如果缺值時間越長，選用的不同的內插方法彼此 間會差異越大，因此應該針對數據選擇恰當的補遺方式。本例使用線性內插(Linear interpolation) ，目的 在希望缺漏段不要有額外的人工訊號產生，影響時頻分析結果。而有些明顯的異常跳動，我們可以使用 RemoveBump 將資料局部驟升或驟降的訊號修正，圖二為資料前處利厚的結果 。 芳草(1) - ToRegular- FillNull-Remove Bump 16 12 2004 2006 2008 2010 Date 圖二 芳草(1)經前處理後之地下水位 接著我們以Selection 選取欲分析的資料區間。選取完 2003 年 1 月 1 日至2011 年 1 月 1 日的資料 後，首先利用接上 EquiphaseStats 元件計算年平均水位變化，來確認一年中水位的變化情形。在 EquiphaseStats 參數中選擇週期為一年，如圖三 ，並於Viewer 中的時日格式中選取MonthOnly ，結果如 圖四所示，可以更直覺的觀察前面所提到的季節變化情形。 圖三 EquiphaseStats 參數設定中，選取時間單位，週期 1 ，時間單位為年 Selection - EquiphaseStat 16 13 一月 四月 七月 十月 Date 圖四 芳草(1)之八年平均地下水位 接著準備將資料轉換到頻率域中呈現。為了觀看時方便，我們先利用XAxisUnit 元件將原始資料的 轉換橫坐標單位從小時轉成天。XAxisUnit 之設定如圖五 。 圖五 XAxisUnit 之參數設定，轉換X 軸時間單位 由Selection 接至XAxisUnit 後再接上 FFT 元件進行傅立葉轉換，FFT 的結果，同樣能以 Channel Viewer 呈現，如圖六芳草之頻譜圖所示，可以看到一天 1 次的頻率特別強，而一天2 次、一天 3 次、一天4 次等頻率，可能為一天 1 次的倍頻。一天