张　丽　于华伟（山东科技大学 资源与土木工程系　泰安 709）（中国石油大学（华东）地球科学与技术学院　青岛 66580）



D-D次生源计数对密度测量的影响

张丽1于华伟2
1（山东科技大学 资源与土木工程系泰安 271019）
2（中国石油大学（华东）地球科学与技术学院青岛 266580）

摘要为研究D-D源密度测量过程中俘获γ计数对密度测量的影响，本文首先分析D-D源密度测量的原理，其次研究D-D源与地层反应产生的γ射线类型，并采用蒙特卡罗模拟方法对比分析了D-D源和D-T源密度测量过程中次生γ强度高低。结果表明，D-D源密度测量过程中次生γ以俘获γ为主；尽管D-D源的中子产额远低于D-T源，但在源距大于50 cm后低孔隙地层中D-D源俘获γ计数高于D-T源非弹γ计数；此外，高孔隙地层中D-D源俘获γ计数也可以通过调整源距大小满足密度测量的要求。因此，使用D-D源测量的次生γ计数能够满足密度测量的统计精度要求。

关键词D-D源，俘获γ，密度测量，蒙特卡罗模拟

山东科技大学人才引进科研启动基金(No.2014RCJJ040)、国家自然科学基金(No.41304095)、山东省自然科学基金(No.ZR2013DQ012)、青岛市科

技计划(No.14-2-4-65-jch)、中国石油大学教学改革项目(No.SY-A201405)资助

第一作者：张丽，女，1984年出生，2013年于中国石油大学（华东）获博士学位，讲师，主要研究方向为核测井蒙特卡罗模拟研究

Supported by Talent Introduction and Scientific Research Initiating Fund of Shandong University of Science and Technology(No.2014RCJJ040),National Natural Science Foundation of China(No.41304095),Shandong Provincial Natural Science Fund(No.ZR2013DQ012),Qingdao City Science and Technology Program(No.14-2-4-65-jch),Teaching Reform Project of China University of Petroleum(No.SY-A201405)

First author:ZHANG Li,female,born in 1984,graduated from China University of Petroleum(East China)with a doctor’s degree in 2013,lecturer,focusing on the simulation of Monte Carlo and nuclear logging

Effect of density measurement on D-D induced gamma counting

ZHANG Li1YU Huawei2
1(Department of Resources and Civil Engineering,Shandong University of Science and Technology,Tai’an 271019,China)2(College of Geo-Resources and Information,China University of Petroleum,Qingdao 266580,China)

AbstractBackground:The controllable neutron source is used to replace the chemical gamma source to measure the formation density in oil exploration.Purpose:The effect of capture gamma counting on density is studied in the process of D-D source density measurement.Methods:Firstly,the principle of D-D source density measurement is analyzed,and the type of gamma-ray source in the process of D-D source and the formation is studied.In addition,the intensity of induced gamma between D-D and D-T source density measurement is compared using Monte Carlo simulation method.Results &Conclusion:The capture gamma takes a large proportion on D-D induced gamma density measurement.Although neutron yield of D-D source is far lower than that of D-T source,the capture gamma counting of D-D source is higher than inelastic gamma counting of D-T source with the distance from the source more than 50 cm in the low porosity formation.In addition,the capture gamma counting can be obtained by adjusting the source size to meet the requirements of density measurement in high porosity formation.Therefore,the counting of D-D induced gamma can meet the statistical accuracy of the density measurement.

Key wordsD-D source,Capture gamma,Density measurement,Monte Carlo simulation

可控中子源替代传统的化学放射源进行密度测井成为核测井发展的必然趋势[1-3]。它克服了传统密度测井中放射源辐射大、操作困难等缺点，成为当前密度测井的关键技术[4-5]。1995年，国外油服公司（如斯伦贝谢）开始研究使用D-T可控中子源的密度测井方法及仪器，并在随后的几十年里将该项技术应用于实际油田作业。但是D-T源寿命短（国内产品大约几十到300 h）、成本高、使用的放射性氚靶存在潜在危险[6-8]。与D-T源相比，D-D可控中子源具有寿命长（超过1000 h）、成本低，而且由于不使用放射性氚气、安全性更高等优点[9]。为此，Badruzzaman等[10-12]提出应用D-D源进行密度测量的方法，但是他们认为由于D-D源的能量较低且中子产额较D-T源低一个量级，而使中子在地层中次生的γ射线强度太低，以致于测量γ射线的统计误差较大，不足以满足测井仪器的计数精度要求，因此没有对D-D源次生γ在地层中的强度分布做进一步的研究。

