[go: up one dir, main page]

TW200928303A - Observation method for sea-level wave fluctuation by using GPS - Google Patents

Observation method for sea-level wave fluctuation by using GPS Download PDF

Info

Publication number
TW200928303A
TW200928303A TW96148575A TW96148575A TW200928303A TW 200928303 A TW200928303 A TW 200928303A TW 96148575 A TW96148575 A TW 96148575A TW 96148575 A TW96148575 A TW 96148575A TW 200928303 A TW200928303 A TW 200928303A
Authority
TW
Taiwan
Prior art keywords
wave
gps
spectrum
axis
observation method
Prior art date
Application number
TW96148575A
Other languages
Chinese (zh)
Other versions
TWI356897B (en
Inventor
Jia-Jiun Gau
Yu-Wei Jang
Yan-Bin Lin
Ching-Ruei Lin
Dung-Jing Dung
Original Assignee
Jia-Jiun Gau
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jia-Jiun Gau filed Critical Jia-Jiun Gau
Priority to TW96148575A priority Critical patent/TW200928303A/en
Publication of TW200928303A publication Critical patent/TW200928303A/en
Application granted granted Critical
Publication of TWI356897B publication Critical patent/TWI356897B/zh

Links

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

This invention relates to an observation method for sea-level wave fluctuation by using a GPS, in which a GPS module is used to detect and obtain the wave movement speed variations on X-axis, Y-axis, and Z-axis, and convert the same to output a wave movement speed time series. An analyzing and processing module filters out the low frequency noise in the movement speed time series and leave only the actual wave movement speed time series. The movement speed time series is then processed by Fourier transform analysis in addition to the wave theory so as to calculate the wave data such as velocity spectrum, level spectrum, significant wave height, average cycle, direction spectrum, etc. Thus, this invention is not limited by constraint of GPS base stations and not interfered by low frequency noise. It can be used to observe sea-level wave fluctuation in any offshore area as long as a buoy used as the carrier is capable of moving with the wave particles. Thus, without limitation to the type of the buoy, this invention is capable of obtaining accurate information of sea-level wave variation.

Description

200928303 九、發明說明: 【發明所屬之技術領域】 本發明係為一種GPS海面波動觀測方法,尤指一種適用於海 上偵測水粒子運動速度,進而取得水面波動資料之方法。 【先前技術】 按’習知觀測並取得海域之相關波浪資訊,係利用浮體設有 觀測裝置而進行觀測。該觀測裝置係運用習知加速度器觀測波 浪’係依據波浪分析理論,分析加速度、東西向與南北向傾角時 序列資料而得到波浪資料。引用此理論時,浮體須能隨著波浪表 面運動’才得以正確量測傾角變化。惟此習知技術所使用之浮體 必須限定為碟型浮體才能使觀測效果達到最佳,而不能使用其它 型態之浮體。 再者’習知觀測技術係運用加速度器觀測波浪。惟在觀測時200928303 IX. Description of the invention: [Technical field to which the invention pertains] The present invention is a GPS sea surface fluctuation observation method, and more particularly to a method for detecting the velocity of water particles on the sea and obtaining water surface fluctuation data. [Prior Art] According to the conventional observation and acquisition of the relevant wave information of the sea area, observation is made by using the observation device of the floating body. The observation device uses the conventional accelerometer to observe the wave wave system based on the wave analysis theory, and analyzes the sequence data of the acceleration, east-west and north-south inclination to obtain the wave data. When citing this theory, the floating body must be able to accurately measure the change in inclination with the wave surface movement. However, the floating body used in the prior art must be limited to a dish-shaped floating body in order to optimize the observation effect, and other types of floating bodies cannot be used. Furthermore, the conventional observation technique uses an accelerometer to observe waves. Only when observing

該加速度器會產生雜訊,使該加速度器所輸出之加速度時序列數 據依據波/良为析理論計算出加速度譜,再經由頻率四次方倒數 轉換為水位譜’ _制波高、週麟料。啊加速度譜雜訊也 同時、.’里過頻率四次方倒數轉換後令低頻處雜訊被放大之倍率較 高頻處為大’進而得到錯誤的水位譜。 習知技術亦有彻RTK Gps量騎面水_化。當利用 RTK-GPS 時,使用 + # 用兩部以上之GPS接收器,並將其中一部設立在已 知位置之固疋點稱為GpS基準站,該GPS基準站以接收到的觀測 資料計算大氣誤# + 、產(電離層、對流層延遲)與執道誤差造成的影 5 200928303 .響’並透過地面域電錢信緖,將誤纽正值傳送u他各 部GPS接收器’即時解算出接收器與基準站位置之相對向量,求 :得高精度之位置資料。細·-哪在制時有其關,即當基 準站與接收器之間的距離增加大氣誤差影響的改變使得接收 器只能距離基準站約㈣公里才能達到公分級精度。因此,利用 RTK GPS測波H、%應用於沿海海域將此技術應用於深海區 域,量測精度將大幅下降。 〇 練於此,吾等發日狀為解決㈣細技賴纽之問題, 在經過不斷設計及改良後而有本發明之問世。 【發明内容】 爰是本發明提供-種GPS海面波動觀測方法,其目的在於使 设於浮體之觀測裝置不受GPS基準站之限制及低頻雜訊之干擾, 於任何海域進行海面波動觀測,且只需作為載具之浮體能隨波浪 水粒子運動,亦不限定浮體之型態,即可獲得準確的海面變動資 ❹訊。 為能實現上述之發明目的,本發明係為一種Gps海面波動觀 測方法,其方法包含:GPS模組偵測取得X軸向、γ轴向及z軸向 之波浪運動速度變化量,並轉換輸出波浪運動速度時序列,由分 析處理模組濾除運動速度時序列内含之低頻雜訊,保留真正屬於 波浪之運動速度時序列,並將該運動速度時序列應用傅立葉轉換 分析法,並配合波浪理論推算取得速度譜、一維水位譜、示性波 高、平均週期及方向波譜等波浪資料。 200928303 • *依據上述之構想,其中該GPS海峨動觀測方法更包含:速度 上曰維水位4、不性波南、平均週期及方向波譜等波浪資料經 及刀析處理模組編碼為資料字串,並將該資料字串儲存於儲存模 組。 ' “雜上叙縣’射_顧_為硬式剌機,或以F㈣ 己隱卡Compact F1 ash記憶卡為記錄媒體之讀卡機。 依據上述之縣,其愤GPS細波她測方法更包含·該分 〇析處雜組輸出該請字串給傳輸模組,由該傳輸單元以無線傳 輸方式將該資料字串傳送給遠端基地台。 依據上述之縣,其中該傳键_為無線電機、GSM數據機、 GPRS數據機或衛星傳輸裝置。 【實施方式】 、關於本發明之技術手段’ _—健實施麻合圖式於下文 進行詳細說明,俾供鈞上深入了解並認同本發明。 ❹ 首先’請參閱第一圖所示,係為本發曰月GPS海面波動觀測方 法較佳實施例示意圖。於本較佳實施例中,該觀測方法必須藉由 設於浮體200之觀測裝置⑽實施。該觀測裝置⑽係包含有奶 模組10、分析處理模組20、儲存模組3〇、傳輸模組4〇以及電力 模組50。其中: 該GPS模組1〇設有全球定位系統(G1〇bal p〇s出⑽加 System ’ GPS)。該GPS模組10與該分析處理模組2〇為電性連接, 用以當該浮體200在預定觀測海域漂浮時,即時偵測該浮體· 7 200928303 之X軸向、Y軸向及Z軸向之水粒子速度變化,並同時對應輸出該 三軸向之運動速度時序列給該分析處理模組20。前述該χ轴向在 本例中是預設為東西向’即GPS模組10偵測東西向之水粒子速度 變化。該Y轴向在本例中是預設為南北向,即GPS模組1〇偵測南 北向之水粒子速度變化。該Z軸向在本例中是預設為垂直於海平 面之方向,即GPS模組10偵測水位垂直速度變化。 ❹ ❹ 該分析處理模組20與該GPS模組1〇、儲存模組3〇、傳輸模 組40為電性連接。該分析處理模組2〇以分析處理軟體分析計算 由該GPS模組10傳送之運動速度時序列作為輸入值,經處理運算 後得出水粒子在X軸向、γ軸向及z軸向之速度譜、水位譜、示性 波高、平均週期及方向波譜等波浪資料,並將該等波浪資料經編 碼轉換為雜料,同_狀_輸触4()及雜存模組3〇。 該儲存模組30係接收該分析處理模組2〇輸出之資料字串, 並儲存該資料字串。該儲存模組3G係為硬式磁碟機,或以❿油 記憶卡、Conpaet Fiash記針為記錄魏之讀卡機。 該傳輸模組40係接收該分析處魏㈣輸出之資料字串。 該傳輸模組40即以無線傳輸方式傳送給遠端基地台。該傳輸模组 4〇係為無線電機、GSM數據機、GpRS數據機或衛星傳輸裝置。 /電力模組40與該GPS模組1Q、分析處理模組2()、儲存模 組3〇及傳輸模組4G為電性連接,用以提供該等模組1G、20、30 所需之電源。 8 200928303 • 如第二圖所示,當該觀測裝置⑽設於該浮體,並被置放 於預定觀測海域時,該浮體200在海面漂浮並隨著波浪水粒子運 :動,同時該分析處理模組20在預設之時間週期定時啟動該GPS模 組10與該傳輸触40 ’使該GPS模組1〇即時偵測該海域之水粒 子在X軸向、Y轴向及Z軸向之速度變化,並同時對應輸出該三轴 之運動速度時序列給該分析處理模組2〇進行分析運算。 如第三圖所示’該分析處理模组2()之分析步驟如下: ❹ 师—湖:齡析處理概20轉三㈣絲子運動速度 時序列。該GPS模組10偵測取得χ軸向、γ轴向及z轴向之水粒 子速度物理量’並傳送運動速度時序列給該分析處理模組2〇。該 分析處理模組20則以該等速度物理量做為運算處理之輸入值。 步驟二S200:該分析處理模組2〇在對三軸向水粒子運動速度 時序列進行運算前,會先行判斷是否有Gps訊號偏移(咖响The accelerometer generates noise, and the acceleration time sequence data output by the accelerometer is calculated according to the wave/good analysis theory, and then converted to the water level spectrum by the frequency quadratic reciprocal count _ _ wave height, Zhou Lin material. Ah, the acceleration spectrum noise is also at the same time. After the frequency is quadruple-reciprocal conversion, the noise at the low-frequency noise is amplified at a higher frequency than at the high frequency, and the wrong water level spectrum is obtained. The conventional technology also has a full-scale RTK Gps ride. When using RTK-GPS, use + # with more than two GPS receivers, and set one of the fixed points in the known position as the GpS base station, which is calculated based on the received observation data. Atmospheric error # + , production (ionosphere, tropospheric delay) and the impact of the obstacle error 5 200928303 . Ringing 'through the ground domain money money letter, will send the wrong value of the positive value of his GPS receivers' instant calculation The relative vector of the position of the device and the base station is obtained: the position data with high precision is obtained. Fine ·-Whether it is in the system, that is, when the distance between the reference station and the receiver increases, the influence of the atmospheric error makes the receiver only about (four) kilometers away from the reference station to achieve the public classification accuracy. Therefore, using RTK GPS to measure H and % in coastal waters, this technology is applied to deep sea areas, and the measurement accuracy will be greatly reduced.练 In this practice, we have solved the problem of (4) fine technology, and after the continuous design and improvement, the invention has been made. SUMMARY OF THE INVENTION The present invention provides a GPS sea surface fluctuation observation method, the purpose of which is to make the observation device provided on the floating body not subject to the limitation of the GPS reference station and the interference of the low frequency noise, and perform sea surface fluctuation observation in any sea area. And only the floating body of the vehicle can move with the wave water particles, and the type of the floating body is not limited, so that an accurate sea surface change information can be obtained. In order to achieve the above object, the present invention is a GPS wave surface observation method, which comprises: the GPS module detects the amount of change of the wave motion speed in the X axis, the γ axis and the z axis, and converts the output. The wave motion speed time sequence is used by the analysis processing module to filter out the low frequency noise contained in the sequence of the motion speed, retaining the sequence of the motion speed truly belonging to the wave, and applying the Fourier transform analysis method to the motion speed time series, and matching the wave Theoretical calculations obtain wave data such as velocity spectrum, one-dimensional water level spectrum, indicative wave height, average period and direction spectrum. 200928303 • *According to the above concept, the GPS sea motion observation method further includes: the wave data such as the water level 4, the irrational Bonan, the average period and the direction spectrum, and the knives are processed into the data words. String and store the data string in the storage module. ' Miscellaneous Syrian County' shoot_Gu_ is a hard-type machine, or F (four) Hidden card Compact F1 ash memory card as a recording medium card reader. According to the above-mentioned county, its infringing GPS fine wave her test method includes The sub-declining processing group outputs the request string to the transmission module, and the transmission unit transmits the data string to the remote base station by wireless transmission. According to the above-mentioned county, the transmission key is a radio GSM data machine, GPRS data machine or satellite transmission device. [Embodiment] The technical means of the present invention is described in detail below, and the present invention is deeply understood and recognized. ❹ First, please refer to the first figure, which is a schematic diagram of a preferred embodiment of the GPS sea surface fluctuation observation method. In the preferred embodiment, the observation method must be provided by the observation device provided on the floating body 200. (10) The observation device (10) comprises a milk module 10, an analysis processing module 20, a storage module 3A, a transmission module 4A, and a power module 50. wherein: the GPS module 1 has global positioning System (G1〇bal p〇s out (10) Adding a system 'GPS'. The GPS module 10 is electrically connected to the analysis processing module 2, and is used to instantly detect the floating body when the floating body 200 floats in a predetermined observation sea area. The water particle velocity changes in the Y direction, the Y axis, and the Z axis, and simultaneously outputs the sequence to the analysis processing module 20 corresponding to the output of the three axial motion speed. The aforementioned axis is preset to something in this example. The GPS module 10 detects the change of the water particle velocity in the east-west direction. The Y-axis is preset to be north-south in this example, that is, the GPS module detects the change of the water particle velocity in the north-south direction. The axial direction is preset to be perpendicular to the sea level in this example, that is, the GPS module 10 detects the vertical speed change of the water level. ❹ ❹ The analysis processing module 20 and the GPS module 1 and the storage module 3 The transmission module 40 is electrically connected. The analysis processing module 2 calculates, by using the analysis processing software, a sequence of motion speeds transmitted by the GPS module 10 as an input value, and after processing, the water particles are obtained in the X axis. , γ-axis and z-axis velocity spectrum, water level spectrum, indicative wave height, average period and Wave data such as a spectrum is encoded, and the wave data is encoded into a miscellaneous material, which is the same as the _ _ _ touch 4 () and the memory module 3 〇. The storage module 30 receives the analysis processing module 2 〇 The output data string is stored, and the data string is stored. The storage module 3G is a hard disk drive, or the oil memory card and the Conpaet Fiash pin are used to record the Wei card reader. Receiving the data string of the output of the analysis (four). The transmission module 40 is transmitted to the remote base station by wireless transmission. The transmission module 4 is a radio, a GSM data machine, a GpRS data machine or a satellite transmission. The power module 40 is electrically connected to the GPS module 1Q, the analysis processing module 2 (), the storage module 3A, and the transmission module 4G for providing the modules 1G, 20, 30 Power supply required. 8 200928303 • As shown in the second figure, when the observation device (10) is placed on the floating body and placed in a predetermined observation sea area, the floating body 200 floats on the sea surface and moves with the wave water particles, and The analysis processing module 20 activates the GPS module 10 and the transmission contact 40' at a preset time period to enable the GPS module 1 to instantly detect the water particles in the sea area in the X axis, the Y axis, and the Z axis. To the speed of the change, and simultaneously output the sequence of the three-axis motion speed to the analysis processing module 2 〇 analysis operation. As shown in the third figure, the analysis steps of the analysis processing module 2 () are as follows: ❹ Division-lake: age analysis processing 20 turns three (four) silk movement speed time sequence. The GPS module 10 detects the physical quantity of water particle velocity obtained in the x-axis, the γ-axis, and the z-axis, and transmits the sequence of motion speed to the analysis processing module 2〇. The analysis processing module 20 uses the velocity physical quantities as input values for the arithmetic processing. Step 2: S200: The analysis processing module 2 first determines whether there is a GPS signal offset before the sequence of the three-axis water particle motion speed is calculated.

Loss)。所謂GPS訊號偏移⑽L〇ss),指當啦衛星發 生偏移或訊號無法接收時,GPS模組1G所輸出之水粒子運動速度Loss). The so-called GPS signal offset (10) L〇ss) refers to the speed of water particle motion output by the GPS module 1G when the satellite is offset or the signal cannot be received.

