測(cè)井資料處理上課

1、會(huì)計(jì)學(xué)1測(cè)井資料處理上課測(cè)井資料處理上課245測(cè)井?dāng)?shù)據(jù)的處理與解釋測(cè)井?dāng)?shù)據(jù)的處理與解釋 地質(zhì)信息地質(zhì)信息（工程狀況）（工程狀況） 測(cè)井信息測(cè)井信息（物理參數(shù)）（物理參數(shù)） 采集采集處理與解釋處理與解釋678910n處理過程的主要工作在于n資料的預(yù)處理、測(cè)井繪圖和利用程序計(jì)算儲(chǔ)層參數(shù)11 從第一條測(cè)井曲線出現(xiàn)，相應(yīng)的測(cè)井解釋從第一條測(cè)井曲線出現(xiàn)，相應(yīng)的測(cè)井解釋技術(shù)誕生。最初是簡單的定性解釋技術(shù)誕生。最初是簡單的定性解釋-“-“相面相面法法”-根據(jù)測(cè)井曲線的形態(tài)判斷油、水層。根據(jù)測(cè)井曲線的形態(tài)判斷油、水層。伴隨著測(cè)井技術(shù)及其它相關(guān)技術(shù)的發(fā)展，測(cè)井伴隨著測(cè)井技術(shù)及其它相關(guān)技術(shù)的發(fā)展，測(cè)井解釋也在發(fā)展

2、：解釋也在發(fā)展：定性定性 定量定量單井單井 多井多井單學(xué)科單學(xué)科 與其它學(xué)科結(jié)合與其它學(xué)科結(jié)合12 模擬計(jì)錄模擬計(jì)錄定性解釋定性解釋 數(shù)字測(cè)井?dāng)?shù)字測(cè)井定量處理：計(jì)算孔、滲、飽，定量處理：計(jì)算孔、滲、飽， 識(shí)別巖性識(shí)別巖性 數(shù)控測(cè)井?dāng)?shù)控測(cè)井采集更多地質(zhì)信息，提高了評(píng)采集更多地質(zhì)信息，提高了評(píng) 價(jià)油氣層、解決地質(zhì)問題的能力價(jià)油氣層、解決地質(zhì)問題的能力 成像測(cè)井成像測(cè)井獲取更豐富地質(zhì)信息，拓展獲取更豐富地質(zhì)信息，拓展 測(cè)井應(yīng)用領(lǐng)域，可開展深層次測(cè)井應(yīng)用領(lǐng)域，可開展深層次 地質(zhì)應(yīng)用研究地質(zhì)應(yīng)用研究1314衛(wèi)星傳送各種磁帶格式轉(zhuǎn)換測(cè)井曲線數(shù)字化表格數(shù)據(jù)錄入數(shù)據(jù)庫數(shù)據(jù)庫資料預(yù)處理單井處理與解釋多井處理與解

3、釋服務(wù)程序：交會(huì)圖、直方圖、數(shù)據(jù)統(tǒng)計(jì)分析、數(shù)據(jù)庫解釋成果的顯示與輸出、成果歸檔、圖形顯示、硬件拷貝、繪圖儀輸出數(shù)據(jù)輸入數(shù)據(jù)輸入15用計(jì)算機(jī)處理和解釋測(cè)井?dāng)?shù)據(jù)所做的各項(xiàng)工作主要包括數(shù)據(jù)管理數(shù)據(jù)管理、測(cè)井?dāng)?shù)據(jù)預(yù)處理測(cè)井?dāng)?shù)據(jù)預(yù)處理和測(cè)井?dāng)?shù)據(jù)地測(cè)井?dāng)?shù)據(jù)地質(zhì)分析質(zhì)分析三個(gè)方面。數(shù)據(jù)管理包括：數(shù)據(jù)輸入、成果輸出其中測(cè)井?dāng)?shù)據(jù)預(yù)處理主要包含以下幾項(xiàng)：其中測(cè)井?dāng)?shù)據(jù)預(yù)處理主要包含以下幾項(xiàng)：深度對(duì)齊、斜井曲線校正、平滑濾波處理、環(huán)境校正、數(shù)據(jù)標(biāo)準(zhǔn)化、解釋模型和解釋參數(shù)的確定。測(cè)井?dāng)?shù)據(jù)分析：測(cè)井?dāng)?shù)據(jù)分析：按照確定的解釋模型，選用相應(yīng)的測(cè)井分析程序，在計(jì)算機(jī)上用測(cè)井?dāng)?shù)據(jù)計(jì)算出各種地質(zhì)或工程參數(shù)，并用形象直觀的圖標(biāo)顯示出

