測井資料管理組織與綜合解釋實驗

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

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

3、為LLD、記錄單位Q m ;淺側(cè)向測井曲線一一曲線符號為LLS、記錄單位Q m ;微球形聚焦測井曲線一一曲線符號為MSFL、記錄單位Q m ;電阻率測井曲線通常為對數(shù)刻度。第四道為反映孔隙度的測井曲線道，包括：密度測井曲線一一曲線符號為DEN或RHOB，記錄單位g/cm 3；中子測井曲線一一曲線符號為CNL或PHIN，記錄單位%，有時為v/v。聲波測井曲線一一曲線符號為AC或DT,記錄單位us/ft，有時為us/m 。中子和密度測井曲線的刻度的特點是保證在含水砂巖層上兩條曲線重迭，在含氣層上，密度孔隙度大于中子孔隙度，在泥巖層上，中子孔隙度大于密度孔隙度；第五道是反映粘土礦物類型的測井曲線道，

4、包括自然伽馬能譜測井中的三條曲線： 放射性釷測井曲線一一曲線符號為Th或THOR，記錄單位是ppm ； 放射性鈾測井曲線一一曲線符號為U或URAN，記錄單位ppm ； 放射性鉀測井曲線一一曲線符號為K或POTA，記錄單位%，有時為v/v。2. 測井曲線特征（1）砂泥巖剖面的測井曲線特征砂泥巖剖面儲集層（砂巖）的典型特征是，一般自然電位有明顯的異常，異常的方向和幅度取決于泥漿濾液電阻率（Rmf ）和地層水的電阻率（Rw ）,或者說與 Rmf與Rw的比值有關(guān)，如果Rmf> Rw，則為負異常，否則為正異常。如果砂層中不含放射性礦物，自 然伽瑪曲線亦顯示低值。微電極曲線一般在砂巖層幅值高，并出現(xiàn)

5、正幅差。而泥巖的幅度和幅差均較低，當井眼條件不好時，可能會出現(xiàn)曲線跳動現(xiàn)象。砂巖中含灰質(zhì)較多的夾層，因為致密電阻率異 常高，幅度差很小或沒有。一般幅度差的大小標明了儲集層滲透性的好壞。普通電阻率測井曲線在泥巖處顯示為低值。砂巖處顯示為高值，含油砂巖幅值就更高，如有兩條探測深度不同的 Ra曲線，幅值的差別顯示著低侵、高侵。通常在油層上為低 侵，水層上為高侵。井徑在泥巖層擴大，砂巖層縮?。孕∮阢@頭直徑）。具體特征總結(jié)見表 1表1砂泥巖剖面測井曲線特征7'、地層測井曲線儲集層-砂巖非儲集層-泥巖自然電位負異常（Rw<Rmf ）正異常（Rw>Rmf ）泥巖基線自然伽馬低高井徑縮徑

6、擴徑深中淺電阻率高阻低阻聲波<300us/m>300us/m釷低高鈾低高鉀低高（2 ）碳酸鹽巖剖面的測井曲線特征碳酸鹽巖剖面的測井解釋任務(wù)， 就是從致密的圍巖中找出孔隙性、 裂縫性的儲集層, 并判斷其含油性。碳酸鹽巖剖面電阻率一般較高， 自然電位效果不好。 為了區(qū)分巖性和劃分 儲層，一般使用自然伽馬測井曲線。儲集層相對于致密的圍巖具有低阻、低自然伽馬以及孔隙度測井反映孔隙度較大的特 點。3. 劃分儲集層的基本方法與原則基本要求：凡一切可能含油氣的地層都要劃出來，要適當?shù)貏澐置黠@的水層。具體要求為：（1）估計為油層、氣層、油水同層和含油水層的儲集層都必須分層解釋。（2）厚度半米以上的