時序列則會㈣-整段往上餘下偏移的現象⑽sigJThe time sequence will be (4) - the whole segment is shifted upwards (10) sigJ

Loss)。當該分析處理模組2G運用移動平均法_判斷運動速度 時序列有GPS訊號偏移時’則進人遽除Gps訊號偏移之程序。反 之’則該分析處賴組20即對該等運動時序列進行波浪資訊之分 析運算。 ° 步驟三S300 :當該分析處理模、组2〇判斷運動速度時序列有 9 200928303 GPS訊號偏移時’則該分析處理模組2〇進一步將其滤除,而保留 真正屬於波浪速度的減’即針對該訊缝行波浪資訊之分析運 算。 步驟四S400 :於此步驟中,該分析處理模組2〇對三抽向之運 動速度時序舰·處理取得速度譜。·波浪巾雜的變動具 定常性(伽i〇nary)、統-性(erg〇dic)及高斯過程(㈣如 process) ’則其時間上與空間上的統計性n將波形魏為 無限多個獨振幅、辭、她及人概向的餘弦函數疊加來表 不,如下式: + _於2約 +〜} (S4-1) 其中’ 為波形函數’此函數與位置及時間有關,乂為 〇 董份波頻率、_應於,的波數,〜為人射波的波向H為 :應於頻率及波向的振幅及相位。當本較佳實施例之浮體200完 、跟隨海面移動標浮,則該觀測裝置1〇〇所觀測Z轴向之速紅可 以由(S4-1)式對時間微分1次而得: “从〇丨巧=。=笔客(2办一11(2約+〜) (S4-2) W)之自相關函數(aut0—correlati〇n functi〇n)定義如下 式所示: \Vz{t)'Vz{t + r)dT (S4-3)Loss). When the analysis processing module 2G uses the moving average method to determine the motion speed, the sequence has a GPS signal offset, then the program for removing the GPS signal offset is entered. In contrast, the analysis group 20 performs the analysis of the wave information for the motion time series. ° Step 3: S300: When the analysis processing mode, group 2〇 determines the motion speed, the sequence has 9 200928303 GPS signal offset, then the analysis processing module 2 further filters it out, and retains the true wave speed reduction. 'The analysis operation of the wave information for the seam. Step 4: S400: In this step, the analysis processing module 2 obtains a velocity spectrum for the three-way moving speed time series ship processing. · The variation of the wave towel has a regularity (gamma〇nary), an erg-dic (erg〇dic), and a Gaussian process ((4) such as process). Then its temporal and spatial statistics n will be infinitely large. The cosine function of the unique amplitude, the word, her and the person's general direction is superimposed, as shown below: + _ at 2 about +~} (S4-1) where 'is a waveform function' This function is related to position and time, 乂For the wave frequency of Dong Mobo, _ should be, the wave direction of the human wave is the amplitude and phase of the frequency and direction. When the floating body 200 of the preferred embodiment is finished and follows the sea surface movement, the speed of the Z-axis observed by the observation device 1 can be differentiated from the time by (S4-1): The autocorrelation function (aut0-correlati〇n functi〇n) from 〇丨巧=.=笔客(2 一一11(2约+〜)(S4-2) W) is defined as follows: \Vz{ t) 'Vz{t + r)dT (S4-3)

