朱林奇，張 沖，魏 旸，郭 聰，周雪晴，陳雨龍

油氣資源與勘探技術教育部重點實驗室，長江大學地球物理與石油資源學院，武漢430100

基于改進雨林模糊神經(jīng)網(wǎng)絡模型的頁巖儲層總有機碳含量評價方法

朱林奇，張 沖*，魏 旸，郭 聰，周雪晴，陳雨龍

油氣資源與勘探技術教育部重點實驗室，長江大學地球物理與石油資源學院，武漢430100

由于采用常規(guī)測井曲線評價頁巖儲層總有機碳含量的精度不高，泛化能力不強，需要大量樣本。針對這些問題，改進了神經(jīng)網(wǎng)絡算法，以增加模型的預測能力。利用模糊系統(tǒng)優(yōu)化細胞神經(jīng)網(wǎng)絡結構，以增強其邏輯推理能力，提高其對模糊數(shù)據(jù)的敏感性；選擇能有效避免“虛擬碰撞”的雨林算法，并針對其存在的缺陷進行改進；利用改進雨林優(yōu)化算法對網(wǎng)絡的初始權值閾值進行優(yōu)化，避免網(wǎng)絡陷入局部極小。分析測井特征曲線的物理意義，選擇密度測井曲線與自然伽馬能譜測井曲線作為網(wǎng)絡的輸入，以總有機碳含量作為輸出，通過70塊巖心樣本網(wǎng)絡學習與26塊巖心樣本預測，證明了新網(wǎng)絡模型的優(yōu)越性。結果表明，新模型回判將相對誤差從23.189%減小到17.185%，預測相對誤差由52.421%減小到15.158%，具有更強的學習能力與泛化能力，更適用于頁巖儲層總有機質(zhì)含量的測井評價。

頁巖；總有機碳含量；模糊神經(jīng)網(wǎng)絡；改進的雨林算法；泛化能力

1 問題的提出

頁巖氣藏是典型的“自生自儲”油氣藏。氣體在頁巖中主要以游離態(tài)氣及吸附態(tài)氣的形式存在，其中游離氣存在于天然裂縫及微孔隙中，吸附氣吸附于干酪根和粘土表面。由于頁巖氣特殊的成藏方式，一些反映地層中碳含量高低的參數(shù)能大致表征地層中含氣量。國內(nèi)研究表明，總有機碳含量與吸附氣、游離氣含量均有較好的正相關關系（董大忠等，2012；姜福杰等，2012）。這是因為，隨著總有機碳含量的增高，干酪根含量增高，地層吸附氣體的也增強，導致吸附氣量增大。同時，頁巖的生排烴行為在一定程度上形成的多微孔以及微裂縫，改善頁巖儲層的儲集性能，而儲集性能的改善有利于游離氣含量的增長?？傊?，總有機碳含量的多少決定了頁巖氣產(chǎn)量的高低，頁巖氣要想達到商業(yè)開采的標準，必須達到總有機碳含量的最低標準。所以，總有機碳含量計算意義重大，其意義類似于常規(guī)油氣藏中的含烴飽和度，是非常重要的儲層參數(shù)。

總有機碳含量（TOC）與頁巖產(chǎn)氣率有著良好的線性正相關關系，在評價頁巖氣產(chǎn)氣率上有著自己獨特的優(yōu)勢（Ross and Bustin,2007;李延鈞等,2011）。國內(nèi)外研究顯示，利用元素俘獲能譜測井（ECS）、核磁共振測井（CMR）等技術輔助評價總有機碳量的精度較為可靠（Jacobi et al., 2002;Pemper etal.,2010;張曉玲等,2013;趙晨陽等,2002）。但由于測量價格昂貴，在所有井中均進行上述新技術的測量是不現(xiàn)實的，精度較高的測井新技術模型應用性不強。

現(xiàn)階段常規(guī)測井曲線評價總有機碳含量的方法主要有3類：