4、來。 162、測(cè)井?dāng)?shù)據(jù)處理系統(tǒng)171819202122232425一、測(cè)井?dāng)?shù)據(jù)格式一、測(cè)井?dāng)?shù)據(jù)格式2627頭記錄塊( General Heading) 用于記錄有關(guān)測(cè)井的輔助性數(shù)據(jù)和信息,其數(shù)據(jù)在頭記錄塊中應(yīng)有固定的位置和固定的格式,頭記錄塊的長度為276 字節(jié)(1 ) 測(cè)井成果代碼( ECC) :4 字節(jié)浮點(diǎn)數(shù)。2829 數(shù)據(jù)記錄塊(Data Record) 用于記錄測(cè)井各曲線的采樣或計(jì)算值。每個(gè)數(shù)據(jù)記錄塊包含有NOLOG條曲線的SPCPR 個(gè)樣點(diǎn),其記錄長度為: 數(shù)據(jù)記錄塊長度= NOLOG x SPCPR x 4(字節(jié)) 其在記錄塊中的記錄順序如圖2 :30第一個(gè)數(shù)據(jù)記錄塊的第一個(gè)采樣點(diǎn)

5、的深度為STDEP, 下一個(gè)采樣點(diǎn)的深度為STDEP + / -RL EV(取決于測(cè)井方向是上提或下放) 。3132333435測(cè)井解釋程序數(shù)據(jù)塊子程序（1個(gè)）主程序調(diào)用用戶子程序（若干）調(diào)用系統(tǒng)服務(wù)子程序系統(tǒng)服務(wù)子程序（若干）3637文件名起始深度1，終止深度1參數(shù)名11=參數(shù)值11，參數(shù)名12=參數(shù)值12，參數(shù)名1n=參數(shù)值1n起始深度n，終止深度n參數(shù)名n1=參數(shù)值n1，參數(shù)名n2=參數(shù)值n2，參數(shù)名nn=參數(shù)值nn解釋井段解釋井段1解釋井段解釋井段238開始系統(tǒng)服務(wù)程序（RDFLNM/IN/OUT/CONST）讀測(cè)井?dāng)?shù)據(jù)/IN處理完？作某些校正處理輸出/OUTNo結(jié)束Yes4041的曲