—T 其中7為時距’則波譜密度函數可由自相關函數的傅立葉轉換 200928303 求得,其為波浪在某頻率範_所含之能量,如下式表示: ^f(/)=2£/?(r)exp(-/2^v)i/r ( S4-4 ) 其中/為鮮,5v<y)^ z ,z 上轉換函數(¥2)可得-維水位譜办),如下式所示: S(f)=Sv(fh(2^ c、 ν〇4-〇) Ο 〇 步驟五S5〇〇 :於此步驟中,該分析處理模組20對三轴向之運 動速度時序列經運算處理取得方向分布函數,進而得到方向波 譜。交錯譜之定義為交錯相關函數(c簡c〇rrelatic)n Fusion ) 之傅立葉轉換’交錯相關函數Α(Γ)的定義如下: R^ = ^J^(^ + T)dt (S5]) 上式中,冰)代表波形函數,Γ代表稽延時間,則交錯譜可由 下式獲得: (S5-2) (S5-3) 。藉由交錯譜—T where 7 is the time interval' and the spectral density function can be obtained from the Fourier transform 200928303 of the autocorrelation function, which is the energy of the wave in a certain frequency range, expressed as follows: ^f(/)=2£/? (r)exp(-/2^v)i/r ( S4-4 ) where / is fresh, 5v<y)^ z , z upconversion function (¥2) is available - dimension water level spectrum), as follows Shown: S(f)=Sv(fh(2^ c, ν〇4-〇) Ο 〇Step 5 S5〇〇: In this step, the analysis processing module 20 pairs the three-axis motion velocity time series The direction distribution function is obtained by the operation processing, and the direction spectrum is obtained. The definition of the interlaced spectrum is the interleaved correlation function (c simple c〇rrelatic) n Fusion ) The Fourier transform 'interlace correlation function Α(Γ) is defined as follows: R^ = ^ J^(^ + T)dt (S5)) In the above formula, ice) represents the waveform function, and Γ represents the delay time, and the interleaved spectrum is obtained by: (S5-2) (S5-3). Interlaced spectrum

