本發(fā)明涉及地質(zhì)勘探,尤其涉及一種應(yīng)力曲線的去噪方法、裝置、電子設(shè)備和存儲介質(zhì)。
背景技術(shù):
1、地震波在地層中傳播過程中發(fā)生速度頻散和能量衰減,這與地層和巖石的結(jié)構(gòu)、流體的性質(zhì)都有關(guān)系。石油與天然氣勘探、煤田勘探都主要依靠地震勘探來探測地下地質(zhì)構(gòu)造和流體分布的細節(jié),進而制定開發(fā)方案。地震波的速度是主要的參數(shù),其衰減信息也逐漸受到重視。地震巖石物理是將野外觀測數(shù)據(jù)轉(zhuǎn)換為巖石物性參數(shù)的橋梁,地震數(shù)字巖石物理是其中的一個重要研究方面,現(xiàn)有的地震數(shù)字巖石物理的方法和技術(shù)采用的彈性參數(shù)計算方法主要包括有限元靜力學(xué)模擬,透射波模擬,動態(tài)應(yīng)力應(yīng)變模擬,準(zhǔn)靜態(tài)應(yīng)力應(yīng)變模擬和阻尼波動方程模擬方法。
2、在現(xiàn)有技術(shù)中,基于數(shù)字巖心的動態(tài)應(yīng)力應(yīng)變模擬計算數(shù)字巖心的頻散和衰減曲線是當(dāng)前數(shù)字巖石物理領(lǐng)域的一個重要的研究內(nèi)容,但是,動態(tài)應(yīng)力應(yīng)變模擬的應(yīng)力曲線中存在高頻噪聲,應(yīng)力曲線的全頻帶的信息都是有效信息,故帶通濾波技術(shù)很難作為一種有效的壓制噪聲的手段。
3、因此,急需提出一種應(yīng)力曲線的去噪方法、裝置、電子設(shè)備和存儲介質(zhì),解決現(xiàn)有技術(shù)中存在的應(yīng)力曲線的全頻帶的信息都是有效信息,無法通過帶通濾波技術(shù)對應(yīng)力曲線進行去噪的技術(shù)問題。
技術(shù)實現(xiàn)思路
1、有鑒于此,有必要提供一種應(yīng)力曲線的去噪方法、裝置、電子設(shè)備和存儲介質(zhì),用以解決現(xiàn)有技術(shù)中存在的應(yīng)力曲線的全頻帶的信息都是有效信息,無法通過帶通濾波技術(shù)對應(yīng)力曲線進行去噪的技術(shù)問題。
2、為了解決上述問題,本發(fā)明提供一種應(yīng)力曲線的去噪方法,包括:
3、獲取應(yīng)力曲線上的采樣點時間序列和應(yīng)力值序列;
4、基于非線性方法對所述采樣點時間序列和所述應(yīng)力值序列進行稀疏重采樣,得到稀疏重采樣時間點序列和稀疏重采樣應(yīng)力值序列;
5、對所述稀疏重采樣時間點序列和所述稀疏重采樣應(yīng)力值序列進行處理,得到去噪采樣時間點序列和去噪應(yīng)力值序列,并根據(jù)所述去噪采樣時間點序列和所述去噪應(yīng)力值序列得到去噪應(yīng)力曲線。
6、在一種可能的實現(xiàn)方式中,所述基于非線性方法對所述采樣點時間序列和所述應(yīng)力值序列進行稀疏重采樣,得到稀疏重采樣時間點序列和稀疏重采樣應(yīng)力值序列,包括:
7、基于非線性方法對所述采樣點時間序列進行稀疏重采樣,得到稀疏重采樣時間點序列;
8、根據(jù)所述采樣點時間序列、所述應(yīng)力值序列和所述稀疏重采樣時間點序列,得到稀疏重采樣應(yīng)力值序列。
9、在一種可能的實現(xiàn)方式中,所述基于非線性方法對所述采樣點時間序列進行稀疏重采樣,得到稀疏重采樣時間點序列,包括:
10、基于非線性方法對所述采樣點時間序列進行稀疏重采樣,得到初始時間點序列;
11、對所述初始時間點序列中的數(shù)據(jù)進行篩選,得到稀疏重采樣時間點序列。
12、在一種可能的實現(xiàn)方式中,所述對所述初始時間點序列中的數(shù)據(jù)進行篩選,得到稀疏重采樣時間點序列,包括:
13、根據(jù)第一預(yù)設(shè)條件對所述初始時間點序列中的數(shù)據(jù)進行篩選,得到篩選時間點序列;所述第一預(yù)設(shè)條件用于對所述初始時間點序列中的每個時間點進行判斷;
14、根據(jù)第二預(yù)設(shè)條件確定所述篩選時間點序列中補充的重采樣時間點,得到稀疏重采樣時間點序列;所述第二預(yù)設(shè)條件用于對所述篩選時間點序列中的最大值進行判斷。
15、在一種可能的實現(xiàn)方式中,所述根據(jù)所述采樣點時間序列、所述應(yīng)力值序列和所述稀疏重采樣時間點序列,得到稀疏重采樣應(yīng)力值序列,包括:
16、根據(jù)所述稀疏重采樣時間點序列中的每個重采樣時間點的大小,確定所述采樣點時間序列中所述每個重采樣時間點對應(yīng)的第一采樣時間點和第二采樣時間點;
17、根據(jù)所述應(yīng)力值序列,確定所述第一采樣時間點對應(yīng)的第一曲線值和所述第二采樣時間點對應(yīng)的第二曲線值;
18、根據(jù)所述每個重采樣時間點對應(yīng)的所述第一采樣時間點、所述第二采樣時間點、所述第一曲線值和所述第二曲線值,得到所述每個重采樣時間點對應(yīng)的重采樣曲線值;
19、根據(jù)所述稀疏重采樣時間點序列中的所有重采樣曲線值,得到稀疏重采樣應(yīng)力值序列。
20、在一種可能的實現(xiàn)方式中,所述對所述稀疏重采樣時間點序列和所述稀疏重采樣應(yīng)力值序列進行處理,得到去噪采樣時間點序列和去噪應(yīng)力值序列,包括:
21、根據(jù)所述稀疏重采樣時間點序列和所述稀疏重采樣應(yīng)力值序列,得到稀疏重采樣曲線;
22、對所述稀疏重采樣曲線進行去噪處理,得到去噪應(yīng)力曲線,并得到所述去噪應(yīng)力曲線的處理采樣時間點序列和處理應(yīng)力值序列;
23、根據(jù)所述采樣點時間序列、所述處理采樣時間點序列對所述處理應(yīng)力值序列進行加密采樣,得到去噪采樣時間點序列和去噪應(yīng)力值序列。
24、在一種可能的實現(xiàn)方式中,所述根據(jù)所述采樣點時間序列、所述處理采樣時間點序列對所述處理應(yīng)力值序列進行加密采樣,得到去噪采樣時間點序列和去噪應(yīng)力值序列,包括:
25、根據(jù)所述采樣點時間序列中的每個采樣時間點的大小,確定所述處理采樣時間點序列中所述每個采樣時間點對應(yīng)的第一處理采樣時間點和第二處理采樣時間點;
26、根據(jù)所述處理應(yīng)力值序列,確定所述第一處理采樣時間點對應(yīng)的第一處理曲線值和所述第二處理采樣時間點對應(yīng)的第二處理曲線值;
27、根據(jù)所述每個采樣時間點對應(yīng)的所述第一處理采樣時間點、所述第二處理采樣時間點、所述第一處理曲線值和所述第二處理曲線值,得到所述每個采樣時間點對應(yīng)的加密采樣曲線值;
28、根據(jù)所述采樣點時間序列中的所有加密采樣曲線值,得到去噪應(yīng)力值序列,并將所述采樣點時間序列確定為去噪采樣時間點序列。
29、另一方面,本發(fā)明還提供了一種應(yīng)力曲線的去噪裝置,包括:
30、序列獲取模塊,用于獲取應(yīng)力曲線上的采樣點時間序列和應(yīng)力值序列;
31、稀疏重采樣模塊,用于基于非線性方法對所述采樣點時間序列和所述應(yīng)力值序列進行稀疏重采樣,得到稀疏重采樣時間點序列和稀疏重采樣應(yīng)力值序列;
32、曲線去噪模塊,用于對所述稀疏重采樣時間點序列和所述稀疏重采樣應(yīng)力值序列進行處理,得到去噪采樣時間點序列和去噪應(yīng)力值序列,并根據(jù)所述去噪采樣時間點序列和所述去噪應(yīng)力值序列得到去噪應(yīng)力曲線。
33、另一方面,本發(fā)明實施例公開了一種電子設(shè)備,包括:處理器、存儲器及存儲在所述存儲器上并能夠在所述處理器上運行的計算機程序,所述計算機程序被所述處理器執(zhí)行時實現(xiàn)上述應(yīng)力曲線的去噪方法實施例的各個步驟。
34、另一方面,本發(fā)明實施例公開了一種計算機可讀存儲介質(zhì),所述計算機可讀存儲介質(zhì)上存儲計算機程序,所述計算機程序被處理器執(zhí)行時實現(xiàn)上述應(yīng)力曲線的去噪方法實施例的各個步驟。
35、本發(fā)明的有益效果是:獲取應(yīng)力曲線上的采樣點時間序列和應(yīng)力值序列;然后可以基于非線性方法對采樣點時間序列和應(yīng)力值序列進行稀疏重采樣,得到稀疏重采樣時間點序列和稀疏重采樣應(yīng)力值序列;并對稀疏重采樣時間點序列和稀疏重采樣應(yīng)力值序列進行去噪處理,得到去噪采樣時間點序列和去噪應(yīng)力值序列,進而可以根據(jù)去噪采樣時間點序列和去噪應(yīng)力值序列得到去噪應(yīng)力曲線,通過非線性方法進行稀疏重采樣之后,再進行去噪處理,可以保留應(yīng)力曲線中的有效信息,提高了去噪后的應(yīng)力曲線的精度。