[發(fā)明專利]一種疊后地震數(shù)據(jù)繞射波屬性提取方法有效
| 申請?zhí)枺?/td> | 201810896115.1 | 申請日： | 2018-08-08 |
| 公開（公告）號： | CN108594303B | 公開（公告）日： | 2019-12-13 |
| 發(fā)明（設(shè)計）人： | 舒夢珵;康洪全;程濤;賈懷存 | 申請（專利權(quán)）人： | 中國海洋石油集團有限公司;中海油研究總院有限責任公司 |
| 主分類號： | G01V1/28 | 分類號： | G01V1/28 |
| 代理公司： | 11245 北京紀凱知識產(chǎn)權(quán)代理有限公司 | 代理人： | 徐寧;孫楠 |
| 地址： | 100010 北*** | 國省代碼： | 北京;11 |
| 權(quán)利要求書： | 查看更多 | 說明書： | 查看更多 |
| 摘要： | |||
| 搜索關(guān)鍵詞： | 繞射 地震數(shù)據(jù) 波場 主成分分析 屬性提取 維度 主成分分析算法 三維地震數(shù)據(jù) 預處理 本征矢量 計算公式 實際數(shù)據(jù) 計算式 應用 改進 保證 | ||
本發(fā)明涉及一種疊后地震數(shù)據(jù)繞射波屬性提取方法，其步驟：地震數(shù)據(jù)預處理；根據(jù)實際數(shù)據(jù)建立正演模型，通過計算，確定主成分分析技術(shù)提取繞射波時，最適參數(shù)及主成分的維度；將本征矢量計算公式改進為適用于三維地震數(shù)據(jù)計算式，并應用最適參數(shù)及主成分維度，提取實際地震數(shù)據(jù)中的繞射波場。本發(fā)明通過應用主成分分析技術(shù)，實現(xiàn)了在疊后地震數(shù)據(jù)中對繞射波場的提取；通過對傳統(tǒng)主成分分析算法的改進，使其適用于在地震數(shù)據(jù)中對繞射波場進行空間提取的目的，也保證了繞射波提取結(jié)果的準確性。
技術(shù)領(lǐng)域
本發(fā)明涉及一種石油勘探領(lǐng)域的地震數(shù)據(jù)波場分離方法，特別是關(guān)于一種疊后地震數(shù)據(jù)繞射波屬性提取方法。
背景技術(shù)
隨著全球油氣田勘探程度的加深,儲藏中小尺度地質(zhì)體的精細描述逐步成為油田勘探的重點問題。碳酸鹽巖儲層縫洞空間的識別，潛山裂縫型油氣藏裂縫系統(tǒng)的識別及小型水道體系空間分布刻畫等問題受到廣泛的關(guān)注。常規(guī)的反射波儲藏描述方法受到波場分辨率的限制,很難對此類地質(zhì)體進行精細刻畫。
前人的研究已經(jīng)證實,繞射波場可以針對小尺度地質(zhì)體進行較好描述,當?shù)刭|(zhì)體的空間尺度接近或小于地震反射數(shù)據(jù)分辨率的情況下(一般小于四分之一波長),從地震數(shù)據(jù)中分離的繞射波數(shù)據(jù)能夠進行小尺度地質(zhì)體的精細刻畫。所以前人根據(jù)繞射波場的動力學和波動學特征,提出不同的繞射波提取方法,由于疊前繞射波提取方法計算效率較低,不利于儲層的快速識別,疊后繞射波提取方法就成為儲層識別的重要手段,現(xiàn)有的波場分離手段一般針對反射波場進行分離，而針對繞射波場有效分離的技術(shù)尚不成熟。
發(fā)明內(nèi)容
針對上述問題，本發(fā)明的目的是提供一種疊后地震數(shù)據(jù)繞射波屬性提取方法，其能夠在疊后地震數(shù)據(jù)中準確地進行繞射波場分離。
為實現(xiàn)上述目的，本發(fā)明采取以下技術(shù)方案：一種疊后地震數(shù)據(jù)繞射波屬性提取方法，其包括以下步驟：1)地震數(shù)據(jù)預處理；2)根據(jù)實際數(shù)據(jù)建立正演模型，通過計算，確定主成分分析技術(shù)提取繞射波時，最適參數(shù)及主成分的維度；3)將本征矢量計算公式改進為適用于三維地震數(shù)據(jù)計算式，并應用步驟2)中的最適參數(shù)及主成分維度，提取實際地震數(shù)據(jù)中的繞射波場。
進一步，所述步驟1)中，預處理包括以下步驟：1.1)定義如下變量：A＝{A1,A2,...,An}是一個m×n的矩陣，n為輸入的地震道數(shù)，每個地震道為m維向量，代表m個采樣點；An表示第n個地震道的地震數(shù)據(jù)采集序列；進行地震數(shù)據(jù)體的標準化預處理：
式中，xij為標準化后的地震數(shù)據(jù)采集序列；為樣本均值：Aij為原始采集的地震數(shù)據(jù)采集序列；Si為樣本標準差：i＝1，…，n；j＝1，…，m；1.2)處理后得到標準化地震數(shù)據(jù)體為X＝{x1，X2，...，Xn}；其中，X是一個m×n的矩陣，n為輸入的地震道數(shù)，每個地震道為m維向量，代表m個采樣點，xn表示第n個地震道標準化后的地震數(shù)據(jù)采集序列。
進一步，所述步驟2)中，最適參數(shù)及主成分的維度確定方法包括以下步驟：2.1)通過下式將輸入數(shù)據(jù)矢量xi變換成新的矢量si：si＝UTxi，其中，U是一個n×n正交矩陣，它的第i列Ui是樣本協(xié)方差矩陣；i＝1，…，n；2.2)計算第i個本征矢量：
λiui＝Cuii＝1，...，n (1)
其中，C為一個n階方陣；λi是C的一個本征值，ui是對應的本征矢量；
該專利技術(shù)資料僅供研究查看技術(shù)是否侵權(quán)等信息，商用須獲得專利權(quán)人授權(quán)。該專利全部權(quán)利屬于中國海洋石油集團有限公司;中海油研究總院有限責任公司，未經(jīng)中國海洋石油集團有限公司;中海油研究總院有限責任公司許可，擅自商用是侵權(quán)行為。如果您想購買此專利、獲得商業(yè)授權(quán)和技術(shù)合作，請聯(lián)系【客服】
本文鏈接：http://www.szxzyx.cn/pat/books/201810896115.1/2.html，轉(zhuǎn)載請聲明來源鉆瓜專利網(wǎng)。