<^(/)= LR^~2mfTdT 將交錯譜以複數形式表示為如下所示: ^(/)=c,(/)-^(/) 其中,d/)稱為同位譜,込⑺稱為轉向譜 與方向波譜具有如下式所示之關係·· _=P.S(/加 (S5-4) 當該GPS模組10朗該浮體·位於海面定點位置之χ轴 向、γ軸向及Ζ軸向之水粒子運動速度,則交錯譜可為下式表示· 200928303<^(/)= LR^~2mfTdT The staggered spectrum is expressed as a complex number as follows: ^(/)=c,(/)-^(/) where d/) is called the isotope spectrum, 込(7) The steering spectrum and the directional spectrum have the relationship shown in the following equation: _=PS (/plus (S5-4) When the GPS module 10 is the floating body, the χ axial direction at the fixed position of the sea surface, γ axis And the velocity of the water particles in the axial direction, the interlaced spectrum can be expressed as follows: 200928303

Ο ΦίΜ)= _π (S5-5) 當該GPS模組10觀測波面之χ轴向、γ轴向及ζ轴向之運動 速度時,係可彳⑽六組交錯譜,並將六組故譜分職人(Μ) 式後,再財限傅㈣級數法,假設方向分布函數可以展開成有 限階數的傅利葉級數,如下列所示: s{f, 〇>a0{f)^an (/) cos ηθ + bn (/) sin n〇 (S5-6) 再以GPS模組1G觀測所得之三轴向運動速度時序列代入 (S5-3)式解至二階,各項之係數分別以下式計算: a(n- Cn(f) (· (S5-7a) a(f)- Qn(f) Αη~{2π[)27ώ (S5-7b) a (f\- ^22(/)- ^33(/) (2φ)2π (S5-7c) (S5-7d) w)--i2C£l , , cosh Λζ/ , 其中hsinhW W為水深4以下式計算: (S5-7e) i / Q? (/)+^33(/) Ί cu(f) (S5-8) 方向譜可視為一維能譜與方向分佈函數之内積. S(/^) = S(/)-D(e\f) (S5-9) 12 200928303 又 其中町)為—輯譜;棒I/)為波絲量的方向分佈函數 町)為如下式之表示: (S5~i〇) (S5-11) &(/)代表 S(f) = i S{Le)de = Cu{f)K27tff=aQ{f)n 又乃州/)為如下式之表示: i{2+ri^Cos^~(^'^(/))]^2(/)cos^-(2<9-<92(/))]| 、中以’)、巧(/)代表波浪方向分佈寬度。巧(乃、 波浪傳遞出去的方向。 步驟六S600 : 昇取得水位譜後, 如下列式子: 於此步驟中,該分析處理模組20利用步驟四運 再由該水位譜推算示性波高Hs與平均週期t, ^ - 4.004Ο ΦίΜ)= _π (S5-5) When the GPS module 10 observes the velocity of the 波, γ, and ζ axes of the wavefront, it can 彳(10) six sets of staggered spectra, and six sets of spectrum After the splitter (Μ), the financial limit (four) series method is assumed, and the direction distribution function can be expanded into a finite order Fourier series, as shown below: s{f, 〇>a0{f)^an (/) cos ηθ + bn (/) sin n〇(S5-6) Then, the three-axis motion velocity obtained by GPS module 1G is substituted into the second-order solution (S5-3), and the coefficients of each are respectively The following formula is calculated: a(n- Cn(f) (· (S5-7a) a(f)- Qn(f) Αη~{2π[)27ώ (S5-7b) a (f\- ^22(/) - ^33(/) (2φ)2π (S5-7c) (S5-7d) w)--i2C£l , , cosh Λζ/ , where hsinhW W is the water depth of 4 or less: (S5-7e) i / Q? (/)+^33(/) Ί cu(f) (S5-8) The directional spectrum can be regarded as the inner product of the one-dimensional energy spectrum and the direction distribution function. S(/^) = S(/)-D(e (f5~i〇) (S5-11) &(/) stands for S(f) = i S{Le)de = Cu{f)K27tff=aQ{f)n /) is expressed as follows: i{2+ri^Cos^~(^'^(/))]^2(/)cos^-(2<9-<92(/))]| Use '), smart (/) to represent the width of the wave direction distribution. Ingeniously, the direction in which the wave is transmitted. Step 6 S600: After obtaining the water level spectrum, as in the following formula: In this step, the analysis processing module 20 uses the step four to calculate the indicative wave height Hs from the water level spectrum. With an average period t, ^ - 4.004

