本發(fā)明屬于地震疊后數(shù)據(jù)補(bǔ)償，具體涉及一種火山巖復(fù)雜區(qū)疊前空間剩余振幅補(bǔ)償方法。
背景技術(shù)：
1、地震波在地下傳播時(shí)，振幅會隨著傳播過程有所衰減?；鹕綆r反射界面強(qiáng)，界面間波阻抗較大，會嚴(yán)重屏蔽上下伏地層的反射能量。而且，由于地質(zhì)結(jié)構(gòu)比較復(fù)雜，地震波會因地質(zhì)結(jié)構(gòu)的不同導(dǎo)致地震能量的變化。尤其是在介質(zhì)分界面處，當(dāng)激發(fā)的地震波在地下傳播到介質(zhì)分界面時(shí)，由于波抗阻的變化引起質(zhì)點(diǎn)震動位移振幅的變化，振幅發(fā)生了畸變，給油氣勘探、地震解釋和儲層反演都造成了很大的困難。
2、目前，對于振幅補(bǔ)償方式主要有基于幾何擴(kuò)散理論的球面擴(kuò)散補(bǔ)償及基于地表非一致性吸收衰減的地表一致性振幅補(bǔ)償，這兩種振幅補(bǔ)償方法無法解決地下存在火山巖復(fù)雜區(qū)所引起的振幅衰減，從而造成成像剖面上火山巖發(fā)育區(qū)和非火山巖發(fā)育區(qū)空間振幅不一致。
3、為了對介質(zhì)分界面處振幅進(jìn)行補(bǔ)償處理，中國發(fā)明專利cn108181649b公開了一種介質(zhì)分界面振幅補(bǔ)償方法及裝置，包括：根據(jù)待分析地震數(shù)據(jù)，獲取地震道的均方根振幅值；拾取介質(zhì)分界面，根據(jù)所述地震道的均方根振幅值，計(jì)算所述介質(zhì)分界面對應(yīng)的分界面反射振幅值；根據(jù)所述分界面反射振幅值和預(yù)設(shè)的地震能量總值，利用a＝ar+at計(jì)算所述介質(zhì)分界面對應(yīng)的分界面透射振幅值，其中，a表示所述地震能量總值，ar表示所述分界面反射振幅值，at表示所述分界面透射振幅值；根據(jù)所述分界面透射振幅值，利用st＝as/at計(jì)算分界面振幅補(bǔ)償因子，其中，st表示所述分界面振幅補(bǔ)償因子，as表示指定的歸一化振幅級別值，at表示所述分界面透射振幅值；根據(jù)所述分界面振幅補(bǔ)償因子和所述地震道的均方根振幅值，進(jìn)行所述介質(zhì)分界面的振幅補(bǔ)償操作。該補(bǔ)償處理減少了介質(zhì)分界面處振幅強(qiáng)弱特征變化不均勻的現(xiàn)象，使得補(bǔ)償后的地震數(shù)據(jù)振幅強(qiáng)弱變化更加自然。
4、中國發(fā)明專利cn115903044a年公開了一種針對強(qiáng)屏蔽下的能量吸收衰減補(bǔ)償方法，包括步驟：s1，輸入去噪后的地震數(shù)據(jù)；s2，求取地震數(shù)據(jù)的均方根振幅；s3，根據(jù)均方根振幅值，求取地震道振幅隨時(shí)間變化而衰減的指數(shù)模型；s4，根據(jù)指數(shù)模型，建立分時(shí)窗補(bǔ)償函數(shù)，對整體數(shù)據(jù)進(jìn)行補(bǔ)償。該發(fā)明是分時(shí)窗求取衰減系數(shù)，所得到的補(bǔ)償函數(shù)在時(shí)間上和空間上都是可變的，可以反映火成巖等強(qiáng)反射界面的特點(diǎn)，能對火成巖上下伏地層衰減的能量進(jìn)行很好地恢復(fù)。
5、但是，上述的補(bǔ)償方法，并未針對火山巖的復(fù)雜空間剩余振幅進(jìn)行補(bǔ)償方法的設(shè)置，而且，并未消除地震疊加數(shù)據(jù)中可能存在的異常均方根振幅屬性，復(fù)雜空間目的層空間振幅的一致性是否得到提高并不清楚。因此，亟需一種可以應(yīng)用于火山巖復(fù)雜區(qū)空間剩余振幅補(bǔ)償方法，實(shí)現(xiàn)火山巖復(fù)雜區(qū)疊前空間剩余振幅補(bǔ)償，提高目的層振幅一致性。
技術(shù)實(shí)現(xiàn)思路
1、本發(fā)明針對現(xiàn)有技術(shù)存在的問題，提供了一種火山巖復(fù)雜區(qū)疊前空間剩余振幅補(bǔ)償方法，通過對火成巖發(fā)育區(qū)目的層均方根振幅統(tǒng)計(jì)反映出目的層振幅變化規(guī)律；經(jīng)過平滑、歸一化后消除由于地震疊加數(shù)據(jù)中可能存在壞道、異常道引起的異常均方根振幅屬性；通過提取振幅補(bǔ)償因子，得到能夠應(yīng)用于疊前cmp道集的計(jì)算因子；實(shí)現(xiàn)了復(fù)雜區(qū)疊前空間剩余振幅補(bǔ)償，尤其是火山巖復(fù)雜區(qū)疊前空間剩余振幅補(bǔ)償，提高了目的層空間振幅一致性。
2、為實(shí)現(xiàn)上述目的，本發(fā)明采用的技術(shù)方案如下：
3、一方面，本發(fā)明提供一種復(fù)雜區(qū)疊前空間剩余振幅補(bǔ)償方法，包括步驟：
4、s1、基于地震數(shù)據(jù)體，統(tǒng)計(jì)目的層均方根振幅屬性；
5、s2、對s1所述的目的層的均方根振幅屬性進(jìn)行平滑處理、歸一化處理；
6、s3、對s2的歸一化處理后的均方根振幅屬性進(jìn)行提取剩余振幅補(bǔ)償因子；
7、s4、利用剩余振幅補(bǔ)償因子對疊前道集進(jìn)行振幅補(bǔ)償，得到補(bǔ)償后的疊前道集；
8、s5、對所得到的補(bǔ)償后的疊前道集進(jìn)行疊加處理。
9、優(yōu)選地，在進(jìn)行步驟s1前，需要對地震波進(jìn)行采集，得到地震數(shù)據(jù)體；地震數(shù)據(jù)體可以疊加形成疊后地震數(shù)據(jù)體或地震疊加數(shù)據(jù)體。
10、優(yōu)選地，步驟s1中，所述發(fā)育區(qū)目的層均方根振幅屬性，為：地震層位解釋約束下某一時(shí)窗范圍內(nèi)疊加剖面振幅均方根值，得到目的層的均方根振幅屬性。
11、進(jìn)一步優(yōu)選地，所述時(shí)窗范圍，其時(shí)窗范圍的確定具體為：基于地震疊加剖面依據(jù)同相軸連續(xù)性，拾取地層界面，通過起點(diǎn)地層界面和終點(diǎn)地層界面，確定出統(tǒng)計(jì)時(shí)的時(shí)窗范圍。
12、優(yōu)選地，步驟s2中，所述平滑處理，平滑網(wǎng)格大小取最大偏移距。
13、進(jìn)一步優(yōu)選地，所述平滑網(wǎng)格大小為7000×7000m。
14、優(yōu)選地，步驟s2中，所述平滑處理，具體為：將目的層的均方根振幅屬性按預(yù)設(shè)平滑網(wǎng)格內(nèi)振幅值相加，而后除以預(yù)設(shè)網(wǎng)格內(nèi)的地震道數(shù)，得到平滑后的振幅值。
15、優(yōu)選地，步驟s2中，所述歸一化處理，具體為：對平滑后的均方根振幅屬性按照下述公式進(jìn)行歸一化處理：其中，ai為歸一化后均方根值；ai為步驟s2所得平滑后的振幅值；amax為步驟s2所得平滑后振幅值中最大振幅值。
16、優(yōu)選地，步驟s3中，所述剩余振幅補(bǔ)償因子為歸一化后的均方根振幅屬性大小。
17、優(yōu)選地，步驟s3中，所述剩余振幅補(bǔ)償因子的計(jì)算公式為pi＝1-ai，其中，pi為提取的剩余振幅補(bǔ)償因子。
18、優(yōu)選地，步驟s4中，所述疊前道集(cmp，common?middle?point)包括共炮點(diǎn)道集、共接收點(diǎn)道集、共反射點(diǎn)道集，對疊前道集中的每個(gè)cmp點(diǎn)進(jìn)行剩余振幅補(bǔ)償因子的振幅補(bǔ)償。
19、優(yōu)選地，步驟s4中，所述振幅補(bǔ)償，補(bǔ)充步驟具體為：疊前道集中的每個(gè)cmp點(diǎn)對應(yīng)一個(gè)剩余振幅補(bǔ)償因子，每個(gè)疊前道集均與對應(yīng)的剩余振幅補(bǔ)償因子相乘，得到補(bǔ)償后的疊前道集。
20、進(jìn)一步優(yōu)選地，所述振幅補(bǔ)償，補(bǔ)償步驟具體為：
21、s4-1：從疊前道集中確定一個(gè)cmp點(diǎn)；
22、s4-2：將步驟s4-1確定的cmp點(diǎn)與對應(yīng)步驟s3所得剩余振幅補(bǔ)償因子相乘，得到補(bǔ)償后的cmp點(diǎn)，具體計(jì)算公式為：mi＝pi*ci，其中，mi為補(bǔ)償后疊前cmp道集，ci為疊前道集；
23、s4-3：遍歷疊前道集，重復(fù)步驟s4-1和s4-2，得到所有補(bǔ)償后的cmp道集。
24、另一方面，本發(fā)明提供上述方法在檢測復(fù)雜區(qū)疊前空間振幅中的應(yīng)用。
25、然后，本發(fā)明還提供了上述方法所建立的剩余振幅補(bǔ)償因子。
26、最后，本發(fā)明還提供了上述剩余振幅補(bǔ)償因子在疊前道集和/或檢測火山巖空間振幅中的應(yīng)用。
27、相對于現(xiàn)有技術(shù)，本發(fā)明具有以下有益效果：
28、本發(fā)明的方法中，通過統(tǒng)計(jì)復(fù)雜區(qū)目的層的均方根振幅屬性，然后對該平面屬性平滑及歸一化之后用于疊前道集的空間反補(bǔ)償，計(jì)算剩余振幅補(bǔ)償因子，應(yīng)用于疊前道集，提高目的層空間振幅一致性。
1.一種復(fù)雜區(qū)疊前空間剩余振幅補(bǔ)償方法，其特征在于，包括步驟：
2.根據(jù)權(quán)利要求1所述的方法，其特征在于，在進(jìn)行步驟s1前，需要對地震波進(jìn)行采集，得到地震數(shù)據(jù)體；地震數(shù)據(jù)體可以疊加形成疊后地震數(shù)據(jù)體或地震疊加數(shù)據(jù)體。
3.根據(jù)權(quán)利要求1所述的方法，其特征在于，步驟s1中，所述發(fā)育區(qū)目的層均方根振幅屬性，為：地震層位解釋約束下某一時(shí)窗范圍內(nèi)疊加剖面振幅均方根值，得到目的層的均方根振幅屬性。
4.根據(jù)權(quán)利要求3所述的方法，其特征在于，所述時(shí)窗范圍，其時(shí)窗范圍的確定具體為：基于地震疊加剖面依據(jù)同相軸連續(xù)性，拾取地層界面，通過起點(diǎn)地層界面和終點(diǎn)地層界面，確定出統(tǒng)計(jì)時(shí)的時(shí)窗范圍。
5.根據(jù)權(quán)利要求1所述的方法，其特征在于，步驟s2中，所述平滑處理，平滑網(wǎng)格大小取最大偏移距；平滑處理，具體為：將目的層的均方根振幅屬性按預(yù)設(shè)平滑網(wǎng)格內(nèi)振幅值相加，而后除以預(yù)設(shè)網(wǎng)格內(nèi)的地震道數(shù)，得到平滑后的振幅值。
6.根據(jù)權(quán)利要求1所述的方法，其特征在于，步驟s2中，所述歸一化處理，具體為：對平滑后的均方根振幅屬性按照下述公式進(jìn)行歸一化處理：其中，ai為歸一化后均方根值；ai為步驟s2所得平滑后的振幅值；amax為步驟s2所得平滑后振幅值中最大振幅值。
7.根據(jù)權(quán)利要求1所述的方法，其特征在于，步驟s3中，所述剩余振幅補(bǔ)償因子為歸一化后的均方根振幅屬性大??；所述剩余振幅補(bǔ)償因子的計(jì)算公式為pi＝1-ai，其中，pi為提取的剩余振幅補(bǔ)償因子。
8.根據(jù)權(quán)利要求1所述的方法，其特征在于，步驟s4中，所述疊前道集包括共炮點(diǎn)道集、共接收點(diǎn)道集、共反射點(diǎn)道集；對疊前道集中的每個(gè)cmp點(diǎn)進(jìn)行剩余振幅補(bǔ)償因子的振幅補(bǔ)償。
9.根據(jù)權(quán)利要求1所述的方法，其特征在于，步驟s4中，所述振幅補(bǔ)償，補(bǔ)充步驟具體為：疊前道集中的每個(gè)cmp點(diǎn)對應(yīng)一個(gè)剩余振幅補(bǔ)償因子，每個(gè)疊前道集均與對應(yīng)的剩余振幅補(bǔ)償因子相乘，得到補(bǔ)償后的疊前道集。
10.根據(jù)權(quán)利要求9所述的方法，其特征在于，所述振幅補(bǔ)償，補(bǔ)償步驟具體為：
11.權(quán)利要求1-10任一項(xiàng)所述的方法在檢測復(fù)雜區(qū)疊前空間振幅中的應(yīng)用。
12.權(quán)利要求1-10任一項(xiàng)所述的方法所建立的剩余振幅補(bǔ)償因子。
13.權(quán)利要求12所述的剩余振幅補(bǔ)償因子在疊前道集和/或檢測火山巖空間振幅中的應(yīng)用。