7、電性（測井曲線）可疑層（即指從測井曲線上看有油氣的地層） 或錄井顯示為微含油級別以上的儲集層必須分出。（3）選擇出作為確定地層水電阻率 Rw的標準水層（厚度大、巖性純、不含油）要 劃分出來。（4）錄井、氣測有大段油氣顯示而測井曲線顯示不好的儲集層，應(yīng)選取一定層位， 尤其是該組儲層的頂部層位，進行分層。（5）當有多套油水系統(tǒng)，油層組包括若干水層時，只解釋最靠近油層的水層。（6） 對于新區(qū)探井，應(yīng)做細致工作，對各個儲層均應(yīng)酌情選層解釋，以使不漏掉可 能有油氣的地層。4. 正確劃分出儲集層的方法（1）砂泥巖剖面通常是自然電位（SP）曲線的異常確定滲透層的位置， 用微電極曲線確定分層界面， 分層前，應(yīng)

8、將井場收集的井壁取芯、 氣測顯示等有關(guān)油氣顯示的資料標注在綜合測井曲線圖 上，并根據(jù)鄰井的測井和試油等資料對本井的油水關(guān)系作出初步估計。分層時應(yīng)注意：確定分層的界面深度時，應(yīng)左右環(huán)顧，照顧到分層線對每條測井曲線的合理性。 分層的深度誤差不應(yīng)大于 0.1m 。滲透層中，凡是 0.5m以上的非滲透性夾層（泥巖或致密層），應(yīng)將夾層上下 的滲透層分兩層解釋。巖性漸變層頂界（頂部漸變層）或底界（底部漸變層）分層深度應(yīng)在巖性漸變結(jié)束處。一個厚度較大的滲透層，如有兩個以上解釋結(jié)論，應(yīng)按解釋結(jié)論分層。在同一解釋井段，如果油氣層與水層巖性、地層結(jié)構(gòu)和孔隙度基本相同，則油 氣層是純水層的電阻率的3-5倍。純水層的

9、自然電位異常最大，油氣層異常明顯偏小，油水同層介于油、水層之間。并且厚度較大的油水同層，自上而下電阻率有明顯減小的趨勢。（2 ）碳酸鹽巖剖面碳酸巖鹽剖面劃分滲透性地層的5. 測井曲線讀數(shù)分層以后，要從有關(guān)的主要測井曲線將代表該儲層的測井曲線讀數(shù)，以便計算孔隙度、飽和度等地質(zhì)參數(shù)， 在厚度較大的儲集層中按測井曲線變化確定幾個取值區(qū)，對每個取值區(qū)對應(yīng)讀數(shù)計算，幾種主要測井曲線取值區(qū)的最小厚度如下：各種孔隙度測井0.6m。側(cè)向測井0.6m感應(yīng)測井，低阻0.6m，高阻層1.5m。每種測井曲線分層和取值要符合其方法特點，例如聲波測井扣除致密夾層，選用與滲透層相對應(yīng)部分的平均值。 電阻率測井曲線則扣除致密

10、夾層， 選用與滲透層相對應(yīng)部分的 極大值的平均值。另外注意孔隙度與電阻率測井曲線對應(yīng)取值的原則。因為要用兩者結(jié)合計算地層的含水飽和度，兩者當然應(yīng)該是對應(yīng)深度上同一地層或同一取值區(qū)的讀數(shù)。巖層含油性的定性判斷， 主要依據(jù)井曲線的測井曲線特征，而電性特征是巖石物性、巖性和含油性的綜合反映。因此在判斷地層的含油性時，一般應(yīng)將測量井段首先按照地層水 礦化度的不同分為不同的解釋井段，然后才有可能對每一個解釋井段在充分考慮其巖性特點 的前提下進行含油性解釋。由于地下地層復(fù)雜性，儀器的局限性，上述原則是一般性的。要做到正確地解釋,方面應(yīng)多收集資料，認真分析曲線，另一方面還要了解區(qū)域性特點和規(guī)律，要積累經(jīng)驗。