Tz =V«v^r (S6'2)Tz =V«v^r (S6'2)

、’卜。與%分別為水位譜能量之零次矩與二次矩。另外, 量密度最大處所對應之方向 驟五運算得財肖賴後,㈣倾 , 即為主波向。 詳虽該分析處理模組2〇經由上述之程序運算出速度譜、水位 °曰、不性波高、平均聊及方向波譜等波浪資料後,該分析處理 =且20即將將前述之波浪資料進行編碼為資料字_,並將該資料 字串同時輸出於該儲存模組3Q儲存備查及該傳輸模組4Q。該傳輸 模、、且40接收到資料字串後,即以無線傳輸方式將資料字串傳送給 13 200928303 综上所述可得知,本_具打_述之優點, 'Bu. And % are the zeroth moment and the second moment of the water level spectrum energy, respectively. In addition, the direction corresponding to the largest volume density is calculated after the fifth calculation, and (4) is the main wave direction. After the analysis processing module 2 calculates the wave data such as the velocity spectrum, the water level °, the irregular wave height, the average chat and the direction spectrum through the above program, the analysis processing = and 20 will encode the aforementioned wave data. It is a data word _, and the data string is simultaneously outputted to the storage module 3Q for storage and the transmission module 4Q. After the transmission mode, and 40 receives the data string, the data string is transmitted to the wireless transmission mode. 13 200928303 In summary, it can be known that the advantages of the present _

