一種核磁共振測井T₂譜100%純水譜的確定方法
【專利摘要】本發(fā)明公開了一種核磁共振測井T2譜100％純水譜的確定方法，包括：首先根據(jù)已知的巖石物理實驗特征建立勘探地區(qū)儲層束縛水飽和度Swi與T2幾何均值T2g的關(guān)系，然后利用球管薄膜模型逐點計算測井T2譜束縛水飽和度Swi及T2幾何均值T2g，并建立核磁共振離心譜，在核磁共振離心譜基礎(chǔ)上構(gòu)建100％純水核磁共振譜。當(dāng)儲層含油氣時核磁信號與100％純水核磁信號存在顯著的差異，利用這種差異對儲層流體進行準(zhǔn)確的識別。本發(fā)明利用核磁共振測井資料構(gòu)造100％純水譜能夠反映儲層真實的孔隙結(jié)構(gòu)，且通過對比實測核磁共振譜與構(gòu)造純水譜的差異有效識別儲層流體性質(zhì)，回避了由于原油性質(zhì)、觀測模式導(dǎo)致的核磁共振測井流體識別困難的難題，提高了核磁共振測井資料的應(yīng)用，增強了解決地質(zhì)問題的能力。
【專利說明】
-種核磁共振測井T2譜100%純水譜的確定方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明設(shè)及巖屯、核磁共振、核磁共振測井及儲層評價領(lǐng)域，特別設(shè)及一種核磁共 振測井T2譜100 %純水譜的確定方法。
【背景技術(shù)】
[0002] 核磁共振是磁場中原子核對電磁波的一種響應(yīng)。核磁共振測井的原始數(shù)據(jù)是幅度 隨時間衰減的自旋回波串，自旋回波串包含了儲層物性、孔隙類型、孔徑大小、流體類型及 其分布等十分豐富的信息，可W表示為多個衰減指數(shù)的疊加：
[0003]
[0004] 其中:M(t)-一核磁共振測井測量得到的隨時間t變化的回波幅度；
[0005] T2,i-第i個組分的橫向弛豫時間，i = l，. . .，N;
[0006] Mi(0)-對應(yīng)于T2,i組分的核磁共振孔隙度，i = l，...，N。
[0007] 通過適當(dāng)?shù)目潭?，根?jù)自旋回波串信號強度可W得到核磁孔隙度，并通過多指數(shù) 擬合自旋回波串，可W得到每一深度點的7'2,廣--機分布曲線，即核磁共振測井T2譜。該測井 T2譜橫坐標(biāo)為T2時間值，縱坐標(biāo)為孔隙度分量。
[000引核磁共振測井在評價儲層品質(zhì)方面具有獨特的優(yōu)勢，通過解譜得到測井T2譜，可 W直觀的反映孔徑大小及其分布。當(dāng)儲層飽含水時，大孔喉對應(yīng)較大T2值，小孔喉對應(yīng)較小 T2值，從而達到評價孔隙結(jié)構(gòu)的目的。但是當(dāng)儲層含油氣時，計算的核磁孔隙度偏低，測井T2 譜的幅度和峰值位置不再反映儲層孔喉分布，若直接使用測井T2譜進行儲層品質(zhì)的評價必 然是不恰當(dāng)?shù)模鵚很有必要對測井T2譜進行含控校正。
[0009] 目前國內(nèi)研究人員從核磁共振實驗及測井T2譜出發(fā)對儲層含控校正做了相關(guān)研 究工作。
[0010] ①優(yōu)化觀測模式。在原始數(shù)據(jù)采集上，可W采用特殊的觀測模式W便提供孔隙度 校正信息。
[0011] ②數(shù)值模擬方法。參見2016年1月第2期《石油勘探與開發(fā)》雜志中，胡法龍、周燦燦 等著作的《核磁共振測井構(gòu)造水譜流體識別方法及應(yīng)用》的文章，其中記載了數(shù)值模擬方法 的應(yīng)用，利用Iw和Te測量下的T2分布，構(gòu)造 Iw長等待時間和Te下完全含水狀態(tài)下T2分布，從而 達到儲層評價的目的。
[0012] ③其他數(shù)學(xué)統(tǒng)計方法。參見2000年6月第3期《測井技術(shù)》雜志中，高敏等著作的《用 核磁共振測井資料評價儲層的孔隙結(jié)構(gòu)》的文章，其中記載了區(qū)域性T2g與各孔隙結(jié)構(gòu)參數(shù) 關(guān)系，進一步消除儲層含控的影響;參見2012年2月02期《石油天然氣學(xué)報》雜志中，肖飛、何 宗斌等著作的《核磁共振測井連續(xù)表征儲層孔隙結(jié)構(gòu)方法研究》的文章，其中記載了通過測 井T2譜計算孔隙結(jié)構(gòu)等相關(guān)參數(shù)，然后與實驗值對比，最后利用含油飽和度作為校正系數(shù) 消除二者之間的誤差;參見2010年石油技術(shù)交流會上，劉忠華、周燦燦等著作的《一種利用 核磁共振測井資料評價儲層孔隙結(jié)構(gòu)的新方法》的文章，其中記載了核磁=孔隙組分百分 含量與孔隙結(jié)構(gòu)類型的關(guān)系，用w準(zhǔn)確識別高產(chǎn)油氣層，運種方法適用于多相流體存在的 地層。
[0013] W上提出的各種含控校正方法具有各自的優(yōu)勢和局限性。第1種方法的優(yōu)點是可 W從原始數(shù)據(jù)采集上獲得高質(zhì)量的油氣水分離信號，但其應(yīng)用上的局限性是需要慎重進行 測前設(shè)計，但目前來說國內(nèi)測前設(shè)計方面研究較少。方法2的優(yōu)點是可W從根本上得到完全 含水的T2分布，但是算法等復(fù)雜，不易推廣。最后一種方法區(qū)域應(yīng)用性較強，不具有普遍規(guī) 律，沒有從根本上解決含控對譜分布的影響。
[0014] 為了解決核磁測井應(yīng)用中T2譜受含控影響不能很好反映儲層品質(zhì)的實際難題，應(yīng) 該采用全新的研究思路與方法，避開巖屯、實驗或其他信息標(biāo)定的傳統(tǒng)思路與方法，提出確 定核磁共振測井含控校正的一種可靠方法和技術(shù)，克服現(xiàn)行方法的技術(shù)局限性。