第1類為ΔlogR法及其改進方法，即利用電阻率與孔隙度曲線疊加預測總有機碳含量（郝建飛等,2012;劉承民,2012;張作清等,2013;鐘光海等2015）。存在的問題一是需要首先確定出成熟度的參數(shù)LOM，制約了該方法的使用，故現(xiàn)在多用其改進算法；二是logR模型的確定受到人為因素的干擾，因為logR的得到需要先確定細粒的非烴源巖作為基線，基線選擇的合適與否決定了logR的預測精度，在刻度確定的情況下，實際情況中有些井難以找到合適的基線，只能選擇較為合適的深度段作為基線，精度難以保證。

第2類方法為利用單一特征曲線或者多特征曲線進行擬合計算，如利用自然伽馬曲線與TOC建立統(tǒng)計關系，利用干酪根密度較小的概念所提出的補償密度曲線與TOC單因素擬合的方法，利用補償密度曲線、補償中子曲線、聲波時差曲線、深電阻率曲線進行多元擬合建模等（Rossand Bustin, 2008;馬林,2013;李軍等,2014,謝慶明等,2014）。這一類算法的優(yōu)勢是建模方法簡單，在精度要求不高的情況下實現(xiàn)較為方便。存在的問題有兩個：首先，在建模時由于缺少理論支撐，無法確定測井曲線與總有機碳含量的具體方程而僅能判斷具有正相關性或者負相關性，這樣的公式精度難以得到保證。其次，在總有機碳含量較小時，由于測井曲線所測得的信息為地層中各類巖石礦物及流體的綜合，總有機碳含量的信息必然會被干擾信息所壓制，導致其與測井曲線的相關性減弱。

第3類是基于上述問題，近幾年提出了神經(jīng)網(wǎng)絡算法預測總有機碳含量的方法，即利用神經(jīng)網(wǎng)絡強逼近函數(shù)的能力提升總有機碳含量計算的精度（熊鐳等,2014;柳筠,2015;孟召平等,2105）。該方法進一步提升了總有機碳含量計算的精度，使用較為廣泛。但傳統(tǒng)的神經(jīng)網(wǎng)絡存在著其固有缺點：第一，需要大量樣本以便統(tǒng)計其規(guī)律，使其具有良好的泛化能力，而因取心昂貴造價高，且實驗價格不菲，在實際情況下并不能進行大量的巖心實驗；第二，神經(jīng)網(wǎng)絡的權值與閾值初始化時存在著隨機性，使得結果容易陷入局部極小，這就使神經(jīng)網(wǎng)絡方法在運用中的可靠性大打折扣。

基于上述分析，筆者嘗試對神經(jīng)網(wǎng)絡進行改進，以解決其存在的問題，增加神經(jīng)網(wǎng)絡模型預測可靠性，提升總有機碳含量的測井評價精度。

2 總有機碳含量的測井響應

常規(guī)測井系列具有成本低的特點，較非常規(guī)測井系列應用更廣泛，故本文僅討論常規(guī)測井曲線與總有機碳含量的關系。利用焦石壩地區(qū)龍馬溪組某頁巖氣井96塊巖心進行相關性分析，發(fā)現(xiàn)總有機碳含量與密度測井曲線、自然伽馬能譜測井曲線（主要為鈾、鉀含量以及去鈾伽馬值）存在著一定的相關性（圖1）。

密度測井方法測量地層的體積密度，在總有機碳含量較高時，地層中的干酪根含量較多，而干酪根的密度范圍一般為1.1～1.4 g/cm3，故在富含干酪根的地層，密度測井值會存在明顯的回落。所以，總有機碳含量應與密度測井值關系密切（圖1a）。除了密度測井外，總有機碳含量還與去鈾伽馬值(KTH)、鈾元素(K)、鉀元素(U)具有良好的相關性（圖1b、c、d），這是因為干酪根中含有放射元素鈾，自然伽馬能譜測井中測量的鈾元素曲線值越高，說明干酪根含量越多，對應的總有機碳含量值應越高。

式(1)為自然伽馬能譜測井反推地層鈾、釷、鉀含量的方程。其中，A1、A2、A3、B1、B2、B3、C1、C2、C3為對應系數(shù)，在儀器出廠時就已經(jīng)刻度完畢，是一個定值。分析該方程組可知，在系數(shù)W1、W2、W3變化不大的情況下，若鈾含量較高，相應的釷、鉀的含量應較低。所以，去鈾伽馬、鉀測井曲線也應與總有機碳含量關系密切，且與鈾元素與總有機碳含量的相關性相反。

