鄒友龍，胡法龍，李長喜，李潮流

（1.中國石油勘探開發(fā)研究院測井與遙感技術(shù)研究所，北京100083；2.中國地質(zhì)大學(xué)地球物理與信息技術(shù)學(xué)院，北京100083）

利用徑向基函數(shù)預(yù)測巖石滲透率及流體分子組分

鄒友龍1，2，胡法龍1，李長喜1，李潮流1

（1.中國石油勘探開發(fā)研究院測井與遙感技術(shù)研究所，北京100083；2.中國地質(zhì)大學(xué)地球物理與信息技術(shù)學(xué)院，北京100083）

測井反問題研究中，經(jīng)常遇到正演模型未知、經(jīng)驗公式的相關(guān)系數(shù)較小的情形，且現(xiàn)有的反問題處理方法難以獲得理想的效果。徑向基函數(shù)插值方法是目前最有效的反問題處理方法之一，該技術(shù)能夠生成平滑連續(xù)的多個變量的多元函數(shù)，逼近未知的正演模型。介紹徑向基函數(shù)插值方法原理，針對核磁共振實驗的T2分布，采用徑向基函數(shù)插值方法預(yù)測了巖石滲透率及流體分子組分，實現(xiàn)了多維到一維映射及多維到多維映射的2種情況，并且預(yù)測結(jié)果具有較高精度，測量值的對比驗證了該方法解決測井反問題的有效性。

測井反問題；正演模型；徑向基函數(shù)插值；滲透率；分子組分

0 引 言

在測井反問題中，精確的正演模型往往未知?，F(xiàn)有計算滲透率的經(jīng)驗?zāi)Ｐ蜑镃oates模型［1］和SDR模型［2］，F(xiàn)reed［3］和Heaton［4］等曾嘗試?yán)贸谠r間和擴散系數(shù)分布建立與分子組分的物理模型。但是經(jīng)驗公式經(jīng)常會遇到相關(guān)系數(shù)低、誤差較大的情況，主要原因是采用的輸入信息太少，受地區(qū)和地質(zhì)條件等多重因素影響較大。經(jīng)驗公式常常是多維輸入到一維輸出的映射，對于多維輸入到多維輸出的映射，經(jīng)驗公式無法實現(xiàn)。2006年Freedman［5］提出了利用徑向基函數(shù)（RBF）插值的方法處理巖心數(shù)據(jù)分析中遇到的一些反問題，并且取得一定應(yīng)用效果。徑向基函數(shù)像軟膜一樣通過各個已知的樣本數(shù)據(jù)點，生成連續(xù)光滑的多變量多元函數(shù)可以逼近正演模型。本文采用徑向基函數(shù)插值方法，利用T2分布分別預(yù)測巖石滲透率與流體分子組分，實現(xiàn)了多維到一維映射及多維到多維映射的2種情況，并且預(yù)測結(jié)果具有較高精度，驗證了利用徑向基函數(shù)解決測井反問題的有效性，為測井解釋中關(guān)鍵參數(shù)獲取提供了一個新思路。

1 徑向基函數(shù)插值原理

逼近函數(shù)形式為

若矩陣的條件數(shù)較大或者實驗數(shù)據(jù)具有誤差時，此時可取正則化解

采用高斯徑向基函數(shù)，即

式中，si為高斯函數(shù)的寬度。s的選取直接影響到預(yù)測結(jié)果的好壞，它的取值與最鄰近距離（NND）成正比，s＝α·NND，α為比例常數(shù)。聯(lián)合式（2）和式（5）有

因此，系數(shù)解的矩陣形式為

式中，N為樣品的個數(shù)；m為樣品預(yù)測的屬性個數(shù)；

考慮到各向異性的問題，φi，j寫為

將式（6）改寫為

如果忽略徑向基函數(shù)之間的重疊，即當(dāng)高斯函數(shù)的寬度足夠小時，式（9）可簡寫為

因此，由式（9）和式（10）有

其次，美國政府具有制裁規(guī)則制定者和執(zhí)行者的雙重身份。制裁政策主要由以海外資產(chǎn)控制辦公室為核心的財政部及相關(guān)部門起草制定，對于政策的解讀和澄清也主要由財政部負責(zé)，海外資產(chǎn)控制辦公室的大量解釋性指導(dǎo)意見和常見問題回答則是制裁政策的重要執(zhí)行依據(jù)，以財政部為核心的政府相關(guān)部門會通力配合實施制裁。因此，美國政府對于制裁誰和如何制裁具有完全的決策權(quán)威。

2 應(yīng)用實例

