一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測方法
【專利摘要】本發(fā)明公開了一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測方法，所述方法包括：(1)測量數(shù)據(jù)，地質(zhì)數(shù)據(jù)包括錄井、試氣以及泥頁巖巖心含氣量測試數(shù)據(jù)；測井數(shù)據(jù)包括聲波時差、密度以及經(jīng)巖心含氣量測試數(shù)據(jù)刻度的含氣量測井解釋曲線；地震數(shù)據(jù)包括常規(guī)疊后或疊前偏移處理的地震資料。(2)利用密度測井曲線與泥頁巖含氣量數(shù)據(jù)之間的高相關(guān)性，建立基于密度測井曲線的含氣量計算模型：VGAS＝f(DEN)(1)式中，VGAS是含氣量，m/t；DEN為密度測井曲線，g/cm；f為含氣量與密度之間的函數(shù)關(guān)系式。(3)開展疊后多屬性反演，獲取含氣量計算模型所需的密度曲線數(shù)據(jù)體。(4)含氣量預測：將步驟(3)得到的密度曲線數(shù)據(jù)體代入到式(1)中，獲得泥頁巖含氣量。
【專利說明】
一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測方法
技術(shù)領(lǐng)域
[0001]本發(fā)明涉及本發(fā)明屬于石油天然氣勘探領(lǐng)域，具體涉及一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測方法。
【背景技術(shù)】
[0002]目前全球頁巖氣資源量已經(jīng)超過全球常規(guī)天然氣資源量，成為非常規(guī)油氣領(lǐng)域的研究熱點。北美的美國和加拿大已經(jīng)擁有數(shù)十年的勘探開發(fā)歷史，泥頁巖的地質(zhì)認識、勘探配套技術(shù)以及工程技術(shù)方面都較為成熟。隨著我國對能源的需求日益擴大，以頁巖氣為代表的非常規(guī)油氣領(lǐng)域的研究在近兩年逐漸成為熱點，也不斷在勘探實踐中取得突破和新的認識。綜合國內(nèi)外近些年勘探開發(fā)經(jīng)驗可知，含氣量是頁巖氣藏評價的關(guān)鍵指標，是頁巖氣藏取得商業(yè)發(fā)現(xiàn)的關(guān)鍵影響因素，這主要是因為含氣量是烴源條件(有機碳含量、熱演化程度)、儲集條件(孔隙及裂縫發(fā)育程度)、保存條件(即頂?shù)装鍡l件、區(qū)域蓋層條件、后期構(gòu)造改造強度和構(gòu)造樣式等)等的綜合體現(xiàn)。目前，已經(jīng)在實驗室測量以及測井評價含氣量方面取得較大進展，如聶海寬等在《頁巖氣聚集條件及含氣量計算一一以四川盆地及其周緣下古生界為例》提出了一種利用多元線性回歸計算頁巖含氣量的公式，但該公式僅能計算井點位置處的含氣量，并不能描述含氣量的橫向變化。

【發(fā)明內(nèi)容】

