【技術實現(xiàn)步驟摘要】
本專利技術屬于地球物理，具體而言，涉及一種基于相似度量的地震資料處理保幅性評價方法。
技術介紹
1、在地震資料處理的過程中需要監(jiān)控每一個步驟的保幅性，通過保幅性來評價處理步驟的優(yōu)缺點。但是，由于缺乏足夠的技術手段，這個過程往往由處理人員憑借經(jīng)驗觀察來完成，沒有直觀的辦法評價保幅性。
2、專利“一種地震資料保幅程度定量評價方法及裝置”(cn117406279a)中，用到了基于小波分析的一種評價方法，該方法是將地震數(shù)據(jù)從時間域變換到頻率域后進行小波分析，再計算其相干性來評價的。專利“保幅處理分析評價方法”(cn102176055a)中，通過用原始segy地震數(shù)據(jù)與保幅處理后segy地震數(shù)據(jù)相減，得到差值地震數(shù)據(jù)，分析差值記錄中是否包含有效信號成分以及差值地震記錄中噪音的成分，判斷處理數(shù)據(jù)的保幅性。現(xiàn)有技術中存在地震資料處理質(zhì)量控制過程中保幅性難以定量準確評估的問題。
技術實現(xiàn)思路
1、為了解決地震資料處理質(zhì)量控制過程中保幅性難以定量準確評估的問題，本申請實施例提供了一種基于相似度量的地震資料處理保幅性評價方法，能夠定量的評價地震資料處理過程中的保幅性情況，為處理人員調(diào)整處理方法及參數(shù)提供直觀依據(jù)。
2、本申請實施例提供了一種基于相似度量的地震資料處理保幅性評價方法，包括：
3、s1，沿著一條主測線方向獲取地震資料處理前和處理后的疊加剖面；
4、s2，分別計算地震資料處理前、處理后相鄰兩條地震道的pearson相似系數(shù)，再求地震資料處理前后
5、s3，分別計算地震資料處理前、處理后相鄰兩條地震道的歐式距離，再求地震資料處理前后歐式距離的相對變化量并歸一化；
6、s4，分別計算地震資料處理前、處理后相鄰兩條地震道的dtw距離，再求地震資料處理前后dtw距離的相對變化量并歸一化；
7、s5，計算所有主測線的地震疊加剖面的平均相似度并加權求和，得到衡量地震資料處理前后地震疊加剖面保幅性的綜合評估值。
8、在一些實施例中，步驟s1中，沿著第j條主測線方向j＝1,2,...,j取出地震資料處理前和處理后的疊加剖面，地震資料處理前的疊加剖面記為s(x,t)，地震資料處理后的疊加剖面記為h(x,t)，其中x表示cmp號，t表示旅行時。
9、在一些實施例中，步驟s2包括：s21，在處理前的疊加剖面第j條主測線中，令si(t)表示第i道的共中心點地震道數(shù)據(jù)，si(t)＝s(xi,t)，計算si(t)與它相鄰的地震道si+1(t)之間的pearson相似系數(shù)；s22，在處理后的疊加剖面第j條主測線中，令hi(t)＝h(xi,t)，表示第i道處理后的地震資料，計算處理后相鄰兩條地震道之間的pearson相似系數(shù)；s23，計算第j條主測線處理前的相似系數(shù)與處理后的相似系數(shù)的相對變化量。
10、在一些實施例中，步驟s21中，si(t)＝s(xi,t)，采用滑動時間窗口算法分段計算si(t)與它相鄰的地震道si+1(t)之間的pearson相似系數(shù)p(si,si+1)(i＝1,2,...,m-1)，設時窗長度為δτ，滑動步長為δτ/2，共設置k個時窗劃分地震道，k的計算公式為：
11、第k個時窗的時間范圍為[τk,τk+δτ]，其中計算第k個時窗內(nèi)的pearson相似系數(shù)公式為：
12、
13、其中e(si,k)和e(si+1,k)分別表示相鄰兩條地震道si(t)和si+1(t)在第k個時窗內(nèi)的振幅平均值；
14、處理前一條主測線s(x,t)所有相鄰兩道的pearson相似系數(shù)p(si,si+1)表示為：
15、
16、步驟s22中，計算處理后相鄰兩條地震道之間的pearson相似系數(shù)p(hi,hi+1)也采用滑動時間窗口算法分段計算，i＝1,2,...,m-1，時窗的劃分與步驟s21一致，第k個時窗內(nèi)處理后相鄰兩條地震道之間的pearson相似系數(shù)為p(hi,k,hi+1,k)，計算公式為：
17、
18、其中e(hi,k)和e(hi+1,k)分別表示處理后相鄰兩條地震道hi(t)和hi+1(t)在第k個時窗內(nèi)的振幅平均值；
19、步驟s23中，得到第j條主測線處理前后第i個相鄰地震道之間相似系數(shù)的變化程度pi,j，公式為：
20、
21、其中i＝1,2,...,m-1。
22、在一些實施例中，步驟s3包括：s31，在處理前的疊加剖面第j條主測線中，計算相鄰兩條地震道si(t)和si+1(t)之間的歐式距離；s32，將相鄰兩條地震道si(t)和si+1(t)之間的歐式距離的度量值變換到[0,1]內(nèi)；s33，計算處理后相鄰兩條地震道hi(t)和hi+1(t)之間的歐式距離；s34，將處理后相鄰兩條地震道hi(t)和hi+1(t)之間的歐式距離的度量值變換到[0,1]內(nèi)；s35，計算第j條主測線處理前的歐氏距離與處理后的歐氏距離的相對變化量。
23、在一些實施例中，步驟s31中，計算相鄰兩條地震道si(t)和si+1(t)之間的歐式距離d(si,si+1)公式為：
24、
25、其中i＝1,2,...,m-1；
26、步驟s32中，相鄰兩條地震道si(t)和si+1(t)之間的歐式距離d(si,si+1)的度量值變換到[0,1]內(nèi)得到d(si,si+1)，變換公式為：
27、
28、其中a為正實數(shù)，用于控制函數(shù)圖像橫坐標伸縮程度；
29、步驟s33中，計算地震資料處理后相鄰兩條地震道hi(t)和hi+1(t)之間的歐式距離d(hi,hi+1)公式為：
30、
31、其中i＝1,2,...,m-1；
32、步驟s34中處理后相鄰兩條地震道hi(t)和hi+1(t)之間的歐式距離d(hi,hi+1)的度量值變換到[0,1]內(nèi)得到d(hi,hi+1)，變換公式為：
33、
34、其中a為正實數(shù)，與步驟s32中的a一致；
35、步驟s35中，得到第j條主測線處理前后第i個相鄰地震道之間歐氏距離的變化程度di,j，公式為：
36、
37、其中i＝1,2,...,m-1。
38、在一些實施例中，步驟s4包括：s41，計算處理前相鄰兩條地震道si(t)和si+1(t)之間的dtw距離，衡量兩條地震道的相似性；s42，將相鄰兩條地震道si(t)和si+1(t)之間的dtw距離的度量值變換到[0,1]內(nèi)；s43，計算處理后相鄰兩條地震道hi(t)和hi+1(t)之間的dtw距離，衡量兩條地震道的相似性；s44，將相鄰兩條地震道hi(t)和hi+1(t)之間的dtw距離的度量值變換到[0,1]內(nèi)；s45，計算第j條主測線處理前的dtw距離與處理后的dtw距離的相對變化量。
本文檔來自技高網(wǎng)...
【技術保護點】
1.一種基于相似度量的地震資料處理保幅性評價方法，其特征在于，包括：
2.根據(jù)權利要求1所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟S1中，沿著第j條主測線方向j＝1,2,...,J取出地震資料處理前和處理后的疊加剖面，地震資料處理前的疊加剖面記為S(x,t)，地震資料處理后的疊加剖面記為H(x,t)，其中x表示CMP號，t表示旅行時。
3.根據(jù)權利要求1或2所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟S2包括：S21，在處理前的疊加剖面第j條主測線中，令si(t)表示第i道的共中心點地震道數(shù)據(jù)，si(t)＝S(xi,t)，計算si(t)與它相鄰的地震道si+1(t)之間的Pearson相似系數(shù)；S22，在處理后的疊加剖面第j條主測線中，令hi(t)＝H(xi,t)，表示第i道處理后的地震資料，計算處理后相鄰兩條地震道之間的Pearson相似系數(shù)；S23，計算第j條主測線處理前的相似系數(shù)與處理后的相似系數(shù)的相對變化量。
4.根據(jù)權利要求3所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟S21中，
5.根據(jù)權利要求1或2所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟S3包括：S31，在處理前的疊加剖面第j條主測線中，計算相鄰兩條地震道si(t)和si+1(t)之間的歐式距離；S32，將相鄰兩條地震道si(t)和si+1(t)之間的歐式距離的度量值變換到[0,1]內(nèi)；S33，計算處理后相鄰兩條地震道hi(t)和hi+1(t)之間的歐式距離；S34，將處理后相鄰兩條地震道hi(t)和hi+1(t)之間的歐式距離的度量值變換到[0,1]內(nèi)；S35，計算第j條主測線處理前的歐氏距離與處理后的歐氏距離的相對變化量。
6.根據(jù)權利要求5所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟S31中，計算相鄰兩條地震道si(t)和si+1(t)之間的歐式距離d(si,si+1)公式為：
7.根據(jù)權利要求1或2所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟S4包括：S41，計算處理前相鄰兩條地震道si(t)和si+1(t)之間的DTW距離，衡量兩條地震道的相似性；S42，將相鄰兩條地震道si(t)和si+1(t)之間的DTW距離的度量值變換到[0,1]內(nèi)；S43，計算處理后相鄰兩條地震道hi(t)和hi+1(t)之間的DTW距離，衡量兩條地震道的相似性；S44，將相鄰兩條地震道hi(t)和hi+1(t)之間的DTW距離的度量值變換到[0,1]內(nèi)；S45，計算第j條主測線處理前的DTW距離與處理后的DTW距離的相對變化量。
8.根據(jù)權利要求7所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟S41中，得到處理前相鄰兩條地震道的DTW距離dtw(si,si+1)，dtw(si,si+1)＝dtw(si,N,si+1,N)；
9.根據(jù)權利要求1或2所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，在步驟S5中得到衡量地震資料處理前后地震疊加剖面保幅性的綜合評估值Comb，計算公式為：
...【技術特征摘要】
1.一種基于相似度量的地震資料處理保幅性評價方法，其特征在于，包括：
2.根據(jù)權利要求1所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟s1中，沿著第j條主測線方向j＝1,2,...,j取出地震資料處理前和處理后的疊加剖面，地震資料處理前的疊加剖面記為s(x,t)，地震資料處理后的疊加剖面記為h(x,t)，其中x表示cmp號，t表示旅行時。
3.根據(jù)權利要求1或2所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟s2包括：s21，在處理前的疊加剖面第j條主測線中，令si(t)表示第i道的共中心點地震道數(shù)據(jù)，si(t)＝s(xi,t)，計算si(t)與它相鄰的地震道si+1(t)之間的pearson相似系數(shù)；s22，在處理后的疊加剖面第j條主測線中，令hi(t)＝h(xi,t)，表示第i道處理后的地震資料，計算處理后相鄰兩條地震道之間的pearson相似系數(shù)；s23，計算第j條主測線處理前的相似系數(shù)與處理后的相似系數(shù)的相對變化量。
4.根據(jù)權利要求3所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟s21中，si(t)＝s(xi,t)，采用滑動時間窗口算法分段計算si(t)與它相鄰的地震道si+1(t)之間的pearson相似系數(shù)p(si,si+1)(i＝1,2,...,m-1)，設時窗長度為δτ，滑動步長為δτ/2，共設置k個時窗劃分地震道，k的計算公式為：
5.根據(jù)權利要求1或2所述基于相似度量的地震資料處理保幅性評價方法，其特征在于，步驟s3包括：s31，在處理前的疊加剖面第j條主測線中，計算相鄰兩條地震道si(t)和si+1(t)之間的歐式距離；s32，將相鄰兩條地震道...
【專利技術屬性】
技術研發(fā)人員：付飛琪，鮑偉，胡國斌，付嘉瑋，
申請(專利權)人：中國石油化工股份有限公司，
類型：發(fā)明
國別省市：
還沒有人留言評論。發(fā)表了對其他瀏覽者有用的留言會獲得科技券。