測(cè)井資料處理與綜合解釋實(shí)驗(yàn)

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

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

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

4、h 或 THOR ，記錄單位是ppm；放射性鈾測(cè)井曲線曲線符號(hào)為U 或 URAN ，記錄單位ppm；放射性鉀測(cè)井曲線曲線符號(hào)為K 或 POTA ，記錄單位 % ，有時(shí)為v/v 。2測(cè)井曲線特征（ 1）砂泥巖剖面的測(cè)井曲線特征砂泥巖剖面儲(chǔ)集層（砂巖）的典型特征是，一般自然電位有明顯的異常，異常的方向和幅度取決于泥漿濾液電阻率（ Rmf ）和地層水的電阻率（ Rw），或者說(shuō)與 Rmf 與 Rw 的比值有關(guān) ,如果 Rmf> Rw ，則為負(fù)異常，否則為正異常。如果砂層中不含放射性礦物，自然伽瑪曲線亦顯示低值。微電極曲線一般在砂巖層幅值高，并出現(xiàn)正幅差。而泥巖的幅度和幅差均較低，當(dāng)井眼條件不好時(shí)

5、， 可能會(huì)出現(xiàn)曲線跳動(dòng)現(xiàn)象。 砂巖中含灰質(zhì)較多的夾層， 因?yàn)橹旅茈娮杪十惓８撸炔詈苄』驔](méi)有。一般幅度差的大小標(biāo)明了儲(chǔ)集層滲透性的好壞。普通電阻率測(cè)井曲線在泥巖處顯示為低值。砂巖處顯示為高值，含油砂巖幅值就更高，如有兩條探測(cè)深度不同的Ra 曲線，幅值的差別顯示著低侵、高侵。通常在油層上為低侵，水層上為高侵。井徑在泥巖層擴(kuò)大，砂巖層縮?。孕∮阢@頭直徑）。具體特征總結(jié)見(jiàn)表1表 1砂泥巖剖面測(cè)井曲線特征地層儲(chǔ)集層砂巖非儲(chǔ)集層泥巖測(cè)井曲線自然電位負(fù)異常（ Rw<Rmf ）泥巖基線正異常（ Rw>Rmf ）自然伽馬低高井徑縮徑擴(kuò)徑深中淺電阻率高阻低阻聲波<300us/m>30

6、0us/m釷低高鈾低高鉀低高（ 2）碳酸鹽巖剖面的測(cè)井曲線特征碳酸鹽巖剖面的測(cè)井解釋任務(wù)， 就是從致密的圍巖中找出孔隙性、 裂縫性的儲(chǔ)集層，并判斷其含油性。 碳酸鹽巖剖面電阻率一般較高， 自然電位效果不好。 為了區(qū)分巖性和劃分儲(chǔ)層，一般使用自然伽馬測(cè)井曲線。儲(chǔ)集層相對(duì)于致密的圍巖具有低阻、低自然伽馬以及孔隙度測(cè)井反映孔隙度較大的特點(diǎn)。3劃分儲(chǔ)集層的基本方法與原則基本要求：凡一切可能含油氣的地層都要?jiǎng)澇鰜?lái)，要適當(dāng)?shù)貏澐置黠@的水層。具體要求為：（ 1）估計(jì)為油層、氣層、油水同層和含油水層的儲(chǔ)集層都必須分層解釋。（ 2）厚度半米以上的電性（測(cè)井曲線）可疑層（即指從測(cè)井曲線上看有油氣的地層）或錄井顯示

