地層元素測井儀器結(jié)構(gòu)參數(shù)的蒙特卡羅數(shù)值模擬

嚴(yán)慧娟，岳愛忠，趙均，張建民，王樹聲，林巖棟

（1.西安交通大學(xué)，陜西西安710049；2.中國石油集團(tuán)測井有限公司技術(shù)中心，陜西西安710077）

地層元素測井儀器結(jié)構(gòu)參數(shù)的蒙特卡羅數(shù)值模擬

嚴(yán)慧娟1，岳愛忠2，趙均1，張建民1，王樹聲2，林巖棟2

（1.西安交通大學(xué)，陜西西安710049；2.中國石油集團(tuán)測井有限公司技術(shù)中心，陜西西安710077）

建立地層元素測井的蒙特卡羅方法數(shù)值模擬模型，分別對儀器的屏蔽體、硼套、源距等儀器參數(shù)對測井響應(yīng)的影響進(jìn)行數(shù)值模擬研究。單一屏蔽體材料中鎢鎳鐵屏蔽率最高，組合材料中鎢鎳鐵＋硼砂和鎢鎳鐵＋聚乙烯的屏蔽率比相同厚度鎢鎳鐵的屏蔽率略高一些，鎢鎳鐵＋聚乙烯的效果更明顯。實際儀器模型中11cm厚的鎢鎳鐵屏蔽率大于0.9；硼套為10B與氟橡膠混合材料時其10B含量達(dá)到10%以上，達(dá)到降低儀器內(nèi)部元素影響的效果；儀器外殼材料均可以采用沉淀不銹鋼或銅，其對地層元素的解譜幾乎沒有影響。這些模擬結(jié)果可以為優(yōu)化儀器結(jié)構(gòu)參數(shù)提供參考。

地層元素測井；蒙特卡羅方法；數(shù)值模擬；儀器結(jié)構(gòu)

0 引 言

地層元素測井通過測量中子在地層中產(chǎn)生的俘獲伽馬能譜和非彈性散射伽馬能譜獲取地層巖石骨架的主要元素種類和含量，確定巖性和骨架參數(shù)。最新的地層元素測井能夠確定Si、Ca、Fe、S、Ti、Gd、 Mg、K、Mn、Al、H、Cl、C和O等14種元素的份額以及其中前10種元素的含量。在國外，30多年前出現(xiàn)自然伽馬能譜和次生伽馬能譜測井的時候，就開始研究地層元素測井技術(shù)，當(dāng)時稱之為地球化學(xué)測井。經(jīng)過多年的發(fā)展，直到20世紀(jì)末才逐漸推出成型儀器。斯倫貝謝公司于1996年推出元素俘獲譜測井儀器ECS［1］；哈里伯頓公司于2009年推出了一款類似儀器GEMTM［2］，這2款儀器都使用同位素中子源；貝克休斯公司于2009年推出了一款使用可控中子源的地層元素測井儀器FLeXTM［3］。

隨著復(fù)雜巖性油氣藏、頁巖氣、致密氣等非常規(guī)油氣藏勘探開發(fā)的深入，地層元素測井的作用得到越來越多重視，應(yīng)用日益廣泛。在頁巖氣和致密氣勘探中，地層元素測井已經(jīng)成為必測項目。地層元素測井可以直接、準(zhǔn)確地確定巖性，這是其他測井方法無法比擬的。

大慶、勝利、華北、青海、克拉瑪依和塔里木等油田利用國外儀器進(jìn)行了測井試驗和應(yīng)用研究，在火山巖等復(fù)雜巖性識別與解釋評價方面進(jìn)行了有益的探討［4－8］。但限于國外儀器應(yīng)用的諸多不便，地層元素測井在國內(nèi)并未得到規(guī)模應(yīng)用。2011年，中國石油天然氣集團(tuán)公司支持立項研制地層元素測井儀器。本文就儀器設(shè)計部分的研究成果作一介紹。