圖1 常規(guī)測井曲線響應與總有機碳含量相關關系Fig.1 Correlation between conventional log curvesand totalorganic carbon content

圖1顯示，雖然上述特征曲線響應與TOC含量之間存在著正相關或者負相關關系，但是這種關系僅能定性的對TOC進行分析，對于定量計算來說精度過低，尤其當TOC含量低時，呈現(xiàn)明顯的離散。這是因為常規(guī)測井曲線測量的是地層的綜合響應，在TOC含量不高時，其所表達的信息必定會被其他信息所壓制而難以顯現(xiàn)，有必要利用神經(jīng)網(wǎng)絡進行特征數(shù)據(jù)信息挖掘。

3 雨林模糊神經(jīng)網(wǎng)絡模型及其驗證

3.1 模糊神經(jīng)網(wǎng)絡

由于模糊系統(tǒng)與神經(jīng)網(wǎng)絡既有著功能上的共性又存在著相輔相成的理論特點，兩者越來越多的被有機的結合在一起使用。模糊神經(jīng)網(wǎng)絡繼承了神經(jīng)網(wǎng)絡的強樣本學習能力，又由于結合了模糊系統(tǒng)，存在著一定的邏輯推理能力，其方法較舊有的BP神經(jīng)網(wǎng)絡啟發(fā)性、魯棒性、透明性均更強，不僅能處理明確的數(shù)據(jù)，對模糊的數(shù)據(jù)同樣敏感（劉立峰等,2014;徐智浩等,2014）。

模糊神經(jīng)網(wǎng)絡通常由4層組成，分別為輸入層、模糊化層、模糊規(guī)則層以及去模糊化層，在網(wǎng)絡學習時以誤差平方和最小為學習目標，利用梯度法與誤差反向傳播法進行權值的更新。在多層神經(jīng)網(wǎng)絡下，層數(shù)越多，其誤差反向傳播速度越慢（嚴鴻和管燕平,2009）。這樣的后果一是使得網(wǎng)絡的學習速度變得異常緩慢，二是容易陷入局部極小值。所以，需要利用優(yōu)化算法優(yōu)化權值及閾值，以加快收斂速度，防止局部極小值的產(chǎn)生。

3.2 改進雨林算法

為了克服模糊神經(jīng)網(wǎng)絡存在的上述問題，通常在進行網(wǎng)絡訓練之前先進行優(yōu)化算法的尋優(yōu)，以確定較為合適的權值與閾值，輸入到模糊神經(jīng)網(wǎng)絡中進行迭代。筆者選用較新的雨林優(yōu)化算法，并對其進行改進，以加強其全局尋優(yōu)能力。

3.2.1 雨林算法（RFA）及其改進思路

雨林算法在2013年被提出（高維尚等，2013；高維尚和邵誠，2014）?？紤]到自然界各種動物群體均存在著“碰撞的現(xiàn)象”，而樹叢的生長并不存在碰撞，而是有機的向外延伸，吸收陽光，具有模擬動物群體算法所不具備的優(yōu)勢。在算法方面，該算法采用全局分區(qū)樣與局部分級采樣相結合、均勻采樣與非均勻采樣相結合的思路，以適應未知目標函數(shù)的復雜分布，有效避免虛擬碰撞，使算法更易尋找到全局最優(yōu)解（圖1）。

雨林算法雖較其他啟發(fā)式算法有著明顯的優(yōu)點，但也存在著一些不合理性。比如，其存在著在迭代后期收斂過快的問題，就有可能未能搜尋到全局最優(yōu)解而可能收斂到局部最優(yōu)解的凸集。所以，筆者嘗試對于雨林算法進行改進，以在最大程度上解決這個問題?？紤]到雨林算法的收斂速度主要決定于學習因子α，通過改進學習因子α使其在后期學習中收斂變慢，以較慢的速度、較多的迭代次數(shù)去更好的找到全局最優(yōu)解。

算法行為規(guī)劃所對應的動態(tài)方程組為：