[0003]本發(fā)明的目的即在于克服現(xiàn)有技術(shù)的不足，提供一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測方法，解決并不能描述含氣量的橫向變化的問題。
[0004]本發(fā)明的通過下述技術(shù)方案實現(xiàn):一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測方法，所述方法包括:
(I)測量數(shù)據(jù)，地質(zhì)數(shù)據(jù)包括錄井、試氣以及泥頁巖巖心含氣量測試數(shù)據(jù)；測井數(shù)據(jù)包括聲波時差、密度以及經(jīng)巖心含氣量測試數(shù)據(jù)刻度的含氣量測井解釋曲線；地震數(shù)據(jù)包括常規(guī)疊后或疊前偏移處理的地震資料。
[0005](2)利用密度測井曲線與泥頁巖含氣量數(shù)據(jù)之間的高相關(guān)性，建立基于密度測井曲線的含氣量計算模型:
VGAS = f(DEN) (I)
式中，VGAS是含氣量，m /t ；DEN為密度測井曲線，g/cm ;f為含氣量與密度之間的函數(shù)關(guān)系式。
[0006](3)開展疊后多屬性反演，獲取含氣量計算模型所需的密度曲線數(shù)據(jù)體。
[0007](4)含氣量預測:將步驟(3)得到的密度曲線數(shù)據(jù)體代入到式(I)中，獲得泥頁巖含氣量。
[0008]進一步的，所述步驟(3)是這樣實現(xiàn)的:以含氣量計算模型中的密度測井曲線作為培訓目標曲線，尋找最佳的地震屬性與波阻抗反演數(shù)據(jù)體組合，使得在井點位置利用這些屬性計算得到的目標曲線與實際密度測井曲線相關(guān)性最高，具體包括:輸入密度測井曲線、井旁地震屬性、疊后反演曲線。
[0009]統(tǒng)計尋找密度測井曲線與井旁地震屬性、疊后波阻抗反演曲線之間的關(guān)系；進行盲井交互驗證，確定最佳地震屬性與波阻抗反演數(shù)據(jù)體組合以及式(2)所述函數(shù)關(guān)系。
[0010]DEN(x，y，t)= f[Attribute l(x,y,t),Attribute2(x,y,t),......,Attributem
(x,y,t)l (2)
式中，DEN(x，y，t)為步驟(2)中含氣量計算模型所需的密度測曲線。
[0011]Attributei(x,y,t)為地震屬性或反演數(shù)據(jù)為DEN(x，y，t)與最佳屬性之間的函數(shù)關(guān)系;m為屬性個數(shù)。
[0012]將統(tǒng)計關(guān)系應用到地震數(shù)據(jù)體，獲得含氣量計算模型所需的密度曲線數(shù)據(jù)體。
[0013]本發(fā)明與現(xiàn)有技術(shù)相比，具有如下的優(yōu)點和有益效果:通過本發(fā)明所述技術(shù)方法的實施可完成對優(yōu)質(zhì)泥頁巖含氣量橫向變化的精細描述，進而減少鉆井風險，提高頁巖氣勘探效?。
【附圖說明】
[0014]此處所說明的附圖用來提供對本發(fā)明實施例的進一步理解，構(gòu)成本申請的一部分，并不構(gòu)成對本發(fā)明實施例的限定。在附圖中:
圖1為本發(fā)明一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測步驟示意圖；
圖2為本發(fā)明含氣量敏感曲線反演流程圖。
【具體實施方式】
[0015]為使本發(fā)明的目的、技術(shù)方案和優(yōu)點更加清楚明白，下面結(jié)合實施例和附圖，對本發(fā)明作進一步的詳細說明，本發(fā)明的示意性實施方式及其說明僅用于解釋本發(fā)明，并不作為對本發(fā)明的限定。
[0016]如圖1所示的，本發(fā)明一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測方法所述方法包括:
(1)測量數(shù)據(jù)，地質(zhì)數(shù)據(jù)包括錄井、試氣以及泥頁巖巖心含氣量測試數(shù)據(jù)；測井數(shù)據(jù)包括聲波時差、密度以及經(jīng)巖心含氣量測試數(shù)據(jù)刻度的含氣量測井解釋曲線；地震數(shù)據(jù)包括常規(guī)疊后或疊前偏移處理的地震資料；
(2)利用密度測井曲線與泥頁巖含氣量數(shù)據(jù)之間的高相關(guān)性，建立基于密度測井曲線的含氣量計算模型:
VGAS = f(DEN) (I)
式中，VGAS是含氣量，m /t ；DEN為密度測井曲線，g/cm;f為含氣量與密度之間的函數(shù)關(guān)系式；
(3)開展疊后多屬性反演，獲取含氣量計算模型所需的密度曲線數(shù)據(jù)體；
(4)含氣量預測:將步驟(3)得到的密度曲線數(shù)據(jù)體代入到式(I)中，獲得泥頁巖含氣量。
[0017]所述步驟(3)是這樣實現(xiàn)的:以含氣量計算模型中的密度測井曲線作為培訓目標曲線，尋找最佳的地震屬性與波阻抗反演數(shù)據(jù)體組合，使得在井點位置利用這些屬性計算得到的目標曲線與實際密度測井曲線相關(guān)性最高，具體包括:輸入密度測井曲線、井旁地震屬性、疊后反演曲線；
統(tǒng)計尋找密度測井曲線與井旁地震屬性、疊后波阻抗反演曲線之間的關(guān)系；進行盲井交互驗證，確定最佳地震屬性與波阻抗反演數(shù)據(jù)體組合以及式(2)所述函數(shù)關(guān)系；
DEN(x，y，t)= f[Attribute l(x,y,t),Attribute2(x,y,t),......,Attributem(x,y,
t)] (2)
式中，DEN(x，y，t)為步驟(2)中含氣量計算模型所需的密度測井曲線；
Attributei(x，y，t)為地震屬性或反演數(shù)據(jù)；f為DEN(x，y，t)與最佳屬性之間的函數(shù)關(guān)系為屬性個數(shù)；
將統(tǒng)計關(guān)系應用到地震數(shù)據(jù)體，獲得含氣量計算模型所需的密度曲線數(shù)據(jù)體。
[0018]圖2為含氣量計算模型所需密度目標測井曲線反演流程，實現(xiàn)方法為:①對每一屬性(瞬時頻率、瞬時振幅、瞬時相位等地震屬性)計算用其預測密度測井曲線時的誤差，誤差最小的即為屬性1，然后將屬性I與其它屬性組成屬性對，再次計算預測誤差，誤差最小的屬性對中的另一屬性為屬性2，依此類推便可選出所需的屬性；②由于步驟①中的屬性對可以無限制的增加，這樣會存在過擬合的問題，因此需要盲井交互驗證，即在每增加一種屬性時依次去掉每一口井后計算其預測誤差，所有井的平均預測誤差最小時的屬性與波阻抗組合即為最佳組合;③將井曲線與屬性組合之間的關(guān)系應用到整個數(shù)據(jù)體就可以計算得到密度數(shù)據(jù)體。
[0019]以上所述的【具體實施方式】，對本發(fā)明的目的、技術(shù)方案和有益效果進行了進一步詳細說明，所應理解的是，以上所述僅為本發(fā)明的【具體實施方式】而已，并不用于限定本發(fā)明的保護范圍，凡在本發(fā)明的精神和原則之內(nèi)，所做的任何修改、等同替換、改進等，均應包含在本發(fā)明的保護范圍之內(nèi)。
【主權(quán)項】
1.一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測方法，其特征在于:所述方法包括: (1)測量數(shù)據(jù)，地質(zhì)數(shù)據(jù)包括錄井、試氣以及泥頁巖巖心含氣量測試數(shù)據(jù)；測井數(shù)據(jù)包括聲波時差、密度以及經(jīng)巖心含氣量測試數(shù)據(jù)刻度的含氣量測井解釋曲線；地震數(shù)據(jù)包括常規(guī)疊后或疊前偏移處理的地震資料； (2)利用密度測井曲線與泥頁巖含氣量數(shù)據(jù)之間的高相關(guān)性，建立基于密度測井曲線的含氣量計算模型: VGAS = f(DEN) (I) 式中，VGAS是含氣量，m /t ;DEN為密度測井曲線，g/cm;f為含氣量與密度之間的函數(shù)關(guān)系式； (3)開展疊后多屬性反演，獲取含氣量計算模型所需的密度曲線數(shù)據(jù)體； (4)含氣量預測:將步驟(3)得到的密度曲線數(shù)據(jù)體代入到式(I)中，獲得泥頁巖含氣量。2.根據(jù)權(quán)利要求1所述的一種利用地震數(shù)據(jù)進行泥頁巖含氣量預測方法，其特征在于:所述步驟(3)是這樣實現(xiàn)的:以含氣量計算模型中的密度測井曲線作為培訓目標曲線，尋找最佳的地震屬性與波阻抗反演數(shù)據(jù)體組合，使得在井點位置利用這些屬性計算得到的目標曲線與實際密度測井曲線相關(guān)性最高，具體包括:輸入密度測井曲線、井旁地震屬性、置后反演曲線； 統(tǒng)計尋找密度測井曲線與井旁地震屬性、疊后波阻抗反演曲線之間的關(guān)系；進行盲井交互驗證，確定最佳地震屬性與波阻抗反演數(shù)據(jù)體組合以及式(2)所述函數(shù)關(guān)系； DEN(x,y,t)= f[Attribute l(x,y,t),Attribute2(x,y,t),......，Attributem(x，y，t)] (2) 式中，DEN(x，y，t)為步驟(2)中含氣量計算模型所需的密度測井曲線； Attributei(x,y,t)為地震屬性或反演數(shù)據(jù)為DEN(x，y，t)與最佳屬性之間的函數(shù)關(guān)系為屬性個數(shù)； 將統(tǒng)計關(guān)系應用到地震數(shù)據(jù)體，獲得含氣量計算模型所需的密度曲線數(shù)據(jù)體。
【文檔編號】G01V1/30GK105842734SQ201610152263
【公開日】2016年8月10日
【申請日】2016年3月17日
【發(fā)明人】李其鑫, 張鑫, 楊波
【申請人】成都創(chuàng)源油氣技術(shù)開發(fā)有限公司