6、線段，并把對(duì)應(yīng)點(diǎn)連校深線進(jìn)行校正。而自動(dòng)校深是應(yīng)用相關(guān)對(duì)比法，自動(dòng)求出校深線從而進(jìn)行校正。 為什么要作深度校正？為什么要作深度校正？421、利用深度控制曲線進(jìn)行深度校正通常是用GR曲線作為深度控制曲線，每次測(cè)量都帶測(cè)一條GR曲線，并以某次測(cè)量的GR曲線為基準(zhǔn)，把其他次測(cè)量的曲線深度向它對(duì)齊。43GR0GR1聲波時(shí)差CNL4445NiNiNixyyiyNxixNyiyxixN12121)(1)(1)()(146X(i)Y(i)474849505152h0hih0-hi標(biāo)準(zhǔn)曲線對(duì)比曲線53C1C2d11d12d21d2254212221111211ddddddddyx)(212211121121d

7、dddddddxy55C1C2d11d12d21d22dxdy56)(11iyiiiiiyddddVVVViiiiiyiyddVVddVV1157C1C2d11d12d21d22dxdyii+1ViVi+15859H0h060H0h0H1h101A61H0h0H1h101A621212hhdhdH0h0H1h101AB h22H2dhdHdhhdh1212cosdhdH )sin(sin1212121221hhdHHHhh636465Ti-1TiTi+1ZZ06621110)(taaTQti00aQ01aQ)(31110iiiiTTTaT即位Ti點(diǎn)的滑動(dòng)平均值。67Ti-1TiTi+1Z0Z=

8、a0+a1t+a2t2Z682222210)(ttitataaTQ00aQ01aQ02aQ同樣利用最小二乘法，即同樣利用最小二乘法，即iiiiiiTTTTTaT17)(12)( 33511122069102030405060701.31.41.51.61.71.81.92.0 Y Axis TitleX Axis Title B102030405060701.31.41.51.61.71.81.92.0 Y Axis TitleX Axis Title B 3 point AA Smoothing of Data2_B線性函數(shù)平滑70102030405060701.31.41.51.61.71

9、.81.92.0 Y Axis TitleX Axis Title B102030405060701.31.41.51.61.71.81.92.0 Y Axis TitleX Axis Title B 5 point AA Smoothing of Data2_B二次函數(shù)平滑71720.00.51.01.52.02.53.03.54.0-0.2-0.10.00.10.20.30.40.50.60.70.8 校正值/g/cm3間隙/cm 密度差值為0.9 密度差值為1.573TNPHO=-0.0333R3+1.4000R2+4.1667CAL+12.0000TNPH5=-0.0008R3+0.0

10、613R2-0.5076CAL+3.6284TNPH10=-0.0029R3+0.1512R2-1.4378CAL+10.4400TNPH15=-0.0014R3+0.0956R2-0.3848CAL+8.2795優(yōu)點(diǎn)：簡單、全面、迅速、有效。7475767778表示在該解釋井段上，滿足N=15，b=2.42的采樣點(diǎn)有2個(gè)7980表示在該解釋井段上，滿足N=15，b=2.42對(duì)應(yīng)的2個(gè)采樣點(diǎn)處的GR的平均級(jí)別是210GMINGMAXGMINLOGZ81828384858687單孔隙度測(cè)井泥質(zhì)砂巖處理程序，單孔隙度測(cè)井泥質(zhì)砂巖處理程序，Atlas公司商用軟件模塊。公司商用軟件模塊。i曲線編號(hào)；S

11、HLGi第i種測(cè)井曲線值；GMINi第i條曲線純砂巖數(shù)值；GMAXi第i條曲線純砂巖數(shù)值；Shi第i條曲線計(jì)算的泥巖相對(duì)值；GCUR經(jīng)驗(yàn)系數(shù)，公式（11-18b）；Vshi第i種曲線計(jì)算的泥質(zhì)含量 Vsh輸出結(jié)果：輸出結(jié)果：Vsh=min(Vshi) Vsh=if (Vsh=SHCT, 1, Vsh)88 shshDfmashmashfmabmaVV shshpsmafmashshpmafmaVCttttVCtttt1 NshshNV 89POR程序采用一種孔隙度測(cè)井方法，按含水泥質(zhì)砂巖體積模型公式計(jì)算POR，進(jìn)行了泥質(zhì)校正，未作油氣校正。用DEN計(jì)算POR = (DEN-DMA)/(DF-D

12、MA)-VSH*(DSH-DMA)/(DF-DMA)當(dāng)DSH = DG時(shí)，不作泥質(zhì)校正。用AC計(jì)算計(jì)算CP：如果CPFG=0，用固定壓實(shí)系數(shù)CP = CP 如果CPFG=1，用TSH計(jì)算壓實(shí)系數(shù)TSH單位為us/ft, CP=TSH/100TSH單位為us/m, CP=0.3048TSH/100 如果CPFG=2，用經(jīng)驗(yàn)公式計(jì)算壓實(shí)系數(shù)CP = BCP-ACP*DEP*0.001CP限制條件：如CP1 or Sw0.7時(shí)，無需迭代（說明地層中無油氣或油氣含量低）；4）當(dāng)Sw=0.7時(shí)，做油氣校正； a) 重新計(jì)算Shr0.5（1Sw）； b）重新計(jì)算 ； c）對(duì)中子、密度做油氣校正 d）當(dāng)相鄰

