遺傳最優化算法在砂礫巖儲層測井評價中的應用

韓雪，潘保芝，張意，王飛

（吉林大學地球探測科學與技術學院，吉林長春130026）

遺傳最優化算法在砂礫巖儲層測井評價中的應用

韓雪，潘保芝，張意，王飛

（吉林大學地球探測科學與技術學院，吉林長春130026）

砂礫巖儲層作為非常規的油氣藏其儲層具有埋藏深、巖性復雜、低孔隙度低滲透率、孔隙結構復雜、非均質性強等特點，難以精確劃分巖性并建立準確的解釋模型，給儲層參數的求取帶來很大的困難。最優化測井解釋可充分利用所有測井信息、測量誤差和響應方程誤差及地質資料和工作經驗，從所有可能的解釋結果中找出最合理的解釋結果。將遺傳算法引入砂礫巖儲層最優化測井解釋，取代經典的數學算法求取儲層孔隙度等參數。對松遼盆地實際測井數據進行了處理，取得了較好的效果。

測井解釋；砂礫巖儲層；最優化測井解釋；遺傳算法；多組分模型

0 引 言

非常規油氣資源主要是指油頁巖、油砂礦、煤層氣、頁巖氣、致密砂巖氣等［1］。砂礫巖儲層測井評價、油氣層解釋、地質和工程應用在國內外石油勘探領域均屬評價難題，需要新的理論和先進的工程技術支撐［2］。描述及評價砂礫巖有效儲層，建立適合于砂礫巖儲層的測井評價方法，對這類油氣儲層的勘探開發具有重要意義。

最優化測井解釋是20世紀80年代發展起來的技術，它首次根據廣義地球物理反演理論，將最優化數學方法和概率統計理論應用于測井資料解釋，為綜合應用所有測井資料評價復雜儲層開辟了新途徑。在尋優過程中傳統的優化方法是基于精確數學方法，對數據的確定性和準確性有嚴格的要求，且注重理論的最優性，而對于復雜的實際問題，適應性更強的元啟發算法，如遺傳算法、粒子群算法、模擬退火以及蟻群算法等在尋找最優解上表現出越來越強的優勢，因而越來越多的學者將這些算法引入到最優化算法中。例如馮國慶等用遺傳算法進行最優化測井解釋，并結合整體模擬退火算法對潛山油藏測井資料進行了解釋取得了可靠的結果［3］；董勇等人將自適應混沌粒子群算法引入最優化算法中，在油水兩相流動分析中的應用取得了較好的效果，得到了油水的日產量［4］。將遺傳算法引入最優化測井解釋很有意義［5］。本文將遺傳算法引入到最優化測井解釋中并應用于砂礫巖儲層的評價，得到了較好的效果。

1 最優化測井解釋基本原理

最優化測井解釋根據地球物理學中的廣義反演理論，以經過環境校正較為真實反映地層及其孔隙流體特性的測井值ai為基礎，以儲層參數和礦物相對體積為自變量x，采用適當的解釋模型和測井響應方程，通過合理選擇區域性解釋參數z與儲層參數初始值x，反演出相應的理論測井值ai＝fi（x，z），并將它們與實際測井值比較，根據非線性加權最小二乘法原理和誤差理論，結合地區地質經驗參數建立最優化測井解釋的數學模型，應用最優化技術不斷調整未知儲層參數x，使計算出的各理論測井值不斷地逼近相應的實際測井值，兩者充分逼近時，目標函數值達到了極小值，此時用以計算各理論測井值的自變量x就是充分反映實際地層的儲層參數及礦物相對含量，即最優化測井的解釋結果［6］。其基本流程見圖1。

圖1 最優化測井解釋基本原理

應用最優化方法進行測井解釋的優化數學模型可歸結為

式中，F（x，a）為測井解釋的目標函數；ai為實際測井值；x為未知儲層參數向量；z為區域性解釋參數向量；m為測井曲線個數；σi為第i種測井值的測量誤差；τi為響應方程誤差；gj（x）≥0為不等式約束；hk（x）＝0為等式約束。

未知參數向量的分量一般為孔隙度、含水飽和度、沖洗帶含水飽和度、泥質含量和n種地層骨架礦物含量

建立好最優化測井解釋的數學模型之后即可用各種最優化方法搜索最優解，這里將遺傳算法引入最優化算法中來尋求最佳儲層參數。

2 遺傳算法基本原理及遺傳最優化解釋流程

