一種識(shí)別高鹽富鉀區(qū)注水開發(fā)油藏碎屑巖巖性的方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于儲(chǔ)層評(píng)價(jià)技術(shù)領(lǐng)域，涉及一種識(shí)別高鹽富鉀區(qū)注水開發(fā)油藏碎肩巖巖 性的方法。
【背景技術(shù)】
[0002] 測(cè)井巖性識(shí)別是準(zhǔn)確求取儲(chǔ)層參數(shù)，進(jìn)行儲(chǔ)層評(píng)價(jià)工作的基礎(chǔ)和關(guān)鍵。注水開發(fā) 油藏碎肩巖巖性一般通過自然電位(SP)、自然伽馬(GR)、聲波時(shí)差(AC)、密度(DEN)和電阻 率(RT)等測(cè)井資料進(jìn)行識(shí)別。
[0003] 常規(guī)的識(shí)別方法主要包括交會(huì)圖法、神經(jīng)網(wǎng)絡(luò)法和曲線重疊法等，但在高鹽富鉀 地區(qū)存在識(shí)別率低，不能連續(xù)識(shí)別地層巖性的缺陷。灘海地區(qū)地質(zhì)條件復(fù)雜，且具有高鹽富 鉀的特點(diǎn)，鹽水泥漿基液使得自然電位反向出現(xiàn)正異常，電阻率值比正常情況下降低從而 無法真實(shí)反映原狀地層特性;地層中存在高鉀礦物，又使得自然伽馬曲線識(shí)別精度降低。因 此單條測(cè)井曲線不能有效反映碎肩巖巖性特征，常規(guī)方法也無法準(zhǔn)確區(qū)分不同類巖性。
[0004] 為了更好地識(shí)別高鹽富鉀區(qū)注水開發(fā)油藏碎肩巖巖性并開展后續(xù)儲(chǔ)層評(píng)價(jià)工作， 放大和提取各類測(cè)井參數(shù)的差異，綜合多參數(shù)多資料進(jìn)行巖性識(shí)別就變得非常有意義。

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