13、兩次迭代的 小于給定的數(shù)值（如 ）時(shí)， 終止迭代；否則轉(zhuǎn)入第2步，重新計(jì)算，重復(fù)上述迭代運(yùn)算。 101103含水泥質(zhì)雙礦物地層處理方法含油氣泥質(zhì)雙礦物地層處理方法1041051、計(jì)算泥質(zhì)含量VshiiiiGMINGMAXGMINSHLGSH106107CNL/%密度b0204060801003.02.01.0水點(diǎn)C1C3C4C2123C1石英石英C2方解石方解石C3白云石白云石C4硬石膏硬石膏108109211ccshVVV2211ccccfshshbVVV2211NccNccNfNshshCNLVVV這是個(gè)3個(gè)方程，4個(gè)未知數(shù)的方程組，怎么解？ 110AC=AC-(0.8*TSH*SH)/(

14、1-0.8*SH)shNshshshshNbNCcbcVtVtt1)/()/()/(DEN=DEN-(0.8*DSH*SH)/(1-0.8*SH)CNL=CNL-(0.8*NSH*SH)/(1-0.8*SH)111則響應(yīng)方程組變?yōu)椋?11ccVV2211ccccfbVV2211NccNccNfCNLVV為了使用統(tǒng)一的計(jì)算方法進(jìn)行求解，用一個(gè)通式表示上述方程組，為：y=V1Y1+V2Y2+V3Y3x=V1X1+V2X2+V3X31=V1+V2+V3112上面方程組解得V1=A1X+B1Y+C1V2=A2X+B2Y+C2V3=1-V1-V2A1、B1、C1和A2、B2、C2成為交會(huì)三角形系數(shù)，得1

15、321DXXB322311XXYYBAC1=-(A1X2+B1Y2)D1=(X2-X3)(Y1-Y2)-(X1-X2)(Y2-Y3)X、Y就是要求的那個(gè)采樣點(diǎn)的縱橫坐標(biāo)1132312DXXB311322XXYYBAC2=-(A2X1+B2Y1)D2=(X1-X3)(Y1-Y2)-(X1-X2)(Y2-Y3)114V1(I)=A1(I)X+B1(I)Y+C1(I)V2(I)=A2(I)X+B2(I)Y+C2(I)V3(I)=100-V1(I)-V2(I)從第1個(gè)三角形開始，對(duì)解釋結(jié)果進(jìn)行挑選：115CNL/%密度b0204060801003.02.01.0水點(diǎn)C1C3C4C2123C1石英石英C

16、2方解石方解石C3白云石白云石C4硬石膏硬石膏P1）當(dāng)V1(1) 、 V2(1) 、V3(1)均大于0，資料點(diǎn)落在什么地方？此時(shí)，此時(shí)， V1(1) 為孔隙度，為孔隙度， V2(1)為石英體積，為石英體積，V3(1)為方解石體積為方解石體積116CNL/%密度b0204060801003.02.01.0水點(diǎn)C1C3C4C2123C1石英石英C2方解石方解石C3白云石白云石C4硬石膏硬石膏2）當(dāng)V1(1) 0，說明資料點(diǎn)落在什么地方？P V1(1) 為孔隙度，為孔隙度， V2(1)為石英體積，為石英體積，V3(1)為方為方解石體積解石體積117CNL/%密度b0204060801003.02.0

17、1.0水點(diǎn)C1C3C4C2123C1石英石英C2方解石方解石C3白云石白云石C4硬石膏硬石膏3）當(dāng)V3(1) 0，說明資料點(diǎn)落在什么地方？PP V1(1) 為孔隙度，為孔隙度， V2(1)為石英體積，為石英體積，V3(1)為方為方解石體積解石體積118CNL/%密度b0204060801003.02.01.0水點(diǎn)C1C3C4C21234）當(dāng)V2(1) 0，說明資料點(diǎn)落在什么地方？P V1(1) 為孔隙度，為孔隙度， V2(1)為石英體積，為石英體積，V3(1)為方為方解石體積解石體積C1石英石英C2方解石方解石C3白云石白云石C4硬石膏硬石膏119CNL/%密度b0204060801003.0

18、2.01.0水點(diǎn)C1C3C4C21235）當(dāng)V2(3) 0.1，m=2.1。m最大值為4。3）仍為Archie公式，但a、n、m需要由解釋人員輸入122原狀地層含水孔隙度：PORw=POR*Sw沖洗帶含水孔隙度：PORxo=POR*Sxo123前面用中子密度交會(huì)圖求出的孔隙度是地層的總有效孔隙度，而聲波測(cè)井資料求出的孔隙度是巖石的原生孔隙度，這樣可由這兩種孔隙度的差值來求取縫洞孔隙度：pmafmafCtttt112415. 01SHfbga巖石的視顆粒密度和視顆粒時(shí)差是包含泥質(zhì)在內(nèi)的全部顆粒的平均密度和平均時(shí)差。其作用是判別巖性。15. 01SHtttfga12524 . 1136. 0wir

19、rSK通過標(biāo)示符XKM，用兩種方式計(jì)算滲透率：23wirrSXKMKXKM=0時(shí)XKM0時(shí)孔隙度和飽和度均為百分?jǐn)?shù)孔隙度為小數(shù)，飽和度為百分?jǐn)?shù)XKM常取25000。126miiRLEVPF1 PF和HF表示從某深度開始，累計(jì)得到的純粹孔隙厚度和純粹油氣厚度。在成果圖中，通常在指定位置用一些粗短線表示，每相鄰兩個(gè)短線之間累計(jì)孔隙厚度或累計(jì)油氣厚度為1m或1ft，短線越多，孔隙越發(fā)育，油氣越多。miwiiRLEVSHF1)1 (128129130一、前言一、前言131二、什么是最優(yōu)化二、什么是最優(yōu)化132axxxaxf2)2()(133134njjmiiba11由由aiai到到bjbj的運(yùn)費(fèi)單價(jià)分

20、別已知為的運(yùn)費(fèi)單價(jià)分別已知為Cij(Cij(元元/ /噸噸)(i=1,2,3,)(i=1,2,3,m;j=1,2,3,m;j=1,2,3,n),n)。問由每個(gè)產(chǎn)地。問由每個(gè)產(chǎn)地到每個(gè)銷售地運(yùn)輸量各為多少時(shí)，既保證需要量，到每個(gè)銷售地運(yùn)輸量各為多少時(shí)，既保證需要量，又使總運(yùn)費(fèi)最少？又使總運(yùn)費(fèi)最少？135136njmixnjbxmiaxtsxCijmijijnjiijminjijij,.,2 , 1;,.,2 , 10,.,2 , 1,.,2 , 1. .min1111137(1-Sxo)(1-Sxo)(Sxo-Sw)(Sxo-Sw)SwSwVclVmax多礦物地層體積物理模型示意圖多礦物地層體積

21、物理模型示意圖138hmfmfmfmfmfnimaimaiclclormfbPKATRPTPVVAS15. 1)15. 011. 1 (07. 1)8277. 6(004898. 05000807 . 01)(18. 11其中其中mf 和和P分別是泥漿濾液密度值和礦化度分別是泥漿濾液密度值和礦化度,h 是油氣是油氣密度值密度值,是經(jīng)驗(yàn)系數(shù)是經(jīng)驗(yàn)系數(shù),通常取通常取1. 139140141142143144145146njjjNiiiiiTxGxfaxF1221222)()()(147148149儲(chǔ)層參數(shù)儲(chǔ)層參數(shù)初始值初始值X0響應(yīng)方程響應(yīng)方程理論測(cè)井值理論測(cè)井值實(shí)際測(cè)井值實(shí)際測(cè)井值比比較較充分逼

22、充分逼近否？近否？輸出最優(yōu)化輸出最優(yōu)化解釋結(jié)果解釋結(jié)果X用最優(yōu)化方法調(diào)整未知量用最優(yōu)化方法調(diào)整未知量X，X0XYN150151iiaiia)(*xfiia152mnNmnNxFRinc)()(154155GEOFRAME軟件包中的P包、G包是一套測(cè)井資料處理軟件。它可以根據(jù)區(qū)域地質(zhì)模型，利用測(cè)井資料準(zhǔn)確的計(jì)算出地層的各種參數(shù)，并作出綜合評(píng)價(jià)。156157158159160勘探測(cè)井解釋平臺(tái)版權(quán)所有版權(quán)所有版權(quán)所有中國石油天然氣集團(tuán)公司油氣勘探部中國石油天然氣集團(tuán)公司油氣勘探部中國石油天然氣集團(tuán)公司油氣勘探部測(cè)井軟件開發(fā)項(xiàng)目組測(cè)井軟件開發(fā)項(xiàng)目組測(cè)井軟件開發(fā)項(xiàng)目組162163164165166167

23、168169170171172173174175176DataInDataIn測(cè)井?dāng)?shù)據(jù)解編測(cè)井?dāng)?shù)據(jù)解編DataOutDataOutgudgud數(shù)據(jù)格式轉(zhuǎn)換其它數(shù)據(jù)格式數(shù)據(jù)格式轉(zhuǎn)換其它數(shù)據(jù)格式DataManDataMangudgud數(shù)據(jù)管理數(shù)據(jù)管理DigiLogDigiLog測(cè)井曲線數(shù)字化測(cè)井曲線數(shù)字化TapeReadTapeRead讀磁帶數(shù)據(jù)讀磁帶數(shù)據(jù)SystemFolderSystemFolder工作路徑設(shè)置工作路徑設(shè)置GeologySymbolsGeologySymbols地質(zhì)符號(hào)庫管理地質(zhì)符號(hào)庫管理ShowingSymbolsShowingSymbols油氣顯示符號(hào)管理油氣顯示符號(hào)管理L