11、6. 計算出孔隙度、飽和度等參數(shù)。讀數(shù)以后，還要做一些定量計算，常用的公式:孔隙度:ttmatmatf含水飽和度：SwaRwmRtma上式中t為當前層的聲波時差，tf為地層水的聲波時差，189us/ft(623us/m)55.5us/ft為固體骨架的聲波時差，對于砂巖骨架，主要礦物為石英，其聲波時差為常數(shù)，對于砂巖地層通常取 1.0 , Rt為當前層的電阻率，m為膠結(jié)指數(shù)。五、實驗結(jié)果處理結(jié)果如圖2所示具體分層如下表所示:頂部深度底部深厚度測井曲線讀數(shù)含水儲層序號(m )度(m)(m )SPGRDEN孔隙度飽和度2.3113.8411122.4114219.6-6.49100.014586.4

12、4103.572.3711.6121154.2116216.8-8.02333586.332115.092.2813.9431218.0512267.95-10.50429883.914104.682.3313.2241236.21252.616.4-14.69560978.8071286.1106.992.3413.9051272.4513.75-17.63688163.182由自然電位的劃分原則可知，在圖2的自然電位的曲線上，有五個明顯的偏離泥基線的負異常，故可以劃分出五個儲集層。然后就可以求出泥質(zhì)含量，進而求出泥質(zhì)校正后的孔隙度, 然后根據(jù)阿爾奇公式，就可求出含水飽和度、滲透率，進而劃分

13、油氣水層。0PE102微球形聚焦20002補償密度30釷30自然電位-3020深度深側(cè)向聲波時差鉀220001404000.5自然伽馬(m)0150淺側(cè)向補償中子扎隙度鈾6井徑16220000.42-0.18020W90井測井曲線圖26802720269027002710273027402750圖1某井的綜合測井曲線圖2770SN183井測井解釋綜合圖MNKt/Qa-JB止MR心圖2SN183井測井解釋綜合圖實驗二利用綜合方法估計地層泥質(zhì)含量一、實驗?zāi)康模和ㄟ^實際計算，鞏固掌握利用多種測井資料確定泥質(zhì)含量的方法。二、實驗要求自編程序，在計算機上運算出地層泥質(zhì)含量。三、實驗場地、用具與設(shè)備計算中

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

15、自然伽馬或自然伽馬能譜或自然電位來求取泥質(zhì)含量效果最好，但自然伽馬要求儲層中除了泥質(zhì)外， 其他物質(zhì)不含放射性礦物。 自然電位要求地層水電阻率保持不變， 且儲層中的泥質(zhì)與相鄰泥巖的的成分相同。用其他方法計算泥質(zhì)含量則要求更為苛刻的條 件：如電阻率方法要求儲層的孔隙度和含水飽和度均要很小。中子和聲波方法則要求孔隙度很小。3.確定Vsh的方法:(1)自然伽瑪法V ' GRGRminGRmaxGRmin2G0URVsh'1VshGcuR/2 1式中，GRmin GRmax分別是砂巖和泥巖層的自然伽馬值，GCUR是與地層有關(guān)的經(jīng)驗系數(shù),新地層（第三系地層）GCUR=3.7，老地層GCUR

16、=2.0.(2)自然電位法VSPSpninSFJnaxSPmin式中，SP是當前層的自然電位讀數(shù)，SPnin和SPmax分別是純地層和泥質(zhì)地層的自然電位讀(3)電阻率Vsh(匹)1小(b=1.5 )Rt(4)中子法VshNNsh式中，n是當前層的視中子孔隙度讀數(shù),Nsh是泥巖層的視中子孔隙度讀數(shù)。（5） 交會圖法以中子一密度測井交會圖為例，通過對圖2所示的石英點（Q ）、水點（W ）和泥巖點（SH ）構(gòu)成的三角形進行分解，依據(jù)資料點所落入三角形中的位置，可以推測出來泥質(zhì)含 量?；蛘呃孟率竭M行計算（依據(jù)點到直線的距離計算方法）V A N B b C sh A NC1 B ci C式中，A n