本文着重论证D-D源测量地层的可行性问题。在分析D-D源与地层相互作用机理、明确D-D源密度测量过程中次生γ类型的基础上，对比D-D源与D-T源在不同地层中次生γ强度的高低，通过探讨基于D-D次生源开展密度测量的可行性，在此过程中对次生γ源的影响因素进行归纳，为开展下一步方法研究提供依据。

1　可控中子源密度测量原理

在可控源密度测量过程中，脉冲中子发生器发出的快中子与地层原子发生反应，其反应过程有弹性散射、非弹性散射、热中子俘获，而在反应时产生γ射线的过程只有非弹性散射和热中子俘获。该γ射线在输运中发生衰减，衰减规律与传统γ-γ密度测井机理相同。因此，可以选用快中子与地层元素原子核发生俘获反应产生的γ射线作为次生γ源来测量地层密度[13-15]。

γ射线从产生之后就会发生衰减，最后在探测器处测量未被吸收的γ射线的响应为[16]：

式中：φ（En,）是中子通量；Ni表示地层中第i种元素的原子密度；σi(E)是第i种元素的非弹性散射截面；Nγ(Eγ)表示中子与第i种元素→发生非弹性散射之后释放出的γ射线数；f（μ（ Eγ）,）是γ扩散函数，此函数与地层的密度相关；μ是γ扩散因子，并且DRF(Eγ)是γ射线探测器响应函数。

由式(1)可知，测量的γ强度不仅由γ射线的衰减（衰减函数与地层密度相关）决定，同时还受到了中子空间分布的影响。可控源密度测量中应用较多的是D-T源，D-D源的性能在一定程度上也满足密度测量的要求，但两者在性能上有一定差别，具体区别如表1所示[17]。

表1　D-D和D-T源性能对比图Table 1　The comparison chart of D-D and D-T source performance.

从表1可以看出，在应用上述两种中子源进行密度测量时，由于其性能有一定的区别，与地层元素的反应也有区别。目前商用的D-D中子源产额大部分为106、D-T源是108。但最近国外文献显示：现在的D-D中子源甚至可以具备108产额了，斯伦贝谢公司最近将测井用的D-T源提高到了109产额[18-19]。在应用可控中子源进行密度测量时，γ射线主要来自于中子与原子核发生非弹性散射和热中子俘获反应。

1.1 D-T源密度测量

氘-氚(D-T)中子发生器发射的能量为14 MeV的高能中子射入地层后，依次与地层中的元素发生非弹性、弹性散射以及俘获作用，其中快中子的非弹性散射和热中子的俘获都会释放出γ射线。这些次生的γ射线可以用来替代传统密度测井中的137Cs放射源来测量地层密度，其中很多研究表明[20-23]，非弹性散射产生的γ射线与地层密度相关性更高一些，因此可以优选快中子与地层元素作用产生的非弹性散射γ射线作为一个“次生的”γ源来测量地层的密度值。

其中，地层快中子非弹性散射γ射线计数，主要来自碳、氧、硅、钙的贡献[24]。表2为以上几种元素的非弹反应阈能。

表2　地层常见元素非弹性散射的阈能Table 2　Inelastic scattering threshold value of common formation elements.

由表2可见，D-T源密度测井中，14 MeV的中子与地层主要元素即12C、16O、28Si、40Ca反生非弹散射产生的γ射线的能量分别为4.78 MeV、6.51 MeV、1.84 MeV、3.82 MeV。

1.2 D-D源密度测量

在应用D-D源进行密度测量时，脉冲源释放出的能量为2.5 MeV的快中子与地层反应，产生γ射线的过程有非弹性散射、热中子俘获。由表2可以看出，当使用D-D中子源进行密度测量，可以发生非弹性散射的主要是Si元素，而C、O、Ca元素的主要非弹反应阈能都要大于2.5 MeV，因此不以非弹性散射为主。