24、ithSymbolsLithSymbols用戶常用巖性符號(hào)管理用戶常用巖性符號(hào)管理HelpHelp幫助幫助EZLogEZLog常規(guī)測(cè)井綜合解釋軟件；常規(guī)測(cè)井綜合解釋軟件；LogDipLogDip地層傾角測(cè)井分析軟件；地層傾角測(cè)井分析軟件；LogViewLogView聲電成像測(cè)井分析軟件；聲電成像測(cè)井分析軟件；PoroDistPoroDist聲電成像孔隙分析軟件；聲電成像孔隙分析軟件；MagresonMagreson核磁共振測(cè)井分析軟件核磁共振測(cè)井分析軟件177操作系統(tǒng)：操作系統(tǒng)： 中英文中英文 WINDOW NT/ 2000 / XPWINDOW NT/ 2000 / XP軟件安裝方式：軟件安裝方式： 單機(jī)版和網(wǎng)絡(luò)版單機(jī)版和網(wǎng)絡(luò)版快速設(shè)置工作路徑快速設(shè)置工作路徑中英文語言界面切換中英文語言界面切換178GUDGUD是作為是作為LogVisionLogVision平臺(tái)內(nèi)部數(shù)據(jù)交換格平臺(tái)內(nèi)部數(shù)據(jù)交換格式，因此任何格式的測(cè)井?dāng)?shù)據(jù)首先必須轉(zhuǎn)式，因此任何格式的測(cè)井?dāng)?shù)據(jù)首先必須轉(zhuǎn)換為換為GUDGUD格式，才能共享格式，才能共享LogVisionLogVision平臺(tái)所平臺(tái)所提供的各種特性。提供的各種特性。 179本平臺(tái)為各軟件產(chǎn)品提供統(tǒng)一的數(shù)據(jù)管理、統(tǒng)一的本平臺(tái)為各軟件產(chǎn)品提供統(tǒng)一的數(shù)據(jù)管理、統(tǒng)一的繪圖模板、統(tǒng)一的繪圖文檔和通用處理框架；也就繪圖模板、統(tǒng)一的繪圖文檔和通用處理框架；