1 蒙特卡羅計算模型的建立

在調(diào)研國外同類儀器后提出了地層元素測井儀器初步的幾何結(jié)構(gòu)和參數(shù)（包括源結(jié)構(gòu)、源距、探測器、儀器外徑和長度等）以及材料的類型等，建立儀器測井的蒙特卡羅方法數(shù)值計算模型（見圖1）。測井模型基本參數(shù)中地層高度為130cm，直徑為150 cm，井眼半徑為10cm；儀器分為2段，2段尺寸分別為Φ12.7cm×70cm和Φ9cm×45cm，總長度為115cm；源為Am－Be中子源，屏蔽體材料為鎢鎳鐵；探測器材料為BGO晶體；硼套材料為45%的硼粉打混氟橡膠。

圖1 地層元素測井儀器蒙特卡羅方法數(shù)值計算模型

2 儀器屏蔽體的屏蔽性能

屏蔽體是為防止中子源產(chǎn)生的中子直接進(jìn)入儀器探測器影響探測結(jié)果而設(shè)置的。其屏蔽性能是屏蔽材料選取和確定厚度的唯一目標(biāo)。屏蔽性能一般用屏蔽率表示，描述了屏蔽體阻擋粒子穿過的能力。

2.1 理論計算模型

為測試屏蔽體的屏蔽效果，首先參照ISO技術(shù)標(biāo)準(zhǔn)［9］建立镅鈹中子源，設(shè)置模型最左側(cè)的垂直面為面源，向右水平發(fā)射中子，經(jīng)過10cm厚的空氣層進(jìn)入屏蔽體，分別記錄中子經(jīng)過屏蔽體前后的中子流，最后獲得屏蔽材料的屏蔽率。屏蔽體材料分別取鎢鎳鐵、鐵、鉛、銅、石墨、硼砂和聚乙烯7種。各種材料屏蔽率隨厚度的變化規(guī)律見圖2。

圖2 不同屏蔽材料的屏蔽率隨厚度的變化規(guī)律

由圖2可以很明顯看出，隨著屏蔽體厚度的增加，各種材料的屏蔽率都逐漸增大，其中鎢鎳鐵是所選幾種材料中屏蔽效果最好的，并且在厚度為14 cm時屏蔽率達(dá)到0.80，達(dá)到較好的屏蔽效果。

考慮到隨著屏蔽層的加厚，熱中子會逐漸增加。如果采用組合材料，即先用一定厚度的鎢鎳鐵，再加一段其他材料，可能會取得更好的屏蔽效果。于是，計算了10cm厚度鎢鎳鐵加3cm及5cm厚度其他材料作屏蔽體的情況，分別與屏蔽體是鎢鎳鐵13 cm和15cm時屏蔽率進(jìn)行比較，得到結(jié)果見表1。

表1 相同厚度情況下不同材料組合的屏蔽率比較

由表1可以看出，在所有組合材料中只有鎢鎳鐵＋硼砂和鎢鎳鐵＋聚乙烯的屏蔽率比相同厚度鎢鎳鐵的屏蔽率略高一些，而其他組合材料反比鎢鎳鐵的屏蔽率要低。

2.2 儀器模型

儀器結(jié)構(gòu)如圖1所示，在源和屏蔽體之間是沉淀不銹鋼，它對中子也起到一定的屏蔽作用。從實際儀器設(shè)計考慮在源和屏蔽體之間也可以是空腔或其他材料。為了分析這2種情況下屏蔽體的屏蔽效果，建立了2種計算模型（見圖3）。圖3（a）所示為源和屏蔽體之間是沉淀不銹鋼的局部模型圖；圖3（b）所示為源和屏蔽體之間是空腔的局部模型圖。模擬計算得到2種情況下屏蔽率隨屏蔽體厚度變化的規(guī)律（見圖4）。

