基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明專利涉及測(cè)量技術(shù)，特別涉及基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方 法。
【背景技術(shù)】
[0002] 在現(xiàn)有的技術(shù)中，只有基于測(cè)井資料的孔隙扁率反演方法，主要是針對(duì)單井的巖 石等效孔隙縱橫比進(jìn)行反演，該方法只能獲得井點(diǎn)處的孔隙結(jié)構(gòu)信息，不能有效地獲得反 映研宄區(qū)孔隙結(jié)構(gòu)信息的三維數(shù)據(jù)。

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

[0003] 本發(fā)明的目的在于提供基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，以有效地 獲得反映研宄區(qū)孔隙結(jié)構(gòu)信息的三維數(shù)據(jù)。
[0004] 本發(fā)明提供一種基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，包括：設(shè)定當(dāng)前 采樣點(diǎn)的線號(hào)、道號(hào)和時(shí)間采樣點(diǎn)號(hào)，提取當(dāng)前采樣點(diǎn)的孔隙度Φ'、基質(zhì)礦物的等效體積 模量K' ^、基質(zhì)礦物的等效剪切模量μ' ^、砂巖的百分比含量v's、泥巖的百分比含量或泥 質(zhì)含量V'。、密度P '和縱波速度t s;設(shè)定當(dāng)前采樣點(diǎn)的砂巖孔隙扁率a ' s的初始值公 式、搜索范圍和搜索間隔；并設(shè)定每個(gè)采樣點(diǎn)的泥巖孔隙的扁率α '。的初始值、搜索范圍 和搜索間隔；設(shè)定完當(dāng)前采樣點(diǎn)的砂巖孔隙扁率〇、和泥巖孔隙扁率α'。的初始值后，計(jì) 算孔隙度為Φ'時(shí)的干巖石剪切模量μ' ty;計(jì)算當(dāng)前采樣點(diǎn)的橫波速度V' s >進(jìn)行當(dāng) 前采樣點(diǎn)的計(jì)算的兩種橫波速度V' s (l和V' s的誤差分析，并保存當(dāng)前砂巖孔隙的扁率 a ' s和泥巖孔隙的扁率α '。的值所對(duì)應(yīng)的誤差值σ。
[0005] 優(yōu)選的，還包括：修改當(dāng)前采樣點(diǎn)的砂巖孔隙扁率α ' 3和泥巖孔隙扁率α '。的 值，直到搜索范圍內(nèi)所有的扁率值計(jì)算完畢；尋找砂巖孔隙的扁率α ' 3和泥巖孔隙的扁 率α'。的值所對(duì)應(yīng)的最小誤差值〇，并輸出對(duì)應(yīng)的砂巖孔隙的扁率a' s、泥巖孔隙的扁 率α '。和φ s'、Φ。'作為當(dāng)前采樣點(diǎn)的計(jì)算結(jié)果。
[0006] 優(yōu)選的，還包括：修改當(dāng)前采樣點(diǎn)的線號(hào)、道號(hào)和時(shí)間采樣點(diǎn)號(hào)，直到所有的采樣 點(diǎn)計(jì)算完畢，并輸出砂巖孔隙的扁率a s及其百分比含量Φ s、泥巖孔隙的扁率α。及其百分 比含量Φ。的三維數(shù)據(jù)體；對(duì)a s、φ3、α。、Φ。進(jìn)行三維中值濾波，去除異常值的影響，并作 為最終的結(jié)果輸出。
[0007] 優(yōu)選的，還包括：輸入砂泥巖研宄區(qū)的測(cè)井參數(shù)：縱波速度\、橫波速度Vs、孔隙度 Φ和密度P ;統(tǒng)計(jì)得到\與P的關(guān)系式，擬合得到系數(shù)α和β ;再利用雙波法統(tǒng)計(jì)φ與 \和V 3的關(guān)系式，擬合得到系數(shù)A、B、C。
[0008] 本發(fā)明的有益效果在于：基于砂泥巖研宄區(qū)的三維地震資料和測(cè)井資料，計(jì)算反 映砂泥巖研宄區(qū)的孔隙類型的三維數(shù)據(jù)體，獲得地下每個(gè)采樣點(diǎn)的孔隙類型信息（包括組 成孔隙的各種孔隙類型的扁率及其百分比含量），從而有效地描述砂泥巖儲(chǔ)層的孔隙類型， 為砂泥巖儲(chǔ)層預(yù)測(cè)和鉆探提供重要的參考數(shù)據(jù)。
【附圖說明】
[0009] 圖1是基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法的流程圖。
【具體實(shí)施方式】
[0010] 為使本發(fā)明的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，下文中將結(jié)合實(shí)施例進(jìn)行詳 細(xì)說明。
[0011] 如圖1所示，本發(fā)明基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法具體包括如下 步驟：
[0012] 步驟（1)，輸入砂泥巖研宄區(qū)的測(cè)井參數(shù)：縱波速度（Vp)、橫波速度（v s)、孔隙度 (Φ)和密度（P)。利用Gardner公式分別統(tǒng)計(jì)得到\與P的關(guān)系式（見公式4-1)，擬合 得到系數(shù)α和β ;再利用雙波法統(tǒng)計(jì)φ與VjPV 3的關(guān)系式（見公式4-2)，擬合得到系數(shù) A、B、C ;
【主權(quán)項(xiàng)】
1. 一種基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，其特征在于，包括： 設(shè)定當(dāng)前采樣點(diǎn)的線號(hào)、道號(hào)和時(shí)間采樣點(diǎn)號(hào)，提取當(dāng)前采樣點(diǎn)的孔隙度巾'、基質(zhì)礦 物的等效體積模量K'^、基質(zhì)礦物的等效剪切模量y'^、砂巖的百分比含量v's、泥巖的百 分比含量或泥質(zhì)含量v'。、密度P'和縱波速度V' s; 設(shè)定當(dāng)前采樣點(diǎn)的砂巖孔隙扁率a's的初始值公式、搜索范圍和搜索間隔；并設(shè)定 每個(gè)采樣點(diǎn)的泥巖孔隙的扁率a'。的初始值、搜索范圍和搜索間隔； 設(shè)定完當(dāng)前采樣點(diǎn)的砂巖孔隙扁率a、和泥巖孔隙扁率a'。的初始值后，計(jì)算孔隙度 為小'時(shí)的干巖石剪切模量U'd,y; 計(jì)算當(dāng)前采樣點(diǎn)的橫波速度Vs' > 進(jìn)行當(dāng)前采樣點(diǎn)的計(jì)算的兩種橫波速度V's(l和's的誤差分析，并保存當(dāng)前砂巖孔隙的 扁率a's和泥巖孔隙的扁率a'。的值及其所對(duì)應(yīng)的誤差值〇。
2. 如權(quán)利要求1所述的基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，其特征在于， 還包括： 修改當(dāng)前采樣點(diǎn)的砂巖孔隙扁率a' 3和泥巖孔隙扁率a'。的值，直到搜索范圍內(nèi) 所有的扁率值計(jì)算完畢； 尋找砂巖孔隙的扁率a's和泥巖孔隙的扁率a'。的值所對(duì)應(yīng)的最小誤差值〇，并 輸出對(duì)應(yīng)的砂巖孔隙的扁率a's、泥巖孔隙的扁率a'。和（})'s、巾'。作為當(dāng)前采樣點(diǎn)的計(jì) 算結(jié)果。
3. 如權(quán)利要求2所述的基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，其特征在于， 還包括： 修改當(dāng)前采樣點(diǎn)的線號(hào)、道號(hào)和時(shí)間采樣點(diǎn)號(hào)，直到所有的采樣點(diǎn)計(jì)算完畢，并輸出砂 巖孔隙的扁率a3及其百分比含量巾s、泥巖孔隙的扁率a。及其百分比含量巾。的三維數(shù) 據(jù)體； 對(duì)as、（})s、a。、(}>。進(jìn)行三維中值濾波，去除異常值的影響，并作為最終的結(jié)果輸出。
4. 如權(quán)利要求3所述的基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，其特征在于， 還包括：輸入砂泥巖研宄區(qū)的測(cè)井參數(shù)：縱波速度Vp、橫波速度Vs、孔隙度巾和密度P;統(tǒng) 計(jì)得到\與P的關(guān)系式，擬合得到系數(shù)a和0 ;再利用雙波法統(tǒng)計(jì)巾與VjPVs的關(guān)系 式，擬合得到系數(shù)A、B、C。
5. 如權(quán)利要求4所述的基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，其特征在于， 還包括：輸入砂泥巖研宄區(qū)的基于疊后反演得到的伽瑪GR三維地震數(shù)據(jù)，并采計(jì)算泥質(zhì)含 量即泥巖的百分比含量V。的三維地震數(shù)據(jù)。
6. 如權(quán)利要求5所述的基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，其特征在于， 還包括：輸入基于疊前彈性反演得到縱波阻抗IP、橫波阻抗Is的三維地震數(shù)據(jù)；計(jì)算Vp、Vs 和P的三維地震數(shù)據(jù)；計(jì)算孔隙度巾的三維地震數(shù)據(jù)。
7. 如權(quán)利要求6所述的基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，其特征在于， 還包括：輸入研宄區(qū)的砂巖和泥巖的剪切模量，采計(jì)算基質(zhì)礦物的等效體積模量Ktl和等效 剪切模量Utl的三維數(shù)據(jù)體。
8. 如權(quán)利要求7所述的基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，其特征在于， 還包括：計(jì)算砂巖孔隙在總孔隙中的百分比含量巾 3和泥巖孔隙在總孔隙中的百分比含量 的三維數(shù)據(jù)體。
【專利摘要】本發(fā)明提供一種基于三維地震數(shù)據(jù)的砂泥巖孔隙類型反演方法，包括：設(shè)定當(dāng)前采樣點(diǎn)的線號(hào)、道號(hào)和時(shí)間采樣點(diǎn)號(hào)，提取當(dāng)前采樣點(diǎn)的孔隙度φ'、基質(zhì)礦物的等效體積模量K0'、基質(zhì)礦物的等效剪切模量μ'0、砂巖的百分比含量v's、泥巖的百分比含量即泥質(zhì)含量v'c、密度ρ'和縱波速度V's；設(shè)定當(dāng)前采樣點(diǎn)的砂巖孔隙扁率α's的初始值公式、搜索范圍和搜索間隔；并設(shè)定每個(gè)采樣點(diǎn)的泥巖孔隙的扁率α'c的初始值、搜索范圍和搜索間隔；設(shè)定完當(dāng)前采樣點(diǎn)的砂巖孔隙扁率α's和泥巖孔隙扁率α'c的初始值后，計(jì)算孔隙度為φ'時(shí)的干巖石剪切模量μ'dry；計(jì)算當(dāng)前采樣點(diǎn)的橫波速度V's_0；進(jìn)行當(dāng)前采樣點(diǎn)的計(jì)算的兩種橫波速度V's_0和V's的誤差分析，并保存當(dāng)前砂巖孔隙的扁率α's和泥巖孔隙的扁率α'c的值及其所對(duì)應(yīng)的誤差值σ。
【IPC分類】G01V1-30, G01N15-08
【公開號(hào)】CN104792684
【申請(qǐng)?zhí)枴緾N201510189404
【發(fā)明人】熊曉軍, 譚峰, 侯秋平, 李良泉, 黃智勇
【申請(qǐng)人】成都理工大學(xué)
【公開日】2015年7月22日
【申請(qǐng)日】2015年4月18日