測井資料處理上課

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

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

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

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

5、的深度為STDEP, 下一個采樣點的深度為STDEP + / -RL EV(取決于測井方向是上提或下放) 。3132333435測井解釋程序數(shù)據(jù)塊子程序（1個）主程序調(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）讀測井?dāng)?shù)據(jù)/IN處理完？作某些校正處理輸出/OUTNo結(jié)束Yes4041的曲

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

7、dddddddxy55C1C2d11d12d21d22dxdy56)(11iyiiiiiyddddVVVViiiiiyiyddVVddVV1157C1C2d11d12d21d22dxdyii+1ViVi+15859H0h060H0h0H1h101A61H0h0H1h101A621212hhdhdH0h0H1h101AB h22H2dhdHdhhdh1212cosdhdH )sin(sin1212121221hhdHHHhh636465Ti-1TiTi+1ZZ06621110)(taaTQti00aQ01aQ)(31110iiiiTTTaT即位Ti點的滑動平均值。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)點：簡單、全面、迅速、有效。7475767778表示在該解釋井段上，滿足N=15，b=2.42的采樣點有2個7980表示在該解釋井段上，滿足N=15，b=2.42對應(yīng)的2個采樣點處的GR的平均級別是210GMINGMAXGMINLOGZ81828384858687單孔隙度測井泥質(zhì)砂巖處理程序，單孔隙度測井泥質(zhì)砂巖處理程序，Atlas公司商用軟件模塊。公司商用軟件模塊。i曲線編號；S

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

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

13、兩次迭代的 小于給定的數(shù)值（如 ）時， 終止迭代；否則轉(zhuǎn)入第2步，重新計算，重復(fù)上述迭代運算。 101103含水泥質(zhì)雙礦物地層處理方法含油氣泥質(zhì)雙礦物地層處理方法1041051、計算泥質(zhì)含量VshiiiiGMINGMAXGMINSHLGSH106107CNL/%密度b0204060801003.02.01.0水點C1C3C4C2123C1石英石英C2方解石方解石C3白云石白云石C4硬石膏硬石膏108109211ccshVVV2211ccccfshshbVVV2211NccNccNfNshshCNLVVV這是個3個方程，4個未知數(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)一的計算方法進行求解，用一個通式表示上述方程組，為：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成為交會三角形系數(shù)，得1

15、321DXXB322311XXYYBAC1=-(A1X2+B1Y2)D1=(X2-X3)(Y1-Y2)-(X1-X2)(Y2-Y3)X、Y就是要求的那個采樣點的縱橫坐標(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個三角形開始，對解釋結(jié)果進行挑選：115CNL/%密度b0204060801003.02.01.0水點C1C3C4C2123C1石英石英C

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

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

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

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

20、別已知為的運費單價分別已知為Cij(Cij(元元/ /噸噸)(i=1,2,3,)(i=1,2,3,m;j=1,2,3,m;j=1,2,3,n),n)。問由每個產(chǎn)地。問由每個產(chǎn)地到每個銷售地運輸量各為多少時，既保證需要量，到每個銷售地運輸量各為多少時，既保證需要量，又使總運費最少？又使總運費最少？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)驗系數(shù)是經(jīng)驗系數(shù),通常取通常取1. 139140141142143144145146njjjNiiiiiTxGxfaxF1221222)()()(147148149儲層參數(shù)儲層參數(shù)初始值初始值X0響應(yīng)方程響應(yīng)方程理論測井值理論測井值實際測井值實際測井值比比較較充分逼

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

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

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