7、為微含油級(jí)別以上的儲(chǔ)集層必須分出。（ 3）選擇出作為確定地層水電阻率Rw 的標(biāo)準(zhǔn)水層（厚度大、巖性純、不含油）要?jiǎng)澐殖鰜?lái)。（ 4）錄井、氣測(cè)有大段油氣顯示而測(cè)井曲線顯示不好的儲(chǔ)集層，應(yīng)選取一定層位，尤其是該組儲(chǔ)層的頂部層位，進(jìn)行分層。（ 5）當(dāng)有多套油水系統(tǒng)，油層組包括若干水層時(shí)，只解釋最靠近油層的水層。（ 6）對(duì)于新區(qū)探井，應(yīng)做細(xì)致工作，對(duì)各個(gè)儲(chǔ)層均應(yīng)酌情選層解釋，以使不漏掉可能有油氣的地層。4正確劃分出儲(chǔ)集層的方法（ 1）砂泥巖剖面通常是自然電位 （ SP）曲線的異常確定滲透層的位置，用微電極曲線確定分層界面，分層前， 應(yīng)將井場(chǎng)收集的井壁取芯、 氣測(cè)顯示等有關(guān)油氣顯示的資料標(biāo)注在綜合測(cè)井曲

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

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

10、與電阻率測(cè)井曲線對(duì)應(yīng)取值的原則。 因?yàn)橐脙烧呓Y(jié)合計(jì)算地層的含水飽和度，兩者當(dāng)然應(yīng)該是對(duì)應(yīng)深度上同一地層或同一取值區(qū)的讀數(shù)。巖層含油性的定性判斷， 主要依據(jù)井曲線的測(cè)井曲線特征， 而電性特征是巖石物性、巖性和含油性的綜合反映。 因此在判斷地層的含油性時(shí)， 一般應(yīng)將測(cè)量井段首先按照地層水礦化度的不同分為不同的解釋井段， 然后才有可能對(duì)每一個(gè)解釋井段在充分考慮其巖性特點(diǎn)的前提下進(jìn)行含油性解釋。由于地下地層復(fù)雜性，儀器的局限性，上述原則是一般性的。要做到正確地解釋，一方面應(yīng)多收集資料，認(rèn)真分析曲線，另一方面還要了解區(qū)域性特點(diǎn)和規(guī)律，要積累經(jīng)驗(yàn)。6計(jì)算出孔隙度、飽和度等參數(shù)。讀數(shù)以后，還要做一些定量計(jì)

11、算，常用的公式：ttma孔隙度：tmat f含水飽和度： SwaRwm Rt上式中 t 為當(dāng)前層的聲波時(shí)差，t f 為地層水的聲波時(shí)差,189us/ft(623us/m) ，tma 為固體骨架的聲波時(shí)差，對(duì)于砂巖骨架，主要礦物為石英，其聲波時(shí)差為55.5us/ft。a 是常數(shù)，對(duì)于砂巖地層通常取 1.0， Rt 為當(dāng)前層的電阻率， m 為膠結(jié)指數(shù)。五、實(shí)驗(yàn)結(jié)果處理結(jié)果如圖2 所示具體分層如下表所示：儲(chǔ)層序號(hào)頂部深度底部深厚度測(cè)井曲線讀數(shù)孔隙度含水（m）度（ m）（ m）飽和度SPGRDEN11122.4114219.6-6.49100.012.31413.84586.4421154.21162

12、16.8-8.023103.5732.37311.61586.33231218.0512267.95-10.504115.0922.28913.94883.91441236.21252.616.4-14.695104.6862.33013.22978.80751272.41286.1513.75-17.636106.9982.34813.90163.182由自然電位的劃分原則可知，在圖 2 的自然電位的曲線上，有五個(gè)明顯的偏離泥基線的負(fù)異常，故可以劃分出五個(gè)儲(chǔ)集層。然后就可以求出泥質(zhì)含量，進(jìn)而求出泥質(zhì)校正后的孔隙度，然后根據(jù)阿爾奇公式，就可求出含水飽和度、滲透率，進(jìn)而劃分油氣水層。0PE10-

13、30自然電位200自然伽馬1506井徑162深度(m) 2226802690270027102720273027402750W903井測(cè)井曲線圖微球形聚焦補(bǔ)償密度釷200023030深側(cè)向聲波時(shí)差鉀20001404000.5淺側(cè)向補(bǔ)償中子孔隙度鈾20000.42-0.18020圖 1某井的綜合測(cè)井曲線圖277012345圖 2SN183 井測(cè)井解釋綜合圖實(shí)驗(yàn)二利用綜合方法估計(jì)地層泥質(zhì)含量一、實(shí)驗(yàn)?zāi)康模和ㄟ^(guò)實(shí)際計(jì)算，鞏固掌握利用多種測(cè)井資料確定泥質(zhì)含量的方法。二、實(shí)驗(yàn)要求自編程序，在計(jì)算機(jī)上運(yùn)算出地層泥質(zhì)含量。三、實(shí)驗(yàn)場(chǎng)地、用具與設(shè)備計(jì)算中心，尺子、像皮和計(jì)算機(jī)；四、實(shí)驗(yàn)容：什么是泥質(zhì)含量：泥質(zhì)