遺傳算法從代表問題可能潛在解集的一個種群開始，計算開始時，隨機產生一定數目的個體，并且計算每個個體的適應度，隨即產生第1代，若不滿足優化準則，則通過選擇、交叉和變異重新計算個體的適應度，子代被插入到種群中將父代取代，構成新的一代。這一過程循環執行，直到滿足優化準則為止。

2.1 遺傳算法編碼

選擇實數編碼對解空間進行編碼，考慮實空間Rm中的向量，假設解向量是n維，并且每個分量都在有限區間上，設解向量s＝（v1，v2，…，vm），則有vi∈［ai，bi］，其中i＝1，2，…，m。

2.2 適應度函數

遺傳算法中以數值方式描述個體優劣程度的指標，由適應度函數f（x）計算得到。適應度函數f（x）一般是求最大值。對于最優化測井解釋，適應度函數可表示為

2.3 遺傳操作算子

（2）交叉算子。選擇單點交叉算子，在染色體隨機地設定一個交叉點，將該點處2個染色體的某個對應的變量進行互換，并生成2個新的個體。

（3）變異算子。選擇基本位變異，依照一定的變異率，隨機選擇染色體變異位置，重新生成自變量范圍內的值。

2.4 遺傳最優化測井解釋流程

遺傳算法求解優化問題的具體步驟為5個部分（流程圖見圖2）。

圖2 遺傳算法解決函數優化問題的流程圖

3 正演模擬

砂礫巖儲層巖性復雜，所含的石英、長石及各種巖屑的含量變化大（見圖3），以至于儲層參數尤其骨架參數不容易選準。多礦物模型分析的最優化解釋能采用多種測井信息和多種模型研究和評價復雜油氣層［7］，因此測井解釋模型宜選擇多礦物模型。由于砂礫巖儲層中巖屑類礦物成分復雜多變，為簡化問題，將各巖屑類的礦物作為一個整體（見圖4），將骨架細分為石英、長石及巖屑等3個部分建立多組分模型，與多礦物模型相同，多組分模型也是把一個巖性礦物成分復雜的地層看成由局部均勻的幾部分組成的：有效孔隙流體、泥質和幾種骨架礦物。地層的有效孔隙度為φ，泥質含量為Vcl，礦物總體積為∑iVmai。本文所選取的儲層未知參數為φ、Vcl、Vmai（i＝1，2，3）（對應石英、長石、巖屑）。

圖3 砂礫巖各組分含量三角圖

圖4 多組分體積模型

根據多組分體積模型建立測井響應方程。由于遺傳算法的引入，方程的個數與未知量個數的關系可以不嚴格要求，模型未知參數為5個。選擇參與計算的測井曲線有4條：GR、CNL、DEN、AC。測井響應方程為

平衡方程為

式中，平衡方程作為懲罰項加入最優化目標函數中控制未知量的和。按常規解方程組的算法，式（4）至式（8）可組成一個適定的方程組并解出一組唯一解，但是用這種方法解出的未知量值不可避免會出現大量的實際意義之外的值，即大于1或小于0，并且由于未知量數目較多，很難人為調整其數值使之符合實際。用遺傳算法求解是在未知量有意義的定義域范圍內搜索最優解，更加便于應用。

為測試程序的可靠性，首先構造了一定深度段各深度點的孔隙度、泥質含量及骨架礦物含量等5個儲層參數。根據式（4）至式（7）帶入構造的儲層參數及適當的解釋參數，構造偽測井曲線GR、CNL、DEN、AC，用遺傳最優化程序反演理論測井值，用最小二乘理論將反演的理論測井值與構造的偽測井值相比，經遺傳算法循環執行尋優過程，最終挑選出最優的儲層參數。圖5為反演的儲層參數與構造值交會圖，反演孔隙度相對構造孔隙度的平均絕對誤差為0.009，平均相對誤差為0.11；反演泥質含量相對構造泥質含量的平均絕對誤差為0.08，平均相對誤差為0.43；反演石英含量相對構造石英含量的平均絕對誤差為0.05，平均相對誤差為0.22。可以用該程序進行計算。

4 實際資料處理

4.1 常規測井方法求孔隙度

密度測井和中子測井確定含流體純巖石孔隙度的關系式分別為

圖5 儲層參數反演結果與構造值交會圖

圖6 常規算法孔隙度與巖心孔隙度交會圖