为验证D-D源密度测量过程中非弹γ和俘获γ的贡献，采用蒙特卡罗的方法模拟得到孔隙度分别为0%、10%、20%的饱含淡水砂岩(SiO2)地层中不同源距(6 cm、18 cm、30 cm、45 cm、60 cm)处非弹γ占总γ贡献的比例，如表3所示。

表3　砂岩地层非弹γ比例Table 3　The inelastic gamma ratio of sandstone formation.

从表3可以看出，纯砂岩地层中距D-D源18 cm处，总γ射线以非弹γ为主，其计数占总γ计数的47%；而在60 cm处以俘获γ为主，非弹γ贡献减弱，只占3%；在高孔隙地层中，非弹γ比例分别降为10%和1%。

产生该结果的原因是纯砂岩地层只有Si元素发生非弹性散射和热中子俘获，O元素基本不参与反应；而在高孔隙砂岩地层，H元素的加入影响俘获γ的空间分布。由于其它元素不发生非弹性散射，因此当地层不是砂岩或石英含量很少时，这3%的砂岩误差不会对结果产生较大影响，但如果是低孔隙的纯砂岩地层，3%的误差已经超过了仪器的允许误差，必须进行校正。同样，纯灰岩(CaCO3)地层中非弹和俘获γ随源距的变化规律如表4所示。

表4　灰岩地层非弹γ比例Table 4　The inelastic gamma ratio of limestone formation.

由表4，在18 cm处非弹γ的计数仅占1%，总γ射线以俘获γ为主，而到了离源距离更大的60 cm处，非弹γ的比例仅有0.1%，因此使用D-D源在灰岩地层产生的γ射线主要是俘获γ的贡献。

综上所述，D-D中子源与地层反应进行密度测量，产生的γ射线主要是俘获γ，该γ射线在输运中发生衰减，衰减规律与传统γ-γ密度测井机理相同。因此，可以选用快中子与地层元素原子核发生反应产生的俘获γ射线作为次生γ源来测量地层密度。其中，地层中热中子俘获截面大，对热中子俘获γ计数贡献较大的核素主要有：1H、28Si、35Cl、40Ca、56Fe，它们的特征峰对应的能量分别为2.23 MeV、3.54 MeV、4.93 MeV、6.11 MeV、6.64 MeV、7.42 MeV、6.42 MeV、7.64 MeV。

2　模型的建立

为考察应用D-D源进行密度测量过程中次生俘获γ计数的影响，采用蒙特卡罗模拟程序建立三维地层模型，模拟不同条件下D-D和D-T源发射的快中子经地层减速到达探测器的计数。模型的具体条件为：井眼直径为7.9 cm，地层径向半径为150cm、高为250 cm，把整个地层划分成环距为5cm、厚度10 cm的相邻栅元；测井仪器外径为7.9cm，脉冲中子源处于仪器前端5-15 cm处，脉冲宽度为40μs；探测器采用5 cm×10 cm的NaI晶体。在中子源和探测器之间放置理想屏蔽体。通过记录不同源距处的γ射线通量，对比分析应用D-T源和D-D源进行密度测量的次生γ源的分布情况。

3　D-D与D-T源次生γ对比

本文设置D-T发生器的产额为1×108n·s-1、D-D发生器的产额为2×107n·s-1，则得到的纯灰岩和纯砂岩地层中次生γ射线随不同源距的分布，结果如图1所示。

图1　纯灰岩(a)和纯砂岩(b)地层次生γ源随源距分布Fig.1　Distribution of the induced gamma change with spacing in pure limestone formation(a)and pure sandstone formation(b).

从图1可以看出，对于D-T源，当源距小于50cm时，非弹γ占主要贡献；50 cm左右时非弹γ和俘获γ对次生γ源的贡献相当，之后俘获γ贡献占优。而对于D-D源，在源距大于40 cm之后，总γ射线以俘获γ为主。并且D-D源与D-T源相比，源距大于50 cm后，D-D源产生的俘获γ的强度要大于D-T源的非弹γ。国外某厂家使用D-T源的随钻密度测井仪器的γ源距约为90 cm，在此距离处的D-D俘获γ计数要远高于D-T非弹γ计数。因此，在纯灰岩和纯砂岩地层的D-D次生γ不会由于计数偏低而造成较大的统计误差。