14、是指顆粒直徑小于0.01mm的碎屑物質(zhì)，泥質(zhì)含量，也叫做泥質(zhì)體積，是指泥質(zhì)的體積占巖石總體積的比：VshV泥100%V巖確定 V sh 的重要性泥質(zhì)含量的確定， 在泥質(zhì)砂巖儲(chǔ)集層的定量解釋中具有重要意義。多年來(lái)人們提出許多計(jì)算泥質(zhì)含量的理論和方法。目前求取泥質(zhì)含量的方法大致可分為兩類，一類是用每種測(cè)點(diǎn)各求出一個(gè)泥質(zhì)含量，然后求出最佳值。當(dāng)巖石含有泥質(zhì)時(shí)，各種測(cè)井曲線均或多或少地受到泥質(zhì)的影響，其影響的程度受V sh的決定，評(píng)價(jià)巖石的特性時(shí)，只有已知V sh ，才知道由于泥質(zhì)帶來(lái)的影響，從而將泥質(zhì)的影響校正掉。一般而言， 用自然伽馬或自然伽馬能譜或自然電位來(lái)求取泥質(zhì)含量效果最好，但自然伽馬要求儲(chǔ)

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

16、min式中， SP是當(dāng)前層的自然電位讀數(shù)，SP和 SP分別是純地層和泥質(zhì)地層的自然電位讀minm ax數(shù)（）電阻率Vsh( Rsh )1 / b（ b=1.5）Rt（）中子法VshNNsh式中，N 是當(dāng)前層的視中子孔隙度讀數(shù)，Nsh 是泥巖層的視中子孔隙度讀數(shù)。（）交會(huì)圖法以中子密度測(cè)井交會(huì)圖為例，通過(guò)對(duì)圖2 所示的石英點(diǎn)（Q）、水點(diǎn)（ W ）和泥巖點(diǎn)（ SH）構(gòu)成的三角形進(jìn)行分解，依據(jù)資料點(diǎn)所落入三角形中的位置，可以推測(cè)出來(lái)泥質(zhì)含量。或者利用下式進(jìn)行計(jì)算（依據(jù)點(diǎn)到直線的距離計(jì)算方法）：VshAANNc1BBbc1CC式中， ANBbC 0 是石英點(diǎn)（ Q）和水點(diǎn)（ W ）連線的直線方程。依據(jù)

17、任意兩WbQSHN點(diǎn)的直線，用石英點(diǎn)(Nma ,ma )和水點(diǎn) ( Nf ,f )兩個(gè)點(diǎn)的參數(shù)可以推出：A(maf ) ， B ( NfNma) ， C( Nma f Nf ma )N (maf )( Vshf )(Nsh (maNfNma )bNmaffmaNfNma )shNmaffma當(dāng)然，也可以用中子聲波、 聲波密度交會(huì)圖的類似方法求Vsh 。如果解釋層段上沒(méi)有純泥巖層時(shí)， 上述交會(huì)圖法所定出的泥巖點(diǎn)位置并不代表實(shí)際的泥巖參數(shù)，導(dǎo)致用交會(huì)圖方法估計(jì)的泥質(zhì)含量比實(shí)際值偏大。上述幾種方法計(jì)算的泥質(zhì)含量往往都有一定的條件，當(dāng)條件滿足時(shí)，泥質(zhì)含量的近似結(jié)果， 當(dāng)條件不滿足時(shí)， 計(jì)算的泥質(zhì)含量均