式中，ρma和ρf分別為巖石骨架和孔隙流體的密度；ρb為密度測井讀數；φNma和φNf分別為巖石骨架和孔隙流體的中子值；φCNL為中子測井讀數。

儲層有效孔隙度的確定，有

所求孔隙度與巖心孔隙度交會圖見圖6。經計算，常規算法求得的孔隙度與巖心孔隙度的平均相對誤差為0.23，平均絕對誤差為0.012。

4.2 遺傳最優化算法求儲層參數

用遺傳最優化測井解釋程序對松遼盆地×井砂礫巖儲層進行資料處理。圖7為儲層參數計算結果，深度道左側3個道為實際測井曲線，右側第1道為計算出的巖性剖面，φcor為巖心孔隙度；第2道至第4道為石英、長石及巖屑含量與薄片分析的礦物含量對比。由于程序中連續性約束條件所設置的誤差較大，計算的儲層參數結果仍有一定程度的毛刺。從圖7可見，所求的孔隙度與巖心孔隙度有較好的吻合度。經統計，平均相對誤差為0.11，平均絕對誤差為0.010，較常規算法取得的效果稍好一些。

圖7 ×井測井資料遺傳最優化處理成果圖

圖8為重構曲線與實際測井曲線對比的質量檢驗圖，CNL0、DEN0、AC0及GR0為重構的理論曲線，由于遺傳算法搜尋最優解要在具有實際意義的解空間內，并要滿足懲罰方程及約束條件等，因而理論曲線會出現個別點的跳躍，屬正常現象。從圖8可見理論與實際曲線吻合較好，表明解釋結果可靠。

圖8 ×井遺傳最優化計算結果的質量檢驗圖

5 結 論

用遺傳算法對砂礫巖儲層進行最優化測井解釋，充分利用了現有的測井信息，一次獲得較全面的儲層參數，且結果較為準確。缺點是在利用大量測井信息的同時，構建測井響應方程所需要確定的參數增多，為解釋帶來一定的難度。遺傳算法結構優化是下一步需要加強的部分。

［1］ 雷群，王紅巖，趙群，等.國內外非常規油氣資源勘探開發現狀及建議［J］.天然氣工業，2008，28（12）：7－10.

［2］ 張孝珍.砂礫巖儲層測井評價方法研究［D］.東營：中國石油大學，2009.

［3］ 馮國慶，陳軍，張烈輝，等.最優化測井解釋的遺傳算法實現［J］.天然氣工業，2002，22（6）：48－51.

［4］ 董勇，郭海敏，張閃.自適應混沌粒子群算法在油水兩相流動分析中的應用［J］.石油天然氣學報：江漢石油學院學報，2011，33（2）：81－84.

［5］ 潘保芝，閆桂京.遺傳最優化測井解釋方法及應用研究［C］∥中國地球物理學會.中國地球放理學會2001年刊：中國地球物理學會第十七屆年會論文集.2001.

［6］ 雍世和.最優化測井解釋［M］.東營：中國石油大學出版社，1996.

［7］ 田云英，夏宏泉.基于多礦物模型分析的最優化測井解釋［J］.西南石油學院學報，2006，28（4）：9－11.

GA－Optimal Log Interpretation Applied in Glutenite Reservoir Evaluation

HAN Xue，PAN Baozhi，ZHANG Yi，WANG Fei（College of Geo－exploration Science and Technology，Jilin University，Changchun，Jilin 130026，China）

As unconventional oil resource，glutenite reservoir gets more and more attention in oil exploration.However it is hard to determine lithology and set up the interpretation model accurately because the reservoir is deep burial，complex lithology，complex pore structures，low porosity and permeability and high heterogeneity.The optimal log interpretation can take full advantage of the logging information，measurement errors，logging responses equation errors，geological information and working experience so to find the most reasonable interpretation result from all the above possible solutions.In this paper，genetic algorithm（GA）is introduced to optimal log interpretation to get the reservoir parameters，such as porosity and saturation，finally the method is applied in Songliao basin and achieves the desired result.

log interpretation，glutenite reservoir，optimal log interpretation，genetic algorithm，multi－component model

P631.84

A

2012－05－03 本文編輯 李總南）

1004－1338（2012）04－0392－05

國家自然科學基金（41174096）和國家重大專項（2011ZX05009、2011ZX05044）聯合資助

韓雪，女，1985年生，碩士研究生，從事測井數據處理與解釋研究。