【發(fā)明內(nèi)容】

[0015] 本發(fā)明的目的是利用測井T2譜構(gòu)造100%純水譜而實現(xiàn)含控校正，利于后續(xù)流體 識別及儲層品質(zhì)評價等。
[0016] 為實現(xiàn)上述目的，本發(fā)明提供了一種核磁共振測井T2譜100%純水譜的確定方法， 包括W下步驟：
[0017] 步驟(1)通過巖屯、核磁共振實驗獲取實驗核磁共振T2飽和譜和離屯、譜數(shù)據(jù)；
[001引步驟（2)根據(jù)步驟（1)得到的巖屯、核磁共振T2飽和譜和離屯、譜實驗數(shù)據(jù)驗數(shù)據(jù)計 算巖屯、束縛水飽和度Swi及核磁共振T2譜幾何均值T2g并建立二者關(guān)系，采用最小二乘法確定 待定系數(shù)a、b;
[0019] 步驟(3)利用球管薄膜模型計算核磁共振T2譜的巖屯、束縛水飽和度Swi，并構(gòu)建離 屯、譜；
[0020] 步驟(4)根據(jù)步驟(2)確定的待定系數(shù)a、b，計算含水飽和度Sw=l時的核磁共振T2 譜幾何均值T2g;
[0021] 步驟(5)借助正態(tài)分布函數(shù)在核磁共振T2譜上構(gòu)建一個自由水譜；
[0022] 步驟(6)，將步驟(3)構(gòu)建的離屯、譜與步驟(5)構(gòu)建的自由水譜相加得到一個總譜， 利用規(guī)劃求解的方法使得構(gòu)建的總譜孔隙度批和幾何均值T2g/與步驟(3)中計算巖屯、束縛 水飽和度Swi產(chǎn)生的中間參數(shù)測井孔隙度與t及步驟(4)中確定的核磁共振T2譜幾何均值T2g 值誤差最小，此時得到的總譜即為100%純水譜。
[0023] 進一步的，所述步驟（1)中，核磁共測井振T2實驗按照《巖樣核磁共振參數(shù)實驗室 測量規(guī)范SY/T 6490-2007》標(biāo)準(zhǔn)規(guī)定的流程進行。
[0024] 進一步的，所述步驟(2)中巖屯、束縛水飽和度Swi及核磁共振實驗T2譜幾何均值T2g， 具體按照如下公式進行計算：
[0025]
[0026]
[0027]其中：與細相為巖屯、離屯、譜孔隙分量;取幼巖屯、飽和譜孔隙分量;T2,功實驗T2譜的 橫向弛豫時間。
[002引進一步的，所述步驟(2)中，建立巖屯、束縛水飽和度Swi及核磁共振T2譜幾何均值T2g 的關(guān)系如下：
[0029]
[0030] 進一步的，在進行步驟(3)操作時，需對核磁測井資料處理得到標(biāo)準(zhǔn)T2譜，處理方 法對于P型核磁參考《DPP核磁處理技術(shù)手冊》。
[0031] 在進行步驟(3)操作之前，還需利用中子-密度曲線進行測井孔隙度祭t的計算，通 過下式實現(xiàn)：
[0032]
[003；3]其中：與DEJV為測井密度孔隙度;機:亂為測井中子孔隙度。
[0034] 所述測井孔隙度取的計算按照石油天然氣行業(yè)標(biāo)準(zhǔn)《生產(chǎn)測井解釋規(guī)程SY/ T5793-1993》規(guī)定的流程進行。
[0035] 進一步的，所述步驟(3)利用球管薄膜模型計算測井T2譜的束縛水飽和度Swi并構(gòu) 建離屯、譜，是根據(jù)球管薄膜模型確定的權(quán)系數(shù)計算得到：
[0036]
[0037]
[00；3 引
[0039]
[0040] 其中：取swU為離屯、譜分量；辦如為束縛水孔隙度;W功權(quán)系數(shù)分量;m為系數(shù);Swi 為束縛水飽和度;取t為測井孔隙度;滬i為測井T2譜孔隙分量;T2,i為測井T2譜的橫向弛豫時 間。
[0041] 進一步的，所述步驟(5)借助正態(tài)分布函數(shù)在測井T2譜上構(gòu)建一個自由水譜：
[0042]
[0043] 其中：A為構(gòu)建自由水譜的幅度值;U為構(gòu)建自由水譜的期望值;0為構(gòu)建自由水譜 的標(biāo)準(zhǔn)差;X為構(gòu)建自由水譜各組分的時間值，與測井T2譜時間一致;f(x)為構(gòu)建自由水譜 的幅值。
[0044] 本發(fā)明將自由水譜與離屯、譜相加得到一個總譜，利用規(guī)劃求解方法使得構(gòu)建的總 譜孔隙度取f和測井T2譜幾何均值T2g/與步驟(3)中的測井孔隙度餐t及步驟(4)中確定的T2g 值誤差最小，此時得到的總譜即為100%純水譜，
[0045] 本發(fā)明實現(xiàn)了多相流體存在情況下利用核磁共振測井資料構(gòu)造100%純水譜，純 水譜能夠反映儲層真實的孔隙結(jié)構(gòu)，并且通過對比實測核磁共振譜與構(gòu)造純水譜的差異能 夠有效識別儲層流體性質(zhì)，在構(gòu)建水譜基礎(chǔ)上可W進行儲層品質(zhì)的評價，回避了由于原油 性質(zhì)、觀測模式等導(dǎo)致的核磁共振測井流體識別困難的難題，并且可W通過對比構(gòu)造水譜 與實測譜之間的差異達到流體識別的目的。提高了核磁共振測井資料的應(yīng)用，增強了解決 地質(zhì)問題的能力。
【附圖說明】
[0046] 圖1為本發(fā)明一種核磁共振測井T2譜100%純水譜的確定方法流程示意圖。
[0047] 圖2為本發(fā)明一種核磁共振測井T2譜100 %純水譜的確定方法實施例T2g確定方法 示意圖。
[004引圖3為本發(fā)明一種核磁共振測井T2譜100%純水譜的確定方法實施例m值確定方法 示意圖。
[0049] 圖4為本發(fā)明一種核磁共振測井T2譜100%純水譜的確定方法實施例單峰型示意 圖。
[0050] 圖5本發(fā)明一種核磁共振測井T2譜100%純水譜的確定方法實施例雙峰型示意圖。
[0051] 圖6本發(fā)明一種核磁共振測井T2譜100%純水譜的確定方法實施例雙峰型示意圖。
[0052] 圖7為本發(fā)明一種核磁共振測井T2譜100%純水譜的確定方法計算機流程圖。
[0053] 圖8為本發(fā)明一種核磁共振測井T2譜100%純水譜的確定方法某地區(qū)某井實際處 理的效果圖。
[0054] 圖9為本發(fā)明一種核磁共振測井T2譜100%純水譜的確定方法某地區(qū)計算的孔隙 度與巖屯、分析孔隙度結(jié)果對比圖。
[0055] 圖10為本發(fā)明一種核磁共振測井T2譜100%純水譜的確定方法某地區(qū)計算的T2g與 巖屯、分析T2g結(jié)果對比圖。
【具體實施方式】
[0056] 下面結(jié)合附圖和實施例對發(fā)明作進一步的詳細說明，但并不作為對發(fā)明做任何限 制的依據(jù)。
[0057] 通過對核磁共振T2譜的形態(tài)研究，在大量實驗數(shù)據(jù)和實測資料的基礎(chǔ)上，提出了 一種核磁共振測井T2譜100%純水譜的確定方法，該方法的流程示意圖如圖1所示。為使本 發(fā)明的目的、技術(shù)方案和優(yōu)點更加清楚，下面將結(jié)合附圖對本發(fā)明實施方式作進一步地詳 細描述。
[005引第一步，按照《巖樣核磁共振參數(shù)實驗室測量規(guī)范SY/T 6490-2007》標(biāo)準(zhǔn)規(guī)定的流 程進行巖屯、核磁T2實驗，獲取實驗核磁共振&譜飽和譜和離屯、譜數(shù)據(jù)。
[0059] 第二步，根據(jù)第一步得到的巖屯、核磁共振T2飽和譜和離屯、譜實驗數(shù)據(jù)計算束縛水 飽和度Swi及核磁共振T2譜幾何均值T2g，并建立二者關(guān)系，如圖2所示，具體確定方法如下：
[0060]
[0061]
[0062]
[0063] 其中：取Swa為巖屯、離屯、譜孔隙分量;斬為巖屯、飽和譜孔隙分量;T2,i為實驗T2譜的 橫向弛豫時間；
[0064] 具體實施過程中，采用最小二乘法確定待定系數(shù)a、b，a值取0.878，b值取1.149。
[0065] 第=步，利用球管薄膜模型計算測井T2譜的束縛水飽和度Swi并構(gòu)建離屯、譜，具體 確定方了-
[0066]
[0067]
[006引
[0069]
[0070]
[007。 其中：從堿i為離屯、譜分量；從wi為束縛水孔隙度;Wi為權(quán)系數(shù)分量;m為系數(shù)；Swi 為束縛水飽和度;輯￡為測井孔隙度;取呪酣為測井密度孔隙度;與為測井中子孔隙度;特左 為測井T2譜孔隙分量;T2,功測井T2譜的橫向弛豫時間。
[0072] 我們發(fā)現(xiàn)，利用球管模型時權(quán)系數(shù)曲線更加平滑，取值更加合理：當(dāng)T2小于某一個 值時，孔隙空間充滿了束縛流體，此時權(quán)系數(shù)為1，隨著孔喉半徑的增大，不可動流體所占空 間比例不斷減小，權(quán)系數(shù)趨近于0,權(quán)系數(shù)是T2值的函數(shù)。對于要準(zhǔn)確構(gòu)造離屯、譜，則上述參 數(shù)m值取值很關(guān)鍵，我們通過大量實驗數(shù)據(jù)分析，發(fā)現(xiàn)m值與T2截止值T2cutDff有很好的線性關(guān) 系，即：
[0073] m = c*T2cut〇ff+d
[0074] 如圖3所示，具體實施過程中，采用最小二乘法確定待定系數(shù)系數(shù)c、d，c值取 0.196，d值取-0.364;當(dāng)缺少地區(qū)經(jīng)驗T2截止值時，默認m取2.46。
[0075] 第四步，根據(jù)第二步中確定的T2g~Swi關(guān)系W及待定系數(shù)a、b，計算細=1時的核磁 共振T2譜幾何均值T2g。
[0076] 第五步，借助正態(tài)分布函數(shù)在測井T2譜上構(gòu)建一個自由水譜：
[0077]
[0078] 其中：A為構(gòu)建自由水譜的幅度值;U為構(gòu)建自由水譜的期望值;0為構(gòu)建自由水譜 的標(biāo)準(zhǔn)差;X為構(gòu)建自由水譜各組分的時間值，與測井T2譜時間一致;f(x)為構(gòu)建自由水譜 的幅值。
[0079] 第六步，將第五步構(gòu)建的自由水譜與第=步構(gòu)建的離屯、譜相加得到一個總譜，具 體實施過程中，我們采用規(guī)劃求解的方法，使得構(gòu)建的總譜孔隙度取寸日核磁共振幾何均值 T2g/與第=步中的測井孔隙度巧t及第四步中確定的T2g值誤差最小，此時得到的總譜即為 100 %純水譜，即：
[0080]
[0081]
[0082] 對于要準(zhǔn)確構(gòu)造純水譜，則上述A、u、〇 =個參數(shù)的確定也是至關(guān)重要，在具體實施 過程中，采用規(guī)劃求解方法確定，但是通過大量數(shù)據(jù)分析發(fā)現(xiàn)，=個參數(shù)都有一個明確的取 值范圍，A值取值在0.35附近，U值取值在1.70附近，0取值在0.44附近。
[0083] 按照W上思路，利用構(gòu)造水譜法我們對某地區(qū)的部分巖屯、做了計算，如圖4-6,巖 屯、數(shù)據(jù)如表1。
[0084] 同時我們在geolog中編寫了該方法的程序，實現(xiàn)了自動處理。流程圖如圖7所示。
[0085] 表 1 「OORAl
[0087]
[0088] 利用構(gòu)造水譜法自動構(gòu)造核磁共振測井T2譜100%純水譜方法W及自主開發(fā)的軟 件包，我們對某油田某井進行了處理，如圖8所示為某油田某井實際處理的效果圖，處理的 效果圖如下：
[0089] 圖中第一道為CAL(井徑曲線）、SP(自然電位曲線)和GR(自然伽馬曲線），第二道為 深度，第=道為=孔隙度曲線，第四道為電阻率曲線，第五道為實測T2譜孔隙分布譜，第六 道和第屯道分別為構(gòu)造離屯、譜與純水譜。從圖上可W看出，利用構(gòu)造水譜法得到的100%純 水譜在含水地層與實測譜形態(tài)基本一致，說明該方法的可信性。
[0090] 和現(xiàn)有的通過含控校正得到純水譜的方法相比，構(gòu)造水譜法可W快速便捷的逐點 計算每個深度點的100%純水譜，受流體類型等的影響較小。利用自主開發(fā)的軟件包，我們 就可W得到一系列連續(xù)的和測井井段相匹配的純水譜，進而利用該孔隙分布更加準(zhǔn)確的對 儲層品質(zhì)及流體類型進行評價。利用構(gòu)造法，我們對某地區(qū)的巖屯、做了計算，圖9中的兩條 邊界線分別是±5%為誤差的分界線，圖10中的兩條邊界線是W±5%為誤差分界線，從圖9 和圖10上可W看出構(gòu)造法的結(jié)果可W很好的替代巖屯、分析的結(jié)果。
[0091] 最后所應(yīng)說明的是，W上實施例僅用W說明本發(fā)明的技術(shù)方案而非限制，盡管參 照較佳實施例對本發(fā)明進行了詳細說明，本領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解，可W對本發(fā)明 的技術(shù)方案進行修改或者等同替換，而不脫離本發(fā)明技術(shù)方案的精神和范圍，其均涵蓋在 本發(fā)明的權(quán)利要求范圍當(dāng)中。
【主權(quán)項】
1. 一種核磁共振測井1'2譜100 %純水譜的確定方法，其特征在于，所述方法包括以下步 驟： 步驟(1 )，對巖心進行核磁共振T2實驗，獲取核磁共振!^飽和譜和離心譜實驗數(shù)據(jù)； 步驟(2)，根據(jù)步驟（1)得到的巖心核磁共振T2飽和譜和離心譜實驗數(shù)據(jù)計算巖心束縛 水飽和度SW1及核磁共振T2譜幾何均值T2g，并建立二者關(guān)系，采用最小二乘法確定待定系數(shù) a、b; 步驟（3)，利用球管薄膜模型計算核磁共振^譜的巖心束縛水飽和度SW1，并構(gòu)建離心 譜； 步驟(4)，根據(jù)步驟(2)確定的待定系數(shù)a、b，計算含水飽和度Sw=l時的核磁共振T2譜幾 何均值T2g; 步驟(5)，借助正態(tài)分布函數(shù)在核磁共振!^譜上構(gòu)建一個自由水譜； 步驟(6)，將步驟(3)構(gòu)建的離心譜與步驟(5)構(gòu)建的自由水譜相加得到一個總譜，利用 規(guī)劃求解的方法使得構(gòu)建的總譜孔隙度#和幾何均值T2f/與步驟(3)中計算巖心束縛水飽 和度SW1產(chǎn)生的中間參數(shù)測井孔隙度爐t及步驟(4)中確定的核磁共振T 2譜幾何均值T2g值誤 差最小，此時得到的總譜即為100 %純水譜。2. 如權(quán)利要求1所述的方法，其特征在于，所述步驟(1)中，巖心核磁共振^實驗按照《巖 樣核磁共振參數(shù)實驗室測量規(guī)范SY/T 6490-2007》標(biāo)準(zhǔn)規(guī)定的流程進行。3. 如權(quán)利要求1所述的方法，其特征在于，所述步驟(2)中，巖心束縛水飽和度SW1及核磁 共振T2譜幾何均值T 2g，具體按照如下公式進行計算：其中：識Swi,i為離心譜孔隙分量；_為飽和譜孔隙分量;T2>1為實驗^譜的橫向弛豫時 間。4. 如權(quán)利要求1所述的方法，其特征在于，所述步驟(2)中，建立巖心束縛水飽和度SW1及 核磁共振!^譜幾何均值T 2g的關(guān)系如下：5. 如權(quán)利要求1所述的方法，其特征在于，在進行步驟(3)操作之前，還需利用中子-密 度曲線進行測井孔隙度街的計算，通過下式實現(xiàn)：其中：識0SW為測井密度孔隙度;識C継為測井中子孔隙度。6. 如權(quán)利要求5所述的方法，其特征在于，所述測井孔隙度識t的計算按照石油天然氣行 業(yè)標(biāo)準(zhǔn)《生產(chǎn)測井解釋規(guī)程SY/T 5793-1993》規(guī)定的流程進行。7. 如權(quán)利要求1所述的方法，其特征在于，所述步驟(3)中計算的束縛水飽和度并構(gòu)建 離心譜，是根據(jù)球管薄膜模型確定的權(quán)系數(shù)計算得到的：其中，爐swi,i力離心譜分量;為束縛水孔隙度;Wi為權(quán)系數(shù)分量;m為系數(shù);Swi為束縛 水飽和度;爐t為測井孔隙度;#i為測井T2譜孔隙分量;T2, i為測井T2譜的橫向弛豫時間。8. 如權(quán)利要求1所述的方法，其特征在于，所述步驟(5)借助正態(tài)分布函數(shù)在核磁共振!^ 譜上構(gòu)建一個自由水譜：其中：A為構(gòu)建自由水譜的幅度值;u為構(gòu)建自由水譜的期望值;〇為構(gòu)建自由水譜的標(biāo) 準(zhǔn)差;x為構(gòu)建自由水譜各組分的時間值，與測井!^譜時間一致;f(x)為構(gòu)建自由水譜的幅 值。9. 如權(quán)利要求1所述的方法，其特征在于，通過所述步驟(6)中建立的100%純水譜進行 儲層品質(zhì)評價，并通過與測井!^譜對比進行儲層流體識別。
【文檔編號】E21B49/08GK106050225SQ201610394292
【公開日】2016年10月26日
【申請日】2016年6月6日
【發(fā)明人】王明方, 萬金彬, 李慧瑩, 程道解, 白松濤, 黃科, 許思勇, 金小慧, 王琪, 張旭, 黃婭, 田揚, 趙建斌, 任龍, 曾靜波, 何羽飛, 袁野
【申請人】中國石油天然氣集團公司, 中國石油集團測井有限公司