18、可能偏高，所以在實(shí)際處理時(shí)，最后選取其中的最小值作為接近實(shí)際的泥質(zhì)含量。五、實(shí)驗(yàn)結(jié)果1.已知條件： GRmin=50GRmax=110SPmin=0 SPmax=100b=1.5 Rsh=5Nsh30 push2.54ma 2.65f 1.0Nf 100Nma42.測(cè)井?dāng)?shù)據(jù)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.0

19、0107.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.00

20、53.4865.949.262690.0227.0055.8145.6645.832.5210.0064.0665.319.163.計(jì)算結(jié)果spmsflllslldrhob-nphinphigr泥質(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.800001.892830.889510.764500.6

21、09270.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.

22、027910.270000.200220.228890.228320.214510.333331.279430.20022附錄：程序源代碼! 功能 :利用測(cè)井?dāng)?shù)據(jù)求泥質(zhì)含量! 時(shí)間 :2016/12/22!-! 參數(shù)說(shuō)明 :!cmdfile: 參數(shù)文件!file_data: 輸入測(cè)井?dāng)?shù)據(jù)文件! file_contant: 常數(shù)參量文件!file_sp: 輸出自然電位估計(jì)的泥質(zhì)含量文件!file_msfl: 輸出微球型聚焦估計(jì)的泥質(zhì)含量文件!file_lls: 輸出淺側(cè)向估計(jì)的泥質(zhì)含量文件!file_lld: 輸出深側(cè)向估計(jì)的泥質(zhì)含量文件!file_nphi: 輸出補(bǔ)償中子估計(jì)的泥質(zhì)含量文件!f

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

24、nonecharacter*80 cmdfilecharacter*80 file_data,file_contantcharacter*80 file_sp,file_msfl,file_lls,file_lld,file_nphi,file_gr,file_n_d integer nreal grmin,grmax,spmin,spmax,b,r_sh,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma real,allocatable:sp(:),msfl(:),lls(:),lld(:),nphi(:),gr(:),rhob(:)cmdfile='cmd.

25、txt'callread_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)callread_contant(file_contan

26、t,grmin,grmax,spmin,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(fil

27、e_n_d,n,nphi,rhob,fai_nsh,p_sh,p_ma,p_f,fai_nf,fai_nma)deallocate(sp,msfl,lls,lld,nphi,gr,rhob)end!-! 功能：讀取參數(shù)文件! 輸入?yún)?shù)說(shuō)明：!cmdfile: 參數(shù)文件! 輸出參數(shù)說(shuō)明：!file_data: 輸入測(cè)井?dāng)?shù)據(jù)文件! file_contant: 常數(shù)參量文件!file_sp: 輸出自然電位估計(jì)的泥質(zhì)含量文件!file_msfl: 輸出微球型聚焦估計(jì)的泥質(zhì)含量文件!file_lls: 輸出淺側(cè)向估計(jì)的泥質(zhì)含量文件!file_lld: 輸出深側(cè)向估計(jì)的泥質(zhì)含量文件!file_nphi:

28、輸出補(bǔ)償中子估計(jì)的泥質(zhì)含量文件!file_gr: 輸出自然伽馬估計(jì)的泥質(zhì)含量文件!file_n_d: 輸出中子 -密度交會(huì)圖估計(jì)的泥質(zhì)含量文件!-subroutineread_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,fil

29、e_lld,file_nphi,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)end subroutine read_cmd!-! 功能：讀取測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)!

30、輸入?yún)?shù)說(shuō)明：!filename: 輸入測(cè)井?dāng)?shù)據(jù)文件! 輸出參數(shù)說(shuō)明：!number: 測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)!-subroutine read_number(filename,number)implicit nonecharacter*(*) filenameinteger numberopen(11,file=filename,status='old')number=0do while(.not.eof(11)read (11,*,end=100,ERR=100)number=number+1100 end doclose(11)number=number-1end subrout

31、ine read_number!-! 功能：讀取測(cè)井?dāng)?shù)據(jù)! 輸入?yún)?shù)說(shuō)明：!file_data: 輸入測(cè)井?dāng)?shù)據(jù)文件!n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)! 輸出參數(shù)說(shuō)明：!sp(:): 自然電位!msfl(:): 微球型聚焦!lls(:): 淺側(cè)向!lld(:): 深側(cè)向!nphi(:): 補(bǔ)償中子!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),l

