CN111638559B - 一种基于脉冲中子测井的快中子散射截面表征方法 - Google Patents
一种基于脉冲中子测井的快中子散射截面表征方法 Download PDFInfo
- Publication number
- CN111638559B CN111638559B CN202010571427.2A CN202010571427A CN111638559B CN 111638559 B CN111638559 B CN 111638559B CN 202010571427 A CN202010571427 A CN 202010571427A CN 111638559 B CN111638559 B CN 111638559B
- Authority
- CN
- China
- Prior art keywords
- section
- scattering cross
- fast neutron
- neutron scattering
- gamma
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001956 neutron scattering Methods 0.000 title claims abstract description 75
- 238000012512 characterization method Methods 0.000 title claims abstract description 29
- 230000004907 flux Effects 0.000 claims abstract description 40
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims abstract description 4
- 238000012417 linear regression Methods 0.000 claims abstract description 4
- 238000004088 simulation Methods 0.000 claims description 6
- 230000005251 gamma ray Effects 0.000 claims description 4
- 101150077194 CAP1 gene Proteins 0.000 claims description 3
- 101150014715 CAP2 gene Proteins 0.000 claims description 3
- 101100438378 Neurospora crassa (strain ATCC 24698 / 74-OR23-1A / CBS 708.71 / DSM 1257 / FGSC 987) fac-1 gene Proteins 0.000 claims description 3
- 101100326803 Neurospora crassa (strain ATCC 24698 / 74-OR23-1A / CBS 708.71 / DSM 1257 / FGSC 987) fac-2 gene Proteins 0.000 claims description 3
- 239000004215 Carbon black (E152) Substances 0.000 claims 1
- 229930195733 hydrocarbon Natural products 0.000 claims 1
- 150000002430 hydrocarbons Chemical class 0.000 claims 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 4
- 238000011161 development Methods 0.000 abstract description 2
- 239000003345 natural gas Substances 0.000 abstract description 2
- 239000003209 petroleum derivative Substances 0.000 abstract description 2
- 238000005755 formation reaction Methods 0.000 description 23
- 239000007789 gas Substances 0.000 description 9
- 238000005259 measurement Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 235000019738 Limestone Nutrition 0.000 description 2
- 229910003271 Ni-Fe Inorganic materials 0.000 description 2
- 230000004323 axial length Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001739 density measurement Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000006028 limestone Substances 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
- G01V5/04—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
- G01V5/08—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays
- G01V5/10—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using neutron sources
- G01V5/104—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using neutron sources and detecting secondary Y-rays as well as reflected or back-scattered neutrons
- G01V5/105—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using neutron sources and detecting secondary Y-rays as well as reflected or back-scattered neutrons the neutron source being of the pulsed type
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种基于脉冲中子测井的快中子散射截面表征方法,涉及石油天然气开发领域。包括以下步骤:通过对非弹伽马场分布的公式变形,获得快中子散射截面的表征形式,得出影响快中子散射截面表征的影响因素包括地层密度和探测器源距;根据选定的脉冲中子测井仪器参数,通过建立地层密度的表征方法,得到利用非弹伽马通量以及非弹、俘获伽马通量比信息表征快中子散射截面的数学表达式;建立三维MCNP数值计算模型,模拟得到多探测器伽马通量信息数据库;利用最小二乘方法对不同地层条件的快中子散射截面与伽马通量和通量比信息做三元线性回归分析,得到快中子散射截面的刻度公式的数学形式,最终形成快中子散射截面的表征方法。
Description
技术领域
本发明涉及石油天然气开发领域,具体涉及一种基于脉冲中子测井的快中子散射截面表征方法。
背景技术
快中子散射截面为一种地层核特性,是衡量地层与快中子相互作用能力的物理参数,主要用于地层含油气饱和度评价。中子测井中,常规饱和度评价方法易受到中子与某些元素的相互作用影响,例如B、Cl和Gd用于Sigma测量,而H用于孔隙率测量。而快中子散射截面仅对气液两相流体变化敏感,对饱和度评价的适用范围更广。同时,当无法进行伽马-伽马密度测量时(例如禁止使用放射性同位素源的情况),快中子散射截面还可用于裸眼井密度测量。
目前,测量快中子散射截面的计数主要是通过中子测井仪器长源距探测器的非弹伽马通量来进行测量。但是由于测井仪器自身的限制以及伽马射线的计数统计性,对快中子散射截面的测量存在一定的误差,难以满足现场油气评价的需求。
发明内容
本发明的目的是针对上述不足,提出了一种利用脉冲中子测井技术中多探测器的伽马信息组合的表征快中子散射截面的方法。
本发明具体采用如下技术方案:
一种基于脉冲中子测井的快中子散射截面表征方法,适用于D-T源脉冲中子测井仪器,包括以下步骤:
步骤1:通过对非弹伽马场分布的公式变形,获得快中子散射截面的表征形式,得出影响快中子散射截面表征的影响因素包括地层密度和探测器源距;
步骤2:根据选定的脉冲中子测井仪器参数,通过建立地层密度的表征方法,代入快中子散射截面的表征公式,得到利用非弹伽马通量以及非弹、俘获伽马通量比信息表征快中子散射截面的数学表达式;
步骤3:建立三维MCNP数值计算模型,模拟得到多探测器伽马通量信息数据库;
步骤4:利用最小二乘方法对不同地层条件的快中子散射截面与伽马通量和通量比信息做三元线性回归分析,得到快中子散射截面的刻度公式的数学形式,最终形成快中子散射截面的表征方法。
优选地,步骤1中对非弹伽马场分布的公式变形得到的快中子散射截面表达式如式(1)所示:
FNXS为快中子散射截面;R为源距;φin(R)为非弹伽马射线通量;ρ为地层密度值;α在给定的脉冲中子仪器参数下为常数;μm为地层衰减系数,φ0为D-T源强度,i为一个中子与地层一个原子核碰撞产生的平均光子数,∑in为快中子的非弹性散射截面,且远小于快中子散射截面,四者均视为常数。
优选地,步骤2中选取脉冲中子仪器包含一个D-T中子源,两个伽马探测器脉冲中子仪器。
优选地,步骤2中对密度的表征表达式如下:
ρ为地层密度值;φin1、φin2为近、远源距的非弹伽马通量;φcap1、φcap2为近、远源距的俘获伽马通量,A、B和C为常数系数。
优选地,步骤3中利用MCNP模拟得到的是不同岩性、孔隙度、含油气饱和度条件下的多探测器非弹、俘获伽马计数通量数据库。
优选地,步骤4中求取快中子散射截面表征公式中的系数采用线性多元回归。
本发明具有如下有益效果:
该方法从中子伽马耦合理论入手,通过对快中子散射截面的数学表征,利用双探测器伽马信息组合,消除地层密度对的快中子散射截面的影响,得到了快中子散射截面的数学表达形式。结合MCNP计算多种岩性、含油含气性的多探测器伽马通量,通过数据分析拟合得到快中子散射截面与伽马信息的刻度公式的常数系数,最终形成快中子散射截面的表征方法。该方法计算精度高、普适性广,受地层岩性、含油气性影响小,弥补了利用单探测器进行快中子散射截面表征的不足。
附图说明
图1为仪器地层模型结构示意图;
图2为砂岩地层快中子散射截面计算值与真实值对比;
图3为灰岩地层快中子散射截面计算值与真实值对比;
图4为白云岩地层快中子散射截面计算值与真实值对比;
图5为快中子散射截面井中测量效果图,各个井段快中子散射截面的真实值与计算值的对比。
具体实施方式
下面结合附图和具体实施例对本发明的具体实施方式做进一步说明:
基于脉冲中子测井的快中子散射截面表征方法,适用于D-T源脉冲中子测井仪器,采用一个D-T中子源和至少两个伽马探测器(轴向阵列探测器同样适用),包括以下步骤:
步骤1:通过对非弹伽马场分布的公式变形,获得快中子散射截面的表征形式,得出影响快中子散射截面表征的影响因素包括地层密度和探测器源距。
非弹伽马场分布公式可以表示为:
FNXS为快中子散射截面;R为源距;φin(R)为非弹伽马射线通量;ρ为地层密度值;α在给定的脉冲中子仪器参数下为常数;μm为地层衰减系数,φ0为D-T源强度,i为一个中子与地层一个原子核碰撞产生的平均光子数,∑in为快中子的非弹性散射截面,且远小于快中子散射截面,四者均视为常数。
对非弹伽马场分布的公式(3)变形得到的快中子散射截面表达式如式(1)所示:
FNXS为快中子散射截面;R为源距;φin(R)为非弹伽马射线通量;ρ为地层密度值;α在给定的脉冲中子仪器参数下为常数;μm为地层衰减系数,φ0为D-T源强度,i为一个中子与地层一个原子核碰撞产生的平均光子数,∑in为快中子的非弹性散射截面,且远小于快中子散射截面,四者均视为常数。
步骤2:根据选定的脉冲中子测井仪器参数,通过建立地层密度的表征方法,代入快中子散射截面的表征公式,得到利用非弹伽马通量以及非弹、俘获伽马通量比信息表征快中子散射截面的数学表达式。
步骤2中选取脉冲中子仪器包含一个D-T中子源,两个伽马探测器脉冲中子仪器,但该方法适用于所有的电缆、随钻及多探测器的脉冲中子仪器。
对密度的表征表达式如下:
ρ为地层密度值;φin1、φin2为近、远源距的非弹伽马通量;φcap1、φcap2为近、远源距的俘获伽马通量,A、B和C为常数系数。
结合式(1)和(2),快中子散射截面可以表征为:
K、L、M和N均为常数系数,且
步骤3:建立三维MCNP数值计算模型,模拟得到多探测器伽马通量信息数据库,利用MCNP模拟得到的是不同岩性、孔隙度、含油气饱和度条件下的多探测器非弹、俘获伽马计数通量数据库。以双源距脉冲中子测井仪器为例,建立MCNP数值计算模型,设置地层为不同岩性、孔隙度和含油气饱和度,将地层划分为径向1.0cm×轴向1.0cm的环形栅元,模拟得到近、远探测器的非弹、俘获伽马通量以及其通量比。
如图1所示,双源距脉冲中子测井仪器的地层模型包含D-T中子源1,钨镍铁屏蔽体2,近探测器3,钨镍铁屏蔽体4,远探测器5,地层6,井眼7及仪器外壳8。模拟过程中,地层划分为径向1.0cm×轴向1.0cm的环形栅元以保证模拟的统计准确性,设置地层为0%-30%间隔5%孔隙度的砂岩、灰岩和白云岩地层,并设置地层饱含油和饱含气的条件,模拟得到近、远探测器的非弹、俘获伽马通量以及其通量比。
步骤4:利用最小二乘方法对不同地层条件的快中子散射截面与伽马通量和通量比信息做三元线性回归分析,得到快中子散射截面的刻度公式的数学形式,最终形成快中子散射截面的表征方法。求取快中子散射截面表征公式中的系数采用线性多元回归,求取公式(4)中的系数K、L、M和N,最终得到利用脉冲中子测井技术多探测器伽马信息表征快中子散射截面的计算式。并对比三种岩性条件下利用探测器公式计算的快中子散射截面值和地层实际的快中子散射截面值,两者的绝对误差在0.2m-1以内,如图2、图3和图4所示。
如图5所示,为本发明的快中子散射截面井中测量效果图,第一道为深度;第二道为地层岩性剖面,第三道为地层孔隙度,第四道为地层含气饱和度,第五道为计算得到的快中子散射截面(虚线)和理论的快中子散射截面(灰度充填)。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (6)
1.一种基于脉冲中子测井的快中子散射截面表征方法,适用于D-T源脉冲中子测井仪器,其特征在于,包括以下步骤:
步骤1:通过对非弹伽马场分布的公式变形,获得快中子散射截面的表征形式,得出影响快中子散射截面表征的影响因素包括地层密度和探测器源距;
步骤2:根据选定的脉冲中子测井仪器参数,通过建立地层密度的表征方法,代入快中子散射截面的表征公式,得到利用非弹伽马通量以及非弹、俘获伽马通量比信息表征快中子散射截面的数学表达式;
步骤3:建立三维MCNP数值计算模型,模拟得到多探测器伽马通量信息数据库;
步骤4:利用最小二乘方法对不同地层条件的快中子散射截面与伽马通量和通量比信息做三元线性回归分析,得到快中子散射截面的刻度公式的数学形式,最终形成快中子散射截面的表征方法。
3.如权利要求1所述的一种基于脉冲中子测井的快中子散射截面表征方法,其特征在于,步骤2中选取脉冲中子仪器为包含一个D-T中子源,两个伽马探测器脉冲中子仪器。
5.如权利要求1所述的一种基于脉冲中子测井的快中子散射截面表征方法,其特征在于,利用MCNP模拟得到的是不同岩性、孔隙度、含油气饱和度条件下的多探测器非弹、俘获伽马计数通量数据库。
6.如权利要求1所述的一种基于脉冲中子测井的快中子散射截面表征方法,其特征在于,步骤4中求取快中子散射截面表征公式中的系数采用线性多元回归。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010571427.2A CN111638559B (zh) | 2020-06-22 | 2020-06-22 | 一种基于脉冲中子测井的快中子散射截面表征方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010571427.2A CN111638559B (zh) | 2020-06-22 | 2020-06-22 | 一种基于脉冲中子测井的快中子散射截面表征方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111638559A CN111638559A (zh) | 2020-09-08 |
CN111638559B true CN111638559B (zh) | 2022-12-13 |
Family
ID=72329310
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010571427.2A Active CN111638559B (zh) | 2020-06-22 | 2020-06-22 | 一种基于脉冲中子测井的快中子散射截面表征方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111638559B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5900627A (en) * | 1997-06-19 | 1999-05-04 | Computalog Research, Inc. | Formation density measurement utilizing pulse neutrons |
CN106597560A (zh) * | 2016-11-30 | 2017-04-26 | 中国石油大学(华东) | 一种利用快中子场分布表征的中子伽马密度测井方法 |
CN107288629A (zh) * | 2017-07-25 | 2017-10-24 | 中国石油大学(华东) | 一种基于新型n‑γ双粒子探测器的中子伽马密度测井方法 |
CN110486002A (zh) * | 2019-08-26 | 2019-11-22 | 中国石油大学(北京) | 中子伽马密度测井中地层体积密度确定方法及设备 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2649474A4 (en) * | 2010-11-11 | 2015-01-21 | Services Petroliers Schlumberger | NEUTRON GAME DRAWING THROUGH A NORMALIZED UNELASTIC RATIO |
-
2020
- 2020-06-22 CN CN202010571427.2A patent/CN111638559B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5900627A (en) * | 1997-06-19 | 1999-05-04 | Computalog Research, Inc. | Formation density measurement utilizing pulse neutrons |
CN106597560A (zh) * | 2016-11-30 | 2017-04-26 | 中国石油大学(华东) | 一种利用快中子场分布表征的中子伽马密度测井方法 |
CN107288629A (zh) * | 2017-07-25 | 2017-10-24 | 中国石油大学(华东) | 一种基于新型n‑γ双粒子探测器的中子伽马密度测井方法 |
CN110486002A (zh) * | 2019-08-26 | 2019-11-22 | 中国石油大学(北京) | 中子伽马密度测井中地层体积密度确定方法及设备 |
Non-Patent Citations (2)
Title |
---|
Measuring the Density of Different Materials by Using the Collimated Fast Neutron Beam;D. Sudac 等;《2015 4th International Conference on Advancements in Nuclear Instrumentation Measurement Methods and their Applications (ANIMMA)》;20151231;第1-4页 * |
随钻脉冲中子_伽马密度测井响应数值模拟;张锋 等;《地球科学(中国地质大学学报)》;20130930;第38卷(第5期);第1116-1120页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111638559A (zh) | 2020-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dong et al. | Research on recognition of gas saturation in sandstone reservoir based on capture mode | |
CN108643890B (zh) | 一种确定地层孔隙度的方法 | |
CN112904436B (zh) | 一种元素产额与热中子计数比组合的孔隙度测量方法 | |
JPH0156397B2 (zh) | ||
US20090108210A1 (en) | Neutron logging tool having source and target with deuterium-tritium gas admixture | |
Evans et al. | Sourceless neutron-gamma density SNGD: a radioisotope-free bulk density measurement: physics principles, environmental effects, and applications | |
CN111638559B (zh) | 一种基于脉冲中子测井的快中子散射截面表征方法 | |
US5808298A (en) | Method for determining formation hydrocarbon saturation and wellbore hydrocarbon holdup from multidetector carbon-oxygen measurements | |
CN110469324B (zh) | 一种基于脉冲中子测井的计算地层密度方法 | |
CN115291288B (zh) | 基于铁中子标记的随钻脉冲中子孔隙度智能处理方法 | |
Fu et al. | A new method of gas reservoir evaluation based on neutron cross section logging | |
CN112523741B (zh) | 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法 | |
CN112523742B (zh) | 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 | |
CN115012920A (zh) | 基于双clyc双粒子探测器的可控中子源多谱测井仪及测井方法 | |
CN112630848B (zh) | 一种基于能谱测井特征谱段的铀矿定量剥离系数求法 | |
CN115267930A (zh) | 一种基于d-t脉冲中子源的高灵敏度中子孔隙度测量方法 | |
CN111894571A (zh) | 一种基于岩性扫描测井资料的流体性质识别方法 | |
Yuan et al. | Can the evaluation accuracy of elemental concentration be further enhanced in geochemical logging?—A break attempt to obtain purer inelastic gamma spectrum | |
Moake | What is the depth of investigation of openhole nuclear tools? | |
Guo et al. | Impact of density tool depth of investigation on dip determination in high angle and horizontal wells | |
Fan et al. | A New Method of Porosity Determination by DT Neutron Generator and Dual CLYC Detector | |
CN112016215B (zh) | 一种基于牛顿迭代法的密度测井六参数反演计算方法 | |
CA1250060A (en) | Method and apparatus for differentiating low porosity limestones from high porosity gas sands | |
Zhang et al. | A new calculating formation porosity method using PNC logging technology in cased wells | |
CN115614022B (zh) | 铀裂变瞬发中子测井的井眼水层原位修正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |