1、 測井資料處理與解釋測井資料處理與綜合解釋實驗專業名稱：勘查技術與工程 學生姓名：張奎濤學生學號：201326020109指導老師：岳崇旺、王飛完成日期：2017-1-2目 錄實驗一 定性劃分儲集層并定量解釋.1實驗二 利用綜合方法估計地層泥質含量.7實驗三 含泥質復雜巖性地層綜合測井處理.18實驗一 定性劃分儲集層并定量解釋一、 實驗目的：通過對測井曲線特征的分析和認識，掌握定性劃分砂泥巖剖面儲集層的基本方法，并應用阿爾奇公式，進行儲層參數的計算，鞏固已經學過的鉆井地球物理課程的主要內容與應用。二、 實驗要求正確劃分出儲集層和非儲集層，對砂泥巖剖面能區分開較明顯的油水層。進行測井曲線讀數，簡

2、單地計算出孔隙度、飽和度等參數。三、 實驗場地、用具與設備測井實驗室或一般的教室，長直尺、鉛筆、像皮和計算器四、 實驗內容：1 測井曲線圖的認識；圖1是某井的綜合測井曲線圖。圖中共有5道，第一道主要為反映巖性的測井曲線道，包括：自然電位測井曲線曲線符號為SP、記錄單位mv；自然伽馬測井曲線曲線符號為GR、記錄單位API；井徑測井曲線曲線符號為CAL，記錄單位in或cm；巖性密度測井曲線（光電吸收界面指數）曲線符號為PE；第二道是深度道；通常的深度比例尺為1：200 或1：500第三道是反映含油性的測井曲線道，包括深中淺三條電阻率測井曲線，分別是：深側向測井曲線曲線符號為LLD、記錄單位m；淺側

3、向測井曲線曲線符號為LLS、記錄單位m；微球形聚焦測井曲線曲線符號為MSFL、記錄單位m；電阻率測井曲線通常為對數刻度。第四道為反映孔隙度的測井曲線道，包括：密度測井曲線曲線符號為DEN或RHOB，記錄單位g/cm3；中子測井曲線曲線符號為CNL或PHIN，記錄單位%，有時為v/v。聲波測井曲線曲線符號為AC或DT，記錄單位us/ft，有時為us/m。中子和密度測井曲線的刻度的特點是保證在含水砂巖層上兩條曲線重迭，在含氣層上，密度孔隙度大于中子孔隙度，在泥巖層上，中子孔隙度大于密度孔隙度；第五道是反映粘土礦物類型的測井曲線道，包括自然伽馬能譜測井中的三條曲線：放射性釷測井曲線曲線符號為Th或T

4、HOR，記錄單位是ppm；放射性鈾測井曲線曲線符號為U或URAN，記錄單位ppm；放射性鉀測井曲線曲線符號為K或POTA，記錄單位%，有時為v/v。2 測井曲線特征（1）砂泥巖剖面的測井曲線特征砂泥巖剖面儲集層（砂巖）的典型特征是，一般自然電位有明顯的異常，異常的方向和幅度取決于泥漿濾液電阻率（Rmf）和地層水的電阻率（Rw），或者說與Rmf與Rw的比值有關,如果Rmf Rw，則為負異常，否則為正異常。如果砂層中不含放射性礦物，自然伽瑪曲線亦顯示低值。微電極曲線一般在砂巖層幅值高，并出現正幅差。而泥巖的幅度和幅差均較低，當井眼條件不好時，可能會出現曲線跳動現象。砂巖中含灰質較多的夾層，因為致密