广哪海面波動觀測方法係運用奶模組分析取得X轴 ^西向)γ轴向(南柏)及2軸向(上下向)之水粒子運 =度時相資料1剌波浪㈣。只科體麟著波浪水粒 GP^l卩,可柯<讀_之限制。再者’本個之方法藉由 、.且並配口以波浪理論為基礎之分析方法而無需架設奶 土準站’即能精確取得波浪資訊,且亦無低頻雜訊之干擾問題。 綜觀上述’本發明就其整體之組合及特徵而言,既未曾見諸 於同類產品中,申請前絲見公開,誠已符合專利法之法定要件, 茇依法提出發明專利申請。 _ X_L所述者,僅為本發明讀佳實細,當不能以此限定本 =月實^之㈣,即大凡依本發明申請專利朗及發明說明書内 谷所作之等效變化與修飾,皆應域本發明專利之涵蓋範圍内。 【圖式簡單說明】 第-圖係本發明GPS海面波峨财法較佳實顧之觀測裝置 架構方塊圖。說明觀測裝置之架構關係。 $二_本發明GPS海面波動觀測方法較佳實施例之觀測裝置 設於浮體’而浮體漂浮贿蚊錢㈣示意圖。 第二圖係本發明GPS海面波動觀測方法較佳實施例之分析處理 200928303 • 模組之分析處理流程圖。 【主要元件符號說明】 100—觀測裝置 10 GPS模組 20 分析處理模組 30--- -一儲存單元 40------傳輸模組The wide sea surface fluctuation observation method uses the milk module analysis to obtain the X-axis ^west direction) γ axis (Nanbai) and 2 axial (up and down) water particle transport = degree phase data 1 剌 wave (4). Only the body body with the wavy water particles GP^l卩, Ke Ke < read _ restrictions. Furthermore, the method of this method is based on the wave theory-based analysis method without the need to set up a milky quasi-station to accurately obtain wave information, and there is no interference problem of low-frequency noise. In view of the above-mentioned invention, the combination and characteristics of the present invention have not been seen in the same kind of products, and the disclosure before the application is open, and the company has already complied with the statutory requirements of the patent law, and has filed an invention patent application according to law. _ X_L is only for the purpose of reading the present invention. If it is not possible to limit the equivalent of the present invention, the equivalent change and modification of the patent application and the invention. Within the scope of the invention patent. [Simple diagram of the diagram] The first diagram is a block diagram of the observation apparatus of the GPS sea surface wave method of the present invention. Explain the architectural relationship of the observing device. $二_ The observation device of the preferred embodiment of the GPS sea surface fluctuation observation method of the present invention is set in the floating body ‘the floating body floats the mosquito money (4). The second figure is the analysis and processing of the preferred embodiment of the GPS sea surface fluctuation observation method of the present invention. 200928303 • The analysis and processing flow chart of the module. [Main component symbol description] 100—observation device 10 GPS module 20 analysis processing module 30----one storage unit 40------transmission module

