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

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

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