其中，i、j、k、s分別為集合數(shù)、迭代次數(shù)、決策空間和播種節(jié)點的示意次數(shù)，n、r、nn為代理的數(shù)量、代理的范圍和保留播種節(jié)點量，α為學習因子。筆者將其改進為一與迭代次數(shù)j與最大迭代次數(shù)Tmax有關的函數(shù)，以達到前期快速確定可行域，后期緩慢在可行域中尋找全局最優(yōu)解的目的。

3.2.2 改進雨林算法實現(xiàn)步驟

（1）初始樹木播種。在優(yōu)化算法中，需要對初始尋優(yōu)點進行初始化采樣。其采樣的合理性決定了算法是否能尋找到全局最優(yōu)點。在雨林算法中，采樣以一種近似均勻的方式進行，來確保有一個以上的采樣點進入關鍵域內(nèi)，方式類似于樹木播種。這種采樣方式會使得可行域中留下的最大空白最小。在算法的可行域中以陣列的方式播種，并使得該種子形成根節(jié)點進行生長。

（2）萌發(fā)。依據(jù)信息熵算法計算播種后的信息熵。雨林算法中信息熵計算方法為：

式中，t代表迭代的數(shù)目；H為采樣的信息熵。

依據(jù)信息熵、節(jié)點適應度計算采樣價值樹信息價值與采樣價值變化率，以確定下一次迭代中每個播種節(jié)點周圍新生采樣點數(shù)量值以及伸展的范圍，信息增益越大，伸展的范圍越大，采樣點數(shù)量越多。

（3）生長。依據(jù)萌發(fā)過程得到的伸展范圍以及采樣點數(shù)，不斷更新進行繼續(xù)采樣，并在更新的同時計算信息增益，以確定伸展范圍及新生采樣點數(shù)量。

（4）競爭。進行節(jié)點的更新，僅保留若干個較為優(yōu)勢的節(jié)點作為下次播種的節(jié)點，以增加優(yōu)勢節(jié)點的生長能力，使較優(yōu)節(jié)點獲得良好的生長的空間。

（5）計壽。若樹木在最近的幾次生長中適應度增加，則認為為較優(yōu)樹，算法將該雨林節(jié)點保存，并跳過播種階段，直接進行萌發(fā)。

（6）繁衍。利用優(yōu)化結果與實際結果的差異以及迭代次數(shù)去衡量程序是否滿足了結束的條件。若未滿足條件，則繼續(xù)第一步的播種階段，在播種時播種更多的樹，以便于尋找到全局最優(yōu)解。

最終，通過雨林算法的優(yōu)化，可得到較為可靠的神經(jīng)網(wǎng)絡初始權值與閾值。對應的改進雨林算法模糊神經(jīng)網(wǎng)絡流程見圖2。

3.2.3 方法驗證

利用前人方法建立模型，以便于進行方法與精度的對比。利用BP神經(jīng)網(wǎng)絡模型，以前文所確定的4種特征曲線為輸入，以總有機碳含量作為輸出，建立4×10××10×10×1的神經(jīng)網(wǎng)絡模型，其中建立1個隱層是因為1隱層結構神經(jīng)網(wǎng)絡即具有逼近一切函數(shù)的能力（范佳妮等，2005）。取96塊巖心中的70塊進行建模，26塊進行模型檢測。最后得到的神經(jīng)網(wǎng)絡模型學習結果以及預測結果（圖3）。

分析圖3可以發(fā)現(xiàn)，神經(jīng)網(wǎng)絡預測TOC的精度較常規(guī)擬合方法精度更高，在精度要求不高的情況下可以達到預測TOC含量的目的。但是，對比學習樣本點回歸精度與未參與建模樣本預測結果，模型預測結果要差于回歸結果，這說明神經(jīng)網(wǎng)絡本身固有的泛化能力較弱，對于未知樣本點的預測穩(wěn)定性不高，不確定性強。值得一提的是，由于初始權值與閾值的隨機性，神經(jīng)網(wǎng)絡極易陷入局部極小，筆者在建立模型時反復運行了27次，以確定最優(yōu)神經(jīng)網(wǎng)絡模型。這種帶有強人為干擾的模型并不利于TOC含量預測，因為并不了解建立的模型是否為最優(yōu)模型。