由圖4可以看出，屏蔽體厚度4cm情況下源和屏蔽體之間是空氣時的屏蔽率為0.63；源和屏蔽體之間是沉淀不銹鋼時，由于鐵對中子的慢化的作用，中子入射面的中子能量會相應(yīng)降低，使得屏蔽率已經(jīng)達(dá)到0.82。所以儀器內(nèi)部的鐵的影響不可忽略。在實際儀器中源和屏蔽體之間會填充輕型的慢化材料，屏蔽效果應(yīng)該介于中間材料是空腔和沉淀不銹鋼之間，厚度11cm的鎢鎳鐵可以達(dá)到屏蔽率0.9以上，完全能夠滿足實際儀器要求。

3 硼套對探測器響應(yīng)的影響

硼套是為了防止儀器外部的熱中子進(jìn)入儀器與儀器材料發(fā)生反應(yīng)影響探測結(jié)果而設(shè)置的。在硼套厚度一定時，將10B與氟橡膠按照不同質(zhì)量比例進(jìn)行混合，使硼的濃度分別為0.3%、2.5%、10%、30%和50%，混合后的硼表面密度見表2。計算沒有硼套情況進(jìn)行比較，得到的能譜曲線見圖5。

表2 不同混合比例下硼的表面密度

圖5 不同硼混合比例能譜比較

由圖5可看出，在無硼套到硼混合比例逐漸增大時，計數(shù)率逐漸降低，而且硼混合比例低時變化幅度大，硼混合比例高時變化幅度小。當(dāng)硼混合比例達(dá)到10%之后隨著硼混合比例增加曲線變化幅度非常小，可以近似認(rèn)為硼混合比例為10%時接近飽和程度，即再增加硼混合比例，計數(shù)率變化不大，鐵峰也沒有顯著的變化。硼混合比例10%的情況下，就可以較好地滿足屏蔽儀器內(nèi)部元素影響的效果。

4 源距對探測器響應(yīng)的影響

根據(jù)儀器的基本結(jié)構(gòu)，模擬計算中源距分別為41、43、45、47、49、51cm和53cm，得到總計數(shù)率隨不同源距的變化曲線及不同源距探測器計數(shù)能譜曲線見圖6。由圖6可知，探測器總計數(shù)率隨源距增加而減??；探測器計數(shù)能譜曲線也隨源距增加而降低，且各主要特征峰幾乎同幅度降低。源距增大，可以減小井眼的影響，但考慮到隨著源距增大計數(shù)率不斷減低，過低的計數(shù)率不利于儀器測量，所以儀器的源距不能過大。另外，儀器源距的選擇還要綜合考慮儀器的探測深度等因素。

圖6 不同源距計數(shù)率變化曲線圖

5 儀器外殼材料對探測器響應(yīng)的影響

儀器外殼材料可以取沉淀不銹鋼或銅。這里分別計算了2種材料對探測器響應(yīng)影響規(guī)律，計算結(jié)果的能譜曲線見圖7。

圖7 外殼材料為沉淀不銹鋼和銅的能譜曲線

由圖7可以看出，外殼為沉淀不銹鋼和銅的差別除在鐵峰處差異稍大外，在其他能量區(qū)間相對偏差均在15%以內(nèi)，對各主要峰值影響很小，對地層元素的解譜幾乎沒有影響。所以，外殼材料的選取采用沉淀不銹鋼或銅均可以。

6 結(jié) 論

（1）在所選的單一屏蔽體材料中，鎢鎳鐵是屏蔽率最高的；在選定的組合材料中，鎢鎳鐵＋硼砂和鎢鎳鐵＋聚乙烯的屏蔽率比相同厚度鎢鎳鐵的屏蔽率略高一些，鎢鎳鐵＋聚乙烯的效果更明顯，所以屏蔽體的材料可以選取鎢鎳鐵和聚乙烯的組合材料；在實際儀器模型中，11cm厚的鎢鎳鐵可以達(dá)到屏蔽率大于0.9的要求。

