本發(fā)明涉及地球物理勘探和開發(fā)，尤其涉及一種基于速度低頻趨勢(shì)建模的地震重構(gòu)方法。

背景技術(shù)：

1、由于深洼帶儲(chǔ)層速度和密度受到壓實(shí)作用的影響，隨著埋深增大、壓實(shí)作用增強(qiáng)，直接影響到三類儲(chǔ)層的速度和密度變化特征，也即影響到波阻抗變化特征。受壓實(shí)影響，平面上不同沉積相砂泥巖波阻抗存在差異性，測(cè)井曲線垂直向上砂泥巖波阻抗差隨著埋深的增大逐漸減小，使得同一類儲(chǔ)層不僅導(dǎo)致測(cè)井巖性波阻抗交叉重疊嚴(yán)重，而且波阻抗差異性減小產(chǎn)生弱反射降低了地震識(shí)別能力，壓實(shí)作用引起的同一類儲(chǔ)層波阻抗差異變化，將直接影響到地震記錄特征的變化，使得從現(xiàn)有地震記錄中提取三類有效儲(chǔ)層對(duì)應(yīng)的地震屬性變化的規(guī)律性不強(qiáng)，難以優(yōu)選出合適的地震屬性來(lái)描述三類儲(chǔ)層。

2、常用的地層去壓實(shí)校正方法有壓實(shí)率恢復(fù)法、壓實(shí)系數(shù)恢復(fù)法、去壓實(shí)恢復(fù)法、地層原始厚度恢復(fù)圖版法、層序地層分析法、鏡質(zhì)組最大反射率恢復(fù)法、地球物理法等。壓實(shí)系數(shù)恢復(fù)法是最常用的方法，其中壓實(shí)系數(shù)是反演盆地埋藏過(guò)程的重要參數(shù)，考慮壓實(shí)作用引起的同一類儲(chǔ)層波阻抗差異變化，將直接影響到地震記錄特征的變化，也是準(zhǔn)確獲取正常壓實(shí)段泥巖壓實(shí)系數(shù)的前提，波阻抗的低頻(時(shí)間域)或低波數(shù)(深度域)分量隨深度變化代表了地層波阻抗的壓實(shí)趨勢(shì)。為此，通過(guò)提取波阻抗的低頻或低波數(shù)分量，可以求得時(shí)間域或深度域地層波阻抗壓實(shí)系數(shù)。

技術(shù)實(shí)現(xiàn)思路

1、鑒于上述問(wèn)題，提出了本發(fā)明以便提供克服上述問(wèn)題或者至少部分地解決上述問(wèn)題的一種基于速度低頻趨勢(shì)建模的地震重構(gòu)方法。

2、根據(jù)本發(fā)明的一個(gè)方面，提供了一種基于速度低頻趨勢(shì)建模的地震重構(gòu)方法包括：

3、步驟s1：將測(cè)井資料轉(zhuǎn)化為時(shí)間域，對(duì)測(cè)井波阻抗進(jìn)行分頻處理，選取與趨勢(shì)相符的低頻分量為測(cè)井低頻分量數(shù)據(jù)體；

4、步驟s2：基于地震速度數(shù)據(jù)體及測(cè)井密度與速度的擬合關(guān)系計(jì)算地震波阻抗低頻分量，利用所述測(cè)井低頻分量數(shù)據(jù)體對(duì)所述地震波阻抗低頻分量進(jìn)行校正，得到校正后的三維空間地震波阻抗低頻分量數(shù)據(jù)體；

5、步驟s3：根據(jù)校正后的三維空間地震波阻抗低頻分量數(shù)據(jù)體確定目的層段對(duì)應(yīng)的壓實(shí)趨勢(shì)線，計(jì)算目的層去壓實(shí)后的波阻抗值，計(jì)算波阻抗去壓實(shí)系數(shù)；

6、步驟s4：建立去壓實(shí)前后的波阻抗關(guān)系、反射系數(shù)數(shù)學(xué)關(guān)系，建立多道地震記錄與地震子波和多道波阻抗差分的數(shù)學(xué)關(guān)系；