圖2 改進雨林模糊神經(jīng)網(wǎng)絡流程Fig.2 The processof im proving the fuzzy neuralnetw ork w ith the im proved rain forestalgorithm

圖3 BP神經(jīng)網(wǎng)絡建模精度與預測精度Fig.3 Modeling precision and prediction accuracy of BPneuralnetw ork

利用筆者所提出的改進雨林模糊神經(jīng)網(wǎng)絡進行TOC含量預測。建模樣本與預測樣本與上述BP神經(jīng)網(wǎng)絡模型建立樣本相同，以便于比較兩模型在建模精度與泛化能力上的優(yōu)劣。預測時進行了雨林算法與改進雨林算法的權值閾值優(yōu)化比較，結果見圖4。從優(yōu)化結果中可以看出，改進的雨林算法（IRFA）雖然在迭代速度上不如原有雨林算法（RFA），但是在尋優(yōu)能力上強于雨林算法，考慮到改進的雨林算法優(yōu)化神經(jīng)網(wǎng)絡也僅僅運行了426.32 s，認為筆者提出的改進的雨林算法應更適應于函數(shù)的尋優(yōu)，其尋優(yōu)效果更明顯。

圖5為改進雨林模糊神經(jīng)網(wǎng)絡建模精度與預測精度圖。對比圖3可以看出，筆者提出的方法在建模精度與預測精度上明顯優(yōu)于BP神經(jīng)網(wǎng)絡，尤其是預測精度有較大改善，體現(xiàn)了良好的泛化能力。計算兩種模型平均相對誤差，運用筆者提出的模型將回歸70個樣本的相對誤差從23.189%降低到17.185%，預測相對誤差由52.421%降低到15.158%，建模效果大為改善。這僅僅是運行了1次的結果，并未通過反復運行以確定最優(yōu)模型。并且，由該相對誤差數(shù)據(jù)可以看出，回歸精度反而大于預測精度，說明該模型未存在傳統(tǒng)BP神經(jīng)網(wǎng)絡極易存在的過擬合現(xiàn)象，利用筆者提出的方法預測TOC含量時，在樣品點達到70個的情況下即可得到較為合理的模型。

圖4 改進雨林算法與雨林算法尋優(yōu)對比Fig.4 Comparison of the improved rain forestalgorithm and rain forestalgorithm in optim ization

圖5 改進雨林模糊神經(jīng)網(wǎng)絡建模精度與預測精度Fig.5 Modeling precision and prediction accuracy of fuzzy neuralnetwork optim ized by the improved rain forestalgorithm

4 實例分析

利用改進模型進行某井段總有機碳含量進行評價，結果見圖6。其中，第6道為上述BP神經(jīng)網(wǎng)絡預測模型所評價出來的總有機碳含量曲線，第7道為筆者提出的改進雨林模糊神經(jīng)網(wǎng)絡模型評價得到的總有機碳含量曲線。

對比顯示，BP神經(jīng)網(wǎng)絡模型預測曲線總體波動較大，并且在2354～2359m深度段，預測結果在精度和趨勢上均嚴重不符合實際情況。這是因為BP神經(jīng)網(wǎng)絡對輸入數(shù)據(jù)過于敏感，U測井曲線在建模樣本點中并不存類似該深度段的較大值。對于這種問題，BP神經(jīng)網(wǎng)絡只有增加樣本點的數(shù)量以增加精度，若利用該模型將會得到不準確的評價認識。

筆者模型預測曲線波動較小，且預測精度與趨勢均較符合實際情況，說明模型具有一定的邏輯思維，尤其在2354～2359m深度段，在建模時未出現(xiàn)類似高鈾值樣本的情況下，其也能得到較好的結果，體現(xiàn)了良好的泛化能力，優(yōu)于原始BP神經(jīng)網(wǎng)絡模型預測結果，對總有機碳含量的評價有所幫助。

5 結論

（1）依據(jù)巖心與常規(guī)測井曲線數(shù)據(jù)，確定了常規(guī)測井曲線與總有機碳含量的相關性及其物理意義，并認為在總有機碳含量較低的情況下，常規(guī)測井中所反映的總有機碳含量信息被地層中其他信息所干擾覆蓋，使得相關性變差，故需要借助神經(jīng)網(wǎng)絡的數(shù)據(jù)挖掘能力進行總有機碳含量的求取。