17、B b C = 0是石英點（Q）和水點（W）連線的直線方程。依據(jù)任意N兩點的直線，用石英點（ma）和水點（Nf，f）兩個點的參數(shù)可以推出:A ( maf) , B ( Nf Nma) , C ( Nma f Nf ma)N ( maf )( Nf Nma ) b Nma ff maNsh (maf )( Nf Nma) sh Nma f f ma當然，也可以用中子一聲波、聲波一密度交會圖的類似方法求Vsh。如果解釋層段上沒有純泥巖層時，上述交會圖法所定出的泥巖點位置并不代表實際的泥巖參數(shù),導致用交會圖方法估計的泥質(zhì)含量比實際值偏大。上述幾種方法計算的泥質(zhì)含量往往都有一定的條件，當條件滿足時,泥

18、質(zhì)含量的近似結(jié)果，當條件不滿足時，計算的泥質(zhì)含量均可能偏高，所以在實際處理時,最后選取其中的最小值作為接近實際的泥質(zhì)含量。五、實驗結(jié)果1.已知條件：GRmi n=50GRmax=110SPmi n=0 SPmax=100b=1.5 Rsh=5Nsh30 push2.54ma2.65f1 .0Nf100Nma42.測井數(shù)據(jù)depthspmsflllslldrhobn phigrdtcal2680.9888.504.809.0910.842.5517.00103.1773.1311.992682.0983.402.026.588.212.0720.00108.2577.7514.192683.01

19、70.001.585.717.171.5924.00107.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.2

20、068.9450.0651.542.5011.0053.4865.949.262690.0227.0055.8145.6645.832.5210.0064.0665.319.163.計算結(jié)果spmsflllslldrhob-n phin phigr泥質(zhì)含量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.80

21、0001.892830.889510.764500.609270.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.2

22、11130.204880.366670.027910.027910.270000.200220.228890.228320.214510.333331.279430.20022附錄：程序源代碼!功能:利用測井數(shù)據(jù)求泥質(zhì)含量!時間：2016/12/22!參數(shù)說明：!cmdfile:參數(shù)文件!file_data:輸入測井數(shù)據(jù)文件! file_co nta nt: 常數(shù)參量文件!file_sp:輸出自然電位估計的泥質(zhì)含量文件file_msfl:輸出微球型聚焦估計的泥質(zhì)含量文件file lls:輸出淺側(cè)向估計的泥質(zhì)含量文件 file_lld:輸出深側(cè)向估計的泥質(zhì)含量文件file_nphi:輸出補償中子

23、估計的泥質(zhì)含量文件file_gr:輸出自然伽馬估計的泥質(zhì)含量文件file_n_d:輸出中子-密度交會圖估計的泥質(zhì)含量文件 n:測井數(shù)據(jù)的個數(shù)grmi n,grmax:自然伽馬的最小值與最大值spmi n, spmax:自然電位的最小值與最大值b,r_sh:已知常數(shù)fai_nsh:已知常數(shù)P_sh,p_ma,p_f:已知常數(shù)fai_n f,fai_ nma:已知常數(shù)sp(:):自然電位!gr(:):自然伽馬!rhob(:):巖性密度!program cjsyimplicit nonecharacter*80 cmdfilecharacter*80 file_data,file_c ontantc

24、haracter*80 file_sp,file_msfl,filels,fileld,file _n phi,file_gr,file_ n_din teger nreal grmi n,grmax,spmi n,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.txt'callread_cmd(cmdfile,file_data,file_c on ta nt,f