7、步驟s5：在波阻抗低頻分量趨勢(shì)約束情形下，建立去壓實(shí)后的波阻抗反演目標(biāo)函數(shù)，采用非線性反演方法得到壓實(shí)校正后的地震記錄。

8、可選的，所述步驟s1：將測(cè)井資料轉(zhuǎn)化為時(shí)間域，對(duì)測(cè)井波阻抗進(jìn)行分頻處理，選取與趨勢(shì)相符的低頻分量為測(cè)井低頻分量數(shù)據(jù)體具體包括：

9、測(cè)井?dāng)?shù)據(jù)計(jì)算時(shí)間域測(cè)井波阻抗數(shù)據(jù)體，采用廣義s變換時(shí)頻譜分析，確定目的層井段代表儲(chǔ)層壓實(shí)趨勢(shì)的波阻抗低頻分量范圍；

10、采用廣義s反變換對(duì)測(cè)井波阻抗低頻分量進(jìn)行重構(gòu)，獲得測(cè)井時(shí)間域波阻抗低頻分量。

11、可選的，所述步驟s2：基于地震速度數(shù)據(jù)體及測(cè)井密度與速度的擬合關(guān)系計(jì)算地震波阻抗低頻分量，利用測(cè)井低頻分量對(duì)地震波阻抗分量進(jìn)行校正，得到校正后的三維空間地震波阻抗低頻分量數(shù)據(jù)體具體包括：

12、基于地震偏移速度數(shù)據(jù)體、解釋的主要目的層數(shù)據(jù)，采用層控多約束層速度反演mcvi方法計(jì)算地震層速度；

13、進(jìn)行目的層測(cè)井速度和測(cè)井密度交匯圖分析，獲得密度與速度的擬合公式，并應(yīng)用于地震層速度換算出地層密度，計(jì)算地震波阻抗體；

14、利用廣義s變換獲取低頻波阻抗體；

15、利用測(cè)井波阻抗低頻分量進(jìn)行校正，得到校正后的三維空間地震波阻抗低頻分量數(shù)據(jù)體。

16、可選的，所述步驟s3：根據(jù)校正后的三維空間地震波阻抗低頻分量數(shù)據(jù)體確定目的層段對(duì)應(yīng)的壓實(shí)趨勢(shì)線，計(jì)算目的層去壓實(shí)后的波阻抗值，計(jì)算波阻抗去壓實(shí)系數(shù)具體包括：

17、根據(jù)低頻波阻抗值是因?yàn)閴簩?shí)作用隨深度變化特征，確定目的層段對(duì)應(yīng)的壓實(shí)趨勢(shì)線，計(jì)算目的層去壓實(shí)后的波阻抗值；

18、利用去壓實(shí)前后的波阻抗值計(jì)算波阻抗壓實(shí)系數(shù)，得到目的層段測(cè)井和三維波阻抗壓實(shí)系數(shù)隨深度變化曲線。

19、可選的，所述步驟s3中計(jì)算井旁地震波阻抗壓實(shí)系數(shù)與測(cè)井波阻抗去壓實(shí)系數(shù)的差異，結(jié)合解釋層位，采用克里金插值校正技術(shù)，建立層控地震波阻抗壓實(shí)系數(shù)誤差數(shù)據(jù)體；

20、對(duì)地震波阻抗壓實(shí)系數(shù)進(jìn)行校正，得到校正后的三維空間地震波阻抗壓實(shí)系數(shù)數(shù)據(jù)體。

21、可選的，所述步驟s4：建立去壓實(shí)前后的波阻抗關(guān)系、反射系數(shù)數(shù)學(xué)關(guān)系，建立多道地震記錄與地震子波和多道波阻抗差分的數(shù)學(xué)關(guān)系具體包括：

22、根據(jù)反射系數(shù)與波阻抗之間的非線性關(guān)系和差分算子，建立反射系數(shù)、多道波阻抗差分和地震記錄的關(guān)系。

23、可選的，所述步驟s5：在波阻抗低頻分量趨勢(shì)約束情形下，建立去壓實(shí)后的波阻抗反演目標(biāo)函數(shù)，采用非線性反演方法得到壓實(shí)校正后的地震記錄具體包括：

