270 likes | 457 Views
GPS 連續觀測資料分析. 指導老師 : 余水倍 學生 : 莊舒雲 林幸瑩. 內容提要. 1.GPS 基本概念 : 組成系統 定位原理 2.GPS 時間序列的雜訊分析 為何要分析時間序列上的雜訊 ? 雜訊的種類、來源與特性 3.GPS 連續觀測資料之分析方法及步驟 選擇雜訊模式 建立模型 模型 參數估計與資料分析 4. 未來目標. GPS基本觀念. A 。組成系統 1. 現有24+3顆衛星。 2. 有六個運行軌道,每一軌道上有4顆衛星在運轉 ﹔ 軌道的平均高度為 20,200km ,運行週期則為 11 小時 58 分。
E N D
GPS連續觀測資料分析 指導老師: 余水倍 學生: 莊舒雲 林幸瑩
內容提要 1.GPS基本概念: 組成系統 定位原理 2.GPS時間序列的雜訊分析 為何要分析時間序列上的雜訊? 雜訊的種類、來源與特性 3.GPS連續觀測資料之分析方法及步驟 選擇雜訊模式 建立模型 模型參數估計與資料分析 4.未來目標
GPS基本觀念 A。組成系統 1.現有24+3顆衛星。 2.有六個運行軌道,每一軌道上有4顆衛星在運轉﹔ 軌道的平均高度為20,200km,運行週期則為11小時58分。 3.地球上任何角落至少可同時接收到4顆衛星, 最多可達11顆,平均為5~6顆。 4.時間系統:UTC國際標準時間。
各種掌上型 GPS信號接收機
GPS基本觀念 B。定位原理 1.絕對定位:或稱單點定位,在地球協議座標系中, 直接確定觀測站相對座標系原點, 也就是地球質心的位置。 2.相對定位﹔於地球協議座標系中,確定觀測站 與某一地面參考點之間的相對位置。
絕對定位 相對定位
GPS連續觀測之目的 1.能夠獲得完整的連續時間資料點座標,進而分析時間序列上的相關性。 2.找出連續觀測產生的干擾誤差類型,進一步消除誤差干擾,得到更精確的時間資料點。 3.提供許多震前、同震及震後變形資料,讓我們更加了解地震的震源特性、地殼的應變累積與能量釋放過程。
為何我們要分析時間序列上的雜訊? 了解GPS時間序列上的雜訊特性,以正確估算各項參數的誤差值。 不同地點受到的干擾程度不同。 在熱帶地區的測站,雜訊會比中緯度地區來的大。
雜訊的種類、來源與特性 1.通常可分為自然雜訊和人為雜訊。 2.自然雜訊來源: 測站的地質特性、隨季節性雨量多寡或地下水位改變和大氣干擾等所造成的座標變化。 人為雜訊來源: 測站天線架不穩定、更換天線造成的系統誤差,或是衛星軌道、大氣折射、天線中心漂移以及參考座標系等無法以模型完全去除的影響。
3.雜訊的特性: GPS雜訊的頻譜特性,可以用以下的式子來表示: P( f )=P0( f / f0 )κPower-law process 如果將兩邊取對數值並作圖可以得到: Power spectrum slopeκ
根據作圖結果我們可以依斜率κ值大小將雜訊模式分為三大類:根據作圖結果我們可以依斜率κ值大小將雜訊模式分為三大類: κ= 0 white noise κ= -1 flicker noise κ= -2 random walk noise (此波動模式又可稱為Brownian motion) 通常在長期觀測的自然情況下發生的干擾多屬低頻干擾,且其κ值在 -3 ~ -1之間,根據κ值的範圍我們可以知道雜訊特性。
GPS連續觀測之分析方法及步驟 • 選擇雜訊模式 • 建立模型 • 模型參數估計與資料分析
雜訊模式的選擇 • 近來研究利用Maximum-likelihood estimator (MLE) 去計算統計上述三種雜訊出現在時間序列上的殘差振幅之後發現最好的雜訊模式為white + flicker noise model ,它會將”機率函數(probability function)”最大化。 • Probability function: or
模型參數 • 各測站在各方向上觀測到的運動位移可以用下式表示: H : Heaviside function (階梯函數) 以下例子用的模型參數:
北部測站 BANC FCWS FIVE HNSN HSIN HUAL ILAN NSHE PEPU PLAN SANI SUAO THAI YMSM S101 南部測站 • CHIA • PKGM • S103 • W029 • YUSN • W030 • CHYN • CHNL • KULN • S011 • S092 • HOKN • S106 • S012 • GAIS • S169 • TSHO • CK01 • S23R 澎湖
P H Z Y φ λ X Y X 模擬操作與分析 1.座標轉換。 • 首先,我們將921地震之後每天觀測的43個測站資料點座標從笛卡兒座標系統XYZ轉成以S01R(澎湖)為參考主站的NEU座標系統。
測站CHIA 尚未去除分散點 outliers 2.將各方向觀測資料點對時間序列作圖。 Position (mm) Years
3.去除分散點/異常點(outliers)。 • 為了得到正確震後位移,我們必須進行消除誤差的動作,這其中包括去除觀測期間同震位移、影響資料點分析的分散點和更換測站天線所造成的偏移量等。 Coseismic offset aftershock
aftershock 殘差值越小越集中表示資料點可信度越高
4.去除同震位移(coseismic offset)得到smooth的曲線得到震後位移(postseismic displacement) Hanging wall Epochs of two major aftershocks
Epochs of two major aftershocks Footwall
未來工作 • 以較完整的模式分區(南北兩區)分析921地震前後多年資料,以了解震前變形、同震位移、震後變形及地震前後地殼運動速度的變化。 • 未來展望-藉由長時間更高精度的GPS連續觀測資料分析,偵測可能的地震前兆,探討地震預測的可行性。