〇 15〇 15

Claims (1)

200928303 十、申請專利範圍: 1. 一種GPS海面波動觀測方法,其方法包含·· GPS模組偵測取得X軸向、γ轴向及%軸向之波浪運動速度變化 量’並轉才奐輪出波浪運動速度時序列,&分析處理模组瀘除運動 速度時序列内含之低頻雜訊,保留真正屬於波浪之運動速度時序 列,並將該運動速度時序列應用傅立葉轉換分析法,並配合波浪 理論推算取得速度譜、一維水位譜、示性波高平均週期及方向 ❹波譜等波浪資料。 2.如申明專利範圍第1項所述之Gps海面波動觀測方法,其中該 GPS海面波動觀測方法更包含: 欠速度谱、一維水位譜、示性波高、平均週期及方向波譜等波 浪身料經該分析處理模組編碼為資料字串,並將該資料字串儲存 於儲存模組。 q :如申明專利範圍第2項所述之GPS海面波動觀測方法,其中該 堵存模組係為硬式磁碟機,或以Flash記憶卡、㈣沾F馳記 憶卡為記錄媒體之讀卡機。 、4.如申請專利範圍第2項所述之GPS海面波動觀測方法,其中該 S海面波動觀測方法更包含: • 該分析處理概輸出該資料字串給傳輸模紐,由該傳輸單元以 …、線傳輪方式將該資料字串傳送給遠端基地台。 申请專利範圍第4項所述之Gp s海面波動觀測方法,其中該 200928303 ’ 傳輸模組係為無線電機、GSM數據機、GPRS數據機或衛星傳輸裝置。 〇200928303 X. Patent application scope: 1. A GPS sea surface wave observation method, the method includes: · GPS module detects the X-axis, γ-axis and %-axis wave motion speed change amount' and turns the wheel When the wave motion speed sequence is used, the & analysis processing module removes the low frequency noise contained in the sequence when the motion speed is removed, retains the sequence of the motion speed truly belonging to the wave, and applies the Fourier transform analysis method to the motion speed time series, and According to the wave theory, the wave data such as velocity spectrum, one-dimensional water level spectrum, high-order periodic wave and directional chopping spectrum are obtained. 2. The Gps sea surface wave observation method described in claim 1 of the patent scope, wherein the GPS sea surface wave observation method further comprises: a wave body material such as an underspeed spectrum, a one-dimensional water level spectrum, an indicative wave height, an average period, and a direction spectrum. The analysis processing module is encoded into a data string, and the data string is stored in the storage module. q: The GPS sea surface fluctuation observation method described in claim 2, wherein the storage module is a hard disk drive, or a card reader with a flash memory card and (4) a memory card as a recording medium. . 4. The GPS sea surface fluctuation observation method described in claim 2, wherein the S sea surface fluctuation observation method further comprises: • the analysis processing outputs the data string to the transmission module, and the transmission unit ... The line-passing mode transmits the data string to the remote base station. The Gp s sea surface fluctuation observation method described in claim 4, wherein the 200928303 ’ transmission module is a radio, a GSM data machine, a GPRS data machine or a satellite transmission device. 〇 1717
TW96148575A 2007-12-19 2007-12-19 Observation method for sea-level wave fluctuation by using GPS TW200928303A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW96148575A TW200928303A (en) 2007-12-19 2007-12-19 Observation method for sea-level wave fluctuation by using GPS

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW96148575A TW200928303A (en) 2007-12-19 2007-12-19 Observation method for sea-level wave fluctuation by using GPS

