本發(fā)明涉及礦產(chǎn)油氣資源勘查，特別涉及一種利用聚類方法識(shí)別多井沉積微相的方法。
背景技術(shù)：
：砂巖儲(chǔ)層是可以儲(chǔ)集和滲濾流體的巖層，即儲(chǔ)集層，能夠儲(chǔ)存和滲濾油氣的巖層。沉積微相對(duì)砂巖儲(chǔ)層尋找和評(píng)價(jià)有著重要指導(dǎo)作用，通過(guò)研究沉積微相空間配置關(guān)系確定砂巖儲(chǔ)層分布特征及其儲(chǔ)集性能評(píng)價(jià)是油氣勘探的主要內(nèi)容之一?，F(xiàn)有沉積微相研究方法一般是通過(guò)單井資料研究某點(diǎn)沉積特征，然后通過(guò)連井分析剖面沉積特征，最后確定沉積微相空間配置關(guān)系?，F(xiàn)有研究方法主要適用于井資料較少地區(qū)，借助研究者既有經(jīng)驗(yàn)，充分挖掘單井所提供的地質(zhì)信息。而對(duì)于勘探程度較高，井資料較多，上述方法工作量大，研究周期長(zhǎng)，不能滿足實(shí)際勘探需要。技術(shù)實(shí)現(xiàn)要素：本發(fā)明所要解決的技術(shù)問(wèn)題是提供一種利用聚類方法識(shí)別多井沉積微相的方法，如何快速判別多井沉積微相，通過(guò)利用數(shù)學(xué)統(tǒng)計(jì)的方法實(shí)現(xiàn)多井沉積微相的定量識(shí)別，從而提高沉積微相評(píng)價(jià)效率，為油氣勘探提供技術(shù)支撐。本發(fā)明解決上述技術(shù)問(wèn)題的技術(shù)方案如下：一種利用聚類方法識(shí)別多井沉積微相的方法，包括以下步驟:s1，收集各井位錄井地質(zhì)數(shù)據(jù)從而獲得巖性組合指數(shù)；s2，從步驟s1中獲得的數(shù)據(jù)中選取能反映不同沉積微相間巖性組合指數(shù)；s3，對(duì)步驟s1中各所述井位進(jìn)行分層，根據(jù)分層數(shù)據(jù)統(tǒng)計(jì)各層內(nèi)步驟s2中選取的所述巖性組合指數(shù)；s4，對(duì)步驟s3中得到的所述巖性指數(shù)的值域控制在存在最小值為0，最大值為1的范圍內(nèi)，對(duì)超過(guò)所述范圍的所述巖性指數(shù)進(jìn)行歸一化處理；s5，通過(guò)步驟s4所得巖性指數(shù)對(duì)多井樣本空間進(jìn)行聚類分析；s6，根據(jù)聚類分析結(jié)果，選擇不同的分類級(jí)別進(jìn)行沉積相組合和沉積微相識(shí)別；以及s7，利用單井沉積微相分析結(jié)果，根據(jù)沉積微相空間配置規(guī)律繪制沉積微相平面分布圖。本發(fā)明的有益效果是：本發(fā)明方法能夠充分利用現(xiàn)有巖性組合數(shù)據(jù)資料，統(tǒng)計(jì)運(yùn)算過(guò)程自動(dòng)化，實(shí)現(xiàn)沉積微相識(shí)別定量分析，在提高分辨率的同時(shí)大幅度減少工作量，縮短研究周期，研究成果可以快速被工業(yè)化應(yīng)用。在上述技術(shù)方案的基礎(chǔ)上，本發(fā)明還可以做如下改進(jìn)。進(jìn)一步，所述一種利用聚類方法識(shí)別多井沉積微相的方法，步驟s2中所述巖性組合指數(shù)包括巖性指數(shù)i1、砂巖粒級(jí)指數(shù)i2、砂巖粒級(jí)分異度i3、砂巖厚度指數(shù)i4、砂巖厚度分異度i5、泥巖顏色指數(shù)i6和泥巖厚度指數(shù)i7；公式(1)中h1為泥巖厚度，h2為泥質(zhì)砂巖厚度，h3為砂巖厚度，h0為地層厚度；公式(2)中h1為粉砂巖厚度，h2為細(xì)砂巖厚度，h3為中砂巖厚度，h4為粗砂巖厚度；公式(3)中α1為粉砂巖在砂巖中所占百分比，α2為細(xì)砂巖在砂巖中所占百分比，α3為中砂巖在砂巖中所占百分比，α4為中砂巖在砂巖中所占百分比；公式(4)中n0為砂巖的總層數(shù)，h1為粉砂巖厚度，h2為細(xì)砂巖厚度，h3為中砂巖厚度，h4為粗砂巖厚度；公式(5)中h1為粉砂巖厚度，h2為細(xì)砂巖厚度，h3為中砂巖厚度，h4為粗砂巖厚度，h0為砂巖平均厚度；公式(6)中h0為泥巖厚度，h1為深灰色泥巖厚度，h2為灰色泥巖厚度，h3為綠色泥巖厚度，h4為棕紅色泥巖厚度；公式(7)中n1為泥巖總層數(shù)；h0為泥巖總厚度。采用上述進(jìn)一步方案的有益效果是采用能夠用于識(shí)別沉積微相的巖性組合指數(shù)定量化表證巖性組合特征，把復(fù)雜的巖性組合特征用數(shù)字形式直觀地反映出來(lái)，通過(guò)定義各指數(shù)所代表的地質(zhì)意義，為進(jìn)一步分析沉積微相奠定基礎(chǔ)。進(jìn)一步，所述一種利用聚類方法識(shí)別多井沉積微相的方法，步驟s4中所述歸一化處理采用公式公式(8)中i＝(1,2,3,…,m)表示m口井的所需歸一化巖性組合指數(shù)，ii為第i口井所需歸一化巖性組合指數(shù)的值，imin為所需歸一化巖性組合指數(shù)的最小值，imax為所需歸一化巖性組合指數(shù)的最大值。采用上述進(jìn)一步方案的有益效果是把各巖性組合指數(shù)值域統(tǒng)一到0～1之間，在聚類分析計(jì)算樣品之間的距離時(shí)，各巖性組合指數(shù)之間互為平等的關(guān)系，所得到的距離不會(huì)因某巖性組合過(guò)大而掩蓋了其它參數(shù)。進(jìn)一步，所述一種利用聚類方法識(shí)別多井沉積微相的方法，步驟s5中所述聚類分析包括以下步驟：s51，將m口某研究層位所代表的沉積微相自劃分一類，共m個(gè)類，然后根據(jù)步驟s4所得巖性組合指數(shù)進(jìn)行聚類分析；s52，使用類平均法計(jì)算沉積微相之間的距離；s53，計(jì)算比較各沉積微相之間的距離，將距離最小的兩個(gè)沉積微相合并為一類；s54，重復(fù)步驟s53，，直到m個(gè)類歸為一個(gè)總類為止。采用上述進(jìn)一步方案的有益效果是在聚類過(guò)程中，選擇合適的距離對(duì)所有井進(jìn)行聚類分析，所選取的距離小，所得到的類越多。在實(shí)際聚類過(guò)程中，最初無(wú)法確定一個(gè)合適的距離，按照距離的遠(yuǎn)近，逐一進(jìn)行聚類分析，后根據(jù)分析結(jié)果，選擇合適的距離，確定合理分類。進(jìn)一步，所述一種利用聚類方法識(shí)別多井沉積微相的方法，步驟s51中，設(shè)有s種沉積微相，m口井，第g口井為ng，每口井某層段共提取p項(xiàng)巖性組合指數(shù)，用xgjk表示第g種沉積微相內(nèi)第j口井的第k項(xiàng)巖性組合指數(shù)；沉積微相1沉積微相2沉積微相ss、m、nm、p為正整數(shù)；計(jì)算需判別樣品屬于第g種沉積微相的概率p{g/y}，其中g(shù)＝1,2，...s；然后比較p{1/y}，p{2/y}，...，p{s/y}的大小；最后將該樣品歸于概率最大的沉積微相；qg為第g組的先驗(yàn)概率，fg(y1,y2,...yp)為樣品屬于沉積微相g的概率密度；所述先驗(yàn)概率采用樣本頻率作為估計(jì)值，式中，ng為第g口井，n為全部井位總數(shù)；求解樣品屬于沉積微相g的概率密度f(wàn)g(y1,y2,...yp)時(shí)設(shè)m個(gè)總體均服從正態(tài)分布，第g個(gè)總體的均值為μg向量，協(xié)方差矩陣為vg，可得概率密度f(wàn)g(y1,y2,...yp)，式(10)中x＝(x1,x2,…,xp)’是p維向量，μg＝(μg1，μg2，…，μgp)’是均值向量，是協(xié)方差矩陣，是協(xié)方差矩陣的逆矩陣；假設(shè)v1＝v2＝…＝vm＝v，此時(shí)在確定沉積微相的過(guò)程中，采用樣本均值來(lái)估計(jì)μg，用樣本的協(xié)方差矩陣s來(lái)估計(jì)v；其中g(shù)＝1,2，…，m；(12)式(12)中其中k＝1,2，…，p；式(13)中，其中sg是第g種沉積微相的離差矩陣，式(13)中n是全部井位，將式(12)、(13)代入式(11)中得出將式(14)代入式(9)得，對(duì)式(15)的分子取自然對(duì)數(shù)，得對(duì)式(16)中與g有關(guān)的項(xiàng)記作fg(y),得其中，g＝1,2,3，…，m；式(17)是判別函數(shù)的矩陣形式，利用矩陣運(yùn)算得式(18)中式中是s的逆矩陣s-1中第k行第t列元素。采用上述進(jìn)一步方案的有益效果是對(duì)聚類分析結(jié)果進(jìn)行檢驗(yàn)和新井沉積微相判別。通過(guò)聚類分析得到的分類結(jié)果，利用已知井位進(jìn)行校驗(yàn)所得模型，驗(yàn)證模型準(zhǔn)確性，對(duì)于新井沉積微相識(shí)別，可單獨(dú)采用此步驟進(jìn)行判別分析。進(jìn)一步，所述一種利用聚類方法識(shí)別多井沉積微相的方法，步驟s52中，沉積微相之間的距離dg(p,q)為式(19)中，d、i和j代表任意兩沉積微相之間的距離，gp和gq用來(lái)代指任意兩個(gè)沉積微相，l和m分別為gp和gq中的樣本數(shù)，d、i和j為gp沉積微相中第i個(gè)樣品與gq沉積微相中第j樣品之間的距離。采用上述進(jìn)一步方案的有益效果是通過(guò)定義并計(jì)算兩沉積微相間的距離，把各巖性組合指數(shù)從p維空間降到一維空間內(nèi)，在實(shí)際分類中，通過(guò)設(shè)定距離的大小，既可實(shí)現(xiàn)沉積微相聚類。進(jìn)一步，所述一種利用聚類方法識(shí)別多井沉積微相的方法，步驟s6中，根據(jù)聚類分析結(jié)果，劃分為不同的類別；巖性指數(shù)i1接近于1時(shí)砂巖所占比例大，砂巖粒級(jí)指數(shù)i2接近于1時(shí)粗砂巖含量高，砂巖粒級(jí)分異度i3接近于0.75時(shí)砂巖粒級(jí)越單一，砂巖厚度指數(shù)i4值越大平均單層砂巖厚度越小，砂巖厚度分異度i5值越大砂層之間厚度差異性越大，泥巖顏色指數(shù)i6接近于1時(shí)氧化色泥巖所占比例大，泥巖厚度指數(shù)i7值越大單層泥巖厚度越小。采用上述進(jìn)一步方案的有益效果是把各巖性組合指數(shù)賦予其地質(zhì)意義，判別各類所代表沉積微相。進(jìn)一步，所述一種利用聚類方法識(shí)別多井沉積微相的方法，步驟s7中，利用單井沉積微相分析結(jié)果，根據(jù)沉積微相空間配置規(guī)律，結(jié)合沉積微相平面配置關(guān)系繪制沉積微相平面分布圖。采用上述進(jìn)一步方案的有益效果是在平面圖上直觀表達(dá)沉積微相聚類分析結(jié)果。附圖說(shuō)明圖1為本發(fā)明方法步驟框圖；圖2為本發(fā)明聚類分析平均聯(lián)接樹狀圖；圖3為本發(fā)明某一研究層沉積微相平面分布圖。具體實(shí)施方式以下結(jié)合附圖對(duì)本發(fā)明的原理和特征進(jìn)行描述，所舉實(shí)例只用于解釋本發(fā)明，并非用于限定本發(fā)明的范圍。本發(fā)明提供一種利用聚類方法識(shí)別多井沉積微相的方法，包括以下步驟:s1，收集各井位錄井地質(zhì)數(shù)據(jù)從而獲得巖性組合指數(shù)；s2，從步驟s1中獲得的數(shù)據(jù)中選取能反映不同沉積微相間巖性組合指數(shù)；s3，對(duì)步驟s1中各所述井位進(jìn)行分層，根據(jù)分層數(shù)據(jù)統(tǒng)計(jì)各層內(nèi)步驟s2中選取的所述巖性組合指數(shù)；s4，對(duì)步驟s3中得到的所述巖性指數(shù)的值域控制在存在最小值為0，最大值為1的范圍內(nèi)，對(duì)超過(guò)所述范圍的所述巖性指數(shù)進(jìn)行歸一化處理；s5，通過(guò)步驟s4所得巖性指數(shù)對(duì)多井樣本空間進(jìn)行聚類分析；s6，根據(jù)聚類分析結(jié)果，選擇不同的分類級(jí)別進(jìn)行沉積相組合和沉積微相識(shí)別；以及s7，利用單井沉積微相分析結(jié)果，根據(jù)沉積微相空間配置規(guī)律繪制沉積微相平面分布圖。步驟s4中巖性組合指數(shù)的值域控制在存在最小值為0，最大值為1的范圍內(nèi)，具體是同一研究層對(duì)應(yīng)的各井的同一指數(shù)中的值域需要最小值為0，最大值為1。在步驟s1中,巖性錄井資料也就是各井位錄井地質(zhì)數(shù)據(jù)是確定巖性組合指數(shù)的主要數(shù)據(jù)來(lái)源,但巖性錄井資料受起下鉆、接單根、鉆井取心、鉆時(shí)、鉆井液循環(huán)流量等因素的影響，降低了錄井?dāng)?shù)據(jù)的分辨率，因此需要通過(guò)測(cè)井和取心資料對(duì)錄井資料進(jìn)行校正。沉積微相所對(duì)應(yīng)的縱向尺度應(yīng)控制在百米級(jí)范圍內(nèi)。在聚類分析之前，根據(jù)旋回關(guān)系，在三級(jí)、四級(jí)層序格架的基礎(chǔ)上，建立五級(jí)層序格架，并完成所有井位的分層工作。在本實(shí)施例中，主要描述某一研究層沉積微相分析及沉積微相平面圖的繪制過(guò)程。某一研究層沉積微相主要受環(huán)湖物源控制，主要發(fā)育三角洲平原沼澤、三角洲平原河道、三角洲前緣水下分流河道、河口壩、遠(yuǎn)砂壩、濱淺湖泥、半深湖泥、深湖泥等沉積微相，碳酸鹽相不發(fā)育。對(duì)此，步驟s2中選取巖性指數(shù)i1、砂巖粒級(jí)指數(shù)i2、砂巖粒級(jí)分異度i3、砂巖厚度指數(shù)i4、砂巖厚度分異度i5、泥巖顏色指數(shù)i6和泥巖厚度指數(shù)i7參數(shù)來(lái)描述研究層段錄井巖性組合指數(shù)。公式(1)中h1為泥巖厚度，h2為泥質(zhì)砂巖厚度，h3為砂巖厚度，h0為地層厚度。公式(2)中h1為粉砂巖厚度，h2為細(xì)砂巖厚度，h3為中砂巖厚度，h4為粗砂巖厚度。公式(3)中α1為粉砂巖在砂巖中所占百分比，α2為細(xì)砂巖在砂巖中所占百分比，α3為中砂巖在砂巖中所占百分比，α4為中砂巖在砂巖中所占百分比。公式(4)中n0為砂巖的總層數(shù)，h1為粉砂巖厚度，h2為細(xì)砂巖厚度，h3為中砂巖厚度，h4為粗砂巖厚度。其中砂巖厚度指數(shù)i4選取100米砂巖的層數(shù)，當(dāng)然也可以選取其他尺寸的砂巖層數(shù)。公式(5)中h1為粉砂巖厚度，h2為細(xì)砂巖厚度，h3為中砂巖厚度，h4為粗砂巖厚度，h0為砂巖平均厚度。公式(6)中h0為泥巖厚度，h1為深灰色泥巖厚度，h2為灰色泥巖厚度，h3為綠色泥巖厚度，h4為棕紅色泥巖厚度。公式(7)中n1為泥巖總層數(shù)；h0為泥巖總厚度。其中泥巖厚度指數(shù)i7選取100米泥巖的層數(shù)。步驟s4中所述歸一化處理采用公式公式(8)中i＝(1,2,3,…,m)表示m口井的所需歸一化巖性組合指數(shù)，ii為第i口井所需歸一化巖性組合指數(shù)的值，imin為所需歸一化巖性組合指數(shù)的最小值，imax為所需歸一化巖性組合指數(shù)的最大值。針對(duì)步驟s1中根據(jù)分層數(shù)據(jù)，統(tǒng)計(jì)各單井某一研究層內(nèi)的巖性組合指數(shù)，在此僅給出部分?jǐn)?shù)據(jù),如表1所示。具體統(tǒng)計(jì)方法可利用相關(guān)軟件程序進(jìn)行。表1序號(hào)表示研究層對(duì)應(yīng)的井號(hào),每口井有一個(gè)編號(hào),共78口井。。i3為砂巖粒級(jí)分異度最小值為0.39，大于0，砂巖厚度指數(shù)i4值域大于1，砂巖厚度分異度i5值域大于1，泥巖厚度指數(shù)i7值域大于1。在聚類分析之前需要對(duì)此四項(xiàng)巖性組合指數(shù)進(jìn)行歸一化處理。在本實(shí)施例中歸一化處理采用公式公式(8)中i＝(1,2,3,…,m)表示m口井的所需歸一化巖性組合指數(shù)，ii為第i口井所需歸一化巖性組合指數(shù)的值，imin為所需歸一化巖性組合指數(shù)的最小值，imax為所需歸一化巖性組合指數(shù)的最大值。通過(guò)歸一化處理，將歸一化后的數(shù)據(jù)進(jìn)行聚類分析。巖性指數(shù)i1、砂巖粒級(jí)指數(shù)i2和泥巖顏色指數(shù)i6值域在0～1之間，且存在最小值為0，最大值為1無(wú)需處理。通過(guò)歸一化處理，得到歸一化后的巖性組合指數(shù)如表3所示。表2歸一化處理是為了在聚類過(guò)程中，避免巖性組合指數(shù)單個(gè)數(shù)據(jù)過(guò)大或過(guò)小，對(duì)聚類結(jié)果產(chǎn)生影響，此數(shù)據(jù)僅用于聚類分析，在沉積微相判別過(guò)程中，使用未歸一化的數(shù)據(jù)。如表3所示。表3巖性組合指數(shù)值域代表意義巖性指數(shù)i10～1接近于1時(shí)砂巖所占比例大砂巖粒級(jí)指數(shù)i20～1接近于1時(shí)粗砂巖含量高砂巖粒級(jí)分異度i30～1接近于1時(shí)砂巖粒級(jí)越單一砂巖厚度指數(shù)i4>1值越大平均單層砂巖厚度越小砂巖厚度分異度i5>0值越大砂層之間厚度差異性越大泥巖顏色指數(shù)i60～1接近于1時(shí)氧化色泥巖所占比例大泥巖厚度指數(shù)i7>1值越大單層泥巖厚度越小把表2中歸一化的數(shù)據(jù)輸入到spss軟件內(nèi)，把每種巖性組合指數(shù)作為變量，在貝葉斯準(zhǔn)則下進(jìn)行系統(tǒng)聚類分析。步驟s5中所述聚類分析包括以下步驟：s51，將m口某研究層位所代表的沉積微相自劃分一類，共m個(gè)類，然后根據(jù)步驟s4所得巖性組合指數(shù)進(jìn)行聚類分析；s52，使用類平均法計(jì)算沉積微相之間的距離；s53，計(jì)算比較各沉積微相之間的距離，將距離最小的兩個(gè)沉積微相合并為一類；s54，重復(fù)步驟s53，每重復(fù)一次減少一個(gè)類，直到n個(gè)類歸為一個(gè)總類為止。在一實(shí)施例中，步驟s51中，設(shè)有s種沉積微相，m口井，第g口井為ng，每口井某層段共提取p項(xiàng)巖性組合指數(shù)，用xgjk表示第g種沉積微相內(nèi)第j口井的第k項(xiàng)巖性組合指數(shù)；沉積微相1沉積微相2沉積微相ss、m、nm、p為正整數(shù)。計(jì)算需判別樣品屬于第g種沉積微相的概率p{g/y}，其中g(shù)＝1,2，...s；然后比較p{1/y}，p{2/y}，...，p{s/y}的大小；最后將該樣品歸于概率最大的沉積微相；qg為第g組的先驗(yàn)概率，fg(y1,y2,...yp)為樣品屬于沉積微相g的概率密度，是關(guān)于y1,y2,...yp的函數(shù)，y1代表礦物成分1屬于g沉積微相的概率密度，y2代表礦物成分2屬于g沉積微相的概率密度，yn代表礦物成分n屬于g沉積微相的概率密度。先驗(yàn)概率采用樣本頻率作為估計(jì)值，式中，ng為第g口井，n為全部井位總數(shù)。求解樣品屬于沉積微相g的概率密度f(wàn)g(y1,y2,...yp)時(shí)設(shè)m個(gè)總體均服從正態(tài)分布，第g個(gè)總體的均值為μg向量，協(xié)方差矩陣為vg，可得概率密度f(wàn)g(y1,y2,...yp)。式(10)中x＝(x1,x2,…,xp)’是p維向量，p是巖性組合指數(shù)的項(xiàng)數(shù)，μg＝(μg1，μg2，…，μgp)’是均值向量，μg1代表g個(gè)總體的第1種輕重礦物成分均值，μg2代表g個(gè)總體的第2種輕重礦物成分均值，μgp代表g個(gè)總體的第p種輕重礦物成分均值。是協(xié)方差矩陣，是協(xié)方差矩陣的逆矩陣。σ(g)用來(lái)表示兩種礦物成分共同產(chǎn)生的方差，用來(lái)衡量?jī)蓚€(gè)中礦物成分的相互影響。例如用來(lái)表示第g個(gè)總體中第1種礦物成分與第2中礦物成分之間的方差，用來(lái)表示第g個(gè)總體中第2種礦物成分與第1中礦物成分之間的方差。假設(shè)v1＝v2＝…＝vm＝v，此時(shí)在確定沉積微相的過(guò)程中，采用樣本均值來(lái)估計(jì)μg，用樣本的協(xié)方差矩陣s來(lái)估計(jì)v。其中g(shù)＝1,2，…，m；(12)式(12)中其中k＝1,2，…，p。xg·p是樣本空間m中屬于某種沉積微相的所有井的巖性組合指數(shù)平均值向量。式(13)中，其中sg是第g種沉積微相的離差矩陣。式(13)中n是全部井位，ng為第g口井。m為井?dāng)?shù)。將式(12)、(13)代入式(11)中得出將式(14)代入式(9)得，對(duì)式(15)的分子取自然對(duì)數(shù)，得對(duì)式(16)中與g有關(guān)的項(xiàng)記作fg(y),得其中，g＝1,2,3，…，m。m是井?dāng)?shù)。fg(y)是y1,y2,...yp的線性函數(shù)，共有m個(gè)；顯然，當(dāng)fg(y)在某沉積微相內(nèi)取得最大值時(shí)，p{g/y}也屬于該沉積微相的概率也達(dá)到最大，樣品y就應(yīng)歸于此沉積微相，因此式(17)是判別函數(shù)的矩陣形式，利用矩陣運(yùn)算得，式(18)中式中是s的逆矩陣s-1中第k行第t列元素。qg為第g組的先驗(yàn)概率。ckg、c0g用來(lái)簡(jiǎn)化物源判別函數(shù)，無(wú)實(shí)際地質(zhì)意義。在一實(shí)施例中，步驟s52中，沉積微相之間的距離dg(p,q)為式(19)中，d、i和j代表任意兩沉積微相之間的距離，gp和gq用來(lái)代指任意兩個(gè)沉積微相，l和m分別為gp和gq中的樣本數(shù)，d、i和j為gp沉積微相中第i個(gè)樣品與gq沉積微相中第j樣品之間的距離。步驟s6中，根據(jù)聚類分析結(jié)果，劃分為不同的類別；巖性指數(shù)i1接近于1時(shí)砂巖所占比例大，砂巖粒級(jí)指數(shù)i2接近于1時(shí)粗砂巖含量高，砂巖粒級(jí)分異度i3接近于0.75時(shí)砂巖粒級(jí)越單一，砂巖厚度指數(shù)i4值越大平均單層砂巖厚度越小，砂巖厚度分異度i5值越大砂層之間厚度差異性越大，泥巖顏色指數(shù)i6接近于1時(shí)氧化色泥巖所占比例大，泥巖厚度指數(shù)i7值越大單層泥巖厚度越小。分析各類中巖性組合指數(shù)所代表的意義，綜合判定各類別沉積微相，進(jìn)而得到各單井研究層段的沉積微相。具體的，聚類過(guò)程中先將所有的樣本都自成一類，然后計(jì)算哪兩者的距離較近，最近的先聚，即先成一類，直到所有樣品都規(guī)為一類位置，最終得到q型樣本聚類樹狀圖如圖2所示，最左側(cè)一列數(shù)字代表井號(hào)，井號(hào)14與井號(hào)58之間的距離為0.9，首先劃歸為一類，然后依次將井號(hào)14～井號(hào)9歸為1類；井號(hào)4在距離為1.4時(shí)，可與井號(hào)14～井號(hào)9歸為一類。通過(guò)分析聚類樹狀圖，最終選擇距離10和距離15將78井共分為7類，分別標(biāo)記為①～⑦。①巖性指數(shù)為0，砂巖粒級(jí)分異度在0.56～1之間，泥巖顏色指數(shù)在0～0.25之間，泥巖厚度指數(shù)在1.4～12.7之間，表現(xiàn)為穩(wěn)定還原條件下泥巖沉積，為半深湖沉積微相；②巖性指數(shù)在0～0.04之間，砂巖粒級(jí)指數(shù)在0～0.5之間，砂巖粒級(jí)分異度在0.51～0.87之間，砂巖厚度指數(shù)在0～200之間，砂巖厚度分異度在0～2之間，泥巖顏色指數(shù)在0～0.26之間，泥巖厚度指數(shù)在7.69～57.69之間，為平原沼澤沉積微相；③巖性指數(shù)在0.08～0.77之間，砂巖粒級(jí)指數(shù)在0.25～0.69之間，砂巖粒級(jí)分異度在0.5～0.87之間，砂巖厚度指數(shù)在11.76～100之間，砂巖厚度分異度在2～20.96之間，泥巖顏色指數(shù)在0.29～1之間，泥巖厚度指數(shù)在6.94～44.83之間，為平原河道沉積微相；④巖性指數(shù)在0.05～0.24之間，砂巖粒級(jí)指數(shù)在0.25～0.75之間，砂巖粒級(jí)分異度在0.74～0.1之間，砂巖厚度指數(shù)在26.09～80之間，砂巖厚度分異度在2.6～8.9之間，泥巖顏色指數(shù)在0～0.36之間，泥巖厚度指數(shù)在7.79～22.22之間，為前三角洲沉積微相；⑤巖性指數(shù)在0.09～0.53之間，砂巖粒級(jí)指數(shù)在0.25～0.75之間，砂巖粒級(jí)分異度在0.87～1.0之間，砂巖厚度指數(shù)在14.29～85.71之間，砂巖厚度分異度在2.08～14之間，泥巖顏色指數(shù)在0～0.25之間，泥巖厚度指數(shù)在20.37～53.66之間，為遠(yuǎn)砂壩沉積微相；⑥巖性指數(shù)在0.33～0.83之間，砂巖粒級(jí)指數(shù)在0.25～0.75之間，砂巖粒級(jí)分異度在0.87～1之間，砂巖厚度指數(shù)在10.71～20.51之間，砂巖厚度分異度在11.11～22.81之間，泥巖顏色指數(shù)在0～0.68之間，泥巖厚度指數(shù)在6.67～47.62之間，為水下分流河道沉積微相；⑦巖性指數(shù)在0.15～0.74之間，砂巖粒級(jí)指數(shù)在0.25～0.75之間，砂巖粒級(jí)分異度在0.39～0.7之間，砂巖厚度指數(shù)在17.95～65.22之間，砂巖厚度分異度在3.7～11.2之間，泥巖顏色指數(shù)在0～0.21之間，泥巖厚度指數(shù)在15.56～54.55之間，為河口壩沉積微相。步驟s7中，利用單井沉積微相分析結(jié)果，根據(jù)沉積微相空間配置規(guī)律，結(jié)合沉積微相平面配置關(guān)系繪制沉積微相平面分布圖。在確定78口單井研究層沉積微相之后，把井位坐標(biāo)投到平面圖中，結(jié)合沉積微相平面配置關(guān)系，繪制沉積微相平面分布圖如圖3所示。以上所述僅為本發(fā)明的較佳實(shí)施例，并不用以限制本發(fā)明，凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。當(dāng)前第1頁(yè)12