5、電阻率異常高，幅度差很小或沒有。一般幅度差的大小標明了儲集層滲透性的好壞。普通電阻率測井曲線在泥巖處顯示為低值。砂巖處顯示為高值，含油砂巖幅值就更高，如有兩條探測深度不同的Ra 曲線，幅值的差別顯示著低侵、高侵。通常在油層上為低侵，水層上為高侵。井徑在泥巖層擴大，砂巖層縮小（略小于鉆頭直徑）。具體特征總結見表1表1 砂泥巖剖面測井曲線特征 地層測井曲線儲集層砂巖非儲集層泥巖自然電位負異常（RwRmf）泥巖基線自然伽馬低高井徑縮徑擴徑深中淺電阻率高阻低阻聲波300us/m釷低高鈾低高鉀低高（2）碳酸鹽巖剖面的測井曲線特征碳酸鹽巖剖面的測井解釋任務，就是從致密的圍巖中找出孔隙性、裂縫性的儲集層，并

6、判斷其含油性。碳酸鹽巖剖面電阻率一般較高，自然電位效果不好。為了區分巖性和劃分儲層，一般使用自然伽馬測井曲線。 儲集層相對于致密的圍巖具有低阻、低自然伽馬以及孔隙度測井反映孔隙度較大的特點。3 劃分儲集層的基本方法與原則基本要求：凡一切可能含油氣的地層都要劃出來，要適當地劃分明顯的水層。具體要求為：（1） 估計為油層、氣層、油水同層和含油水層的儲集層都必須分層解釋。（2） 厚度半米以上的電性（測井曲線）可疑層（即指從測井曲線上看有油氣的地層）或錄井顯示為微含油級別以上的儲集層必須分出。（3） 選擇出作為確定地層水電阻率Rw 的標準水層（厚度大、巖性純、不含油）要劃分出來。（4） 錄井、氣測有大

7、段油氣顯示而測井曲線顯示不好的儲集層，應選取一定層位，尤其是該組儲層的頂部層位，進行分層。（5） 當有多套油水系統，油層組包括若干水層時，只解釋最靠近油層的水層。（6） 對于新區探井，應做細致工作，對各個儲層均應酌情選層解釋，以使不漏掉可能有油氣的地層。4 正確劃分出儲集層的方法（1）砂泥巖剖面通常是自然電位（SP）曲線的異常確定滲透層的位置，用微電極曲線確定分層界面，分層前，應將井場收集的井壁取芯、氣測顯示等有關油氣顯示的資料標注在綜合測井曲線圖上，并根據鄰井的測井和試油等資料對本井的油水關系作出初步估計。分層時應注意：l 確定分層的界面深度時，應左右環顧，照顧到分層線對每條測井曲線的合理性

8、。l 分層的深度誤差不應大于0.1m。l 滲透層中，凡是0.5m以上的非滲透性夾層（泥巖或致密層），應將夾層上下的滲透層分兩層解釋。l 巖性漸變層頂界（頂部漸變層）或底界（底部漸變層）分層深度應在巖性漸變結束處。l 一個厚度較大的滲透層，如有兩個以上解釋結論，應按解釋結論分層。l 在同一解釋井段，如果油氣層與水層巖性、地層結構和孔隙度基本相同，則油氣層是純水層的電阻率的3-5倍。純水層的自然電位異常最大，油氣層異常明顯偏小，油水同層介于油、水層之間。并且厚度較大的油水同層，自上而下電阻率有明顯減小的趨勢。（2）碳酸鹽巖剖面碳酸巖鹽剖面劃分滲透性地層的5 測井曲線讀數分層以后，要從有關的主要測井

9、曲線將代表該儲層的測井曲線讀數，以便計算孔隙度、飽和度等地質參數，在厚度較大的儲集層中按測井曲線變化確定幾個取值區，對每個取值區對應讀數計算，幾種主要測井曲線取值區的最小厚度如下：各種孔隙度測井0.6m。側向測井0.6m感應測井，低阻0.6m，高阻層1.5m。每種測井曲線分層和取值要符合其方法特點，例如聲波測井扣除致密夾層，選用與滲透層相對應部分的平均值。電阻率測井曲線則扣除致密夾層，選用與滲透層相對應部分的極大值的平均值。另外注意孔隙度與電阻率測井曲線對應取值的原則。因為要用兩者結合計算地層的含水飽和度，兩者當然應該是對應深度上同一地層或同一取值區的讀數。巖層含油性的定性判斷，主要依據井曲線

10、的測井曲線特征，而電性特征是巖石物性、巖性和含油性的綜合反映。因此在判斷地層的含油性時，一般應將測量井段首先按照地層水礦化度的不同分為不同的解釋井段，然后才有可能對每一個解釋井段在充分考慮其巖性特點的前提下進行含油性解釋。由于地下地層復雜性，儀器的局限性，上述原則是一般性的。要做到正確地解釋，一方面應多收集資料，認真分析曲線，另一方面還要了解區域性特點和規律，要積累經驗。6 計算出孔隙度、飽和度等參數。讀數以后，還要做一些定量計算，常用的公式：孔隙度：含水飽和度：上式中為當前層的聲波時差，為地層水的聲波時差,189us/ft(623us/m)，為固體骨架的聲波時差，對于砂巖骨架，主要礦物為石英

11、，其聲波時差為55.5us/ft。a是常數，對于砂巖地層通常取1.0，為當前層的電阻率，m為膠結指數。五、 實驗結果處理結果如圖2所示 具體分層如下表所示：儲層序號頂部深度（m）底部深度（m）厚度（m）測井曲線讀數孔隙度含水飽和度SPGRDEN11122.4114219.6-6.49100.012.31413.84586.4421154.2116216.8-8.023103.5732.37311.61586.33231218.0512267.95-10.504115.0922.28913.94883.91441236.21252.616.4-14.695104.6862.33013.22978

12、.80751272.41286.1513.75-17.636106.9982.34813.90163.182由自然電位的劃分原則可知，在圖2的自然電位的曲線上，有五個明顯的偏離泥基線的負異常，故可以劃分出五個儲集層。然后就可以求出泥質含量，進而求出泥質校正后的孔隙度，然后根據阿爾奇公式，就可求出含水飽和度、滲透率，進而劃分油氣水層。 圖1 某井的綜合測井曲線圖23451圖2 SN183井測井解釋綜合圖實驗二利用綜合方法估計地層泥質含量一、 實驗目的：通過實際計算，鞏固掌握利用多種測井資料確定泥質含量的方法。二、 實驗要求自編程序，在計算機上運算出地層泥質含量。三、 實驗場地、用具與設備計算中心

13、，尺子、像皮和計算機；四、 實驗內容： 什么是泥質含量：泥質是指顆粒直徑小于0.01mm的碎屑物質，泥質含量，也叫做泥質體積，是指泥質的體積占巖石總體積的比：確定Vsh的重要性泥質含量的確定，在泥質砂巖儲集層的定量解釋中具有重要意義。多年來人們提出許多計算泥質含量的理論和方法。目前求取泥質含量的方法大致可分為兩類，一類是用每種測點各求出一個泥質含量，然后求出最佳值。當巖石含有泥質時，各種測井曲線均或多或少地受到泥質的影響，其影響的程度受Vsh的決定，評價巖石的特性時，只有已知Vsh，才知道由于泥質帶來的影響，從而將泥質的影響校正掉。一般而言，用自然伽馬或自然伽馬能譜或自然電位來求取泥質含量效果

14、最好，但自然伽馬要求儲層中除了泥質外，其他物質不含放射性礦物。自然電位要求地層水電阻率保持不變，且儲層中的泥質與相鄰泥巖的的成分相同。用其他方法計算泥質含量則要求更為苛刻的條件：如電阻率方法要求儲層的孔隙度和含水飽和度均要很小。中子和聲波方法則要求孔隙度很小。3. 確定Vsh的方法：（） 自然伽瑪法式中，分別是砂巖和泥巖層的自然伽馬值，GCUR是與地層有關的經驗系數，新地層（第三系地層）GCUR=3.7，老地層GCUR=2.0.（） 自然電位法式中，是當前層的自然電位讀數，和分別是純地層和泥質地層的自然電位讀數（）電阻率（b=1.5）（）中子法式中，是當前層的視中子孔隙度讀數，是泥巖層的視中子

15、孔隙度讀數。（） 交會圖法以中子密度測井交會圖為例，通過對圖2所示的石英點（Q）、水點（W）和泥巖點（SH）構成的三角形進行分解，依據資料點所落入三角形中的位置，可以推測出來泥質含量。或者利用下式進行計算（依據點到直線的距離計算方法）：式中，SHWQ0是石英點（Q）和水點（W）連線的直線方程。依據任意兩點的直線，用石英點()和水點()兩個點的參數可以推出：，當然，也可以用中子聲波、聲波密度交會圖的類似方法求Vsh 。如果解釋層段上沒有純泥巖層時，上述交會圖法所定出的泥巖點位置并不代表實際的泥巖參數，導致用交會圖方法估計的泥質含量比實際值偏大。上述幾種方法計算的泥質含量往往都有一定的條件，當條件

16、滿足時， 泥質含量的近似結果，當條件不滿足時，計算的泥質含量均可能偏高，所以在實際處理時，最后選取其中的最小值作為接近實際的泥質含量。五、 實驗結果1.已知條件：GRmin=50 GRmax=110 SPmin=0 SPmax=100 b=1.5 Rsh=5 2.測井數據depthspmsflllslldrhobnphigrdtcal2680.9888.504.809.0910.842.5517.00103.1773.1311.992682.0983.402.026.588.212.0720.00108.2577.7514.192683.0170.001.585.717.171.5924.00

17、107.8181.8115.592684.0280.001.925.967.481.9623.00108.7581.2114.792685.0484.302.986.818.282.3020.00105.0878.3812.032686.0578.509.607.158.222.5919.00112.1977.1310.592687.0736.0014.1810.9311.732.6417.0098.4472.409.822688.0926.5019.3818.1319.142.5817.0085.5470.739.492689.0025.2068.9450.0651.542.5011.005

18、3.4865.949.262690.0227.0055.8145.6645.832.5210.0064.0665.319.163.計算結果spmsflllslldrhob-nphinphigr泥質含量0.885001.027590.671330.596980.542990.566670.805350.542990.834001.829840.832710.718490.463950.666670.947160.463950.700002.155470.915280.786381.433950.800000.934210.700000.800001.892830.889510.764500.60

19、9270.766670.962040.609270.843001.412000.813860.714430.071650.666670.856730.071650.785000.647340.787850.717900.710030.633331.069200.633330.360000.499110.593700.566390.752580.566670.687470.360000.265000.405270.423690.408650.612850.566670.424370.265000.252000.173910.215270.211130.204880.366670.027910.0

20、27910.270000.200220.228890.228320.214510.333331.279430.20022附錄：程序源代碼! 功能:利用測井數據求泥質含量! 時間:2016/12/22!-! 參數說明:! cmdfile:參數文件! file_data:輸入測井數據文件! file_contant:常數參量文件! file_sp:輸出自然電位估計的泥質含量文件! file_msfl:輸出微球型聚焦估計的泥質含量文件! file_lls:輸出淺側向估計的泥質含量文件! file_lld:輸出深側向估計的泥質含量文件! file_nphi:輸出補償中子估計的泥質含量文件! file_

21、gr:輸出自然伽馬估計的泥質含量文件! file_n_d:輸出中子-密度交會圖估計的泥質含量文件! n:測井數據的個數! grmin,grmax:自然伽馬的最小值與最大值! spmin,spmax:自然電位的最小值與最大值! b,r_sh:已知常數! fai_nsh:已知常數! p_sh,p_ma,p_f:已知常數! fai_nf,fai_nma:已知常數! sp(:):自然電位! msfl(:):微球型聚焦! lls(:):淺側向! lld(:):深側向! nphi(:):補償中子! gr(:):自然伽馬! rhob(:):巖性密度!-program cjsyimplicit nonecha

22、racter*80 cmdfilecharacter*80 file_data,file_contantcharacter*80 file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_dinteger nreal grmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nmareal,allocatable:sp(:),msfl(:),lls(:),lld(:),nphi(:),gr(:),rhob(:) cmdfile=cmd.txtcall read_

23、cmd(cmdfile,file_data,file_contant,file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_d)call read_number(file_data,n)allocate(sp(1:n),msfl(1:n),lls(1:n),lld(1:n),nphi(1:n),gr(1:n),rhob(1:n)call read_data(file_data,n,sp,msfl,lls,lld,nphi,gr,rhob)call read_contant(file_contant,grmin,grmax,sp

24、min,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma)call get_sp(file_sp,n,sp,spmin,spmax)call get_r(file_msfl,n,msfl,b,r_sh)call get_r(file_lls,n,lls,b,r_sh)call get_r(file_lld,n,lld,b,r_sh)call get_nphi(file_nphi,n,nphi,fai_nsh)call get_gr(file_gr,n,gr,grmin,grmax)call get_n_d(file_n_d,n,nphi,rho

25、b,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma)deallocate(sp,msfl,lls,lld,nphi,gr,rhob)end!-! ! 功能：讀取參數文件! 輸入參數說明：! cmdfile:參數文件! 輸出參數說明：! file_data:輸入測井數據文件! file_contant:常數參量文件! file_sp:輸出自然電位估計的泥質含量文件! file_msfl:輸出微球型聚焦估計的泥質含量文件! file_lls:輸出淺側向估計的泥質含量文件! file_lld:輸出深側向估計的泥質含量文件! file_nphi:輸出補償中子估計的泥質含量文件

26、! file_gr:輸出自然伽馬估計的泥質含量文件! file_n_d:輸出中子-密度交會圖估計的泥質含量文件!-subroutine read_cmd(cmdfile,file_data,file_contant,file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_d)implicit nonecharacter*(*) cmdfilecharacter*(*) file_data,file_contantcharacter*(*) file_sp,file_msfl,file_lls,file_lld,file_nphi

27、,file_gr,file_n_dopen(10,file=cmdfile,status=old)read(10,*) file_dataread(10,*) file_contantread(10,*) file_spread(10,*) file_msflread(10,*) file_llsread(10,*) file_lldread(10,*) file_nphiread(10,*) file_grread(10,*) file_n_dclose(10)endsubroutine read_cmd!-! ! 功能：讀取測井數據的個數! 輸入參數說明：! filename:輸入測井數據

28、文件! 輸出參數說明：! number:測井數據的個數!-subroutine read_number(filename,number)implicit nonecharacter*(*) filenameinteger numberopen(11,file=filename,status=old)number=0 do while(.not.eof(11) read (11,*,end=100,ERR=100) number=number+1100 end doclose(11)number=number-1endsubroutine read_number!-! ! 功能：讀取測井數據!

29、輸入參數說明：! file_data:輸入測井數據文件! n:測井數據的個數! 輸出參數說明：! sp(:):自然電位! msfl(:):微球型聚焦! lls(:):淺側向! lld(:):深側向! nphi(:):補償中子! gr(:):自然伽馬! rhob(:):巖性密度!-subroutine read_data(file_data,n,sp,msfl,lls,lld,nphi,gr,rhob)implicit nonecharacter*(*) file_datainteger nreal sp(1:n),msfl(1:n),lls(1:n),lld(1:n),nphi(1:n),gr

30、(1:n),rhob(1:n)integer ireal hopen(10,file=file_data,status=old)read(10,*)do i=1,n,1 read(10,*) h,sp(i),msfl(i),lls(i),lld(i),rhob(i),nphi(i),gr(i)enddoclose(10)endsubroutine read_data!-! ! 功能：讀取常數參量! 輸入參數說明：! file_contant:常數參量文件! 輸出參數說明：! grmin,grmax:自然伽馬的最小值與最大值! spmin,spmax:自然電位的最小值與最大值! b,r_sh:已

31、知常數! fai_nsh:已知常數! p_sh,p_ma,p_f:已知常數! fai_nf,fai_nma:已知常數!-subroutine read_contant(file_contant,grmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma)implicit nonecharacter*(*) file_contantreal grmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nmaopen(10,file=file_contan

32、t,status=old)read(10,*) grmin,grmaxread(10,*) spmin,spmaxread(10,*) bread(10,*) r_shread(10,*) fai_nshread(10,*) p_sh,p_ma,p_fread(10,*) fai_nf,fai_nmaclose(10)end subroutine read_contant!-! ! 功能：用自然電位的估計泥質含量! 輸入參數說明：! n:測井數據的個數! sp(:):自然電位! spmin,spmax:自然電位的最小值與最大值! 輸出參數說明：! filename:輸出自然電位估計的泥質含量文

33、件!-subroutine get_sp(filename,n,sp,spmin,spmax)implicit nonecharacter*(*) filenameinteger nreal sp(1:n)real spmin,spmaxinteger ireal v_shopen(10,file=filename,status=unknown)do i=1,n,1 v_sh=(sp(i)-spmin)/(spmax-spmin) write(10,*) v_shenddoclose(10)end subroutine get_sp!-! ! 功能：用電阻率的估計泥質含量! 輸入參數說明：!

34、n:測井數據的個數! r(:):電阻率! b,r_sh:已知常數! 輸出參數說明：! filename:輸出電阻率估計的泥質含量文件!-subroutine get_r(filename,n,r,b,r_sh)implicit nonecharacter*(*) filenameinteger nreal r(1:n)real b,r_shinteger ireal v_shopen(10,file=filename,status=unknown)do i=1,n,1 v_sh=(r_sh/r(i)*(1.0/b) write(10,*) v_shenddoclose(10)end subro

35、utine get_r!-! ! 功能：用中子的估計泥質含量! 輸入參數說明：! n:測井數據的個數! nphi(:):補償中子! fai_nsh:已知常數! 輸出參數說明：! filename:輸出中子估計的泥質含量文件!-subroutine get_nphi(filename,n,nphi,fai_nsh)implicit nonecharacter*(*) filenameinteger nreal nphi(1:n)real fai_nshinteger ireal v_shopen(10,file=filename,status=unknown)do i=1,n,1 v_sh=np

36、hi(i)/fai_nsh write(10,*) v_shenddoclose(10)end subroutine get_nphi!-! ! 功能：用中子的估計泥質含量! 輸入參數說明：! n:測井數據的個數! gr(:):自然伽馬! grmin,grmax:自然伽馬的最小值與最大值! 輸出參數說明：! filename:輸出自然伽馬估計的泥質含量文件!-subroutine get_gr(filename,n,gr,grmin,grmax)implicit nonecharacter*(*) filenameinteger nreal gr(1:n)real grmin,grmaxint

37、eger ireal v_sh,gcurgcur=2.0open(10,file=filename,status=unknown)do i=1,n,1 v_sh=(gr(i)-grmin)/(grmax-grmin) v_sh=(2*(gcur*v_sh)-1)/(2*gcur)-1) write(10,*) v_shenddoclose(10)end subroutine get_gr!-! ! 功能：用中子-密度交會圖的估計泥質含量! 輸入參數說明：! n:測井數據的個數! nphi(:):補償中子! rhob(:):巖性密度! fai_nsh:已知常數! p_sh,p_ma,p_f:已知

38、常數! fai_nf,fai_nma:已知常數! 輸出參數說明：! filename:輸出中子-密度交會圖估計的泥質含量文件!-subroutine get_n_d(filename,n,nphi,rhob,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma)implicit nonecharacter*(*) filenameinteger nreal nphi(1:n),rhob(1:n)real fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nmainteger ireal v_sh,a,b,ca=p_ma-p_fb=fai_nf-fai_nma

39、c=fai_nma*p_f-fai_nf*p_maopen(10,file=filename,status=unknown)do i=1,n,1 v_sh=(abs(a*nphi(i)+b*rhob(i)+c)/(abs(a*fai_nsh+b*p_sh+c) write(10,*) v_shenddoclose(10)end subroutine get_n_d實驗三 含泥質復雜巖性地層綜合測井處理一、實驗目的：綜合所學習的測井解釋與數據處理知識，掌握一般含泥質復雜巖性地層測井數據處理的步驟和方法。二、實驗要求用C語言編程，調試通過，并能計算出正確的結果。三、實驗場地、用具與設備計算中心，計

40、算機；四、實驗內容：（一） 原理：測井數據處理與解釋的主要任務就是依據“四性”關系，用測井所獲得的地層的各種“物理參數”，定量地計算出地下地層的巖性（）、儲集性(孔隙度和滲透率)以及含油性（飽和度Sw）。當儲層含有泥質時，還要求取泥質含量，對測井數據進行泥質校正。對于含有泥質骨架由兩種礦物構成,且礦物含量隨深度變化的復雜巖性地層，其體積模型如圖3所示。這種泥質儲層測井處理的過程見圖4所示。其基本思路就是用幾種方法分別計算出獨立的泥質含量Vsh ，然后取最極小值作為該深度的泥質含量。對測井數據進行泥質含量校正，兩種孔隙度組合求出孔隙度和兩種礦物的含量，在恢復為含泥質地層，計算出實際地層的孔隙度和

41、礦物含量。在計算過程中所設計的參量和模型的變化見圖3所示。VshV1V2恢復為泥質巖石求Sw K，VshNwWV1V2 VshV1V2Vma1Vma2純巖 石 圖3 體積模型與計算過程中的變化 實際上求上述四種體積含量的方法有兩種：（） 三種孔隙度測井聯立解方程（矩陣方法，計算的核心是求逆陣）接出，Vsh , V1 , V2 (2) 獨立方法計算泥質含量，二種孔隙度測井組合，解出純巖石模型的各個部分的體積。讀入一個深度的測井數據孔隙度，礦物1和礦物2構成的體積模型泥質含量Vsh GR 取最小值作為泥質含量Vsh SPRtTHCNL-DEN泥質含量Vsh泥質含量Vsh泥質含量Vsh泥質含量Vsh

42、 中子、密度測井曲線泥質校正Vsh計算含水飽和度Sw矩陣方法計算純地層孔隙度，礦物1和礦物2的含量計算滲透率K讀入測井解釋參數最后一個深度點結束是否圖3 計算步驟圖先求出泥質含量，再校正泥質的影響后的測井數據聯立求解, Vma1 , Vma2 ，相當按純巖石計算，然后再根據泥質的多少換算出實際巖石的孔隙度和礦物含量。， （二）計算過程。（） 根據GR，SP求Vsh 當Vsh 50%時，按泥巖處理 =0，Sw =1.0, Vsh =100, V1 =0, V2 =0.K=0當Vsh 50%時，按下述方法處理。（2）進行泥質校正 （3）求純巖石的孔隙度和礦物含量用矩陣表示為：CX=L其中：C= 則孔隙度和兩種礦物的含量可以用解常數陣的逆的方法得出：X=C-1L（3）恢復成泥質巖石（） 計算含水飽和度（混合泥質公式）利用1963年Simandoux的（混合泥質公式）計算含水飽和度。可以看成是關于Sw的一元二次方程。對于A、B、C分別上述一元二次方程的系數。當a=1.0,m=2.0時