24、建立低頻約束目標(biāo)函數(shù)，采用非線性反演方法得到壓實(shí)校正后的波阻抗，得到去壓實(shí)校正后的反射系數(shù)；

25、利用井旁資料提取的地震子波，通過(guò)褶積得到去壓實(shí)校正后的地震記錄。

26、可選的，所述步驟s5中的非線性反演方法中采用了譜投影梯度算法算法求解到壓實(shí)校正后的波阻抗。

27、可選的，所述步驟s5中建立去壓實(shí)后的波阻抗反演目標(biāo)函數(shù)具體包括：多道波阻抗反演問(wèn)題表示為

28、

29、其中，w為子波時(shí)間序列組成的行向量，xlow為波阻抗低頻分量，從測(cè)井波阻抗或地震速度譜獲得；l為低通濾波算子；y為已知的多道地震記錄，為所反演的反射系數(shù)時(shí)間序列具有稀疏性和橫向光滑性；為合成地震記錄與實(shí)際地震記錄吻合約束；xlow＝lx為波阻抗低頻分量趨勢(shì)約束；

30、可選的，所述步驟s5中建立去壓實(shí)后的波阻抗反演目標(biāo)函數(shù)還包括：去壓實(shí)后的波阻抗反演問(wèn)題表示為

31、

32、式中，c為壓實(shí)系數(shù)組成的矩陣；

33、反演求得去壓實(shí)后的波阻抗半對(duì)數(shù)xc；

34、利用差分算子計(jì)算去壓實(shí)后的多道反射系數(shù)rc；

35、利用井旁地震紀(jì)錄提取的子波計(jì)算去壓實(shí)后的地震記錄yc。

36、本發(fā)明提供的一種基于速度低頻趨勢(shì)建模的地震重構(gòu)方法包括：步驟s1：將測(cè)井資料轉(zhuǎn)化為時(shí)間域，對(duì)測(cè)井波阻抗進(jìn)行分頻處理，選取與趨勢(shì)相符的低頻分量為測(cè)井低頻分量數(shù)據(jù)體；步驟s2：基于地震速度數(shù)據(jù)體及測(cè)井密度與速度的擬合關(guān)系計(jì)算地震波阻抗低頻分量，利用所述測(cè)井低頻分量數(shù)據(jù)體對(duì)所述地震波阻抗低頻分量進(jìn)行校正，得到校正后的三維空間地震波阻抗低頻分量數(shù)據(jù)體；步驟s3：根據(jù)校正后的三維空間地震波阻抗低頻分量數(shù)據(jù)體確定目的層段對(duì)應(yīng)的壓實(shí)趨勢(shì)線，計(jì)算目的層去壓實(shí)后的波阻抗值，計(jì)算波阻抗去壓實(shí)系數(shù)；步驟s4：建立去壓實(shí)前后的波阻抗關(guān)系、反射系數(shù)數(shù)學(xué)關(guān)系，建立多道地震記錄與地震子波和多道波阻抗差分的數(shù)學(xué)關(guān)系；步驟s5：在波阻抗低頻分量趨勢(shì)約束情形下，建立去壓實(shí)后的波阻抗反演目標(biāo)函數(shù)，采用非線性反演方法得到壓實(shí)校正后的地震記錄。用于校正深洼帶儲(chǔ)層壓實(shí)作用，從而提高深洼帶儲(chǔ)層的薄儲(chǔ)層、儲(chǔ)層物性、儲(chǔ)層含油性的預(yù)測(cè)準(zhǔn)確率。廣泛適用于西部探區(qū)準(zhǔn)中深層乃至中石化探區(qū)的深層油氣藏勘探。

37、上述說(shuō)明僅是本發(fā)明技術(shù)方案的概述，為了能夠更清楚了解本發(fā)明的技術(shù)手段，而可依照說(shuō)明書的內(nèi)容予以實(shí)施，并且為了讓本發(fā)明的上述和其它目的、特征和優(yōu)點(diǎn)能夠更明顯易懂，以下特舉本發(fā)明的具體實(shí)施方式。