核磁共振測井獲得的數(shù)據(jù)T2分布是一個二維數(shù)據(jù)，而現(xiàn)有滲透率模型均采用基于T2分布的某些特定參數(shù)，如T2幾何均值、自由流體指數(shù)和束縛流體指數(shù)等，無法將T2分布中所有數(shù)據(jù)信息應(yīng)用到滲透率計算中，其滲透率求取效果并不理想。2001年Freedman等人發(fā)現(xiàn)弛豫時間和擴散系數(shù)可以提供分子組分的信息［7］，并試圖推導(dǎo)出相互之間的物理模型［3－4］，然而原油是含烴分子的復(fù)雜混合物，推導(dǎo)出準(zhǔn)確的物理模型非常困難［5］。徑向基函數(shù)在解決多維數(shù)據(jù)反問題的研究中優(yōu)勢明顯，以下利用徑向基函數(shù)技術(shù)分別實現(xiàn)T2分布求取滲透率與分子組分，準(zhǔn)確獲得了測井解釋中的關(guān)鍵參數(shù)。

本文采用的是leave－one－out方法［5］，即從所有的樣本中取出1顆待預(yù)測的巖心，然后利用剩下的所有巖心數(shù)據(jù)求解得到系數(shù)矩陣C，得到構(gòu)造函數(shù)），最后計算得到取出巖心的預(yù)測值，這樣依次循環(huán)對每顆巖心進行預(yù)測。通過對預(yù)測結(jié)果與測量值的比較，發(fā)現(xiàn)該方法運用效果較好。

2.1 T2分布預(yù)測滲透率

試驗樣品來自渤海灣盆地東營組巖心資料，主要包括基礎(chǔ)孔隙度、滲透率和核磁共振巖心分析數(shù)據(jù)。圖1所示為東營組12顆巖心經(jīng)過核磁共振實驗得到的T2分布及其對應(yīng)的滲透率，12顆巖心的核磁共振實驗和滲透率數(shù)據(jù)均在相同條件下測定。

圖1 東營組12顆巖心的T2分布及其滲透率

由徑向基函數(shù)插值原理可知，預(yù)測滲透率的徑向基映射函數(shù)可寫為

為了對比RBF插值方法的效果，利用SDR模型，得到滲透率計算公式為

式中，φ為孔隙度；T2，gm為T2幾何均值。圖2為滲透率K與×φ4的交會圖，從圖2中可以看到，兩者之間的擬合關(guān)系較差。

圖2 滲透率K與×φ4交會圖

式中，log vm為測量值取自然對數(shù)；log vp為預(yù)測值取自然對數(shù)。

圖3為經(jīng)驗公式計算的滲透率與測量滲透率的對比圖，對應(yīng)的誤差χ2為1.743 8。由于系數(shù)α的選取直接影響預(yù)測效果的好壞，通過考察α對預(yù)測結(jié)果的影響，選取α的最佳值。取α以0.05為步長從0.1～2.0變化。圖4為預(yù)測值與測量值的誤差χ2隨α的變化關(guān)系。從圖4中可以看出α的最優(yōu)值在0.4～0.7左右。圖5為α＝0.65時預(yù)測結(jié)果與測量值的對比，預(yù)測滲透率與測量滲透率的誤差在1倍范圍之內(nèi)，相應(yīng)的預(yù)測值與測量值誤差χ2為 0.513 8小于SDR模型預(yù)測結(jié)果的誤差，預(yù)測結(jié)果的精度較高。

圖3 經(jīng)驗公式計算的滲透率與測量結(jié)果對比

圖4 預(yù)測滲透率與測量結(jié)果的誤差隨α的變化

圖5 α＝0.65時預(yù)測滲透率與測量結(jié)果的對比

2.2 T2分布預(yù)測流體分子組分

原油的分子組分是評價儲層質(zhì)量，產(chǎn)能和劃分的一個關(guān)鍵參數(shù)［8］。原油是含烴分子的復(fù)雜混合物，而核磁共振測井能為井下含氫流體的定性及定量分析提供重要資料，因此核磁共振測井資料可用來分析原油的組分。利用徑向基函數(shù)插值的方法運用T2分布預(yù)測流體分子組成，其中T2分布與分子組分?jǐn)?shù)據(jù)均采用Freedman文獻中的T2分布與分子組分?jǐn)?shù)據(jù)。預(yù)測分子組分的徑向基映射函數(shù)形式與式（12）一致，輸出參數(shù)為分子組分向量，輸入?yún)?shù)為經(jīng)過最大值標(biāo)準(zhǔn)化后的T2分布。

該部分的預(yù)測值與測量值的誤差χ2定義為

式中，vm為測量值；vp為預(yù)測值。