[0005] 本發(fā)明的目的是提供一種識(shí)別高鹽富鉀區(qū)注水開發(fā)油藏碎肩巖巖性的方法，基于 數(shù)理統(tǒng)計(jì)原理，建立巖性識(shí)別方程，構(gòu)建差值識(shí)別函數(shù)，擬合巖性識(shí)別綜合曲線，達(dá)到準(zhǔn)確 識(shí)別高鹽富鉀區(qū)注水開發(fā)油藏碎肩巖巖性的目的。
[0006] 本發(fā)明所采用的技術(shù)方案是，一種識(shí)別高鹽富鉀區(qū)注水開發(fā)油藏碎肩巖巖性的方 法，按照以下步驟進(jìn)行：
[0007] 步驟1，統(tǒng)計(jì)已知巖性的巖心分析化驗(yàn)資料并將巖性分類;采用數(shù)理統(tǒng)計(jì)方法，剔 除非相關(guān)變量，選取對(duì)巖性反映比較敏感的η個(gè)測(cè)井參數(shù)；
[0008] 步驟2,基于數(shù)理統(tǒng)計(jì)原理，建立巖性識(shí)別權(quán)重函數(shù)F，計(jì)算出的權(quán)重函數(shù)根據(jù)η個(gè) 測(cè)井參數(shù)對(duì)巖性的貢獻(xiàn)率賦予不同的函數(shù)特征值，提取與各組有最大可能多重相關(guān)的、函 數(shù)特征值最大的兩個(gè)權(quán)重函數(shù)，即第一權(quán)重函數(shù)和第二權(quán)重函數(shù)：
[0009] Fi = a〇〇+a〇i · χιΗ-----i~a〇n · xn (1)
[0010] F2 = ai〇+aii · xi+."+ain · xn (2)
[0011] 其中，F(xiàn)i、F2為第一和第二權(quán)重函數(shù)；χ^χπη為選取的n個(gè)測(cè)井參數(shù)變量;a00， a〇i…a〇n和ai〇，a2〇…ain分別對(duì)應(yīng)兩個(gè)權(quán)重函數(shù)的判別系數(shù)；
[0012] 步驟3,將給出的第一權(quán)重函數(shù)和第二權(quán)重函數(shù)用于計(jì)算各樣品在降維后空間坐 標(biāo)系中的坐標(biāo)值，然后計(jì)算出各類別在低維空間中的中心坐標(biāo)，根據(jù)各樣品點(diǎn)距離各類別 重心的距離構(gòu)造出每種巖性的識(shí)別函數(shù)，即：
[0013] fi(X) =bi〇+bii · xi+---+bin · xn (3)
[0014]其中，h為第i類巖性識(shí)別函數(shù);b^bn···!^為第i類巖性識(shí)別函數(shù)的判別系數(shù)；
[0015] 步驟4,對(duì)于兩種巖性的判別，將砂巖巖性識(shí)別函數(shù)減去泥巖巖性識(shí)別函數(shù)，得出 差值識(shí)別函數(shù)，該函數(shù)即為兩種巖性的綜合識(shí)別曲線，記為P(x);即：
[0016] p⑴=A f = f砂-f泥=Δ bo+Δ bi · χι+···+Δ bn · Xn (4)
[0017] 其中，P(x)為兩種巖性的綜合識(shí)別曲線;f砂為砂巖巖性識(shí)別函數(shù)，f泥為泥巖巖性識(shí) 另IJ函數(shù);Χ1，χ 2···χΑ選取的η個(gè)測(cè)井參數(shù)變量；△ bo，Δ br·· Δ bn為兩種巖性綜合識(shí)別曲線的 判別系數(shù)；
[0018] 把待判測(cè)井?dāng)?shù)據(jù)體的觀測(cè)值Χ=(χ(1)χ(2)···χ(η))代入上述差值識(shí)別函數(shù)計(jì)算出巖 性綜合識(shí)別曲線值Poo，若Poo 2 〇則為砂巖巖性，相反綜合曲線值Poo < 〇則為泥巖巖性； [0019]步驟5,對(duì)于Μ類巖性的判別，把樣品的觀測(cè)值Χ=(χ (1)χ(2)···χ(η))代入式(3)中，分 別求得fi( x)，f2(x)"_fM(x)，取最大值記為fg( x)=max{fi(x)}，1 < i SM;由判別準(zhǔn)則知，巖性識(shí)別 函數(shù)值大為相應(yīng)巖性的概率高，則認(rèn)為待判測(cè)井?dāng)?shù)據(jù)體XefgW對(duì)應(yīng)的巖性所屬類型，同時(shí) 巖性綜合識(shí)別曲線Q(x)賦值為m;對(duì)Μ類不同巖性識(shí)別結(jié)果賦值，從而擬合出Μ類巖性的綜合 識(shí)別曲線Q(x);即：
[0021]其中，Quo為Μ類巖性的綜合識(shí)別曲線;f1(x)，f2(x)…fM(x)為Μ類巖性識(shí)別函數(shù);f g(x) 為代入待判數(shù)據(jù)體后Μ類巖性識(shí)別函數(shù)的最大值;?Μ,πινΠΗ為Μ類巖性綜合識(shí)別曲線賦值， nu的賦值區(qū)間為-Μ < mi < Μ。
[0022]進(jìn)一步的，所述步驟1中，η個(gè)測(cè)井參數(shù)包括自然電位SP、自然伽馬GR、聲波時(shí)差A(yù)C、 巖性密度DEN、電阻率RT。
[0023] 本發(fā)明的有益效果是克服常規(guī)巖性識(shí)別方法的不足，提供一種識(shí)別高鹽富鉀區(qū)注 水開發(fā)油藏碎肩巖巖性的方法，即充分利用測(cè)井、巖心分析化驗(yàn)等資料，以數(shù)理統(tǒng)計(jì)原理為 基礎(chǔ)，建立巖性識(shí)別函數(shù)，并擬合單井巖性識(shí)別綜合曲線，利用該曲線識(shí)別巖性效果較好， 自檢符合率高達(dá)90%。本發(fā)明能很好地考慮多種測(cè)井參數(shù)，進(jìn)而擬合出一條準(zhǔn)確的巖性綜 合識(shí)別曲線，計(jì)算結(jié)果合理且準(zhǔn)確率高，實(shí)現(xiàn)對(duì)高鹽富鉀區(qū)注水開發(fā)油藏全井段碎肩巖巖 性的連續(xù)判別，可為類似地區(qū)的巖性識(shí)別提供借鑒。
【附圖說明】
[0024] 圖1是本發(fā)明實(shí)施例中研究區(qū)兩種巖性的綜合識(shí)別曲線。
[0025]圖2是本發(fā)明實(shí)施例中研究區(qū)四種巖性的綜合識(shí)別曲線。
【具體實(shí)施方式】
[0026] 為了使本發(fā)明的目的、技術(shù)方案及優(yōu)點(diǎn)更加清楚明白，以下結(jié)合附圖及實(shí)施例，對(duì) 本發(fā)明進(jìn)行進(jìn)一步詳細(xì)說明。應(yīng)當(dāng)理解，此處所描述的具體實(shí)施例僅僅用以解釋本發(fā)明，并 不用于限定本發(fā)明。
[0027] -種識(shí)別高鹽富鉀區(qū)注水開發(fā)油藏碎肩巖巖性的方法，具體按照以下步驟進(jìn)行： [0028]步驟1，統(tǒng)計(jì)已知巖性的巖心分析化驗(yàn)資料并將巖性分類。采用數(shù)理統(tǒng)計(jì)方法，剔 除非相關(guān)變量，選取對(duì)巖性反映比較敏感的η個(gè)測(cè)井參數(shù)（自然電位(SP)、自然伽馬(GR)、聲 波時(shí)差(AC)、巖性密度(DEN)、電阻率(RT)…）；
[0029]步驟2，基于數(shù)理統(tǒng)計(jì)原理，建立巖性識(shí)別權(quán)重函數(shù)F(包括Fi、F2、Fr··)，計(jì)算出的 權(quán)重函數(shù)可以根據(jù)η個(gè)測(cè)井參數(shù)對(duì)巖性的貢獻(xiàn)率賦予不同的函數(shù)特征值，提取與各組有最 大可能多重相關(guān)的、函數(shù)特征值最大的兩個(gè)權(quán)重函數(shù)，即第一權(quán)重函數(shù)和第二權(quán)重函數(shù)： [0030] Fi = aoo+aoi · xi+---+a〇n · xn (1)
[0031] F2 = aio+aii · xi+."+ain · Xn (2)
[0032] 其中，F(xiàn)i、F2為第一和第二權(quán)重函數(shù)；X1，x 2…xn*選取的n個(gè)測(cè)井參數(shù)變量；aoo， aoi…a〇n和aio，a2〇…ain分別對(duì)應(yīng)兩個(gè)權(quán)重函數(shù)的判別系數(shù)。
[0033] 步驟3,將給出的第一權(quán)重函數(shù)和第二權(quán)重函數(shù)用于計(jì)算各樣品在降維后空間坐 標(biāo)系中的坐標(biāo)值，然后計(jì)算出各類別在低維空間中的中心坐標(biāo)，根據(jù)各樣品點(diǎn)距離各類別 重心的距離構(gòu)造出每種巖性的識(shí)別函數(shù)，即：
[0034] fi(X) =bi〇+bii · χι+...+bin · Xn (3)
[0035] 其中，fi為第i類巖性識(shí)別函數(shù);幻，^_心為選取的11個(gè)測(cè)井參數(shù)變量;1^，1^1···!^ 為第i類巖性識(shí)別函數(shù)的判別系數(shù)。
[0036] 步驟4，對(duì)于兩種巖性的判別，將砂巖巖性識(shí)別函數(shù)減去泥巖巖性識(shí)別函數(shù)，得
[0037] 出差值識(shí)別函數(shù)，該函數(shù)即為兩種巖性的綜合識(shí)別曲線，記為P〇〇;即：
[0038] p⑴=A f = f砂-f泥=Δ bo+Δ bi · χι+···+Δ bn · Xn (4)
[0039] 其中，P(x)為兩種巖性的綜合識(shí)別曲線;f砂為砂巖巖性識(shí)別函數(shù)，f泥為泥巖巖性識(shí) 另IJ函數(shù);Χ1，χ 2···χΑ選取的η個(gè)測(cè)井參數(shù)變量；△ bo，Δ br·· Δ bn為兩種巖性綜合識(shí)別曲線的 判別系數(shù)。
[0040] 把待判測(cè)井?dāng)?shù)據(jù)體的觀測(cè)值X=(x(1)x(2)…x(n))代入上述差值識(shí)