（2）考慮到神經(jīng)網(wǎng)絡固有的問題，借助模糊系統(tǒng)的邏輯推理的能力，應用模糊系統(tǒng)優(yōu)化神經(jīng)網(wǎng)絡結構，提高其對模糊數(shù)據(jù)的敏感程度；嘗試使用避免“虛擬碰撞”的雨林算法并嘗試將其改進，以優(yōu)化神經(jīng)網(wǎng)絡的初始權值閾值，避免其陷入局部極小。

（3）利用70個樣本點對未改進的BP神經(jīng)網(wǎng)絡與改進雨林模糊神經(jīng)網(wǎng)絡分別建模，并利用26個樣本點進行預測。結果顯示，筆者提出的改進雨林模糊神經(jīng)網(wǎng)絡模型精度較傳統(tǒng)的BP神經(jīng)網(wǎng)絡精度更高，并發(fā)現(xiàn)建模樣本在70個以上時就可得到未過擬合模型。實例分析結果顯示，筆者提出的模型更合理，適用于頁巖儲層總有機質(zhì)含量的測井評價。

圖6 S井某井段總有機碳含量預測對比Fig.6 Comparison of the predicted TOC in a section of Swell

（References）：

董大忠，鄒才能，楊樺，等.2012.中國頁巖氣勘探開發(fā)進展與發(fā)展前景[J].石油學報，33(S1):107-114.

范佳妮，王振雷，錢鋒.2005.BP人工神經(jīng)網(wǎng)絡隱層結構設計的研究進展[J].控制工程，12(S0):105-109.

高維尚，邵誠.2014.復雜非凸約束優(yōu)化難題與迭代動態(tài)多樣進化算法[J].自動化學報，40(11):2469-2479.

高維尚，邵誠，高琴.2013.群體智能優(yōu)化中的虛擬碰撞：雨林算法[J].物理學報，62(19):28-43.

郝建飛，周燦燦，李霞，等.2012.頁巖氣地球物理測井評價綜述[J].地球物理學進展，27(4):1624-1632.

姜福杰，龐雄奇，歐陽學成，等.2012.世界頁巖氣研究概況及中國頁巖氣資源潛力分析[J].地學前緣，19(2):198-211.

李軍，路箐，李爭，等.2014.頁巖氣儲層“四孔隙”模型建立及測井定量表征方法[J].石油與天然氣地質(zhì)，35(2):266-271.

李延鈞，劉歡，劉家霞，等.2011.頁巖氣地質(zhì)選區(qū)及資源潛力評價標準[J].西南石油大學學報（自然科學版），33(2):28-34.

劉承民.2012.頁巖氣測井評價方法及應用[J].中國煤炭地質(zhì)，24(8): 77-79.

柳筠.2015.涪陵地區(qū)頁巖氣測井評價參數(shù)研究[J].江漢石油科技，25 (3):34-41.

劉立峰，孫贊東，韓劍發(fā)，等.2014.量子粒子群神經(jīng)網(wǎng)絡碳酸鹽巖流體識別方法研究[J].地球物理學報，57(3):991-1000.

孟召平，郭彥省，劉尉.2015.頁巖氣儲層有機碳含量與測井參數(shù)的關系及預測模型[J].煤炭學報，40(2):247-253.

馬林.2013.頁巖儲層關鍵參數(shù)測井評價方法研究[J].油氣藏評價與開發(fā)，3(6):66-71.

謝慶明，程禮軍，劉俊峰，等.2014.渝東南黔江地區(qū)龍馬溪組頁巖氣儲層測井解釋評價研究[J].地球物理學進展，29(3):1312-1318.

熊鐳，張超謨，張沖，等.2014.A地區(qū)頁巖氣儲層總有機碳含量測井評價方法研究[J].巖性油氣藏，26(3):74-83.

徐智浩，李勝，張瑞雷，等.2014.基于LuGre摩擦模型的機械臂模糊神經(jīng)網(wǎng)絡控制[J].控制與決策，29(6):1097-1102.