Publications (2)

Publication Number Publication Date
TW200928303A true TW200928303A (en) 2009-07-01
TWI356897B TWI356897B (en) 2012-01-21

Family

ID=44863989

Family Applications (1)

Application Number Title Priority Date Filing Date
TW96148575A TW200928303A (en) 2007-12-19 2007-12-19 Observation method for sea-level wave fluctuation by using GPS

Country Status (1)

Country Link
TW (1) TW200928303A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013003081A1 (en) * 2011-06-26 2013-01-03 The Regents Of The University Of California An inexpensive instrument for measuring wave exposure and water velocity
TWI486557B (en) * 2014-03-28 2015-06-01 Cing Jhe Huang Methods and systems for monitoring the water surface elevations of the ocean
CN106597010A (en) * 2017-02-14 2017-04-26 珠江水利委员会珠江水利科学研究院 GPS-based wave direction measurement device and method

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013003081A1 (en) * 2011-06-26 2013-01-03 The Regents Of The University Of California An inexpensive instrument for measuring wave exposure and water velocity
TWI486557B (en) * 2014-03-28 2015-06-01 Cing Jhe Huang Methods and systems for monitoring the water surface elevations of the ocean
CN106597010A (en) * 2017-02-14 2017-04-26 珠江水利委员会珠江水利科学研究院 GPS-based wave direction measurement device and method
CN106597010B (en) * 2017-02-14 2023-02-14 珠江水利委员会珠江水利科学研究院 Wave direction measuring device and method based on GPS

Also Published As

Publication number Publication date
TWI356897B (en) 2012-01-21

Similar Documents

Publication Publication Date Title
CN111896984A (en) A real-time high-precision wave measurement method and device based on GNSS
Lipphardt Jr et al. Blending HF radar and model velocities in Monterey Bay through normal mode analysis
US7613072B2 (en) System and method for determining directional and non-directional fluid wave and current measurements
CN108469620B (en) Underwater terrain measurement method suitable for shallow water sea area of radiation sand ridge group
JP2003302221A (en) Gps wave height, current direction and current speed measuring device and gps wave height, current direction and current speed measuring system
Trifunac et al. A note on the accuracy of computed ground displacements from strong-motion accelerograms
Fong et al. Evaluation of AUV‐based ADCP measurements
CN115451921B (en) A method, device, medium and terminal for synchronous observation of waves, tides and currents in the entire ocean area
CN107063269A (en) Build the inverse distance weighted interpolation algorithm of subaqueous gravity assisting navigation ambient field
Elhoushi et al. Online motion mode recognition for portable navigation using low‐cost sensors
CN103605167A (en) Mallat algorithm-based marine gravity measurement error eliminating method
CN108415096A (en) Subaqueous gravity gradient object detection method based on Newton iteration method
CN106908044A (en) A kind of full marine site tidal observation system and tidal observation method based on the Big Dipper
TW200928303A (en) Observation method for sea-level wave fluctuation by using GPS
CN103529451B (en) Method for calibrating coordinate position of seabed transponder of water-surface mother ship
Beardsley A comparison of the vector‐averaging current meter and new Edgerton, Germeshausen, and Grier, Inc., vector‐measuring current meter on a surface mooring in Coastal Ocean Dynamics Experiment 1
US7352651B2 (en) System and method for determining directional and non-directional fluid wave and current measurements
Sweeney et al. Centimeter-level positioning of seafloor acoustic transponders from a deeply-towed interrogator
Yoerger et al. Drag forces and flow-induced vibrations of a long vertical tow cable—Part I: Steady-state towing conditions
CN112557060B (en) Dynamic comfort-oriented human body three-dimensional sitting posture change detection system, method and terminal
Fu Ocean circulation and variability from satellite altimetry
CN119620120B (en) Real-time tide measuring method and system for single Beidou altimeter buoy without base station
Das et al. Crustal thickness estimation using high frequency Rayleigh waves
CN110887506A (en) Motion amplitude detection method and system of inertial sensor influenced by sea waves
Huang et al. Monitoring High-frequency Ocean signals using low-cost GNSS/IMU buoys

Legal Events

Date Code Title Description
MM4A Annulment or lapse of patent due to non-payment of fees