DT脉冲中子发生器随钻中子孔隙度测井的蒙特卡罗模拟.doc
《DT脉冲中子发生器随钻中子孔隙度测井的蒙特卡罗模拟.doc》由会员分享,可在线阅读,更多相关《DT脉冲中子发生器随钻中子孔隙度测井的蒙特卡罗模拟.doc(7页珍藏版)》请在三一办公上搜索。
1、D 2T 脉冲中子发生器随钻中子孔隙度测井的蒙特卡罗模拟锋1 ,靳秀云2 ,侯爽1张( 1 . 中国石油大学 ( 华东) 地球资源与信息学院 ,山东 青岛 266555 ;2 . 中石化胜利油田 东辛采油厂 ,山东 东营 257094)摘要 : 利用蒙特卡罗方法模拟研究了 D2T 脉冲中子发生器和241 A m2Be 中子源产生的中子与地层的作用过程 ,以探讨 D2T 脉冲中子发生器在随钻中子孔隙度测井中的应用价值 。模拟结果显示 ,使用这两种中子源 , 热中子计数均随源距增加而呈指数下降 ;孔隙度较小时 ,两者的计数差异较小 ,当地层孔隙度达到 40 %时 , D2T 脉冲中子发生器产生的热
2、中子和超热中子计数均比241 A m2Be 中子源高很多 ,其分布范围也更宽 ,近探测器的源距选择 2030 cm ,远探测器的源距选择约 6070 cm ;D2T 脉冲中子发生器用于中子孔隙度测井时对地层孔隙度的灵敏度降低 ,而相同源距条件下探测深度几乎不变 。以上结果提示 ,利用 D2T 脉冲中子 发生器可以进行补偿中子孔隙度测井 ,在增加源距的同时既可以保证计数统计性 ,又可以提高灵敏度和探 测深度 ,在随钻测井仪器设计中可以取代241 A m2Be 中子源 。关键词 : D2T 脉冲中子发生器 ;241 A m2Be 中子源 ;随钻测井 ;中子孔隙度测井 ;蒙特卡罗模拟中图分类号 :
3、P631 . 8 ; TL 816文献标志码 : A文章编号 : 100027512 (2010) 0120015207Monte Carlo Simulat ion on Compensated Neutron Porosity Loggingin L WD With D2T Pul sed Neutron GeneratorZ H A N G Fe ng1 , J IN Xi u2yun2 , HO U Sh ua ng1( 1 . Fac ul t y o f Geo2resou rces a n d I n f o rm at i on , Chi n a U ni ve rs i t
4、 y o f Pet role um , Qi n g d ao 266555 , Chi n a;2 . D on g x i n Oi l R ecove r y Pl a nt , S I N O P EC S hen g l i Oi l Fiel d Com p a n y , D on g y i n g 257094 , Chi n a)Abstract : The p roce ss of ne ut ro n i nt e ractio n i nduce d by D2T p ul sed ne ut ro n ge nerato r a nd241 A m2Be so u
5、rce wa s si mulat ed by u si ng Mo nt e Ca rlo met ho d . It i s co ncl ude d t hat t he t her2 mal ne ut ro n co u nt de sce nd e xpo ne ntiall y a s t he sp aci ng i ncrea si ng. The smaller po ro sit y wa s , t he smaller t he diff ere nce s bet wee n t he t wo so urce s we re . Whe n t he po ro
6、sit y reached40 % , t he ratio of t her mal ne ut ro n co unt ge ne rat e d by D2T p ul sed ne ut ro n so urce wa s muchla r ger t ha n t hat ge nerat e d by 241 A m2Be ne ut ro n so urce , a nd it s di st ri b utio n ra nge wa s wi2der . The nea r sp aci ng select e d wa s 20230 c m , a nd t hat of
7、 f a r sp aci ng wa s a bo ut 60270 c m.收稿日期 : 2009208205 ; 修回日期 : 2009210229基金项目 : 国家自然科学基金资助项目 ( 40874065)作者简介 : 张锋 ( 1970 ) ,男 ,山东莒南人 , 副教授 ( 博士) ,从事核测井方法基础研究 、核测井资料处理及蒙特卡罗模拟等工作The det ectio n dep t h by u si ng D2T p ul se d ne ut ro n so urce wa s al mo st uncha nged under co ndi2tio n of t he
8、sa me sap ci ng , a nd t he se n sitivit y of mea sure me nt to t he fo r matio n po ro sit y de2crea se s. The re sult s sho we d t hat it ca n no t o nl y gua ra nt ee t he st ati stic of co unt , but al so i mp ro ve det ectio n se n sitivit y a nd dep t h at t he sa me ti me of i ncrea si ng sp
9、aci ng. Therefo re ,241 A m2Be ne ut ro n so urce ca n be rep lace d by D2T ne ut ro n t ube i n L WD tool .Key words : D2T p ul sed ne ut ro n ge nerato r ; 241 A m2Be ne ut ro n so urce ; L WD ; co mp e n sat ed ne ut ro n po ro sit y ; Mo nt e Ca rlo si mulatio n过去的 20 年中 ,随钻核测井中普遍采用放射性化学源 。从第一代随
10、钻核测井仪中的补偿 密度中子孔隙度测井 ( CDN ) 到第二代的方位密度中子孔 隙 度 测 井 ( A ND ) 都 采 用 同 位 素 中 子源 122 。目前 ,补偿中子孔隙度测井仍采用电缆 测井时所用的241 A m2Be 中子源和两个3 He 管中子探测器 ,它们是通过记录两个探测器的热中子 计数比获取中子孔隙度 3 。同位素中子源会对 操作人员造成人体辐射 ,在随钻测井过程中利用 放射源比在电缆测井中具有更大的风险 ,利用脉冲中子发生器代替化学源 , 既可以减小辐射风 险 ,又可以增加岩石物理测量的可靠性 。脉冲中 子发生器的无源中子测量在裸眼井测量中已应用了近 50 年 ,主要用
11、于碳氧比能谱测井和中子寿命测井 。已有的应用和研究结果均表明 ,无论 是电缆测井还是随钻测井 ,采用脉冲中子发生器是未来测井发展的方向 4 。考虑到辐射安全以及其他参数参量 ,斯伦贝谢公司 2005 年开始研 究在随钻过程中利用 D2T 脉冲中子发生器进行补偿热中子孔隙度测井 527 。国内还未见到相关的研究 成 果 。本 工 作 拟 利 用 蒙 特 卡 罗 ( Mo nt eCa rlo ,M C) 方法模拟中子与地层物质的作用过 程 ,对比研究随钻过程中241 A m2Be 中子源和 D2T 脉冲中子发生器的补偿中子孔隙度测井的响应特性 ,为实现国产的随钻测井仪器选用中子源 提供理论依据
12、。( s B q) ,是中子 测 井 中 常 用 的 中 子 源 , 实 际 测井时常用241 A m2Be 中子源的中子产额约为 4 107 / s 。由于 D2T 脉冲中子发生器比241 A m2Be 中子 源产生的中子能量高 ,进入地层的减速长度大 , 热中子在地层中的分布范围更宽 ;另外从中子产 额来看 , D2T 脉冲中子发生器 比241 A m2Be 中子源高 12 个数量级 ,在增加源距的条件下既可 以保证探测的热中子计数率统计性 ,又可以提高 随钻补偿中子孔隙度测井的探测深度 ,因此从计 数统计性和探测深度方面看 ,在随钻测井中利用 D2T 脉冲中子发生器替代241 A m2B
13、e 中子源是可行的 。随钻中子孔隙度测井的蒙特卡罗模拟22 . 1蒙特卡罗模拟方法及计算模型2 . 1 . 1蒙特卡罗模拟方法蒙特卡罗方法 ,又称随机抽样技巧或统计实 验方法 。该方法是建立一个概率模型或随机过程 ,通过逐一模拟和记录单个粒子的历程来计算所求参 数 的 统 计 特 征 , 给 出 所 求 解 的 近 似 值 。M CN P 程序是用来模拟中子和光子联合输运过程 ,通过逐一记录单个粒子的历程 ,对中子或光子与原子核发生碰撞时的位置 、能量 、运动方向 、反应类型 、源分布等多方面进行抽样 ,其平均结果反映中子和光子在物质中的输运 10 。利用 M CN P24C 11 程序对 D
14、2T 脉冲中子发 生器产生的快中子与地层元素原子核的作用过 程进行模拟 ,分别记录不同探测器位置处的热中 子 。为 了 研 究 问 题 的 方 便 , 本 研 究 过 程 不 考虑3 He 管探测器的响应特性 ;模拟过程中所选取 的截面数据为 END F/ B2V I. 0 ,模拟时追踪的历 史中子计数为 5 107 , 记录的热中子能量范围 为 00 . 1 eV ,热中子计数为相对一个源中子的 计数 ,计数误差小于 3 % ,计数时间为 150 mi n 。1 D2T 脉冲中子发生器和 Am2Be 中子源的中子特性D2T 脉 冲 中 子 发 生 器 产 生 中 子 的 核 反 应 为 82
15、9 : d + 3 H 4 He + n 。其产生的快中子能量为 14 MeV ,中子产额可以达到 108 / s ,已经在电缆测井中得到广泛应用 。而电缆测井中的补偿 中子孔隙度测井仪采用的241 A m2Be 中子源活度为 0 . 666 TBq ( 18 Ci ) , 平 均 中 子 能 量 为 4 . 2 5 MeV ,中 子 发 射 率 为 ( 2 . 22 2 . 74 ) 10 - 4 /2 . 1 . 2 蒙特卡罗模拟计算模型利用蒙特卡罗方法建立裸眼井条件下的计 算模型 ,井眼直径为 20 cm ,井眼内充满淡水 ;把 地层设 为 内 外 半 径 分 别 为 10 和 70 c
16、m 、高 为147 c m的圆筒状几何模型 ,如图 1 所示 。地层分别填充不同岩性和流体物质 。把整个地层划分 成高 3 c m 、环距 2 cm 的相邻栅元 ,共有 1 470 个 栅元 。测井仪直径为 85 mm ,将 D2T 脉冲中子子源的热中子计数都呈指数衰减 。地层孔隙度对其的影响表现为 : 孔隙度较小时 ,两者的计数 差异较小 ,当地层孔隙度达到 40 %时 ,D2T 脉冲 中子发生器产生的热中子计数比241 A m2Be 中子 源高很多 。这主要是由于 D2T 脉冲中子发生器 产生的中子能量高 ,地层孔隙度越大 ,中子减速 长度和241 A m2Be 中子 源差 异 越大 ,
17、其 热中 子分 布范围越 宽 。在 井 眼 和 地 层 条 件 下 , 利 用 D2T 脉冲中子发生器和241 A m2Be 中子源补偿孔隙度 测井时零源距相差不大 ,大约为 12 cm 。此外 ,从图 2 的模拟结果还可以看出 ,同一 源距处 D2T 脉冲中子发生器记录的热中子计数 和241 A m2Be 中子源不同 。利用图 2 得到相同地 层条件下同一源距处241 A m2Be 中子源和 D2T 脉 冲中子发生器的热中子计数比 N Am / N D ,其随源 距的变化示于图 3 。由 图 3 可 以 看 出 , 源 距 小 于 30 c m 时 ,241 A m2Be 中子源热中子计数高
18、于 D2T 脉冲中子 发生器 ,原因是 D2T 脉冲中子发生器产生的中 子能量高 ,快中子慢化成热中子走过的路程长 ,源距越小的区域慢化的热中子越少 ;当源距大于30 cm 时 ,两种中子源产生的中子慢化后热中子 的数量与地层的含氢指数有关 ;当地层孔隙度较 小时 ,同一位置处 D2T 脉冲中子发生器产生的 热中子数低 于241 A m2Be 中 子 源 ; 当 地 层孔 隙度较大时 ,正好相反 ,且地层孔隙度越大 ,其热中子 计数相差也越大 。在常规地层条件下 ,D2T 补偿中子孔隙度测 井所设计的短源距处热中子计数大约是241 A m2Be 中子源的 90 % ,而长源距处的热中子计数可以
19、达到241 A m2Be 中子源的 1 倍多 ; 再加上中子发 生 器 和24125 cm处 。A m2Be中 子 源 置 于 仪 器 下 端图 1计算模型3结果分析3 . 1D2T 补偿中子孔隙度测井源距的选择3 . 1 . 1热中子计数随源距变化的规律利用图 1 的计算模型 ,井眼为水 、地层是孔 隙度分别为 10 % 、20 % 、30 %和 40 %的饱含水 砂岩 。测井仪的结构和尺寸相同 ,中子源分别为 D2T 脉冲中子发生器和241 A m2Be 中子 源 , 记录 不同源距 L 处的热中子计数 N , 模拟得到热中 子计数随源距的变化规律 ,结果示于图 2 。从图 2 可 以 看
20、出 , D2T 脉 冲 中子 发 生 器 和241 A m2Be 中 子 源 的 热 中 子 计 数 除 与 源 距 有 关 外 ,还与地层的孔隙度有关 。源距增大 ,两种中图 2 热中子计数与源距的关系= 10 % ; = 20 % ; = 30 % ; = 40 %考虑实际补偿中子孔隙度测井 ,在常规地层条件下选取近探测器源距为 2030 cm ,为了选 择最佳长源距 ,固定短源距为 26 c m , 且热中子计数为 N S ,依次增加 3 c m 改变长源距 ,相应的 热中子计数为 N L ,模拟得到孔隙度分别为 10 %和 40 %饱含水砂岩地层的近远探测器热中子计数 比 R = N
21、S / N L 及 相 对 孔 隙 度 的 动 态 变 化( R2 - R1 ) / ( (2 - 1 ) 100) ,结果列于表 1 。 从表 1 数据可以看出 , 短源距一定的条件下 ,近远探测器热中子计数比随着源距的增加而 增大 ;不同孔隙度地层的热中子计数比差值随源距增加而增大 ,动态范围增加 。因此单纯从孔隙度测井响应来看 ,长源距越大 ,对孔隙度的变化图 3 两种中子源的热中子计数比与源距的关系= 10 % ; = 30 % ; = 40 %产额大约是241 A m2Be 中子源的 35 倍 ,因此从计数统计性来看 ,可以增加长源距 ,以便增加随 钻 D2T 脉冲中子发生器补偿中子
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- DT 脉冲 中子 发生器 孔隙 测井 蒙特卡罗 模拟

链接地址:https://www.31ppt.com/p-4077507.html