本發明專利技術公開了一種識別高鹽富鉀區注水開發油藏碎屑巖巖性的方法，統計已知巖性的巖心分析化驗資料并將巖性分類，選取對巖性反映比較敏感的n個測井參數；基于數理統計原理，建立不同巖性識別方程，再將單井未取心段測井數據體分別代入各類巖性識別方程，根據識別準則，識別方程中函數值最大者則為相應巖性的概率最高，該巖性則確定為相應的巖性。通過擬合綜合判別曲線識別高鹽富鉀區注水開發油藏碎屑巖巖性，實現對多類巖性的準確判別，可建立連續的地層巖性剖面，自檢符合率達90％，識別效果較好，對類似地區巖性識別具有一定的借鑒意義。
【技術實現步驟摘要】
本專利技術屬于儲層評價
，涉及一種識別高鹽富鉀區注水開發油藏碎肩巖巖 性的方法。
技術介紹
測井巖性識別是準確求取儲層參數，進行儲層評價工作的基礎和關鍵。注水開發 油藏碎肩巖巖性一般通過自然電位(SP)、自然伽馬(GR)、聲波時差(AC)、密度(DEN)和電阻 率(RT)等測井資料進行識別。 常規的識別方法主要包括交會圖法、神經網絡法和曲線重疊法等，但在高鹽富鉀 地區存在識別率低，不能連續識別地層巖性的缺陷。灘海地區地質條件復雜，且具有高鹽富 鉀的特點，鹽水泥漿基液使得自然電位反向出現正異常，電阻率值比正常情況下降低從而 無法真實反映原狀地層特性;地層中存在高鉀礦物，又使得自然伽馬曲線識別精度降低。因 此單條測井曲線不能有效反映碎肩巖巖性特征，常規方法也無法準確區分不同類巖性。 為了更好地識別高鹽富鉀區注水開發油藏碎肩巖巖性并開展后續儲層評價工作， 放大和提取各類測井參數的差異，綜合多參數多資料進行巖性識別就變得非常有意義。
技術實現思路
本專利技術的目的是提供一種識別高鹽富鉀區注水開發油藏碎肩巖巖性的方法，基于 數理統計原理，建立巖性識別方程，構建差值識別函數，擬合巖性識別綜合曲線，達到準確 識別高鹽富鉀區注水開發油藏碎肩巖巖性的目的。 本專利技術所采用的技術方案是，一種識別高鹽富鉀區注水開發油藏碎肩巖巖性的方 法，按照以下步驟進行： 步驟1，統計已知巖性的巖心分析化驗資料并將巖性分類;采用數理統計方法，剔 除非相關變量，選取對巖性反映比較敏感的η個測井參數； 步驟2,基于數理統計原理，建立巖性識別權重函數F，計算出的權重函數根據η個 測井參數對巖性的貢獻率賦予不同的函數特征值，提取與各組有最大可能多重相關的、函 數特征值最大的兩個權重函數，即第一權重函數和第二權重函數： Fi = a〇〇+a〇i · χιΗ-----i~a〇n · xn (1) F2 = ai〇+aii · xi+."+ain · xn (2) 其中，Fi、F2為第一和第二權重函數；χ^χπη為選取的n個測井參數變量;a00， a〇i…a〇n和ai〇，a2〇…ain分別對應兩個權重函數的判別系數； 步驟3,將給出的第一權重函數和第二權重函數用于計算各樣品在降維后空間坐 標系中的坐標值，然后計算出各類別在低維空間中的中心坐標，根據各樣品點距離各類別 重心的距離構造出每種巖性的識別函數，即： fi(X) =bi〇+bii · xi+---+bin · xn (3)其中，h為第i類巖性識別函數;b^bn···!^為第i類巖性識別函數的判別系數； 步驟4,對于兩種巖性的判別，將砂巖巖性識別函數減去泥巖巖性識別函數，得出 差值識別函數，該函數即為兩種巖性的綜合識別曲線，記為P(x);即： p⑴=A f = f砂-f泥=Δ bo+Δ bi · χι+···+Δ bn · Xn (4) 其中，P(x)為兩種巖性的綜合識別曲線;f砂為砂巖巖性識別函數，f泥為泥巖巖性識 另IJ函數;Χ1，χ 2···χΑ選取的η個測井參數變量；△ bo，Δ br·· Δ bn為兩種巖性綜合識別曲線的 判別系數； 把待判測井數據體的觀測值Χ=(χ(1)χ(2)···χ(η))代入上述差值識別函數計算出巖 性綜合識別曲線值Poo，若Poo 2 〇則為砂巖巖性，相反綜合曲線值Poo < 〇則為泥巖巖性； 步驟5,對于Μ類巖性的判別，把樣品的觀測值Χ=(χ (1)χ(2)···χ(η))代入式(3)中，分 別求得fi( x)，f2(x)"_fM(x)，取最大值記為fg( x)=max{fi(x)}，1 < i SM;由判別準則知，巖性識別 函數值大為相應巖性的概率高，則認為待判測井數據體XefgW對應的巖性所屬類型，同時 巖性綜合識別曲線Q(x)賦值為m;對Μ類不同巖性識別結果賦值，從而擬合出Μ類巖性的綜合 識別曲線Q(x);即：其中，Quo為Μ類巖性的綜合識別曲線;f1(x)，f2(x)…fM(x)為Μ類巖性識別函數;f g(x) 為代入待判數據體后Μ類巖性識別函數的最大值;?Μ,πινΠΗ為Μ類巖性綜合識別曲線賦值， nu的賦值區間為-Μ < mi < Μ。進一步的，所述步驟1中，η個測井參數包括自然電位SP、自然伽馬GR、聲波時差AC、 巖性密度DEN、電阻率RT。 本專利技術的有益效果是克服常規巖性識別方法的不足，提供一種識別高鹽富鉀區注 水開發油藏碎肩巖巖性的方法，即充分利用測井、巖心分析化驗等資料，以數理統計原理為 基礎，建立巖性識別函數，并擬合單井巖性識別綜合曲線，利用該曲線識別巖性效果較好， 自檢符合率高達90%。本專利技術能很好地考慮多種測井參數，進而擬合出一條準確的巖性綜 合識別曲線，計算結果合理且準確率高，實現對高鹽富鉀區注水開發油藏全井段碎肩巖巖 性的連續判別，可為類似地區的巖性識別提供借鑒。【附圖說明】 圖1是本專利技術實施例中研究區兩種巖性的綜合識別曲線。圖2是本專利技術實施例中研究區四種巖性的綜合識別曲線?！揪唧w實施方式】 為了使本專利技術的目的、技術方案及優點更加清楚明白，以下結合附圖及實施例，對 本專利技術進行進一步詳細說明。應當理解，此處所描述的具體實施例僅僅用以解釋本專利技術，并 不用于限定本專利技術。 -種識別高鹽富鉀區注水開發油藏碎肩巖巖性的方法，具體按照以下步驟進行： 步驟1，統計已知巖性的巖心分析化驗資料并將巖性分類。采用數理統計方法，剔 除非相關變量，選取對巖性反映比較敏感的η個測井參數（自然電位(SP)、自然伽馬(GR)、聲 波時差(AC)、巖性密度(DEN)、電阻率(RT)…）；步驟2，基于數理統計原理，建立巖性識別權重函數F(包括Fi、F2、Fr··)，計算出的 權重函數可以根據η個測井參數對巖性的貢獻率賦予不同的函數特征值，提取與各組有最 大可能多重相關的、函數特征值最大的兩個權重函數，即第一權重函數和第二權重函數： Fi = aoo+aoi · xi+---+a〇n · xn (1) F2 = aio+aii · xi+."+ain · Xn (2) 其中，Fi、F2為第一和第二權重函數；X1，x 2…xn*選取的n個測井參數變量；aoo， aoi…a〇n和aio，a2〇…ain分別對應兩個權重函數的判別系數。 步驟3,將給出的第一權重函數和第二權重函數用于計算各樣品在降維后空間坐 標系中的坐標值，然后計算出各類別在低維空間中的中心坐標，根據各樣品點距離各類別 重心的距離構造出每種巖性的識別函數，即： fi(X) =bi〇+bii · χι+...+bin · Xn (3) 其中，fi為第i類巖性識別函數;幻，^_心為選取的11個測井參數變量;1^，1^1···!^ 為第i類巖性識別函數的判別系數。 步驟4，對于兩種巖性的判別，將砂巖巖性識別函數減去泥巖巖性識別函數，得 出差值識別函數，該函數即為兩種巖性的綜合識別曲線，記為P〇〇;即： p⑴=A f = f砂-f泥=Δ bo+Δ bi · χι+···+Δ bn · Xn (4) 其中，P(x)為兩種巖性的綜合識別曲線;f砂為砂巖巖性識別函數，f泥為泥巖巖性識 另IJ函數;Χ1，χ 2···χΑ選取的η個測井參數變量；△ bo，Δ br·· Δ b本文檔來自技高網...
【技術保護點】
一種識別高鹽富鉀區注水開發油藏碎屑巖巖性的方法，其特征在于，按照以下步驟進行：步驟1，統計已知巖性的巖心分析化驗資料并將巖性分類；采用數理統計方法，剔除非相關變量，選取對巖性反映比較敏感的n個測井參數；步驟2，基于數理統計原理，建立巖性識別權重函數F，計算出的權重函數根據n個測井參數對巖性的貢獻率賦予不同的函數特征值，提取與各組有最大可能多重相關的、函數特征值最大的兩個權重函數，即第一權重函數和第二權重函數：F1＝a00+a01·x1+…+a0n·xn????????????????(1)F2＝a10+a11·x1+…+a1n·xn????????????????(2)其中，F1、F2為第一和第二權重函數；x1，x2…xn為選取的n個測井參數變量；a00，a01…a0n和a10，a20…a1n分別對應兩個權重函數的判別系數；步驟3，將給出的第一權重函數和第二權重函數用于計算各樣品在降維后空間坐標系中的坐標值，然后計算出各類別在低維空間中的中心坐標，根據各樣品點距離各類別重心的距離構造出每種巖性的識別函數，即：fi(X)＝bi0+bi1·x1+…+bin·xn????????????????(3)其中，fi為第i類巖性識別函數；bi0，bi1…bin為第i類巖性識別函數的判別系數；步驟4，對于兩種巖性的判別，將砂巖巖性識別函數減去泥巖巖性識別函數，得出差值識別函數，該函數即為兩種巖性的綜合識別曲線，記為P(X)；即：P(X)＝Δf＝f砂?f泥＝Δb0+Δb1·x1+…+Δbn·xn???????(4)其中，P(X)為兩種巖性的綜合識別曲線；f砂為砂巖巖性識別函數，f泥為泥巖巖性識別函數；x1，x2…xn為選取的n個測井參數變量；Δb0，Δb1…Δbn為兩種巖性綜合識別曲線的判別系數；把待判測井數據體的觀測值X＝(x(1)x(2)…x(n))代入上述差值識別函數計算出巖性綜合識別曲線值P(X)，若P(X)≥0則為砂巖巖性，相反綜合曲線值P(X)≤0則為泥巖巖性；步驟5，對于M類巖性的判別，把樣品的觀測值X＝(x(1)x(2)…x(n))代入式(3)中，分別求得f1(x),f2(x)…fM(x)，取最大值記為fg(x)＝max{fi(x)}，1≤i≤M；由判別準則知，巖性識別函數值大為相應巖性的概率高，則認為待判測井數據體X∈fg(x)對應的巖性所屬類型，同時巖性綜合識別曲線Q(X)賦值為m；對M類不同巖性識別結果賦值，從而擬合出M類巖性的綜合識別曲線Q(X)；即：Q(X)=m1····················fg(x)=f1(x)m2···················fg(x)=f2(x)m3····················fg(x)=f3(x)·....................................··.....................................··.....................................·mi·······&Cen...
【技術特征摘要】
【專利技術屬性】
技術研發人員：楊少春，牛海瑞，溫雅茹，
申請(專利權)人：中國石油大學華東，
類型：發明
國別省市：山東;37
還沒有人留言評論。發表了對其他瀏覽者有用的留言會獲得科技券。