[發明專利]一種全井段巖性識別方法

有效

• 申請號: 201410037682.3
• 申請日: 2014-01-26
• 公開（公告）號: CN103775075A
• 公開（公告）日: 2014-05-07
• 發明（設計）人: 崔維平;楊玉卿;于本志;朱衛星
• 申請（專利權）人: 中國海洋石油總公司;中海油田服務股份有限公司
• 主分類號: E21B49/00
• 分類號: E21B49/00
• 代理公司: 北京安信方達知識產權代理有限公司 11262
• 代理人: 惠磊;曲鵬

技術領域

本發明涉及一種巖性識別方法，特別是涉及一種采用交互式多維直方圖對巖性進行全井段巖性類型定量識別的方法。

背景技術

目前利用測井信息識別巖性的方法主要有常規測井資料解釋法、判別分析法、神經網絡算法等。常規測井資料解釋法目的主要是確定地層泥質含量，解釋結果僅為砂巖和泥巖兩種巖性,解釋結果太粗,無法反映地質非均質性情況；利用判別分析法識別地層巖性較為普遍，但對原始樣本數據的代表性和影響巖性類別的指標的全面性要求較高，針對性較差，由于受到流體性質等因素的影響，可靠性也較差，無法滿足儲層物性的精細評價和沉積學研究的需要；神經網絡算法，該方法泛化能力不強，識別結果準確性較差，且需要人工先驗知識，在巖心資料較少的情況下，很難用常規方法進行巖性識別。上述三種方法巖性識別精度不夠，均無法滿足測井精細評價和地質研究的需要。

發明內容

本發明提供一種全井段巖性識別方法，其目的是提供一種利用巖心資料和測井資料通過數理統計學建立多維直方圖模型進行巖性的分類和定量識別，最終建立精細連續的巖性剖面。

為了實現上述發明目的，本發明采取的技術方案如下：

一種全井段巖性識別方法，包括：

從多個深度位置取出多個巖心樣本，確定每一巖心樣本的粒度數據及其巖性類型；

確定所述多個深度位置處的多種測井參數的值，選擇隨著巖性樣本粒度增大而增大或減小的N種測井參數作為特征參數，N≥2；

根據所述多個深度位置處巖心樣本點的巖性類型及所述N種特征參數的值，確定以所述特征參數坐標的N維空間中，各N維單元格與巖性類型之間的對應關系；

根據全井段各深度位置處所述特征參數的值及所述對應關系，確定各深度位置處的代表巖性類型。

進一步地，在所述N維空間中確定所述特征參數與巖性類型之間的對應關系的方法為：

對所述多個深度位置中的每一深度位置，將該深度位置處的巖性類型確定為該深度位置處N種特征參數值所指向的N維單元格對應的一份巖性類型，由此得到所述N維空間中，一個或多個N維單元格對應的一種或多種巖性類型及每一種巖性類型的份數；

對所述N維空間中沒有對應的巖性類型的N維單元格，通過就近插值法插值生成對應的一種或多種巖性類型及其份數。

進一步地，各深度位置處的代表巖性類型為此深度位置所有巖性類型中份數最多的巖性類型。

進一步地，所述巖心樣本點粒度數據容量不少于五個。

進一步地，所述測井參數為自然伽馬參數、密度參數、聲波參數和中子參數。

另一種可替代的技術方案為：

一種交互式多維直方圖全井段巖性識別方法，包括：

步驟一，從多個深度位置取出多個巖心樣本，確定每一巖心樣本的粒度數據，并根據巖心樣本粒度數據確定巖性類型；

步驟二，用同一深度位置處的巖性類型刻度多種測井參數的值，并建立測井參數與巖性類型之間的交會圖；

步驟三，選擇交會圖中，隨著巖性樣本粒度增大而增大或減小的測井參數作為特征參數，分別統計每一個特征參數中各類巖性類型出現的份數，并建立多個特征參數-巖性類型-份數之間的直方圖模型；

步驟四，將步驟三中得到所有直方圖模型統一表示在一個N維坐標空間中，N≥2；

步驟五，對于N維坐標空間的單元格中沒有巖性類型的，查找距擬生成單個格空間距離最近的樣本點單元格，采用就近插值法插值生成無樣本點單元格中的各類巖性類型，并更新N維坐標空間，最終在N維坐標空間形成包括全井段巖性類型分布圖；

步驟六，分別統計N維坐標空間每個單元格中各類巖性類型的數目，并計算各類巖性類型在每個單元格內出現的份數；

步驟七，將全井段的每一個深度區段所對應的特征參數值逐個在N維坐標空間中單元格尋找對應位置，并以此單元格中出現份數最多的巖性作為對應深度區段的巖性類型代表，然后對連續且巖性類型相同的特征參數值合并，形成全井段的連續巖性剖面。

優選地，所述步驟一中，巖心樣本點粒度數據容量不少于五個。

優選地，所述步驟二中，測井參數為自然伽馬參數、密度參數、聲波參數或中子參數。

優選地，所述步驟三，直方圖中每個特征參數軸根據其所代表的特征參數值范圍等間距劃分，形成空間單元格。

優選地，所述步驟四，N維坐標空間中每個維度間距和數值范圍的確定同其所代表特征參數在直方圖中范圍和間距。

本發明和現有技術相比，具有如下有益效果：