嚴鴻，管燕平.2009.BP神經(jīng)網(wǎng)絡隱層單元數(shù)的確定方法及實例[J].控制工程，16(S1):100-102.

張曉玲，肖立志，謝然紅，等.2013.頁巖氣藏評價中的巖石物理方法[J].地球物理學進展，28(4):1962-1974.

張作清，鄭煬，孫建孟.2013.頁巖氣評價“六性關系”研究[J].油氣井測試，22(1):66-74.

趙晨陽，杜禹，蔡振東，等.2015.國外頁巖氣儲層測井評價技術綜述[J].遼寧化工，44(4):473-478.

鐘光海，謝冰，周肖.2015.頁巖氣測井評價方法研究—以四川盆地蜀南地區(qū)為例[J].巖性油氣藏，27(4):96-102.

Jacobi D,Breig J,LeCompte B,et al.2002.Effective geochemical and geomechanical characterization of shale gas reservoirs from the Wellboreenvironment:Caney and theWoodford Shale[A].SPE124231.

Pemper R，Han X G，Mendez F，et al.2010.The direct measurement of carbon in wells containing oil and natural gas using a pulsed neutron mineralogy Ttool[A].SPE133128.

Ross D JK and Bustin R M.2007.Shale gas potential of the Lower Jurassic Gordondale Member，northeastem British Colunbia，Canada[J].AAPG Bulletin，55(1):51-75.

Ross D JK and Bustin R M.2008.Characterizing the shale gas resource potential of Devonian-Mississippian strata in the Western Canada sedimentary basin:Application ofan integrated formation evaluation[J]. AAPGBulletin，92(1):87-125.

TheM ethod for TOC Content Evaluation in Shale Reservoirs Based on Im p roved Rain Forest Fuzzy NeuralNetwork M odel

ZHU Linqi,ZHANG Chong*,WEIYang,GUO Cong,ZHOU Xueqing,CHEN Yulong
Key Laboratory of Exploration Technologies for Oiland Gas Resources,Ministry of Education,Geophysics and Oil Resource Institute,Yangtze University,Wuhan 430100,China

Theaccuracy ofevaluating totalorganic carbon in shale reservoirs is limited by using conventional logging curvesbecauseof their insufficient generalization ability and requirement of a large number of samples.In view of these problems,neural network algorithm was improved to improve the prediction ability of themodel.The cellular neural network structure was optim ized by using a fuzzy system to enhance its logical reasoning ability and to improve its sensitivity to fuzzy data.The rain forest algorithm,which can effectively avoid the virtual collision,was selected,and its defect of slow convergence in the late learning was overcome.The initial weight value and threshold value of the network were optimized by the improved rain forest optimization algorithm to prevent thenetwork from resulting in localminimum,which can improve theaccuracy and generalization ability of themodel.Based on the analysis of the physicalmeaning of the characteristic curve,the density log curves and the natural gamma ray spectrum logging curves were chosen as the input to the network and the total organic carbon contentwas used as the output.Through the networklearning of 70 core samples and the prediction of 26 core samples,the role of the improved rain forest algorithm and fuzzy logic is proved.The superiority of thenew networkmodel isdemonstrated.The resultshows that the relative regressionalerror of thenewmodel is reduced from 23.189%to 17.185%,and the relative prediction error is reduced from 52.421%to 15.158%,whichmeans that the prediction is in accordance with the real situation of formation.From the above,we learn that the newmodel has better learning ability and generalization ability.Thenewmodel ismore suitable for loggingevaluation of totalorganicmatter content in shale reservoirs.

shalegas；TOC；fuzzy neuralnetwork；improved rainforestalgorithm；generalization

ZHANGChong，Associate Professor；E-mail:yzlogging@163.com

P631.82

A文獻標識碼：1006-7493（2016）04-0716-08

10.16108/j.issn1006-7493.2016079

2016-05-31；

2016-06-04

國家自然科學基金重點項目（41404084）；湖北省自然科學基金（2013CFB396）聯(lián)合資助

朱林奇，男，1993年生，碩士，地球探測與信息技術專業(yè);E-mail:445364694@qq.com

*通訊作者：張沖，副教授;E-mail:yzlogging@163.com