25、ile_sp,file_msfl,filels,fileld,file _n phi,file_gr,file_ n_d)call read_ nu mber(file_data ,n)allocate(sp(1: n),msfl(1: n),lls(1: n),lld(1: n),n phi(1: n),gr(1: n),rhob(1: n)call read_data(file_data ,n, sp,msfl,lls,lld ,n phi,gr,rhob)callread_c ontan t(file_c ontan t,grm in, grmax,spm in, spmax,b,r_s

26、h,fai_ nsh,p_sh,p_ma,p_f,fai_n f,fai_ nma) call get_sp(file_sp ,n, sp,spm in, 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,n phi,fai_ nsh)call get_gr(file_gr, n,gr,grmi n,grmax)call get_ n_d(file_ n_d,n

27、,n phi,rhob,fai_ nsh,p_sh,p_ma,p_f,fai_ nf,fai_ nma)deallocate(sp,msfl,lls,lld, nphi,gr,rhob)end輸入?yún)?shù)說明：cmdfile:參數(shù)文件輸出參數(shù)說明：file_data:輸入測井數(shù)據(jù)文件file_co nta nt: 常數(shù)參量文件file_sp:輸出自然電位估計的泥質(zhì)含量文件file_msfl:輸出微球型聚焦估計的泥質(zhì)含量文件file_lls:輸出淺側(cè)向估計的泥質(zhì)含量文件file_lld:輸出深側(cè)向估計的泥質(zhì)含量文件file_nphi:輸出補償中子估計的泥質(zhì)含量文件file_gr:輸出自然伽馬估計的泥

28、質(zhì)含量文件file_n_d:輸出中子-密度交會圖估計的泥質(zhì)含量文件subrouti neread_cmd(cmdfile,file_data,file_c on ta nt,file_sp,file_msfl,filels,fileld,file _n phi,file_ gr,file_ n_d) implicit nonecharacter*) cmdfilecharacter*) file_data,file_contantcharacter*(*) file_sp,file_msfl,filels,fileld,file _n phi,file_gr,file _n_d ope n(1

29、0,file=cmdfile,status='old')read(10,*) file_dataread(10,*) file_co nta ntread(10,*) file_spread(10,*) file_msflread(10,*) file_llsread(10,*) file_lld read(10,*) file_ nphi read(10,*) file_grread(10,*) file_n_dclose(10)end subrouti ne read_cmd!輸入?yún)?shù)說明：!file name:輸入測井數(shù)據(jù)文件!輸出參數(shù)說明：!number:測井數(shù)據(jù)的個數(shù)

30、!subrouti ne read_ nu mber(file name ,nu mber) implicit none character*(*) file name in teger nu mberope n(11,file=file name,status='old')nu mber=0do while( .n ot.eof(11)read (11,*,e nd=100,ERR=100)nu mber =nu mber+1100 end doclose(11)nu mber =nu mber-1end subrouti ne read_ nu mber!輸入?yún)?shù)說明：!f

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

32、(1: n)in teger ireal h ope n(10,file=file_data,status='old') read(10,*)do i=1, n,1read(10,*) h,sp(i),msfl(i),lls(i),lld(i),rhob(i), nphi(i),gr(i) en ddo close(10)end subrout ine read_data輸入?yún)?shù)說明：file_co nta nt: 常數(shù)參量文件!輸出參數(shù)說明：!grmi n,grmax:自然伽馬的最小值與最大值!spmi n, spmax:自然電位的最小值與最大值!b,r_sh:已知常數(shù)!fa

33、i_nsh:已知常數(shù)! p_sh,p_ma,p_f:已知常數(shù)!fai_n f,fai_ nma:已知常數(shù)!subrouti neread_c ontan t(file_c ontan t,grm in, grmax,spm in, spmax,b,r_sh,fai_ nsh,p_sh,p_ma,p_f,fa i_n f,fai_ nma)implicit nonecharacter*) file_c ontantreal grmi n,grmax,spmi n,spmax,b,r_sh,fai_ nsh,p_sh,p_ma,p_f,fai_ nf,fai_ nma ope n(10,file=

34、file_contan t,status='old')read(10,*) grmi n,grmaxread(10,*) spm in, spmaxread(10,*) bread(10,*) r_sh read(10,*) fai_nshread(10,*) p_sh,p_ma,p_fread(10,*) fai_n f,fai_ nmaclose(10) end subrout ine read_c ontant!功能：用自然電位的估計泥質(zhì)含量!輸入?yún)?shù)說明：!n:測井數(shù)據(jù)的個數(shù)!sp(:):自然電位spmi n, spmax: 自然電位的最小值與最大值輸出參數(shù)說明:fil

35、ename:輸出自然電位估計的泥質(zhì)含量文件subrouti ne get_sp(file name ,n, sp,spm in, spmax) implicit nonecharacter*(*) file name in teger nreal sp(1: n) real spm in, spmaxin teger ireal v_shope n(10,file=file name,status='unknown') do i=1, n,1v_sh=(sp(i)-spm in)/(spmax-spm in) write(10,*) v_shen ddoclose(10)end

36、 subrouti ne get_sp功能：用電阻率的估計泥質(zhì)含量輸入?yún)?shù)說明：n:測井數(shù)據(jù)的個數(shù)r(:):電阻率b,r_sh:已知常數(shù)輸出參數(shù)說明：file name:輸出電阻率估計的泥質(zhì)含量文件!subroutine get_r(filename,n,r,b,r_sh)implicit nonecharacter*) file namein teger nreal r(1: n)real b,r_sh in teger ireal v_shope n(10,file=file name,status='unknown') do i=1, n,1v_sh=(r_sh/r(i)

37、*(1.0/b)write(10,*) v_shen ddo close(10)end subrout ine get_r功能：用中子的估計泥質(zhì)含量輸入?yún)?shù)說明:!n:測井數(shù)據(jù)的個數(shù)!nphi(:):補償中子!fai_nsh:已知常數(shù)!輸出參數(shù)說明：!file name:輸出中子估計的泥質(zhì)含量文件!subroutine get_nphi(filename,n,nphi,fai_nsh) implicit none character*) file namein teger nreal n phi(1: n)real fai_nshin teger ireal v_shope n(10,file

38、=file name,status='unknown')do i=1, n,1v_sh=n phi(i)/fai_ nshwrite(10,*) v_shen ddoclose(10) end subrout ine get_ nphi!輸入?yún)?shù)說明：!n:測井數(shù)據(jù)的個數(shù)!gr(:):自然伽馬!grmi n,grmax:自然伽馬的最小值與最大值!輸出參數(shù)說明：!file name:輸出自然伽馬估計的泥質(zhì)含量文件!subroutine get_gr(filename,n,gr,grmin,grmax)implicit nonecharacter*) file namein teg

39、er nreal gr(1: n)real grmi n,grmaxin teger ireal v_sh,gcurgcur=2.0 ope n(10,file=file name,status='unknown') do i=1, n,1v_sh=(gr(i)-grmi n)/(grmax-grmi n) v_sh=(2*(gcur*v_sh)-1)/(2*gcur)-1) write(10,*) v_shen ddoclose(10)end subrouti ne get_gr!功能：用中子-密度交會圖的估計泥質(zhì)含量!輸入?yún)?shù)說明：!n:測井數(shù)據(jù)的個數(shù)!nphi(:):補償

40、中子!rhob(:):巖性密度!fai_nsh:已知常數(shù)! p_sh,p_ma,p_f:已知常數(shù)!fai_n f,fai_ nma:已知常數(shù)!輸出參數(shù)說明：!file name:輸出中子-密度交會圖估計的泥質(zhì)含量文件!subrouti ne get_ n_d(file name ,n,n phi,rhob,fai_ nsh,p_sh,p_ma,p_f,fai_ nf,fai_ nma) implicit nonecharacter*) file namein teger nreal n phi(1: n),rhob(1: n)real fai_n sh,p_sh,p_ma,p_f,fai_ n

