版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、 測(cè)井資料處理與解釋測(cè)井資料處理與綜合解釋實(shí)驗(yàn)專業(yè)名稱:勘查技術(shù)與工程 學(xué)生姓名:張奎濤學(xué)生學(xué)號(hào):201326020109指導(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)康模和ㄟ^對(duì)測(cè)井曲線特征的分析和認(rèn)識(shí),掌握定性劃分砂泥巖剖面儲(chǔ)集層的基本方法,并應(yīng)用阿爾奇公式,進(jìn)行儲(chǔ)層參數(shù)的計(jì)算,鞏固已經(jīng)學(xué)過的鉆井地球物理課程的主要內(nèi)容與應(yīng)用。二、 實(shí)驗(yàn)要求正確劃分出儲(chǔ)集層和非儲(chǔ)集層,對(duì)砂泥巖剖面能區(qū)分開較明顯的油水層。進(jìn)行測(cè)井曲線讀數(shù),簡(jiǎn)
2、單地計(jì)算出孔隙度、飽和度等參數(shù)。三、 實(shí)驗(yàn)場(chǎng)地、用具與設(shè)備測(cè)井實(shí)驗(yàn)室或一般的教室,長(zhǎng)直尺、鉛筆、像皮和計(jì)算器四、 實(shí)驗(yàn)內(nèi)容: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è)
3、向測(cè)井曲線曲線符號(hào)為L(zhǎng)LS、記錄單位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)為Th或T
4、HOR,記錄單位是ppm;放射性鈾測(cè)井曲線曲線符號(hào)為U或URAN,記錄單位ppm;放射性鉀測(cè)井曲線曲線符號(hào)為K或POTA,記錄單位%,有時(shí)為v/v。2 測(cè)井曲線特征(1)砂泥巖剖面的測(cè)井曲線特征砂泥巖剖面儲(chǔ)集層(砂巖)的典型特征是,一般自然電位有明顯的異常,異常的方向和幅度取決于泥漿濾液電阻率(Rmf)和地層水的電阻率(Rw),或者說與Rmf與Rw的比值有關(guān),如果Rmf Rw,則為負(fù)異常,否則為正異常。如果砂層中不含放射性礦物,自然伽瑪曲線亦顯示低值。微電極曲線一般在砂巖層幅值高,并出現(xiàn)正幅差。而泥巖的幅度和幅差均較低,當(dāng)井眼條件不好時(shí),可能會(huì)出現(xiàn)曲線跳動(dòng)現(xiàn)象。砂巖中含灰質(zhì)較多的夾層,因?yàn)橹旅?/p>
5、電阻率異常高,幅度差很小或沒有。一般幅度差的大小標(biāo)明了儲(chǔ)集層滲透性的好壞。普通電阻率測(cè)井曲線在泥巖處顯示為低值。砂巖處顯示為高值,含油砂巖幅值就更高,如有兩條探測(cè)深度不同的Ra 曲線,幅值的差別顯示著低侵、高侵。通常在油層上為低侵,水層上為高侵。井徑在泥巖層擴(kuò)大,砂巖層縮?。孕∮阢@頭直徑)。具體特征總結(jié)見表1表1 砂泥巖剖面測(cè)井曲線特征 地層測(cè)井曲線儲(chǔ)集層砂巖非儲(chǔ)集層泥巖自然電位負(fù)異常(RwRmf)泥巖基線自然伽馬低高井徑縮徑擴(kuò)徑深中淺電阻率高阻低阻聲波300us/m釷低高鈾低高鉀低高(2)碳酸鹽巖剖面的測(cè)井曲線特征碳酸鹽巖剖面的測(cè)井解釋任務(wù),就是從致密的圍巖中找出孔隙性、裂縫性的儲(chǔ)集層,并
6、判斷其含油性。碳酸鹽巖剖面電阻率一般較高,自然電位效果不好。為了區(qū)分巖性和劃分儲(chǔ)層,一般使用自然伽馬測(cè)井曲線。 儲(chǔ)集層相對(duì)于致密的圍巖具有低阻、低自然伽馬以及孔隙度測(cè)井反映孔隙度較大的特點(diǎn)。3 劃分儲(chǔ)集層的基本方法與原則基本要求:凡一切可能含油氣的地層都要?jiǎng)澇鰜?,要適當(dāng)?shù)貏澐置黠@的水層。具體要求為:(1) 估計(jì)為油層、氣層、油水同層和含油水層的儲(chǔ)集層都必須分層解釋。(2) 厚度半米以上的電性(測(cè)井曲線)可疑層(即指從測(cè)井曲線上看有油氣的地層)或錄井顯示為微含油級(jí)別以上的儲(chǔ)集層必須分出。(3) 選擇出作為確定地層水電阻率Rw 的標(biāo)準(zhǔn)水層(厚度大、巖性純、不含油)要?jiǎng)澐殖鰜?。?) 錄井、氣測(cè)有大
7、段油氣顯示而測(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è)井曲線圖上,并根據(jù)鄰井的測(cè)井和試油等資料對(duì)本井的油水關(guān)系作出初步估計(jì)。分層時(shí)應(yīng)注意:l 確定分層的界面深度時(shí),應(yīng)左右環(huán)顧,照顧到分層線對(duì)每條測(cè)井曲線的合理性
8、。l 分層的深度誤差不應(yīng)大于0.1m。l 滲透層中,凡是0.5m以上的非滲透性夾層(泥巖或致密層),應(yīng)將夾層上下的滲透層分兩層解釋。l 巖性漸變層頂界(頂部漸變層)或底界(底部漸變層)分層深度應(yīng)在巖性漸變結(jié)束處。l 一個(gè)厚度較大的滲透層,如有兩個(gè)以上解釋結(jié)論,應(yīng)按解釋結(jié)論分層。l 在同一解釋井段,如果油氣層與水層巖性、地層結(jié)構(gòu)和孔隙度基本相同,則油氣層是純水層的電阻率的3-5倍。純水層的自然電位異常最大,油氣層異常明顯偏小,油水同層介于油、水層之間。并且厚度較大的油水同層,自上而下電阻率有明顯減小的趨勢(shì)。(2)碳酸鹽巖剖面碳酸巖鹽剖面劃分滲透性地層的5 測(cè)井曲線讀數(shù)分層以后,要從有關(guān)的主要測(cè)井
9、曲線將代表該儲(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)部分的極大值的平均值。另外注意孔隙度與電阻率測(cè)井曲線對(duì)應(yīng)取值的原則。因?yàn)橐脙烧呓Y(jié)合計(jì)算地層的含水飽和度,兩者當(dāng)然應(yīng)該是對(duì)應(yīng)深度上同一地層或同一取值區(qū)的讀數(shù)。巖層含油性的定性判斷,主要依據(jù)井曲線
10、的測(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ì)算,常用的公式:孔隙度:含水飽和度:上式中為當(dāng)前層的聲波時(shí)差,為地層水的聲波時(shí)差,189us/ft(623us/m),為固體骨架的聲波時(shí)差,對(duì)于砂巖骨架,主要礦物為石英
11、,其聲波時(shí)差為55.5us/ft。a是常數(shù),對(duì)于砂巖地層通常取1.0,為當(dāng)前層的電阻率,m為膠結(jié)指數(shù)。五、 實(shí)驗(yàn)結(jié)果處理結(jié)果如圖2所示 具體分層如下表所示:儲(chǔ)層序號(hào)頂部深度(m)底部深度(m)厚度(m)測(cè)井曲線讀數(shù)孔隙度含水飽和度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的自然電位的曲線上,有五個(gè)明顯的偏離泥基線的負(fù)異常,故可以劃分出五個(gè)儲(chǔ)集層。然后就可以求出泥質(zhì)含量,進(jìn)而求出泥質(zhì)校正后的孔隙度,然后根據(jù)阿爾奇公式,就可求出含水飽和度、滲透率,進(jìn)而劃分油氣水層。 圖1 某井的綜合測(cè)井曲線圖23451圖2 SN183井測(cè)井解釋綜合圖實(shí)驗(yàn)二利用綜合方法估計(jì)地層泥質(zhì)含量一、 實(shí)驗(yàn)?zāi)康模和ㄟ^實(shí)際計(jì)算,鞏固掌握利用多種測(cè)井資料確定泥質(zhì)含量的方法。二、 實(shí)驗(yàn)要求自編程序,在計(jì)算機(jī)上運(yùn)算出地層泥質(zhì)含量。三、 實(shí)驗(yàn)場(chǎng)地、用具與設(shè)備計(jì)算中心
13、,尺子、像皮和計(jì)算機(jī);四、 實(shí)驗(yàn)內(nèi)容: 什么是泥質(zhì)含量:泥質(zhì)是指顆粒直徑小于0.01mm的碎屑物質(zhì),泥質(zhì)含量,也叫做泥質(zhì)體積,是指泥質(zhì)的體積占巖石總體積的比:確定Vsh的重要性泥質(zhì)含量的確定,在泥質(zhì)砂巖儲(chǔ)集層的定量解釋中具有重要意義。多年來人們提出許多計(jì)算泥質(zhì)含量的理論和方法。目前求取泥質(zhì)含量的方法大致可分為兩類,一類是用每種測(cè)點(diǎn)各求出一個(gè)泥質(zhì)含量,然后求出最佳值。當(dāng)巖石含有泥質(zhì)時(shí),各種測(cè)井曲線均或多或少地受到泥質(zhì)的影響,其影響的程度受Vsh的決定,評(píng)價(jià)巖石的特性時(shí),只有已知Vsh,才知道由于泥質(zhì)帶來的影響,從而將泥質(zhì)的影響校正掉。一般而言,用自然伽馬或自然伽馬能譜或自然電位來求取泥質(zhì)含量效果
14、最好,但自然伽馬要求儲(chǔ)層中除了泥質(zhì)外,其他物質(zhì)不含放射性礦物。自然電位要求地層水電阻率保持不變,且儲(chǔ)層中的泥質(zhì)與相鄰泥巖的的成分相同。用其他方法計(jì)算泥質(zhì)含量則要求更為苛刻的條件:如電阻率方法要求儲(chǔ)層的孔隙度和含水飽和度均要很小。中子和聲波方法則要求孔隙度很小。3. 確定Vsh的方法:() 自然伽瑪法式中,分別是砂巖和泥巖層的自然伽馬值,GCUR是與地層有關(guān)的經(jīng)驗(yàn)系數(shù),新地層(第三系地層)GCUR=3.7,老地層GCUR=2.0.() 自然電位法式中,是當(dāng)前層的自然電位讀數(shù),和分別是純地層和泥質(zhì)地層的自然電位讀數(shù)()電阻率(b=1.5)()中子法式中,是當(dāng)前層的視中子孔隙度讀數(shù),是泥巖層的視中子
15、孔隙度讀數(shù)。() 交會(huì)圖法以中子密度測(cè)井交會(huì)圖為例,通過對(duì)圖2所示的石英點(diǎn)(Q)、水點(diǎn)(W)和泥巖點(diǎn)(SH)構(gòu)成的三角形進(jìn)行分解,依據(jù)資料點(diǎn)所落入三角形中的位置,可以推測(cè)出來泥質(zhì)含量?;蛘呃孟率竭M(jìn)行計(jì)算(依據(jù)點(diǎn)到直線的距離計(jì)算方法):式中,SHWQ0是石英點(diǎn)(Q)和水點(diǎn)(W)連線的直線方程。依據(jù)任意兩點(diǎn)的直線,用石英點(diǎn)()和水點(diǎn)()兩個(gè)點(diǎn)的參數(shù)可以推出:,當(dāng)然,也可以用中子聲波、聲波密度交會(huì)圖的類似方法求Vsh 。如果解釋層段上沒有純泥巖層時(shí),上述交會(huì)圖法所定出的泥巖點(diǎn)位置并不代表實(shí)際的泥巖參數(shù),導(dǎo)致用交會(huì)圖方法估計(jì)的泥質(zhì)含量比實(shí)際值偏大。上述幾種方法計(jì)算的泥質(zhì)含量往往都有一定的條件,當(dāng)條件
16、滿足時(shí), 泥質(zhì)含量的近似結(jié)果,當(dāng)條件不滿足時(shí),計(jì)算的泥質(zhì)含量均可能偏高,所以在實(shí)際處理時(shí),最后選取其中的最小值作為接近實(shí)際的泥質(zhì)含量。五、 實(shí)驗(yàn)結(jié)果1.已知條件:GRmin=50 GRmax=110 SPmin=0 SPmax=100 b=1.5 Rsh=5 2.測(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.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.計(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.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附錄:程序源代碼! 功能:利用測(cè)井?dāng)?shù)據(jù)求泥質(zhì)含量! 時(shí)間:2016/12/22!-! 參數(shù)說明:! 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ì)含量文件! file_
21、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 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!-! ! 功能:讀取參數(shù)文件! 輸入?yún)?shù)說明:! cmdfile:參數(shù)文件! 輸出參數(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ì)含量文件
26、! file_gr:輸出自然伽馬估計(jì)的泥質(zhì)含量文件! file_n_d:輸出中子-密度交會(huì)圖估計(jì)的泥質(zhì)含量文件!-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!-! ! 功能:讀取測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)! 輸入?yún)?shù)說明:! filename:輸入測(cè)井?dāng)?shù)據(jù)
28、文件! 輸出參數(shù)說明:! number:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)!-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!-! ! 功能:讀取測(cè)井?dāng)?shù)據(jù)!
29、輸入?yún)?shù)說明:! file_data:輸入測(cè)井?dāng)?shù)據(jù)文件! n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)! 輸出參數(shù)說明:! 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),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!-! ! 功能:讀取常數(shù)參量! 輸入?yún)?shù)說明:! file_contant:常數(shù)參量文件! 輸出參數(shù)說明:! grmin,grmax:自然伽馬的最小值與最大值! spmin,spmax:自然電位的最小值與最大值! b,r_sh:已
31、知常數(shù)! fai_nsh:已知常數(shù)! p_sh,p_ma,p_f:已知常數(shù)! fai_nf,fai_nma:已知常數(shù)!-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!-! ! 功能:用自然電位的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說明:! n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)! sp(:):自然電位! spmin,spmax:自然電位的最小值與最大值! 輸出參數(shù)說明:! filename:輸出自然電位估計(jì)的泥質(zhì)含量文
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!-! ! 功能:用電阻率的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說明:!
34、n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)! r(:):電阻率! b,r_sh:已知常數(shù)! 輸出參數(shù)說明:! 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,1 v_sh=(r_sh/r(i)*(1.0/b) write(10,*) v_shenddoclose(10)end subro
35、utine get_r!-! ! 功能:用中子的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說明:! n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)! nphi(:):補(bǔ)償中子! fai_nsh:已知常數(shù)! 輸出參數(shù)說明:! filename:輸出中子估計(jì)的泥質(zhì)含量文件!-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!-! ! 功能:用中子的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說明:! n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)! gr(:):自然伽馬! grmin,grmax:自然伽馬的最小值與最大值! 輸出參數(shù)說明:! filename:輸出自然伽馬估計(jì)的泥質(zhì)含量文件!-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!-! ! 功能:用中子-密度交會(huì)圖的估計(jì)泥質(zhì)含量! 輸入?yún)?shù)說明:! n:測(cè)井?dāng)?shù)據(jù)的個(gè)數(shù)! nphi(:):補(bǔ)償中子! rhob(:):巖性密度! fai_nsh:已知常數(shù)! p_sh,p_ma,p_f:已知
38、常數(shù)! fai_nf,fai_nma:已知常數(shù)! 輸出參數(shù)說明:! 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_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實(shí)驗(yàn)三 含泥質(zhì)復(fù)雜巖性地層綜合測(cè)井處理一、實(shí)驗(yàn)?zāi)康模壕C合所學(xué)習(xí)的測(cè)井解釋與數(shù)據(jù)處理知識(shí),掌握一般含泥質(zhì)復(fù)雜巖性地層測(cè)井?dāng)?shù)據(jù)處理的步驟和方法。二、實(shí)驗(yàn)要求用C語言編程,調(diào)試通過,并能計(jì)算出正確的結(jié)果。三、實(shí)驗(yàn)場(chǎng)地、用具與設(shè)備計(jì)算中心,計(jì)
40、算機(jī);四、實(shí)驗(yàn)內(nèi)容:(一) 原理:測(cè)井?dāng)?shù)據(jù)處理與解釋的主要任務(wù)就是依據(jù)“四性”關(guān)系,用測(cè)井所獲得的地層的各種“物理參數(shù)”,定量地計(jì)算出地下地層的巖性()、儲(chǔ)集性(孔隙度和滲透率)以及含油性(飽和度Sw)。當(dāng)儲(chǔ)層含有泥質(zhì)時(shí),還要求取泥質(zhì)含量,對(duì)測(cè)井?dāng)?shù)據(jù)進(jìn)行泥質(zhì)校正。對(duì)于含有泥質(zhì)骨架由兩種礦物構(gòu)成,且礦物含量隨深度變化的復(fù)雜巖性地層,其體積模型如圖3所示。這種泥質(zhì)儲(chǔ)層測(cè)井處理的過程見圖4所示。其基本思路就是用幾種方法分別計(jì)算出獨(dú)立的泥質(zhì)含量Vsh ,然后取最極小值作為該深度的泥質(zhì)含量。對(duì)測(cè)井?dāng)?shù)據(jù)進(jìn)行泥質(zhì)含量校正,兩種孔隙度組合求出孔隙度和兩種礦物的含量,在恢復(fù)為含泥質(zhì)地層,計(jì)算出實(shí)際地層的孔隙度和
41、礦物含量。在計(jì)算過程中所設(shè)計(jì)的參量和模型的變化見圖3所示。VshV1V2恢復(fù)為泥質(zhì)巖石求Sw K,VshNwWV1V2 VshV1V2Vma1Vma2純巖 石 圖3 體積模型與計(jì)算過程中的變化 實(shí)際上求上述四種體積含量的方法有兩種:() 三種孔隙度測(cè)井聯(lián)立解方程(矩陣方法,計(jì)算的核心是求逆陣)接出,Vsh , V1 , V2 (2) 獨(dú)立方法計(jì)算泥質(zhì)含量,二種孔隙度測(cè)井組合,解出純巖石模型的各個(gè)部分的體積。讀入一個(gè)深度的測(cè)井?dāng)?shù)據(jù)孔隙度,礦物1和礦物2構(gòu)成的體積模型泥質(zhì)含量Vsh GR 取最小值作為泥質(zhì)含量Vsh SPRtTHCNL-DEN泥質(zhì)含量Vsh泥質(zhì)含量Vsh泥質(zhì)含量Vsh泥質(zhì)含量Vsh
42、 中子、密度測(cè)井曲線泥質(zhì)校正Vsh計(jì)算含水飽和度Sw矩陣方法計(jì)算純地層孔隙度,礦物1和礦物2的含量計(jì)算滲透率K讀入測(cè)井解釋參數(shù)最后一個(gè)深度點(diǎn)結(jié)束是否圖3 計(jì)算步驟圖先求出泥質(zhì)含量,再校正泥質(zhì)的影響后的測(cè)井?dāng)?shù)據(jù)聯(lián)立求解, Vma1 , Vma2 ,相當(dāng)按純巖石計(jì)算,然后再根據(jù)泥質(zhì)的多少換算出實(shí)際巖石的孔隙度和礦物含量。, (二)計(jì)算過程。() 根據(jù)GR,SP求Vsh 當(dāng)Vsh 50%時(shí),按泥巖處理 =0,Sw =1.0, Vsh =100, V1 =0, V2 =0.K=0當(dāng)Vsh 50%時(shí),按下述方法處理。(2)進(jìn)行泥質(zhì)校正 (3)求純巖石的孔隙度和礦物含量用矩陣表示為:CX=L其中:C= 則孔隙度和兩種礦物的含量可以用解常數(shù)陣的逆的方法得出:X=C-1L(3)恢復(fù)成泥質(zhì)巖石() 計(jì)算含水飽和度(混合泥質(zhì)公式)利用1963年Simandoux的(混合泥質(zhì)公式)計(jì)算含水飽和度??梢钥闯墒顷P(guān)于Sw的一元二次方程。對(duì)于A、B、C分別上述一元二次方程的系數(shù)。當(dāng)a=1.0,m=2.0時(shí)
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 《基于POE的歷史街區(qū)保護(hù)與更新研究》
- 家庭聚餐的健康菜譜設(shè)計(jì)
- 《基于CDIO的SCS創(chuàng)客教學(xué)法在高中信息技術(shù)課程中應(yīng)用研究》
- 知識(shí)付費(fèi)時(shí)代的匯報(bào)制作新趨勢(shì)
- 藝術(shù)行業(yè)創(chuàng)業(yè)策劃書模板
- 《《“藝術(shù)東岸”上海東岸雕塑設(shè)計(jì)方案全球征集》翻譯實(shí)習(xí)報(bào)告》
- 小班生活體驗(yàn)室課程設(shè)計(jì)
- 零售業(yè)資金管理優(yōu)化方案
- 《AP飼料公司存貨內(nèi)部控制問題研究》
- 工程數(shù)學(xué)實(shí)驗(yàn)課程設(shè)計(jì)
- DTU配網(wǎng)自動(dòng)化測(cè)控終端精講
- 道路運(yùn)輸達(dá)標(biāo)車輛客車貨車核查記錄表
- 兒童詩兒童詩的欣賞和創(chuàng)作(課件)
- 人力資源管理工作思路(共3頁)
- 五筆常用字根表3746
- 新生兒肺氣漏
- 氣管切開(一次性氣切導(dǎo)管)護(hù)理評(píng)分標(biāo)準(zhǔn)
- 保安工作日志表
- 姜太公釣魚的歷史故事
- 數(shù)控車床實(shí)訓(xùn)圖紙國(guó)際象棋圖紙全套
- 電子政務(wù)概論教案
評(píng)論
0/150
提交評(píng)論