自適應截斷奇異值結合有效集的T2譜反演方法

李鵬舉， 王梓丞， 任 莉

(1.東北石油大學 地球科學學院， 黑龍江 大慶163318； 2.東北石油大學 非常規油氣成藏與開發省部共建國家重點實驗室培育基地， 黑龍江 大慶163318； 3.中國石油集團測井有限公司大慶分公司解釋評價中心， 黑龍江 大慶163000)

0 引 言

核磁共振測井技術能夠提供地層巖石物理性質和流體性質參數，如孔隙度、滲透率、可動水飽和度等[1,2]。這些參數能否計算準確與T2解譜技術息息相關。解譜技術是核磁共振測井技術的關鍵。目前，廣泛使用的解譜技術的實質是多指數反演，反演方法有奇異值分解(SVD)[1,3-5]和最優化及其改進方法[6]。例如：肖立志等[1,5]基于奇異值分解的反演方法；王為民等[7]提出基于解的自由度選取最優反演模型的反演變換法；K.J.Kunn等[8]，王才志等[9]采用奇異值分解與正則化方法結合截去小的非零奇異值的方法。學者們普遍采用截斷法或正則化方法確保解的穩定與平滑。截斷奇異值的選擇和正則化因子的選取與噪聲有關，一般采用信噪比作為噪聲的統計特征參數。而CPMG回波信號為多指數衰減信號，與穩恒信號不同，信噪比并不是常數[3]。因此，以信噪比作為衰減信號的統計特征參數不適合。以噪聲的方差作為噪聲的統計特征參數，來確定截止奇異值，能夠最大限度提取自旋回波中的有用信息。根據解的物理意義，T2譜的幅值應為非負實數。為了保證解的非負性，廣泛采用改進迭代法隨機搜索非負解，而這種解法可能會喪失解的連續性[4,10]。……