41、f,fai_ nmain teger ireal v_sh,a,b,ca=p_ma-p_fb=fai_ nf-fai_ nmaope n(10,file=file name,status='unknown')do i=1, n,1v_sh=(abs(a* nphi(i)+b*rhob(i)+c)/(abs(a*fai_ nsh+b*p_sh+c)write(10,*) v_shen ddoclose(10)實驗三含泥質(zhì)復(fù)雜巖性地層綜合測井處理一、實驗?zāi)康模壕C合所學習的測井解釋與數(shù)據(jù)處理知識，掌握一般含泥質(zhì)復(fù)雜巖性地層測井 數(shù)據(jù)處理的步驟和方法。二、實驗要求用C語言編程，調(diào)試通過

42、，并能計算出正確的結(jié)果。三、實驗場地、用具與設(shè)備計算中心，計算機；四、實驗內(nèi)容：（一）原理：測井數(shù)據(jù)處理與解釋的主要任務(wù)就是依據(jù)“四性”關(guān)系，用測井所獲得的地層的各種“物理參數(shù)”，定量地計算出地下地層的巖性（ Vma1.VmaN ）、儲集性（孔隙度 和滲透率 ） 以及含油性（飽和度 Sw ）。當儲層含有泥質(zhì)時，還要求取泥質(zhì)含量，對測井數(shù)據(jù)進行泥質(zhì) 校正。對于含有泥質(zhì)骨架由兩種礦物構(gòu)成，且礦物含量隨深度變化的復(fù)雜巖性地層，其體積模型如圖3所示。這種泥質(zhì)儲層測井處理的過程見圖4所示。其基本思路就是用幾種方法分別計算出獨立的泥質(zhì)含量 Vsh，然后取最極小值作為該深度的泥質(zhì)含量。對測井數(shù)據(jù)進行泥質(zhì)含量

43、校正，兩種孔隙度組合求出孔隙度和兩種礦物的含量，在恢復(fù)為含泥質(zhì)地層， 計算3所示。出實際地層的孔隙度和礦物含量。在計算過程中所設(shè)計的參量和模型的變化見圖VshViV2求Sw K圖3體積模型與計算過程中的變化實際上求上述四種體積含量的方法有兩種：(1)三種孔隙度測井聯(lián)立解方程(矩陣方法，計算的核心是求逆陣)接出0, Vsh , Vi , V2(2)獨立方法計算泥質(zhì)含量，二種孔隙度測井組合， 解出純巖石模型的各個部分的體積。圖3計算步驟圖先求出泥質(zhì)含量，再校正泥質(zhì)的影響后的測井數(shù)據(jù)聯(lián)立求解,Vmal , Vma2 ,相當按純巖石計算，然后再根據(jù)泥質(zhì)的多少換算出實際巖石的孔隙度和礦物含量。，(二)計

44、算過程。(1)根據(jù) GR, SP 求 VshGR GRminGRGRmaxGRminVspSP Spninmin (Vgr,Vsp)SFJnax SPmin當 Vsh >50% 時，按泥巖處理 0=0 , Sw =1.0, V sh =100, V 1 =0, V 2 =0.K=0當Vsh <50%時，按下述方法處理。(2)進行泥質(zhì)校正(NNCbC ( bVsh)Vsh)(3 )求純巖石的孔隙度和礦物含量NCNma1Nma2Vma2NfbCma1Vma1ma2Vma21 V ma1Vma2用矩陣表示為:CX=L其中：NmalNma2NfC=ma1ma2f111Vma1V ma2NCbC1則孔隙度和兩種礦物的含量可以用解常數(shù)陣的逆的方法得出:X=C -1L(3 )恢復(fù)成泥質(zhì)巖石ViVmai(1Vsh)V2Vma2(1 Vsh)'(1Vsh)(4) 計算含水飽和度(混合泥質(zhì)公式)利用1963年Simandoux 的(混合泥質(zhì)公式)計算含水飽和度。1RtRshm可以看