除了纯岩石骨架地层，俘获γ还对地层中的流体比较敏感，因此还分析了不同孔隙度的砂岩地层情况，图2(a)、(b)分别是孔隙度20%和40%、饱含淡水的砂岩地层中次生γ射线通量随源距的分布。

对比图1(b)、图2(a)、(b)，虽然砂岩的孔隙度增加了，但对于D-T源，其非弹γ的贡献基本不变，而对于D-D源，俘获γ的计数变化率增大了。当源距大于50 cm后，D-D俘获γ的变化率和D-T非弹γ的变化率基本一致，但是D-D俘获γ的计数略低。分析原因主要是由于地层中氢元素的影响，不仅影响俘获γ的空间分布，对其强度大小也有一定影响。

但是由于距D-D源80 cm处的俘获γ计数与距D-T源90 cm处的非弹γ计数基本相当，所以即使在孔隙度或岩性不同的地层中，也可以通过调整源距的大小，使D-D源次生γ计数达到D-T源密度测井的计数量级。因此，如果使用D-D源进行密度测量，并不会存在次生γ计数偏低而引起统计精度差的问题。

4　结语

本文通过分析D-D源与地层相互作用过程，明确了D-D源密度测量过程中次生γ类型，以及次生γ强度高低对密度测量可行性的影响，得到了以下结论，并提出下一步研究计划。

1)D-D源密度测量时，应用的次生γ主要为俘获γ，特别是当源距大于40 cm时，非弹γ的贡献很小，可忽略。

2)虽然D-D源的中子产额较D-T源低一个数量级，但当地层含氢量较低且源距大于50 cm时，D-D俘获γ计数高于D-T非弹γ计数；即使地层含氢量增加时，D-D俘获γ计数略低，但也满足密度测量的统计精度要求，因此D-D源的俘获γ计数的高低不会对密度测量产生较大影响。

3)D-D源密度测量中俘获γ作为次生γ源，影响因素多且复杂，如地层含氢量、元素种类和含量等，要想提高密度测量的精度，需要进一步分析其变化规律，针对其影响因素提出合理的校正方案。

参考文献

1彭琥.2000-2008年放射性测井技术进展述评[J].测井技术,2009,33(1):1-8 PENG Hu.Review on progress of radioactive well logging technology in 2000-2008[J].Well Logging Technology,2009,33(1):1-8

2Quirein J,Smith J H,Chen D D,et al.Formation density prediction using pulsed neutron capture tools[A].SPWLA 46thannua1 logging symposium[C].New Orleans,June 26-29,2005,Paper QQ

3Rodriguez R,Weller G,Evans M.A new approach for identifying gas response using multiple LWD density measurements[A].SPWLA 50thannual logging symposium[C].Woodlands,Texas,June 21-24,2009,Paper YY

4张锋,袁超,刘军涛.基于可控中子源的多探测器随钻中子测井方法数值模拟研究[A].第十七届测井年会论文集[C].2011:75-84 ZHANG Feng,YUAN Chao,LIU Juntao.Numerical simulation of neutron logging-while-drilling method based on multi-detector controllable neutron source[A].17thlogging symposium[C].2011:75-84

5于华伟,Gardner R,李勇华.随钻脉冲中子伽马密度测井的数值模拟研究[A].第一届随钻测控技术发展研讨会论文集[C].2012 YU Huawei,Gardner R,LI Yonghua.Numerical simulation of pulsed-neutron density logging while drilling[A].The first drilling control technology development symposium[C].2012

6Nicole R,Evans M,Allioli F,et al.Neutron-gamma density(NGD):principles,field test results and log quality control of a radioisotope-free bulk density measurement[A].SPWLA 53thannual logging symposium[C].Cartagena,Colombia,June 16-20,2012

7于华伟,张丽,侯博然.电子对效应对D-T次生γ密度测量的影响[J].核技术,2015,38(8):080502.DOI:10.11889/j.0253-3219.2015.hjs.38.080502 YU Huawei,ZHANG Li,HOU Boran.Effect of pair production on D-T induced gamma density measurement[J].Nuclear Techniques,2015,38(8):080502.DOI:10.11889/j.0253-3219.2015.hjs.38.080502

