【技術實現步驟摘要】
本專利技術涉及地震觀測領域，尤其涉及一種基于eemd的井水位固體潮特征提取方法。
技術介紹
1、地球上任意一點受到來自日、月引力以及地球繞地日(地月)公共質心旋轉而產生的慣性離心力的合力即為引潮力。在引潮力作用下，地球隨之產生周期性形變。
2、井-含水層系統對固體潮動態響應現象指地殼巖體在引潮力作用下，產生周期性彈性變形(膨脹和壓縮)，巖體孔隙中的地下水產生周期性變化，這種現象稱為固體潮效應。井水位變化形態與固體潮變化形態相似，呈同向變化，時間上滯后于固體潮，滯后時間因井而異。井水位對固體潮的響應實質上反映了含水層的水力特征及地下水流動運動狀態，故通過分析井水位的潮汐特征可以達到推求含水層水力參數、判定地下水流運動特征的目的，從而為解釋地震應力對井孔-含水層系統的影響提供依據。
3、潮汐因子用來表征井-含水層的介質彈性參數，可用于研究地區地殼彈性變異。目前，潮汐因子通常通過調和分析和小波分析等手段對其進行計算。常規的潮汐因子計算方法在處理數據前需要進行預處理，去除明顯錯誤數據，對缺失數據進行插值處理，然后采樣低通濾波提取趨勢變化，用源數據減去低通濾波的趨勢變化得到僅保留潮汐成分的水位，利用一元線性回歸計算得到潮汐因子。然而，傳統計算方法僅適用于井水觀測數據質量較好的觀測井，當井水位數據缺失或達不到觀測要求時，獲得的效果較差。
技術實現思路
1、本專利技術針對現有技術的缺陷，目的是提出一種基于eemd的井水位固體潮特征提取方法，通過集合經驗模態分解(ensem
2、為了實現上述目的，本專利技術提供如下技術方案：
3、一種基于eemd的井水位固體潮特征提取方法，包括以下步驟：
4、s1、對小時值井水位數據進行eemd分解，獲得8～12個imf分量和1個趨勢項r，固體潮引起井水位動態變化特征包含在井水位高頻分量imf1～3中；
5、s2、對井水位高頻分量imf1～3相加重構，即得到固體潮引起的井水位動態變化特征；
6、s3、對井水位高頻分量imf1～3相加重構信號與固體潮理論值做回歸分析，即得到潮汐因子。
7、進一步地，步驟s1中，對小時值井水位數據進行eemd分解時，引入高斯噪聲標準差設置為0.1，迭代次數為1000次。
8、進一步地，步驟s1中，引起井水位動態變化的固體潮頻率分布在2.3×10-5～7.9×10-6hz，主要由月球和太陽引起，全日波周期分別是25.82h和23.93h，半日波周期分別是12.42h和12.00h；引起的承壓井水位變化最大為20cm，正常情況下，井水位呈雙峰雙谷變化。
9、與現有技術相比，本專利技術的有益效果為：
10、本專利技術提供的一種基于eemd的井水位固體潮特征提取方法，能夠對非平穩、非線性的井水位信號進行自適應的分解，可用于水位數據缺失和其它原因不滿足進行調和分析或準調和分析的要求時提取固體潮引起的井水位變化信息，與傳統的調和分析方法相比，區別于傅里葉變換集合和小波變換以數學函數為基，eemd方法是基于自身信號獲得的基，能更客觀、更真實反映非線性、非平穩的物理過程，本專利技術通過對井水位觀測數據時值eemd分解后進行高頻分量重構，再進行回歸分析即可獲得潮汐因子，方便快捷，且可用于觀測質量較差的數據。
本文檔來自技高網...【技術保護點】
1.一種基于EEMD的井水位固體潮特征提取方法，其特征在于，包括以下步驟：
2.根據權利要求1所述的基于EEMD的井水位固體潮特征提取方法，其特征在于，步驟S1中，對小時值井水位數據進行EEMD分解時，引入高斯噪聲標準差設置為0.1，迭代次數為1000次。
3.根據權利要求1所述的基于EEMD的井水位固體潮特征提取方法，其特征在于，步驟S1中，引起井水位動態變化的固體潮頻率分布在2.3×10-5～7.9×10-6Hz，正常情況下，井水位對固體潮動態響應特征表現為每日呈雙峰雙谷變化。
【技術特征摘要】
1.一種基于eemd的井水位固體潮特征提取方法，其特征在于，包括以下步驟：
2.根據權利要求1所述的基于eemd的井水位固體潮特征提取方法，其特征在于，步驟s1中，對小時值井水位數據進行eemd分解時，引入高斯噪聲標準差設置為0.1，迭代次數為...
【專利技術屬性】
技術研發人員：許英霞，丁俊柯，馬傳璧，尹寶軍，王震坤，郭建芳，董祎瑋，王金生，張雨，
申請(專利權)人：山西工程技術學院，
類型：發明
國別省市：
還沒有人留言評論。發表了對其他瀏覽者有用的留言會獲得科技券。