（2）硼套為10B與氟橡膠混合材料時，其中的10B含量達(dá)到10%以上，即可以達(dá)到屏蔽儀器內(nèi)部元素影響的效果。

（3）隨著源距的增大，各能量段的計數(shù)率及總計數(shù)率都隨之減小。源距的選取要同時綜合考慮其他因素。

（4）儀器外殼材料采用沉淀不銹鋼或銅均可以，對地層元素的解譜幾乎沒有影響。

［1］ Barson D，Christensen R，Decoster E，et al.Spectroscopy：The Key to Rapid，Reliable Petrophysical Answers［J］.Oilfield Review，2005：14－33.

［2］ Galford J，Quirein J，Shannon S，et al.Field Test Results of a New Neutron－Induced Gamma－Ray Spectroscopy Geochemical Logging Tool［C］∥New Orleans，Louisiana.2009：22.

［3］ Xiaogang Han，Richard Pemper，Teresa Tutt，et al.Environmental Corrections and System Calibration for a New Pulsed－Neutron Mineralogy Instrument［C］∥SPWLA 50th Annual Logging Symposium，June 21－24，2009.

［4］ 袁祖貴，楚澤涵.一種新的測井方法（ECS）在王莊稠油油藏中的應(yīng)用［J］.核電子學(xué)與探測技術(shù)，2003，23（5）：417－423.

［5］ 魏國，趙佐安.元素俘獲譜（ECS）測井在碳酸鹽巖中的應(yīng)用探討［J］.測井技術(shù)，2008，32（3）：285－288.

［6］ 劉緒綱，孫建孟，郭云峰.元素俘獲譜測井在儲層綜合評價中的應(yīng)用［J］.測井技術(shù)，2005，29（3）：236－239.

［7］ 譚鋒奇，李洪奇，姚振華，等.元素俘獲譜測井在火山巖儲層孔隙度計算中的應(yīng)用 ［J］.國外測井技術(shù)，2008，168：27－30.

［8］ 劉應(yīng)華，趙軍，胡洪，等.ECS在塔里木復(fù)雜地層中的應(yīng)用［J］.國外測井技術(shù)，2008，168：21－23.

［9］ ISO.Reference Neutron Radiations——Part 1：Characteristics and Methods of Production［S］.Geneva，Schweiz.2001.

Monte Carlo Numerical Simulation of Structure Parameters in Formation Elements Logging Tool

YAN Huijuan1，YUE Aizhong2，ZHAO Jun1，ZHANG Jianmin1，WANG Shusheng2，Lin Yandong2
（1.Xi’an Jiaotong Univeristy，Xi’an，Shaanxi 710049，China；2.Technical Center，China Petroleum Logging CO.LTD.，Xi’an，Shaanxi 710077，China）

The Monte Carlo simulation model of formation elements logging tool is established.With the model the effects of tool parameters including the shield，boron coating and source spacing are studied on logging responses.In the single shield materials，the shielding effect of Tungsten（W），Nickel（Ni）and Iron（Fe）is the best.In the combined shield materials，the shielding effect of W－Ni－Fe＋Boron（B）and W－Ni－Fe＋polyethylene is better then that of W－Ni－Fe with the same thickness；the shielding effect of W－Ni－Fe＋polyethylene is much better.In practical tool model，the shielding effect of 11cm thick of W－Ni－Fe material is above 0.9；if the Boron coating is the mixed material of10B and fluoro－elastomer，its10B content is more than 10%，which may reduce the tool’s internal elements influencing effect.The tool shells may use settling stainless steel or Copper（Cu），because they don’t have any influence on spectrum deconvolution of formation elements.The above results may be used as a reference for optimizing the tool structure design.

formation elements logging，Monte Carlo method，numerical simulation，tool structure

P631.84 文獻(xiàn)標(biāo)識碼：A

2011－12－28 本文編輯 李總南）

齊德鵬，男，工程師，從事射孔技術(shù)管理工作。