8潘保芝,张瑞,刘坤,等.脉冲中子源密度测井俘获伽马射线强度的数值模拟[J].石油物探,2014,53(6):642-648 PAN Baozhi,ZHANG Rui,LIU Kun,et al.Numerical simulation of pulsed neutron source density logging’s secondary captured gamma ray strength[J].Geophysical Prospecting for Petroleum,2014,53(6):642-648

9严岩,李炳营,郑世平,等.D-D中子源孔隙度测井的蒙特卡罗模拟[J].兰州大学学报(自然科学版),2012,3(48):123-127 YAN Yan,LI Bingying,ZHENG Shiping,et al.Monte Carlo simulation ofthe neutron porosity logging using D-D neutron generator[J].Journal of Lanzhou University(Natural Sciences),2012,3(48):123-127

10 Badruzzaman A,Bames S,Bair F,et al.Radioactive sources in petroleum industry:applications,concerns and alternatives[C].SPE123593,2009.DOI:10.2118/123593-ms

11 Bond L J,Denslow K M,Griffin J W,et al.Evaluation of non-nuclear techniques for well logging technology evaluation[R].Washington:PNNL,2010.DOI:10.2172/1006309

12 Griffin J W,Moran T L,Bond L J.Radiation source replacement workshop[R].Washington:The Pacific Northwest National Laboratory,2010.DOI:10.2172/1062523

13 于华伟.随钻环境下脉冲中子测量地层密度的理论基础研究[D].山东:中国石油大学,2011 YU Huawei.The fundamental research of the pulsed-neutron density logging while drilling[D].Shandong:China University of Petroleum(East China),2011

14 Yu H W,Sun J M,Gardner R P.Accuracy and borehole influences in pulsed neutron gamma density logging while drilling[J].Applied Radiation and Isotopes,2011,69(9):1313-1317.DOI:10.1016/j.apradiso.2011.04.023

15 张丽.随钻方位密度成像测井基础研究[D].山东:中国石油大学,2013 ZHANG Li.The fundamental research of azimuthal density imaging logging while drilling[D].Shandong:China University of Petroleum(East China),2013

16 Han X G,Pemper R,Tutt T,et al.Environmental corrections and system calibration for a new pulsed-neutron mineralogy instrument[A].SPWLA 50thannual logging symposium[C].June 21-24,2009

17 Committee on Radiation Source and Replacement.Radiation source use and replacement abbreviated version[M].Washington:National Academies Press,2008:147-158.DOI:10.17226/11976

18 Chen A X,Antolak A J,Leung K N.Electronic neutron sources for compensated porosity well logging[J].Nuclear Instruments and Methods in Physics Research A,2012,684:52-56.DOI:10.1016/j.nima.2012.04.053

19 Stoller C.Thermal neutron porosity measurements using a pulsed D-T generator in LWD and Wireline[C].USA:Pacific Northwest National Laboratory,2012

20 Wilson R D.Bulk density logging with high-energy gammas produced by fast neutron reactions with formation oxygen atoms[A].Nuclear science symposium and medical imaging conference[C].USA,1995:209-213.DOI:10.1109/nssmic.1995.504211

21 Badruzzaman A,Calif A.Multidetector pulsed-neutron through-tubing cased-hole density measurement sonde[P].US:5825024,October 20,1998

22 Badruzzaman A,Badruzzaman T,Zalan T A,et al.Multi-sensor through-casing density and saturation measurement concepts with a pulsed neutron source:a modeling assessment[A].SPE international petroleum conference in Mexico[C].SPE89884,Puebla,Mexico,November 8-9,2004.DOI:10.2523/89884-ms

23 Odom R O,Bailey S M,Wilson R D,et al.Pulsed neutron density measurements:modeling the depth of investigation and cased-hole wellbore uncertainties[A].SPWLA 40thannual logging symposium[C].May 30-June 3,1999,Paper JJ

24 黄隆基.核测井原理[M].东营:石油大学出版社,2000:89-90 HUANG Longji.Principle of nuclear logging[M].Dongying:University of Petroleum Press,2000:89-90

收稿日期：2016-01-11，修回日期：2016-02-15

Corresponding author:YU Huawei,E-mail:yhwlog@163.com

通信作者：于华伟，E-mail:yhwlog@163.com

DOI:10.11889/j.0253-3219.2016.hjs.39.030502

中图分类号TL99，TE19