一种地层元素测井解谱方法的数值模拟研究

合集下载

利用元素伽马能谱测井识别气层的数值模拟

收稿日期:2018-10-09基金项目:国家自然科学基金项目(41374125,41574119);国家科技重大专项(2017ZX05019005-005);中央高校基本科研业务费专项(15CX06008A)作者简介:张锋(1970-),男,教授,博士,研究方向为核测井方法㊁核测井数据处理及蒙特卡罗模拟㊂E-mail:zhfxy_cn@㊂文章编号:1673⁃5005(2019)04⁃0060⁃09 doi:10.3969/j.issn.1673⁃5005.2019.04.007利用元素伽马能谱测井识别气层的数值模拟张　锋1,2,赵　靓1,张泉滢1,2,李向辉3,遆永周3(1.中国石油大学(华东)地球科学与技术学院,山东青岛266580;2.海洋国家实验室海洋矿产资源评价与探测技术功能实验室,山东青岛266071;3.河南省科学院同位素研究所有限责任公司,河南郑州450000)摘要:利用蒙特卡罗方法建立不同井眼和地层条件下的计算模型,模拟研究非弹俘获计数比R 与含气饱和度的测井响应规律及井筒和地层环境对R 的影响;利用研究成果处理新疆某油田实测井资料,验证元素伽马能谱测井技术在实际生产应用进行地层含气性评价的可行性㊂结果表明:当地层孔隙度一定时,非弹和俘获伽马计数随着含气饱和度的增加而增加,而比值R 减小;当孔隙度增大时,含气饱和度小于某一确定值(约35%)时非弹伽马计数减小,但大于该值时俘获伽马计数减少,比值R 增加;岩性㊁地层水矿化度㊁井眼流体和尺寸以及钻井液侵入等因素都会对地层非弹和俘获计数比产生影响;不同井眼流体对非弹俘获计数比R 的影响较小,在同一地层孔隙度条件下,地层水矿化度越高㊁井径越大,非弹俘获计数比R 越大,且石灰岩地层的非弹俘获计数比R 高于砂岩地层;钻井液侵入较浅时,非弹俘获计数比R 随钻井液侵入的加深而增大,当达到探测深度后,侵入加深但非弹俘获计数比R 保持不变;地层水矿化度越高㊁井径越大㊁钻井液侵入越小越有利于含气饱和度的评价,但在实际处理时需要对岩性㊁地层水矿化度㊁井径进行校正㊂关键词:含气饱和度;元素伽马能谱测井;非弹和俘获计数比;蒙特卡罗模拟中图分类号:P 631 文献标志码:A引用格式:张锋,赵靓,张泉滢,等.利用元素伽马能谱测井识别气层的数值模拟[J].中国石油大学学报(自然科学版),2019,43(4):60⁃68.ZHANG Feng,ZHAO Liang,ZHANG Quanying,et al.Numerical simulation on gas reservoir identification using element gamma spectrum well logging [J].Journal of China University of Petroleum(Edition of Natural Science),2019,43(4):60⁃68.Numerical simulation on gas reservoir identification using elementgamma spectrum well loggingZHANG Feng 1,2,ZHAO Liang 1,ZHANG Quanying 1,2,LI Xianghui 3,DI Yongzhou 3(1.School of Geosciences in China University of Petroleum (East China ),Qingdao 266580,China ;boratory for Marine Mineral Resources ,Qingdao National Laboratory for Marine Science and Technology ,Qingdao 266071,China ;3.Isotope Research Institute of Henan Academy of Sciences Company Limited ,Zhengzhou 450000,China )Abstract :Using Monte Carlo method,we build simulation models under different borehole and formation conditions to simu⁃late the response of R (the ratio of inelastic gamma counts and capture gamma counts)and gas saturation.Meanwhile,the in⁃fluence of R under different borehole and formation environment conditions are also simulated.This method is further appliedto a gas well in Xinjiang Oil field to evaluate gas saturation.The results validate the effectiveness of this method.The simula⁃tion results show that both inelastic and capture gamma counts increase in the same porosity,while the value of R decreases.When formation porosity increases,inelastic gamma counts decrease while gas saturation is less than a certain value (about2019年　第43卷中国石油大学学报(自然科学版) Vol.43　No.4　第4期 Journal of China University of Petroleum Aug.201935%).The trendis opposite when gas saturation is greater,capture gamma counts decrease while the value of R increases. Various factors,such as lithology,formation water salinity,borehole fluid types,borehole size and drilling mud invasion can affect the accuracy of R.Therefore,these factors should be corrected.The effect of the borehole fluid types is not significant. When formation porosity is confirmed,the higher of formation water salinity and the bigger of the borehole diameter are,the greater value of R becomes.In addition,R is larger in limestone formation than in sandstone formation.When the depth of mud invasion is less than the detection depth,the value of R increases with the increasing mud invasion depth.High salinity formation water,big borehole size and deep mud invasion are beneficial for gas saturation evaluation.Keywords:gas saturation;element gamma spectrum well logging;counts ratio of inelastic and capture gamma ray;Monte Carlo simulation 随着油气田勘探开发的深入,地层含气饱和度的评价和定量计算变得尤为重要[1⁃2]㊂传统评价地层含气的方法有脉冲中子曲线叠加技术㊁岩性密度测井㊁三孔隙度重叠法㊁热中子衰减时间测井㊁电阻率测井㊁成像测井[3⁃5];近年来随着脉冲中子测井技术的发展,利用仪器探测到的中子及伽马时间谱信息进行地层含气饱和度评价在油气识别方面发挥了重要作用[6⁃7]㊂国内外专家对定量评价含气饱和度进行了相关研究,Trcka等[8]利用不同探测器的非弹伽马计数比值动态监测储集层含气情况㊂黄导武等[9]提出了在PNC测井中利用长短源距计数率比值计算低矿化度水地层含气饱和度的方法;Guo 等[10]通过利用不同探测器的俘获伽马计数组合来评价地层含气响应㊂张锋等[11]利用三探测器中子测井仪近㊁远探测器和中㊁远探测器热中子计数比值的差值,定量评价含气储层㊂Guo等[12]根据非弹窗对岩性密度敏感的特性将俘获窗分为快㊁慢两个部分,利用在致密储层中非弹窗和慢俘获窗的比值不受地层水矿化度影响的特性识别气层㊂2015年, Gerardo Cedillo等[13]提出了利用多探测器脉冲中子测井仪中不同源距探测器记录的非弹伽马计数比值量化计算地层含气饱和度㊂现有使用脉冲中子源的元素伽马能谱仪器主要包括贝克休斯公司的FLeX 仪器[14]和斯伦贝谢公司的Litho Scanner仪器[15],两种仪器通过利用D-T中子源和单探测器信息在实际勘探开发中完成地层元素评价及有机碳探测,但能否利用元素伽马能谱仪器探测得到的伽马信息进行地层含气饱和度评价仍是亟待解决的问题㊂由于气层和水层在组成和快中子核参数存在差异,笔者利用元素能谱测井仪器记录的伽马能谱信息,通过研究不同井眼和地层条件下非弹及俘获伽马计数比值与地层含气饱和度的关系,提出利用非弹伽马计数与俘获伽马计数的比值确定地层饱和度的方法,为元素伽马能谱测井技术进行气层识别和定量评价提供依据㊂1　元素伽马能谱测井识别气层物理基础相比于常规储集层,煤岩㊁页岩等非常规具有储层岩性复杂㊁储层物性差(孔隙度小于10%㊁渗透率小于1.0×10-3μm2㊁埋藏较深等特点,油气勘探开发难度较大[16⁃17]㊂元素伽马能谱测井技术通过利用中子与地层元素原子核发生作用放出伽马射线信息确定元素含量㊁划分岩性,在非常规储层评价中具有广泛的应用前景㊂脉冲中子源按照一定的频率和脉宽向地层发射快中子,快中子进入地层后发生非弹性散射和弹性散射,产生非弹伽马射线并减速为热中子,热中子在地层中扩散被地层原子核俘获,产生俘获伽马射线㊂不同地层流体对快中子的减速能力和对热中子的俘获能力存在差异,当地层含气时,由于天然气的密度和含氢指数均远小于水,地层对快中子的减速能力降低且热中子的扩散长度增加,因此利用元素伽马能谱测井仪器记录的非弹㊁俘获伽马射线信息能够反映地层流体的变化㊂根据中子扩散理论,假设源强S0的中子源向地层发射快中子,快中子与一个原子核发生非弹性散射放出的伽马射线强度为i,则单位时间内地层体积元d V在探测器处产生的非弹伽马射线强度为d C in=14πX2iΣinφf e-μX d V=iS016π2L2s X2r e-r/L s e-μX d V.(1)C in=iS016π2L2s∭1r e-r/L s e-μX sinθd X dφdθ.(2)式中,C in为非弹伽马计数;Σin为非弹散射截面;φf 为快中子的通量分布;μ为伽马射线的线性吸收系数;L s为快中子的减速长度;r㊁X分别为体积元d V到中子源和探测器的距离㊂同理,热中子与原子核发生俘获反应产生的俘获伽马射线为㊃16㊃第43卷　第4期张　锋,等:利用元素伽马能谱测井识别气层的数值模拟C cap=i′S016π2(L2s-L2t)∭1r e-r/L s e-μX sinθd X dφdθ.(3)式中,C cap为俘获伽马计数;i′为一个热中子与原子核发生俘获反应放出俘获伽马射线强度;L t为热中子的扩散长度㊂因此探测器记录的非弹伽马计数与俘获伽马计数比值R为R=C in Ccap =f(i,i′,L s,L t).(4)式中,比值R只与i㊁i′㊁快中子减速长度L s以及热中子扩散长度L t相关㊂其中i和i′与地层中物质的原子核相关,考虑到i和i′在R中是以比值(i/i′)形式存在的,两者随物质原子核的变化在一定程度上相互抵消,其影响相对L s和L t造成的影响可以忽略,因此不考虑地层物质原子核与中子作用产生非弹和俘获伽马射线强度的比值(i/i′)的影响㊂显然非弹俘获计数比R主要受快中子减速长度L s和热中子扩散长度L t的影响,其主要取决于地层孔隙度φ和含气饱和度S g㊂如图1为砂岩地层(孔隙度为10%和30%)快中子减速长度和热中子扩散长度随地层含气饱和度S g的变化关系[18]㊂图1　减速长度㊁扩散长度与地层含气饱和度关系Fig.1　Relationship between L s and L t withdifferent gas saturation从图1中可以看出,当地层孔隙度一定时,快中子减速长度和热中子扩散长度均随地层含气饱和度的增加而增加,但快中子减速长度受含气饱和度的影响变化幅度更大㊂由于地层中产生的非弹㊁俘获伽马射线分布与快中子减速长度和热中子的扩散长度有关,不同含气饱和度S g地层的非弹俘获计数比R不同,故可利用同一探测器记录的R值来确定地层含气饱和度㊂2　蒙特卡罗模拟及结果分析为了研究非弹和俘获伽马计数比与地层含气饱和度响应关系,利用蒙特卡罗方法模拟不同地层条件下的非弹和俘获伽马场分布㊂建立计算模型参数如下:井眼直径为20cm,井眼内充满淡水,地层半径为10~70cm,高为140cm,地层为砂岩地层,地层孔隙度中分别填充不同的流体;仪器外径为4.3cm,仪器外壳厚0.4cm,仪器贴井壁测量㊂仪器采用由1个D-T脉冲中子发生器和1个BGO探测器组成的测量系统,脉冲宽度为40μs,周期为1000μs,BGO探测器源距为65cm,中子源和探测器之间放置厚度为5cm的钨镍铁屏蔽体㊂为了简化模拟的过程,没有考虑探测器的响应特性,计算模型如图2所示㊂图2　蒙特卡罗计算模型Fig.2　Monte Carlo simulation model模拟中子数目1×108,能量14MeV;模拟时采用F4栅元计数方式,记录伽马时间谱时间范围为0~1000μs,道宽20μs,计数误差小于5%㊂数据处理时,非弹伽马和俘获伽马采用能窗范围均为0~ 8.5MeV,记录非弹伽马时间窗为0~40μs,俘获伽马时间窗为200~600μs㊂利用蒙特卡罗模拟得到同一源距处的非弹伽马和俘获伽马计数,研究利用非弹俘获计数比的方法来定量评价气层㊂2.1　非弹㊁俘获伽马计数随含气饱和度的变化规律利用图2所示模型,井眼内充满淡水,地层孔隙含气(主要成分为CH4,密度为0.2g/cm3),地层含气饱和度为0㊁30%㊁50%㊁70%㊁100%时,改变地层孔隙度为0到40%,间隔为5%,模拟得到不同地层孔隙度条件下,探测器的非弹㊁俘获伽马计数(N I和N C)与地层含气饱和度的关系,如图3所示㊂㊃26㊃中国石油大学学报(自然科学版) 2019年8月图3　不同孔隙度地层伽马计数与含气饱和度关系Fig.3　Relationship between gamma counts and gas saturation with different porosity 由图3可知,探测器记录的非弹㊁俘获伽马计数与地层含气饱和度存在函数关系㊂当地层孔隙度一定时,随着含气饱和度的增加,地层密度下降㊁含氢指数减小,非弹伽马计数和俘获伽马计数均增加;当地层孔隙度改变时,含气饱和度小于某一确定值(约35%)时,非弹伽马计数随着地层孔隙度的增大而减小,而当地层含气饱和度大于此确定值时规律相反㊂当地层孔隙度变化时,俘获伽马计数随着地层孔隙度的增加而减小㊂利用探测器记录到的非弹伽马计数和俘获伽马计数,计算井眼流体为水,地层含气饱和度为0㊁30%㊁50%㊁70%㊁100%,地层孔隙度由0变化为40%时,不同地层条件下的非弹俘获计数比R 随地层含气饱和度的关系如图4所示㊂图4　非弹俘获计数比和地层含气饱和度的响应关系Fig.4　Relationship of R and porosity with differentformation gas saturation由图4可知,非弹和俘获伽马计数比R 的变化能够反映地层含气饱和度的大小㊂当地层孔隙度一定时,非弹和俘获伽马计数均随着地层含气饱和度的增大而增大,计算得到的非弹俘获计数比R 随着地层含气饱和度的增加呈线性减小,地层孔隙度越大,不同含气饱和度引起的R 的差异越大㊂当地层含气饱和度一定时,R 随着地层孔隙度的增大而增大,且地层含气饱和度越小,R 的变化越剧烈㊂根据图4中不同孔隙度条件下的非弹俘获比R 与地层含气饱和度S g 的响应关系,建立地层含气饱和度评价模型S g =(2145-176.6φ+5.93φ2-0.066φ3)+(-2635+233.7φ-7.85φ2+0.088φ3)R.(5)由式(5)可得,地层含气饱和度S g 是非弹俘获计数比R 和地层孔隙度φ的函数,利用探测器探测到的R 结合地层孔隙度参数,能够实现地层含气饱和度评价㊂2.2　影响因素分析及校正由于实际测量中R 受井眼和地层环境的影响㊂利用图2所示模型,模拟得到不同影响因素对R 的影响,为进行实际地层含气饱和度解释评价提供了依据㊂2.2.1　地层水矿化度的影响其他条件不变,改变地层孔隙度由0到40%,间隔为5%,地层完全含水,地层水矿化度C w 为0㊁1%㊁2%㊁5%㊁10%,计算得到地层流体不同时非弹俘获计数比R 随孔隙度的变化规律,如图5所示㊂由图5(a)可知,井眼流体相同,地层孔隙度一定时,R 随着地层水矿化度的增大而增大,这主要是因为地层水矿化度越高,地层对热中子的俘获能力越强,非弹俘获计数比R 越大㊂当地层水矿化度一定时,非弹俘获计数比R 随着地层孔隙度的增大而增大,且地层水矿化度越大,非弹俘获计数比R 随孔隙度的变化越剧烈,差异越明显㊂因此高地层水㊃36㊃第43卷　第4期张　锋,等:利用元素伽马能谱测井识别气层的数值模拟矿化度有利于地层含气饱和度评价㊂图5　非弹俘获计数比和地层水矿化度的响应关系及校正图版Fig.5　Relationship of R and porosity under different formation water salinity and its correction chart of water salinity 尽管地层水矿化度越高越有利于地层含气评价,但矿化度变化会影响地层含气饱和度计算的准确性,由于高地层水矿化度条件下计算得到的R 偏大,从而导致利用公式计算得到的地层含气饱和度较真实值偏小㊂因此在实际利用R 确定地层含气饱和度时还需要结合电阻率曲线计算储层的地层水矿化度,并根据如图5(b)所示的地层水矿化度校正图版对其进行校正;同时,根据俘获伽马计数受地层水矿化度影响明显的特点,还可以通过分析俘获伽马能谱和时间谱,优化俘获伽马计数的能量范围和时间窗提高非弹俘获计数比R 的准确性,从而降低地层水矿化度的影响㊂2.2.2　井眼流体的影响采用图2所示模型,其他条件不变,地层分别饱含气和纯水时,改变地层孔隙度由0改变到40%,间隔为5%,模拟井眼流体为淡水㊁油和矿化度为2%盐水时的测井响应规律,计算得到不同井眼和地层条件下非弹俘获计数比R 随孔隙度的变化规律,如图6所示㊂由图6可得,井眼流体的变化对非弹俘获比R 的影响较小,不同井眼流体条件下计算得到的非弹俘获计数比R 差异较小㊂井眼流体一定时,R 随着地层孔隙度的增大而增大,且井眼流体为水和油条件下计算得到的R 基本相同,对地层含气评价几乎没有影响㊂当井眼流体为矿化度2%盐水时,由于氯对热中子的俘获能力强,相同孔隙度条件下,井眼中盐水的矿化度会影响R 的测量结果,探测器记录的俘获伽马计数随着井眼流体矿化度的增大而减小,导致R 增大㊂图6　不同井眼流体条件下的非弹俘获计数比Fig.6　Rationship of R and porosity withdifferent borehole fluid2.2.3　岩性的影响采用图2所示模型,其他条件不变,地层孔隙度由0改变到40%,间隔为5%,分别模拟地层岩性为砂岩和石灰岩时,地层饱含淡水和气条件下的测井响应规律,计算得到两种岩性条件下非弹俘获计数比R 随孔隙度的变化规律,如图7所示㊂由图7可以看出,岩性对R 计算的绝对值存在影响㊂在不同岩性地层R 均随地层孔隙度的增大而增大,但地层孔隙度增大,岩石骨架不同引起R 的差异减小;且在同一地层孔隙度条件下,石灰岩地层对热中子的俘获能力强,使得其俘获伽马计数比砂岩低,R 更大㊂由于不同岩性计算得到的R 随孔隙度的变化幅度相同,但是R 的绝对值不同,因此在利用R 确定地层含气饱和度时,需要对岩性进行校正㊂㊃46㊃中国石油大学学报(自然科学版) 2019年8月图7　不同岩性条件下的非弹俘获计数比Fig.7　Relationship of R with different lithology2.2.4　井眼尺寸的影响其他条件不变,改变井径为14㊁20㊁24.1cm,地层孔隙中分别饱含淡水和气,地层孔隙度由0改变到40%,间隔为5%,模拟得到3种井径条件下的非弹俘获比,如图8所示㊂由图8(a)可知,非弹俘获计数比R 受井径变化的影响㊂当地层孔隙度一定时,R 随着井径的增大而增大;当井径一定时,R 随着孔隙度的增大而增大㊂通过对比发现,在地层流体相同时,井径的变化仅仅影响R 绝对值的大小,不影响R 随孔隙度变化曲线的斜率,即不同井径条件下,R 随孔隙度变化的相对变化量相同㊂由于井径变化会影响R ,因此可以利用如图8(b)所示校正图版对井径进行校正,从而提高含气饱和度评价的准确性㊂图8　不同井眼尺寸的非弹俘获计数比及校正图版Fig.8　Relationship of R and porosity in different borehole size and its correction chart2.2.5　钻井液侵入的影响地层为孔隙度25%的砂岩地层,地层含气饱和度为80%时,利用图2所建模型,从井壁开始沿着地层径向依次用钻井液驱替一定厚度地层中的气体,使得驱替后的地层含水饱和度为100%,经过若干次驱替,井眼周围地层变为饱含水的砂岩地层㊂通过模拟得到每次驱替后探测器的非弹㊁俘获伽马计数,计算得到非弹俘获计数比R 随钻井液侵入深度的变化关系,如图9所示㊂由图9可以看出,仪器探测深度约为35cm㊂当钻井液侵入较浅(小于35cm)时,R 随着侵入深度的增大而增大,这主要是因为钻井液侵入地层时,相比于含气地层条件下,俘获伽马计数衰减剧烈,从而使得此时计算得到的R 比真实值大;当钻井液侵入深度大于探测深度(大于35cm)时,随着侵入的增加R 基本不变,此时无法利用非弹俘获计数比确定地层含气饱和度㊂图9　不同钻井液侵入深度的非弹俘获计数比Fig.9　Relationship of R with different drillingfluid invasion depth3　应用实例利用元素能谱测井仪器在套管井中实际测量得㊃56㊃第43卷　第4期张　锋,等:利用元素伽马能谱测井识别气层的数值模拟到的非弹性散射伽马计数和俘获伽马计数,计算实测资料的非弹俘获计数比R ,结合模拟研究得到的含气饱和度响应模型对实测数据进行含气饱和度评价㊂新疆某井中5300~5360m 井段利用元素能谱测井资料对地层进行含气性评价,实测解释如图10所示㊂图10　XX 井测井解释Fig.10　Logging interpretation of well XX 结合孔隙度曲线和岩性剖面可知,该井段主要为常规孔隙度砂岩层,夹杂部分泥岩层(泥岩层默认不含气,含气饱和度判定为0),井眼流体为淡水,因此不需要进行岩性和井眼流体校正;根据井径和电法测井曲线得到对应的井眼尺寸和地层水矿化度,利用相应的校正图版对非弹俘获计数比进行了井眼环境校正,得到标准地层情况下的非弹俘获计数比R ;并利用图4得到含气饱和度计算模型,得到了该井段的含气饱和度㊂图10中第5道为探测器测量得到的非弹伽马计数(粉色曲线)和俘获伽马计数(黑色曲线),第6道为非弹俘获计数比R ㊂由图10中可知,当深度为5320~5327m 时,声波时差曲线出现周波跳跃,补偿密度曲线指示地层密度值降低,此时利用图中第5道中的非弹伽马㊃66㊃中国石油大学学报(自然科学版) 2019年8月计数和俘获伽马计数计算得到非弹俘获计数比R,此时R显示低值,综合含气饱和度响应模型求取含气饱和度,结果显示该层为高含气地层;当深度为5327.5~5329m时,地层伽马值为高伽马,岩性显示该层为泥岩层,非弹俘获计数比相对于含气储层高很多,主要由于泥岩层含氢指数高,非弹和俘获伽马计数都降低,但俘获伽马计数降低更剧烈㊂可见,通过结合实际测井资料,利用元素伽马能谱测井技术测量得到的非弹俘获比可以快速划分含气地层,当地层含气时,非弹俘获计数比R显示低值,且含气饱和度越高,R越小;当地层为泥岩层时, R显示为明显的高值㊂同时,结合孔隙度数据,根据模拟得到的含气饱和度响应模型能够对地层含气饱和度定量分析,实现地层含气性评价㊂4　结　论(1)利用元素伽马能谱测井数据信息识别气层,根据伽马时间谱得到非弹俘获计数比R确定地层含气饱和度,可以实现地层含气性评价㊂(2)当地层孔隙度一定时,随着含气饱和度的增加,非弹和俘获伽马计数均增加,而R减小;当地层孔隙度增加时,含气饱和度小于35%时,非弹伽马计数减少,大于该值时规律正好相反,而俘获伽马计数都减少,R都增加㊂(3)井眼流体对非弹俘获计数比R的影响较小,非弹俘获计数比R主要受岩性㊁地层水矿化度㊁井眼尺寸㊁钻井液侵入的影响㊂相同孔隙条件下,石灰岩地层的非弹俘获计数比R高于砂岩地层;当地层条件相同时,地层水矿化度越高㊁井径越大㊁钻井液侵入越小越有利于含气饱和度的评价㊂由于地层环境因素会影响非弹俘获计数比R的计算结果,因此在实际含气评价中需要对相关参数进行校正㊂(4)利用提出的含气饱和度评价模型,对实测井数据进行地层含气饱和度定量评价,解释结果与其他评价方法吻合良好,验证了元素伽马能谱测井技术评价地层含气饱和度的可行性㊂参考文献:[1]　沈慧萍,史振勤,乔国安,等.深层气层识别方法研究[J].天然气工业,2005,25(10):47⁃48.SHEN Huiping,SHI Zhenqin,QIAO Guoan,et al.Stud⁃y on the method of recognizing deep gas zone[J].NaturalGas Industry,2005,25(10):47⁃48.[2]　于世元.复杂气层识别技术研究[J].天然气工业,1999,19(1):67⁃70.YU Shiyuan.Research on the recognition technology ofcomplex gas reservoir[J].Natural Gas Industry,1999,19(1):67⁃70.[3]　邹友龙,谢然红,郭江峰,等.致密储层数字岩心重构及核磁共振响应模拟[J].中国石油大学学报(自然科学版),2015,39(6):63⁃71.ZOU Youlong,XIE Ranhong,GUO Jiangfeng,et al.Re⁃construction of digital core of tight reservoir and simula⁃tion of NMR response[J].Journal of China University ofPetroleum(Edition of Natural Science),2015,39(6):63⁃71.[4]　李曦宁,沈金松,李振苓,等.用多尺度形态学方法实现成像测井电导率图像的缝洞参数表征[J].中国石油大学学报(自然科学版),2017,41(1):69⁃77.LI Xining,SHEN Jinsong,LI Zhenling,et al.Character⁃ization of reservoir fracture and vug parameters by con⁃ductivity image of FMI based on multi⁃scale mathematicalmorphology method[J].Journal of China University of Pe⁃troleum(Edition of Natural Science),2017,41(1):69⁃77.[5]　邵才瑞,张鹏飞,张福明,等.用J函数提高致密砂岩气层饱和度测井评价精度[J].中国石油大学学报(自然科学版),2016,40(4):57⁃65.SHAO Cairui,ZHANG Pengfei,ZHANG Fuming,et al.Improving well log evaluation accuracy of tight sandstonegas saturation using J function[J].Journal of China Uni⁃versity of Petroleum(Edition of Natural Science),2016,40(4):57⁃65.[6]　刘军涛,张锋,张泉滢,等.低孔隙度储层碳氧比测井灵敏度提高方法[J].中国石油大学学报(自然科学版),2016,40(6):57⁃62.LIU Juntao,ZHANG Feng,ZHANG Quanying,et al.Amethod of improving sensitivity of carbon/oxygen well log⁃ging for low porosity formation[J].Journal of China Uni⁃versity of Petroleum(Edition of Natural Science),2016,40(6):57⁃62.[7]　BADRUZZAMAN A,LOGAN J P,BEAN C,et al.Isaccurate gas/steam determination behind pipe feasiblewith pulsed neutron measurements?[R].SPE110098,2007.[8]　TRCKA D E,GILCHRIST W A,RILEY S,et al.Fieldtrials of a new method for the measurement of formation⁃gas saturation using pulsed⁃neutron instruments[R].SPE102350,2006.[9]　黄导武,宋春华,沈晓红,等.脉冲中子俘获测井求取低矿化度水地层含气饱和度方法研究及应用[J].中国海上油气,2007,19(3):169⁃172.㊃76㊃第43卷　第4期张　锋,等:利用元素伽马能谱测井识别气层的数值模拟。

_用地层元素测井(ECS)资料评价复杂地层岩性变化

第25卷　第3期核电子学与探测技术V ol .25　No .3 2005年　5月Nuclear Electro nics&Detectio n Techno log yM ay 2005 用地层元素测井(ECS )资料评价复杂地层岩性变化程华国,袁祖贵(胜利石油管理局石油工程技术管理处,山东东营257001) 摘要:在石油测井中,地层的岩性判别是评价储层参数的首要条件,地层元素测井(ECS )能测出地层中Si 、Ca 、Fe 、S 、Ti 、Cl 、Cr 、Gd 等元素的含量,结合地质录井等资料可准确确定储层的岩性。

通过实例分析,正确评价了复杂地层的岩性变化,为石油勘探开发提供了可靠的地质参数。

关键词:地层元素测井(ECS);岩性;测井解释;γ能谱;元素分析中图分类号:　T E 151 文献标识码:　A 文章编号:　0258-0934(2005)03-0233-06收稿日期:2004-10-15作者简介:程华国(1954-),男,安徽安庆人,高级工程师,从事石油工程技术研究工作随着油气田勘探开发的不断深入,相对简单和整装的油气藏越来越少,非常规储集层如火成岩、变质岩等的研究与评价越来越受到重视,而火成岩、变质岩等复杂储层的岩性识别是石油测井解释中的难题之一。

斯伦贝谢公司在本世纪向中国市场推出了一种新型的测井仪器—地层元素测井(ECS :Elemental Capture Spectro sco py ),并在中国的东北、西北和东部等油田和地区进行了测井,在岩性识别上取得了令人满意的效果。

1　ECS 测井的核物理基础利用快中子和地层中的原子核发生非弹性碰撞,碰撞的同时会发射非弹性散射γ射线,γ射线的能量和被碰撞核的核结构有关,它表征了原子核的性质。

同一种原子核在同快中子发生(n,n ′)反应中,所放出的非弹性散射γ射线的能量和数量都是一定的。

对不同的核在(n,n′)反应中放出的γ射线的能谱分析,可确定在地层中存在哪些原子核,它们的含量是多少。

地层元素测井解释方法研究

地层元素测井解释方法研究喻益明【摘要】配合脉冲中子地层元素测井仪的应用,开展地层元素测井解释方法研究,获取地层元素标准谱库,编制地层元素解释软件,解释得到元素产额、干重(不含结晶水的质量百分含量)、矿物组分含量.取芯井对比结果显示,解释结果与取芯化验结果一致性较好,研究结果可为复杂岩性评价、准确计算孔隙度提供数据.%In order to apply pulse neutron formation elements logging tool,the research on the interpretation method of the formation elements logging was performed.A standard spectrum library of the formation elements was obtained,and an interpretation software was compiled,through which,the element yield,dry weight (mass percent content without crystalline water) and mineralogical composition could be paring the interpretation results with the data from the cored well showed that the interpretation results had a good consistency with the core test results,which provided data bases for the evaluation of complex lithology and accurate calculation of porosity.【期刊名称】《同位素》【年(卷),期】2017(030)003【总页数】7页(P187-193)【关键词】地层元素;能谱测井;标准谱;矿物组分【作者】喻益明【作者单位】大庆油田测试技术服务分公司,黑龙江大庆 163453【正文语种】中文【中图分类】TL99;TE19地层元素能谱测井是利用热中子与地层各种元素的原子核发生辐射俘获核反应，根据瞬发的伽马射线能量不同，探测和记录俘获伽马能谱，从而定性、定量确定地层中的元素成分和含量。

地层元素测井仪器结构参数的蒙特卡罗数值模拟

地层元素测井仪器结构参数的蒙特卡罗数值模拟严慧娟;岳爱忠;赵均;张建民;王树声;林岩栋【摘要】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 ofTungsten(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 11 cm thick of W-Ni-Fe material is above 0. 9; if the Boron coating is the mixed material of 10B and fluoro-elastomer, its 10B 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.%建立地层元素测井的蒙特卡罗方法数值模拟模型,分别对仪器的屏蔽体、硼套、源距等仪器参数对测井响应的影响进行数值模拟研究.单一屏蔽体材料中钨镍铁屏蔽率最高,组合材料中钨镍铁十硼砂和钨镍铁十聚乙烯的屏蔽率比相同厚度钨镍铁的屏蔽率略高一些,钨镍铁十聚乙烯的效果更明显.实际仪器模型中11cm厚的钨镍铁屏蔽率大于0.9；硼套为10B与氟橡胶混合材料时其10B含量达到10％以上,达到降低仪器内部元素影响的效果；仪器外壳材料均可以采用沉淀不锈钢或铜,其对地层元素的解谱几乎没有影响.这些模拟结果可以为优化仪器结构参数提供参考.【期刊名称】《测井技术》【年(卷),期】2012(036)003【总页数】4页(P282-285)【关键词】地层元素测井;蒙特卡罗方法;数值模拟;仪器结构【作者】严慧娟;岳爱忠;赵均;张建民;王树声;林岩栋【作者单位】西安交通大学,陕西西安710049;中国石油集团测井有限公司技术中心,陕西西安710077;西安交通大学,陕西西安710049;西安交通大学,陕西西安710049;中国石油集团测井有限公司技术中心,陕西西安710077;中国石油集团测井有限公司技术中心,陕西西安710077【正文语种】中文【中图分类】P631.84地层元素测井通过测量中子在地层中产生的俘获伽马能谱和非弹性散射伽马能谱获取地层岩石骨架的主要元素种类和含量，确定岩性和骨架参数。

地层元素测井技术的发展及其应用

谱进行进一步的分析才能获得地层中各元素的信息。

这就需要对获得的数据进行一定的处理，否则这些数据依然无法被正确的运用。

2.1 建立标准的伽马能谱利用地层元素测井技术对地层元素进行探测，首先就要建立地层中常见元素的标准伽马能谱，只有建立了这个标准伽马能谱才能进一步对仪器采集的数据进行详尽的分析。

标准的伽马能谱需要通过相关的实验，然后利用数值模拟方法进行数据处理后获得。

在进行地层元素测井的过程中，需要利用地层中常见元素的标准伽马能谱，需要以此为依据对测井仪器进行校正，为测量数据的处理提供依据。

通常情况下，在利用数值模拟方法对数据进行处理，获得标准伽马能谱时，科研人员要根据实际的情况进行计算模型的制定，这样可以更加准确的获得标准的伽马能谱[1]。

2.2 地层元素产额相对产额在底层元素测井过程有着非常重要的作用。

这个数据主要表达了单个元素发出的伽马光子在发射的总伽马光子中的贡献。

如果可以获得某个元素独自存在时，伽马光子的实际数据，然后再获得混合存在时的伽马光子能谱，这样就可以分别获得各个元素发射的伽马光子在伽马光子总数中的贡献。

各个元素独自存在时的伽马光子谱线被称为单原子标准谱，而与之对应的是各元素混合存在的伽马光子谱，也就是混合谱。

在进行解谱的过程中，研究人员通常会将这两者进行归一处理，这样的化相对产额就会拥有百分比的含义。

相对产额的求解，需要单元素标准谱，实际的混合谱以及解谱的算法。

这三者缺一不可。

首先单元素标准谱需要根据实际情况选择地层中含量比较大且对相关工作影响较大的元素，这就需要通过科学的方式获取氮元素标准谱，通常情况下需要进行相关的实验然后通过数学方法获取。

其次就是确定解谱算法，由于测得的数据是通过一定的方式表现在仪器上，从中找出对应的关系。

2.3 灵敏度因子灵敏度因子同样也是这项技术中非常重要的数据之一，它充分表现了每种元素对于热中子的吸收程度。

对于每种不同的元素来说其灵敏度因子也存在着十分巨大的差异，这个参数只与探测仪器以及元素本身有关，地层对这个参数的大小没有任何影响。

基于地层元素测井资料的沉积岩、岩浆岩识别方法

基于地层元素测井资料的沉积岩、岩浆岩识别方法王淑真;李蓓【摘要】地层元素测井技术对沉积岩、岩浆岩岩性识别具有独特意义.通过分析收集的沉积岩、岩浆岩数据,提出利用Si、Ca、Al三种元素建立岩性识别图版,并利用Si、Ca、Al、Fe四种元素,建立神经网络模型.利用元素含量交会法建立岩性图版能准确划分沉积岩与岩浆岩的线性部分,BP神经网络对于解决元素含量与沉积岩、岩浆岩岩性之间的非线性关系具有独特优势.本文基于地层元素测井资料,利用岩性识别图版与BP神经网络,建立沉积岩、岩浆岩识别模型,用未参与建立模型的数据进行验证,证明模型的可靠性.本文的研究对于掌握元素含量识别岩性技术,拓展地层元素测井的应用能力具有重要的意义.【期刊名称】《地下水》【年(卷),期】2017(039)006【总页数】3页(P110-112)【关键词】沉积岩;岩浆岩;地层元素测井;元素图版;BP神经网络【作者】王淑真;李蓓【作者单位】西北大学大陆动力学国家重点实验室,陕西西安710069;西北大学大陆动力学国家重点实验室,陕西西安710069【正文语种】中文【中图分类】P588地层中元素含量对于判断岩石类型具有重要意义。

目前地层中已发现的元素虽然有100多种，但研究证实，各元素在地壳中的分布是极不均匀的，只相对集中于少数几种元素。

其中O(46.13%)、Si(26.00%)、Al(7.45%)、Fe(4.2%)、Ca(3.25%)、Na(2.4%)、Mg(2.35%)、K(2.35%)、H(1.00%)等9种元素占地壳总质量的98.13%，其余元素仅占1.87%。

因此，只要精确测量到这些主要元素的含量，进而研究所测元素与岩性对应关系，就可以通过元素含量鉴别地壳的岩石类型。

地层元素测井，是通过测量中子与井眼周围地层原子核发生非弹性散射、辐射俘获等核反应后产生的次生伽马能谱，得到地层岩石骨架的Si，Ca，Fe，S，Ti，Gd，Mg，K，Mn，Al等十余种元素的含量，进而确定矿物含量和识别岩性[1～4]。

应用奇异值分解算法的核磁共振测井解谱方法

应用奇异值分解算法的核磁共振测井解谱方法
王才志;尚卫忠
【期刊名称】《石油地球物理勘探》
【年(卷),期】2003(038)001
【摘要】由于核磁共振测井数据的信噪比很低,能否很好地降低噪声的影响成为核磁共振测井解谱方法有效与否的关键因素.本文介绍了在解谱中将特征矩阵进行奇异值分解,然后截去小的非零奇异值的方法,以增加特征矩阵的稳定性,从而保证数据在低信噪比时仍然可以得到稳定的弛豫频谱.本方法在实际应用中取得了很好的效果.
【总页数】4页(P91-94)
【作者】王才志;尚卫忠
【作者单位】中国石油勘探开发研究院;中油国际工程有限责任公司
【正文语种】中文
【中图分类】TE8
【相关文献】
1.一种地层元素测井解谱方法的数值模拟研究 [J], 何彪;李艳婷;邱益香;张悦
2.基于奇异值分解的核磁共振测井T2谱反演方法的改进 [J], 李鹏举;葛成;孙国平;陈新;王彦凯
3.Cif2000平台下的核磁共振测井解谱方法研究 [J], 王才志;李宁
4.地层元素测井标准谱建立和解谱方法研究 [J], 徐琳
5.核磁共振共轭梯度解谱方法研究 [J], 李鹏举;宋延杰;俞军;田宏娟;仲从存
因版权原因，仅展示原文概要，查看原文内容请购买。

地层元素测井_ECS_一种全面评价储层的测井新技术_袁祖贵

第38卷增刊原子能科学技术Vol .38,Suppl .　2004年7月Atomic Energy Science and TechnologyJuly 2004地层元素测井(ECS )———一种全面评价储层的测井新技术袁祖贵1,2,成晓宁3,孙娟4(1.石油大学资源与信息学院,北京　102249;2.中国石化集团胜利石油管理局石油工程技术管理处,山东东营　257001;3.胜利石油管理局档案管理中心,山东东营　257001;4.吉林石油集团公司测井公司,吉林松原　138003)摘要:以胜利油田王庄稠油油藏为例,简要介绍地层元素测井的原理及方法,分析地层元素与ECS (ele -mental capture spectroscopy )测井的关系,论述ECS 数据在地质与石油工程上的诸多应用,为全面评价储层提供了一种新的技术与方法。

地层元素测井(ECS )具有广泛的应用前景。

关键词:地层元素测井(ECS );测井解释;王庄地区;储层评价中图分类号:T L816.3;T L151 文献标识码:A 文章编号:1000-6931(2004)S0-0208-06ECS ———New Well -logging Techniquefor Completely Evaluating the ReserviorYUAN Zu -gui 1,2,CHENG Xiao -ning 3,SUN Juan 4(1.Resources and I nformation College ,University of Petroleum ,Beijing 102249,China ;2.Petroleum Engineering &T echnology Management Department ,SIN OPEC Shengli Petroleum Administrative Bureau ,Dongying 257001,China ;3.Archives Management Center of Shengli Petroleum Administrative Bureau ,Dongy ing 257001,China ;4.J ilin Petroleum Group Well -logging Company ,Songyuan 138003,China )A bstract :　Elemental capture spectroscopy (ECS )is a new w ell -logging technique w ith the ex tensive applied foreground .The paper takes Shengli oil -field with the viscous oil reservoir in Wangzhuang for ex ample ,introduces the principle and method of ECS ,analyzes the rela -tio n betw een of the stratum chemical element and ECS ,and discusses the applications of ECS data on the geology and petroleum engineering ,in o rder to provide a new technique and method for completely evaluating the reservoir .Key words :element capture spectroscopy (ECS );w ell -logging interpretation ;Wangzhuang area ;reservoir evaluation收稿日期:2004-04-15;修回日期:2004-05-08作者简介:袁祖贵(1957—),男,湖北钟祥人,高级工程师,博士研究生,地球物理专业王庄油田是我国大型的稠油油藏之一,为了充分认识和整体解剖王庄大型稠油油藏,在该地区应用了先进的测井技术,诸如核磁共振、电阻率成象、地层元素测井ECS(elemental cap-ture spectroscopy)等。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

能量／ＶＭｅ
图２地层模拟谱
３解谱
早期的地层元素测井仪＿元素测量采用的脉冲１］
材料。为提高计算效率，地层划分成网格状；将井眼中充满清水；层尺寸满足无限空间要求。地考虑实际刻度过程的实施，地层模型以中子及
１２３４５６７８９
ｌ模型设计
图ｌ是蒙特卡罗数值模拟计算的几何模型图。探测器是ＢＧＯ晶体；测器与源之间选用屏蔽材探
料以阻止源发射的中子及伴随伽马射线直接照射；测器对应仪器壳外部包裹吸收热中子的硼探
ｒｔｍ．ＴｈｌｏｉｍｓｓＦ８ｃｕｔｇｆｃｏ，ａｄＧａｓｉｎｓｒｔｈｉｍａｅａｃｒｉｇｔＧＯｉｈｅａｇｒｔｈｕｅｏｎｉａｔｒｎｕｓａｔｅｃｓｎｄｃｏｄｎｏＢｄｔｃｏ ’ ｎｒｙｒｓｌｔｎｅｅｔｒＳｅｅｇｅｏｕｉ．Ｆｉｈｗｓｔａｈａｔｒｏｇ２ｓｏｈｔｔｅｃｐｕｅＧＲｐｃｒｍｅｋｆＣａｅｅｎｎｓｅｔｕｐａｓｏｌｍｅｔｉｌｅｔｎｎｉｉａｄｔｎｒｅｙｏｖｏｓＳｅｔｕｄｃｎｏｕｉｎＷｅｇｔｐｒｅｔｇｆｉｓｏｅａｄＳｎｓｎｓｏｅａｅｖｒｂｉｕ．ｐｃｒｍｅｏｖｌｔ．ｍｏｉｈｅｃｎａｅｏ
密度测井仪器刻度中常用的砂岩、石灰岩、白云岩地层为３主要岩性，在这些模型井中参入适量的种再ＡｌＫ、ｅＳＴ、、Ｆ、、ｉＧｄ等元素。
中子发生器中子源发射１Ｖ中子，４Ｍｅ诱发地层的俘获伽马谱，通过俘获伽马谱测量地层的俘获伽马
Ｈ，Ｓ。ＣａａｄＦｅｈｖｕｈｂｔｅｏｒｌｔｎｗｉｈｈｏｅｉａａｕｓｉｎａｅｍｃｅｔｒｃｒｅａｉｔｔｅｔｅｒｔｃｌｌｅ；ＣｏｔｎｆＭｇｉｏ — ｏｈｖｎｅｔｏｄｌｎｏｔａｅｔｒｃｒｅａｉｎｍｉｅｈｓｂｔｅｏｒｌｔ，ｂｔｏｈｒｅｅｎｏ ’ ａｅｇｏｏｒｌｔｎＴｈｔｏｓｎｏｕｔｅｌｍｅｔｄｎｔｈｖｏｄｃｒｅａｉ．ｏｅｍｅｈｄｕｉｇｍｕｔｌｅｒｓｉｎｍｅｈｄｔｉｕｅｏｔｔｅｅｅｅｔｙｄｒｃｌｃｕｒｈｂｏｕｅｃｎｅｔｏｌｐｅｒｇｅｓｏｔｏｏｆｒｕｈｌｍｎｓｍａｉｅｔｙａｑｉｅｔｅａｓｌｔｏｔｎｆｉｇ
何彪，李艳婷，邱益香，张悦
（中国石油集团测井有限公司技术中心，陕西西安７０７）１０７
摘要：针对基于ＡｍＢ－ｅ中子源测量中子俘获谱确定地层元素的测井方法，用蒙特卡罗数值模拟方法模拟测量响应谱，究一种简化的解谱方法，研直接求解地层元素的重量百分含量。模型参数设计包含砂岩、灰岩和白云岩等石３种典型的地层，再参入适量的ＡｌＫ、ｅＳＴｉ、Ｆ、、、Ｇｄ等元素。解谱计算设计了含有Ｈ、ｉｃ、ｅＡｌＳＫ、、、Ｓ、ａＦ、、、ＭｇＴｉＧｄ等１Ｏ种元素的地层；地层孔隙度从０１～３。Ｈ、ｉａＦ模拟解析结果与理论值有很好的相关性，ａ．４Ｓ、、ｅＣｃ的
第３卷５
增刊
测
井
技
术
Ｙ
ＶｏＩ５ＮｏＳｌ３．Ｏ
Ｄｅ２ｃ０１１
２１年ｌ０１２月文章编号：０４１３（０１Ｓ一６２０１０ —３８２１）０４ —４Ｏ
பைடு நூலகம்
ＷＥＬＬＬＯＧＧＩＮＧＴＥＣＨＮ０Ｉ
一
种地层元素测井解谱方法的数值模拟研究
谱；然后利用氧化物闭合模型求得地层元素Ｓ、ａｉＣ、
Ｆ、ｉＴｉＧｄ等元素的含量［ｅＳ、、。
斯伦贝谢、里伯顿等公司的地层元素测井哈
仪都采用了Ａｍ－ｉＢ同位素中子源替代加速器中
子源，比而言，器测量方式简单，相仪更加稳定，仪器长度缩短。中子源发射中子产额稳定，射出中子发
ＨＥｉｏ，ＬＩＹａｔｎＢａｎｉｇ，ＱＩＹｉｉｎＵｘａｇ，ＺＨＡＮＧｅＹｕ
（ｃｎｃｌＣｅｔｒＴｅｈｉａｎｅ，ＣｈｎｔｏｅｍｇｉｇＣＯ．ＬＴＤ．，ＸｉａｉａＰｅｒｌｕＬｏｇｎ ’ ｎ，Ｓｈａｘ１０７ａｎｉ７０７，Ｃｈｎｉａ）
ｅｃｌｍｅｎｔｔｎｄｒｐｅｔｕ．ａｈｅｅｎｔａｄｉｓｓａａｄｓｃｒｍ
Ｋｅｒｓｌｇｉｇｍｅｈｄｌｇ，ｆｒａｉｎｅｅｅｔａｔｒｙｗｏｄ：ｏｇｎｔｏｏｏｙｏｍｔｏｌｍｎ，ｃｐｕｅＧＲ，ＭｏｔｒｏｎｍｅｉａｉｕｎｅＣａｌｕｒｃｌｍ — ｓ
ｔｎｐｏｅｗｅｎｅｃｉｌｔｄｅｅｅｔｃｎｅｔｎｈｏｅｉａａｕ，ｗｅｃｎｓｅｔｅｃｎｅｔｆｉｌｔｂｔｅａｈｓｍｕａｅｌｍｎｏｔｎｄｔｅｒｔｃｌｌｅｏａｖａｅｈｏｔｎｓｏ
量谱中非弹性伽马谱所占的比例比较小。按照常规
作者简介：何彪，，９３男１６年生，高级工程师，从事核测井方法及试验。
第３卷５
增刊
何彪，：等一种地层元素测井解谱方法的数值模拟研究
・６３・４
的解谱方法，要解谱求取各元素的相对产额，通先再过氧化物闭合模型求取各元素的绝对含量。但求取各元素的相对产额就要获取其标准谱，个的解谱整过程相对要求较高。本文通过蒙特卡罗数值模拟计算得出一些典型地层的混合谱，直接解谱获取各元素的绝对含量。
ｅｃｌｍｅｔｉｄｒｖｄｆｏｅｕｔｏｉｅｈａｅ．Ｔｅｏｅｅｅｅｔｕｈａ，Ｓ，ＣａａｈｅｅｎｓｅｉｅｒｍｑａｉｎｇｖｎｉｔｅｐｐｒｎｎｍｒｌｍｎｓｓｃｓＨｉ，Ｆｅ，Ａ１，Ｓ，Ｋ，Ｍｇ，Ｔｉｎｒｉｅｎｔｅｓｅｔｕｄｃｖｌｔｏｌｏｉｍ．Ｆｒｍｏｒｌ — ｄＧｄａｅｇｖｎｉｈｐｃｒｍｅｏｏｕｉｎａｇｒｔａｈｏｃｒｅａ
ＡｂｔａｔＭｏｔｒｏｎｍｅｉａｉｕａｉｎｍｅｈｄｉｕｅｏｃｌｕａｅｓｍｅｔｐｃｌｆｒａｉｎｓｒｃ：ｎｅＣａｌｕｒｃｌｍｌｔｏｔｏｓｓｄｔａｃｌｔｏｙｉａｏｍｔｓｓｏｍｉｅｐｃｒｍ，ａｄｔｅ，ｔｅａｓｌｔｏｔｎｆｅｃｏｍａｉｎｅｅｎＳｏｔｉｅｈｏｇｘｄｓｅｔｕｎｈｎｈｂｏｕｅｃｎｅｔｏａｈｆｒｔｌｍｅｔｉｂａｎｄｔｒｕｈｏｓｅｔｕｄｃｎｏｕｉｎＭｏｅｅｉｎＭｏｔｒｏｓｍｕａｉｎｇｏｔｉｄｌｉｇｖｎ（ｅｐｃｒｍｅｏｖｌｔ．ｏｄ１ｄｓｇ．ｎｅＣａｌｉｌｔｏｅｍｅｒｃｍｏｅｓｉｅｓｅＦｉ１ｇ）．ｉｉｈｔｅｄｔｃｏｓＢｎｗｈｃｈｅｅｔｒｉＧＯｒｓａ，ａｄｔｅｍｏｅ１ｃｎｉｔｆｓｎｓｏｅ，Ｉｅｃｙｔｌｎｈｄｌｗｅｌｏｓｓｓｏａｄｔｎｓｉ — ｍｓｏｅｎｏｏｔｓｗｉｈａｐｏｅｕｎｉｆ１ｔｎｓａｄｄｌｍｉｅｔｒｐｒｑａｔｙｏ，Ｋ，ＦｅｔＡ，Ｓ，Ｔｉｎ，ｅｃｉｕａｉｎａｇ — ，ａｄＧｄｔ．Ｓｍｌｔｏｌｏ
能量较低，主要是４５Ｍｅ能量分布的中子。虽～Ｖ
ｌｔｎ，ｓｅｔｕｄｃｎｏｕｉｎａｉｏｐｃｒｍｅｏｖｌｔｏ
０引言
地层元素测井采用同位素中子源Ａｍ－ｅ通过Ｂ
中子与地层的作用产生非弹性伽马谱及俘获伽马谱
获取地层信息。同位素中子源的能量相对较低，测