同樣考察系數(shù)α對預(yù)測結(jié)果的影響，取α以0.1為步長從0.1～2.5變化。圖6為預(yù)測值與測量值的誤差χ2隨α的變化。從圖6中可以看出，系數(shù)α的最優(yōu)值在1.4左右，相應(yīng)的預(yù)測值與測量值的誤差χ2為1 129.3。圖7為α＝1.4時預(yù)測的分子組分與測量分子組分的對比，預(yù)測值與測量值基本重合，預(yù)測效果比較滿意。

圖6 預(yù)測分子組分與測量結(jié)果的誤差隨α的變化

3 結(jié)論與建議

（1）徑向基函數(shù)法為滲透率和分子組分的預(yù)測提供了一種新思路，并可推廣到其他測井反問題中。

（2）徑向基函數(shù)法可解決測井中的反問題，且具有以下優(yōu)點：① 當(dāng)正演模型未知時，該方法是解決測井反問題的一個很好的工具；② 其不受維數(shù)的限制，輸入與輸出參數(shù)都可為多維；③ 相對于經(jīng)驗公式，可考慮的信息量更多，使得逼近函數(shù)更加趨近于正演模型。

（3）該方法適用于多維數(shù)據(jù)輸入與多維數(shù)據(jù)輸出的情況，且計算量小運算速度快，這是經(jīng)驗公式和其他方法無法比擬的。

（4）該方法預(yù)測的效果受α的影響，α的一般取值范圍在0.5～1.5左右，α取值不宜過大，否則預(yù)測值接近于實驗數(shù)據(jù)的平均值。

圖7 預(yù)測分子組分與測量結(jié)果對比

（5）雖然該方法具有許多優(yōu)點，但也存在一些需要解決的問題，如實驗數(shù)據(jù)的誤差對預(yù)測結(jié)果的影響，以及如何自動選取合適的高斯函數(shù)寬度。

［1］ Coates G R，Marschall D，Mardon D，et al.A New Characterization of Bulk－volume Irreducible Using Magnetic Resonance［C］∥SPWLA 38th Annual Logging Symposium Transactions：Society of Petrophysicists and Well Log Analysts，Paper QQ，1997.

［2］ Kenyon W E，Day P I，Straley C，et al.A Three－part Study of NMR Longitudinal Relaxation Properties of Water－saturated Sandstones［J］.SPE Formation Evaluation，1988，3：622－636.

［3］ Freed D.Molecular Composition from Diffusion or Relaxation Measurements：US 2004／0253743A1［P］.2004－12－16.

［4］ Heaton N，F(xiàn)reedman R.Method for Determining Molecular Properties of Hydrocarbon Mixtures from NMR Data：U.S.Patent 6 859 032［P］.2005.

［5］ Freedman R.New Approach for Solving Inverse Problems Encountered in Well－logging and Geophysical Applications［J］.Petrophysics，2006，47：93－111.

［6］ Lukaszyk S.A New Concept of Probability Metric and Its Applications in Approximation of Scattered Data Sets［J］.Computational Mechanics，2004，33：299－3004.

［7］ Freedman R，Lo S，F(xiàn)laum M，Hirasaki G J，et al.A New NMR Method of Fluid Characterization in Reservoir Rocks：Experimental Confirmation and Simulation Results［J］.SPE Journal，2001，4：52－464.

［8］ Anand V，F(xiàn)reedman R.New Methods for Predicting Properties of Live Oils from NMR［C］∥Paper DD in 50th SPWLA Annual Logging Symposium Transactions：Society of Petrophysicists and Well Log Analysts，June 21－24，Paper AAAA，2009.

Prediction of Rock Permeability and Fluid Molecular Composition Using Radial Basis Function

ZOU Youlong1，2，HU Falong1，LI Changxi1，LI Chaoliu1
（1.Research Institute of Petroleum Exploration ＆Development，Beijing 100083，China；2.School of Geophysics ＆Information Technology，China University of Geosciences，Beijing 100083，China）

In logging inverse problem research，the forward models are often unknown and correlation coefficients of empirical equations are small，and the existing inverse problem approaches are difficult to obtain the desired results.Radial basis function（RBF）interpolation is one of the most effective approaches to deal with inverse problems，which can generate smooth and continuous multivariate functions of many variables to approximate the forward mode.Introduced is the theory of RBF interpolation.According to T2distributions of the NMR experiments，predicted are rock permeability and fluid molecular compositions using RBF interpolation method.Achieved are mapping from multi－dimensional to one－dimensional and from multi－dimensional to multi－dimensional，and prediction accuracy is high.Verified is the validity of RBF interpolation method to solve the logging inverse problems by comparing the predicted results with measurements.

logging inverse problem，forward model，radial basis function（RBF）interpolation，permeability，molecular composition

P631.84

A

2011－12－12 本文編輯 余迎）

1004－1338（2012）03－0225－05

鄒友龍，男，1988年生，碩士研究生，從事測井資料反演與巖石物理實驗分析研究。