32、ls(1:n),lld(1:n),nphi(1:n),gr(1:n),rhob(1:n)integer ireal hopen(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)enddoclose(10)end subroutine read_data!-! 功能：讀取常數(shù)參量! 輸入?yún)?shù)說(shuō)明：! file_contant: 常數(shù)參量文件! 輸出參數(shù)說(shuō)明：! grmin,grmax: 自然伽馬的最小值與

33、最大值! spmin,spmax: 自然電位的最小值與最大值!b,r_sh:已知常數(shù)!fai_nsh:已知常數(shù)! p_sh,p_ma,p_f: 已知常數(shù)! fai_nf,fai_nma: 已知常數(shù)!-subroutineread_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

34、,fai_nf,fai_nmaopen(10,file=file_contant,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!-! 功能：用自然電位的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說(shuō)明：!n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)!sp(:): 自然電位! spmin,spmax:

35、 自然電位的最小值與最大值! 輸出參數(shù)說(shuō)明：!filename: 輸出自然電位估計(jì)的泥質(zhì)含量文件!-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,1v_sh=(sp(i)-spmin)/(spmax-spmin)write(10,*) v_shenddoclose(

36、10)end subroutine get_sp!-! 功能：用電阻率的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說(shuō)明：!n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)!r(:): 電阻率!b,r_sh:已知常數(shù)! 輸出參數(shù)說(shuō)明：!filename: 輸出電阻率估計(jì)的泥質(zhì)含量文件!-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,1v_s

37、h=(r_sh/r(i)*(1.0/b)write(10,*) v_shenddoclose(10)end subroutine get_r!-! 功能：用中子的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說(shuō)明：!n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)!nphi(:): 補(bǔ)償中子fai_nsh:已知常數(shù)! 輸出參數(shù)說(shuō)明：!filename: 輸出中子估計(jì)的泥質(zhì)含量文件!-subroutine get_nphi(filename,n,nphi,fai_nsh)implicit nonecharacter*(*) filenameinteger nreal nphi(1:n)real fai_nshinteger ireal v_sho

38、pen(10,file=filename,status='unknown')do i=1,n,1v_sh=nphi(i)/fai_nshwrite(10,*) v_shenddoclose(10)end subroutine get_nphi!-! 功能：用中子的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說(shuō)明：!n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)!gr(:): 自然伽馬! grmin,grmax: 自然伽馬的最小值與最大值! 輸出參數(shù)說(shuō)明：!filename: 輸出自然伽馬估計(jì)的泥質(zhì)含量文件!-subroutine get_gr(filename,n,gr,grmin,grmax)implicit nonech

39、aracter*(*) filenameinteger nreal gr(1:n)real grmin,grmaxinteger ireal v_sh,gcurgcur=2.0open(10,file=filename,status='unknown')do i=1,n,1v_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!-! 功能：用中子 -密度交會(huì)圖的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說(shuō)明：!n:測(cè)井

40、數(shù)據(jù)的個(gè)數(shù)!nphi(:): 補(bǔ)償中子!rhob(:): 巖性密度!fai_nsh:已知常數(shù)! p_sh,p_ma,p_f: 已知常數(shù)! fai_nf,fai_nma: 已知常數(shù)! 輸出參數(shù)說(shuō)明：!filename: 輸出中子 -密度交會(huì)圖估計(jì)的泥質(zhì)含量文件!-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

41、_f,fai_nf,fai_nmainteger ireal v_sh,a,b,ca=p_ma-p_fb=fai_nf-fai_nmac=fai_nma*p_f-fai_nf*p_maopen(10,file=filename,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_shenddoclose(10)end subroutine get_n_d實(shí)驗(yàn)三含泥質(zhì)復(fù)雜巖性地層綜合測(cè)井處理一、實(shí)驗(yàn)?zāi)康模壕C合所學(xué)習(xí)的測(cè)井解釋與數(shù)據(jù)處理知識(shí)，掌握一般含泥