1、测 井 技 术WELL LOGGING TECHNOLOGYVol.46 No.6 Dec 2022第46卷 第6期 2022年12月文章编号：10041338（2022）06068306基于直接解调法的可控源地层元素 测井联合解谱方法研究王树声1，2，岳爱忠1，2，田文新3，李晓1，2，何绪新1，2，王春林1，2（1.中国石油集团测井有限公司测井技术研究院，陕西 西安 710077；2.中国石油天然气集团有限公司测井技术试验基地，陕西 西安 710077；3.中国石油集团测井有限公司物资装备公司，河北 廊坊 065007）摘要：可控源地层元素测井可以对地层中主要元素含量和矿物含量进行定量测量
2、，获得地层中的碳、氧元素的含量和总有机碳含量（TOC）。矿物含量和TOC对复杂储层和页岩油气等非常规油气的测井评价具有重要意义，而碳、氧元素的定量测量可以在碳捕获、利用与封存中识别和量化近井区中的二氧化碳含量。提出一种基于直接解调法的可控源地层元素测井解谱方法，使用该方法对非弹性能谱和俘获能谱联合解谱获得元素的产额，使用氧化物闭合模型和架桥法获得地层中元素的含量。可控源地层元素测井仪在中国石油测井计量站的地层元素标准井群中进行了一系列的方法试验，结果表明，使用基于直接解调法的可控源地层元素测井联合解谱方法可以有效地提高能谱解析的精度。通过在吉林油田的现场试验验证，该方法计算的元素含量和矿物组分
3、含量与区域地质特征相吻合，可准确反映地层岩性变化规律。关键词：直接解调法；可控中子源；元素标准谱；氧化物闭合模型；元素含量中图分类号：P631.84 文献标识码：ADoi:10.16489/j.issn.10041338.2022.06.007Study on the Method of Combined Spectral Processing of Controlled Source Formation Element Logging Based on Direct Demodulation MethodWANG Shusheng1,2,YUE Aizhong1,2,TIAN Wenxin3
4、,LI Xiao1,2,HE Xuxin1,2,WANG Chunlin1,2(1.Logging Technology Research Institute,China National Logging Corporation,Xian,Shaanxi 710077,China;2.Well Logging Technology Pilot Test Center,China National Petroleum Corporation,Xian,Shaanxi 710077,China;3.Material and Equipment Company,China National Logg
5、ing Corporation,Langfang,Hebei 065007,China)Abstract:Controlled source formation element logging can quantitatively measure the content of main elements and minerals in the formation,obtain the content of carbon,oxygen and total organic carbon(TOC)in the formation.The mineral content and TOC are of
6、great significance for log evaluation of unconventional oil and gas reservoirs such as complex reservoirs and shale oil and gas,while the quantitative measurement of carbon and oxygen elements can identify and quantify the carbon dioxide content in the near wellbore area in the carbon capture,utiliz
7、ation and storage(CCUS).We present a spectral processing method of element logging in controllable source formation element based on direct demodulation method in this paper.Using this method,the inelastic energy spectrum and capture energy spectrum are solved together to obtain the element yield.Th
8、e oxide closure model and bridging method are used to obtain the element content in the formation.A large number of method tests have been carried out on the controllable source formation element logging tool in the formation element standard well group of the oil logging metering station.The result
9、s show that the precision of energy spectrum analysis can be effectively improved by using the combined spectral analysis method of the 基金项目：中国石油天然气集团有限公司项目“175/200 h 测井高性能中子管研制”（2020B-3714）；中国石油集团测井有限公司项目“全域测井装备研制与配套”（CNLC2022-2A）第一作者：王树声，男，1980 年生，高级工程师，博士，从事核测井方法与仪器研制工作。E-mail：通信作者：岳爱忠，男，1968 年生，
10、高级工程师，首席技术专家，从事中子管与核测井仪器研制工作。E-mail：2022年测 井 技 术684 controllable source formation element logging based on the direct demodulation method.The field test in Jilin oilfield shows that the element content and mineral component content are consistent with the regional geological characteristics,which ca
11、n accurately reflect the change law of formation lithology.Keywords:direct demodulation method;controlled neutron source;element standard spectrum;oxide closure model;element content0 引 言可控源地层元素测井仪使用可控中子发生器和溴化镧晶体探测器。可控中子发生器向井眼周围地层发射能量为14 MeV的高能中子，中子与周围物质的原子核发生非弹性散射、弹性散射、辐射俘获以及中子活化等核反应。其中在非弹性散射和辐射俘获核
12、反应中会产生次生伽马射线，通过对非弹伽马能谱和俘获伽马能谱的数据进行联合解析，可以获得地层中包括硅、钙、碳、铝、铁等18种地层主要元素的含量1。可控源地层元素测井不但可以对地层中主要元素和矿物进行定量测量，同时还可以获得地层中的碳、氧元素的含量和总有机碳含量（TOC）。矿物含量和TOC对复杂储层和页岩油气等非常规油气的测井评价具有重要意义，因此，可控源地层元素、高分辨率核磁共振和介电测井被称为非常规油气的新三组合测井技术2-5。可控源地层元素测井可以对碳、氧元素进行定量测量，可以在碳捕获、利用与封存（Carbon Capture，Utilization and Storage）中识别和量化近井
13、区中二氧化碳含量6。可控源地层元素测井常用的解谱方法有加权最小二乘法、极大似然估计、有效集法等，本文把直接解调法（Direct Demodulation）用于可控源地层元素测井的能谱解析过程。直接解调法最初主要用于天文物理成像中的二维图像重建，已经成功地用于多种空间高能天体观测7-12。通过在模型井群及现场试验显示该方法可以有效提高能谱解析的精度。1 直接解调法的基本原理可控源地层元素测井仪用来进行能谱解析的谱主要有净俘获谱和净非弹谱。这2种谱可以看作是不同元素标准谱的线性组合，使用矩阵的表达形式为dPXE=+（1）式中，d为能谱计数率组成的向量；P为由元素标准谱组成的矩阵；X为元素产额组成的
14、向量；E 为误差向量。基于最小二乘法，式 （1）的最优解Q（X）就是误差向量E的L2范数|dPX2最小Q XdPXX P PXd PXd d()=+=-+|22TTTT（2）对X求偏导数，可以得到式 （2）的正则方程P PXP dTT=（3）用直接解调法得到式 （3）的第l 次迭代的近似解见式 （4），解式 （4）时第 l 次迭代得到的解向量第i个元素，pl =PTP，c=PTd。xip i ic ip i j xjp i jlljillj iNl()=-()=+()=()()-()()-(),111x xjxill-()-()()+-()()111（4）式中，p i jl,()为系数矩阵pl
15、中第i行、第j列元素；c i()为常数矩阵c的第i个元素；为松弛因子；x(l)、x(l 1)分别为第l次、l1次迭代的近似解。直接解调法实质上是对式 （3）使用逐次超松弛迭代法（Succes-sive Over Relaxation Method，SOR），逐次超松弛迭代法是高斯-塞德尔方法的一种加速方法，是解大型稀疏矩阵方程组的有效方法，其具有计算公式简单、占用计算机内存少等优点，但需要选择好的加速因子。2 可控源地层元素联合解谱方法可控源地层元素测井处理方法的步骤：对仪器获得净非弹谱和净俘获谱进行数据预处理；通过蒙特卡罗数值模拟结合模型井试验的方式，获得元素的非弹标准谱和俘获标准谱；对净非
16、弹谱和净俘获谱分别使用直接解调法进行解谱获得元素的相对产额；通过对俘获谱使用氧化物闭合模型、对非弹谱使用“伪俘获谱法”，进而获得元素的含量。2.1 元素标准谱的模拟元素标准谱是地层中单一关键元素的原子核与中子发生非弹性散射或辐射俘获形成的伽马能谱，其包括俘获元素标准谱和非弹元素标准谱。每一种核素都会产生1个或者若干个具有特定能量的伽马，如H产生2.23 MeV的伽马，O产生6.13 MeV的伽马，由于闪烁伽马探测器有限的能量分辨率和统计涨落，每种元素标准谱都有一个或者多个展宽的能峰。这王树声，等：基于直接解调法的可控源地层元素测井联合解谱方法研究第46卷 第6期685 些能峰的形状近似以未展宽
17、的伽马能量为中心，对元素的特征能量按高斯分布展宽形成全能峰，类似地形成第1、第2逃逸峰。元素标准谱反映了元素各特征能量对应的伽马计数，是通过解谱反演元素相对产额的基础。元素的标准谱常可以通过实验室测量和蒙特卡罗数值模拟方法获得。实验室测量需建立大量的模型井，不但耗资巨大，而且需要通过复杂的数据处理剔除测量中其他元素的影响。通过蒙特卡罗数值模拟方法可以获得元素非弹标准谱和俘获标准谱（见图1）。COOOSiCaSHSiCaS能量能量相对计数率相对计数率AlAlMgFe图1 可控源地层元素标准谱（a）俘获标准谱（b）非弹标准谱2.2 可控源元素联合解谱方法元素的相对产额反映了元素对测量能谱的贡献，不
18、能直接用来进行岩石物理评价，对于俘获能谱、要通过氧化物闭合模型将元素的产额转化为元素的含量。氧化物闭合模型的基本思想：所有元素的重量百分含量之和为1；地层所有矿物都可以认为是由氧化物或碳酸盐组成，组成矿物的氧化物、碳酸盐含量百分数之和为1。该方法的核心就是用独立的方式对通过热中子辐射俘获核反应测得的每种元素的相对产额重新归一化，从而求得每种元素的百分含量。此模型的优点在于克服了难以定量描述骨架中C、O这2种元素含量的问题，能够直接计算岩石骨架主要元素含量。氧化物闭合模型见式（5），元素的含量Wi计算见式（6）。FAYSiii()=/1（5）WFYSiii=（6）式中，F 为随深度变化的标准化因
19、子，也称为归一化因子；Ai为第i 种元素的氧化物指数；Yi和Si为第i种元素的相对产额和灵敏度因子13。在计算C、O等非弹谱元素含量时，使用“伪俘获谱”或者“架桥法”14-17。由于非弹谱中的元素总数比较少，所以不能构成闭合模型。利用既有俘获反应也发生非弹性散射的Si、Ca、Fe、S等元素，假定这些元素通过俘获反应和非弹性散射得到元素百分含量是相等的，可得到非弹性散射元素的含量计算式WWSYYSIEZIZIZIEIE=()()/（7）式中，W、Y、S 分别为元素的含量、产额和相对灵敏度；I 为对应的非弹性散射的值；Z 为既能发生非弹性散射又能发生辐射俘获反应的元素，如Si、Ca、Fe、S、Mg
20、、Al、Ba等；E 为只发生非弹性散射的元素，如C和O。由式 （7）可得，C、O等元素的含量可以用Si、Ca等元素作为桥梁，通过这些元素的含量、非弹谱产额和非弹灵敏度，以及C、O等元素的非弹产额和非弹灵敏度获得。2.3 可控源地层元素测井数据处理流程可控源地层元素仪器所测的仪器谱包括非弹门测得的总谱、俘获门测得的俘获谱和总谱。其中非弹门测得的总谱是这一时间门内非弹谱、俘获谱和本底谱的叠加；俘获门测得的俘获谱是俘获谱和本底谱的叠加。而进行联合解谱所需的谱必须是净俘获谱和净非弹谱，故在解谱前要获得净俘获谱和净非弹谱。其中净俘获谱等于俘获门中得到的俘获谱扣除本底谱的贡献，净非弹谱等于非弹门中得到的总
21、谱减去俘获谱和本底谱的贡献。此外，由于测量的仪器谱和解谱使用的标准谱可能在能量刻度、能量转换为道址时有差异，所以需要对仪器谱和标准谱进行校正。特别重要的是要进行能量分辨率的匹配，即标准谱的能量分辨率要调整匹配2022年测 井 技 术686 到仪器测量谱的能量分辨率。解谱时首先要确定地层岩石骨架中各主要元素在测量谱中的相对产额，也就是确定各种不同元素标准谱在测量谱中的相对贡献。这个过程要对俘获谱和非弹谱同时进行解析，主要算法采用直接解调法。可控源地层元素测井俘获能谱的数据处理流程见图2，非弹谱的处理流程见图3。图4 可控源地层元素测井仪在地层元素模型井群的测量能谱1098765432101001
22、010.10.01能量/MeV计数率/cps近非弹谱 近俘获谱近总谱Fe 7.64Fe 6.03铁井1098765432101001010.10.01能量/MeV计数率/cps近非弹谱 近俘获谱近总谱1098765432101001010.10.01能量/MeV计数率/cps近非弹谱 近俘获谱近总谱Si 3.54Si 4.93Si 1.78 非弹砂岩 13 号灰岩 26 号白云岩 10 号1098765432101001010.10.011E-3能量/MeV计数率/cpsO 近非弹谱 H 近俘获谱近总谱碳井1098765432101001010.10.011E-3能量/MeV计数率/cpsO
23、近非弹谱 H 近俘获谱近总谱O 6.13C 4.43水井1098765432101001010.10.011E-3能量/MeV计数率/cps近非弹谱 近俘获谱近总谱Ca 4.42 Mg 3.92Ca 6.42Ca 1.94 Ca 6.42Ca 4.42Ca 1.94（a）砂岩井（b）灰岩井（c）白云岩井（d）水井（e）碳井（f）铁井仪器俘获谱 获得净俘获谱 解谱获得俘获谱产额 仪器本底谱 数据预处理 俘获灵敏度 氧化物闭合模型 俘获谱元素含量 图2 可控源地层元素俘获能谱的数据处理流程仪器总谱 获得净非弹谱 解谱获得非弹谱产额仪器俘获谱 数据预处理 非弹灵敏度 架桥法 架桥元素信息 非弹谱元素
24、含量 图3 可控源地层元素非弹能谱的数据处理流程3 测井试验3.1 模型井群试验可控源地层元素测井仪在中国石油测井计量站的地层元素刻度模型井群进行了试验，在灰岩井、砂岩井、白云岩井、水井、铁井、铝井、镁井、碳井等进行了实测谱的测量，部分结果见图4。可控源地层元素测井仪在花岗岩、大理石等天然岩石井进行了元素含量计算方法的试验。元素含量的计算使用了基于直接解调法的可控源地层元素测井联合解谱方法（见表1），结果显示该方法可以有效地提高能谱解析的精度。表 1 可控源地层元素测井仪在花岗岩模型井实测元素含量元素含量/%SiAlFeMgCaNaKSTiC标称值32.926.711.100.290.862.
25、022.070.010.180.23测量值33.667.960.750.002.361.123.150.210.100.12测量偏差0.741.250.350.291.500.901.080.200.080.11王树声，等：基于直接解调法的可控源地层元素测井联合解谱方法研究第46卷 第6期687 3.2 现场试验可控源地层元素测井仪于2020年7月在吉林油田松原井进行了现场试验。图5是测井曲线重复性对比分析，418道为地层中主要元素Si、Ca、Fe、Al、S、Ti、K、Mg、Gd、Mn、Na、Ni、Cu、C的含量曲线。分析表明，可控源地层元素含量曲线重复性较好，基本达到设计指标，层位特征准确，
26、符合地层响应规律。基于元素含量特征，进行了目的层段地层岩性精细解释，通过最优化算法完成矿物组分精细及TOC定量解释（见图6）。地层特征以砂泥岩为主，矿物组分主要包括伊利石、石英、钾长石、钠长石、方解石、铁白云石、黄铁矿。以泥岩夹泥质粉砂岩为主，泥质含量平均41%，石英含量平均26%，长石类型以钾长石为主，含少量斜长石，方解石含量5%10%，个别层段含铁白云石和黄铁矿，TOC平均1.2%3.2%。基于元素含量处理得到的矿物组分含量与区域地质特征相吻合，能够准确反映地层岩性变化规律，TOC有助于进行有机质甜点识别，对于提高综合储层评价效果有重要意义。图6 可控源地层元素测井仪在吉林油田松原井测井元
27、素含量曲线及矿物剖面自然伽马 深度密度中子孔隙度硅含量钙含量铁含量铝含量硫含量钾含量锰含量钠含量镁含量钛含量钆含量碳含量总有机碳矿物组分谱GR/APIDEN/（gcm-3）CNL/p.u.深度/mSi/（ww-1）Ca/（ww-1）Fe/（ww-1）Al/（ww-1）S/（ww-1）K/（ww-1）Mn/（ww-1）Na/（ww-1）Mg/（ww-1）Ti/（ww-1）Gd/（mgL-1）C/（ww-1）TOC/（vv-1）黄铁矿铁白云石方解石钠长石钾长石石英伊利石020013010.5000.4 00.100.200.2 00.02 00.1 00.200.2 00.502500.0500.
28、152 2702 2802 2902 300图5 可控源地层元素测井仪在吉林油田松原井测井曲线重复性对比分析深度/m2001005050500500500200300250200200250502.651.65-1545200CCuNiNaMnGdMgKTiSAl FeCaSiDENCNLGR（mgL-1）%（gcm-3）API2001005050500500500200300250200200250502.651.65-15452002 2702 2802 2902 3002 310CCuNiNaMnGdMgKTiSAl FeCaSiDENCNLGR（mgL-1）%（gcm-3）API202
29、2年测 井 技 术688 4 结 论（1）本文采用直接解调法进行可控源元素测井的解谱，其实质上是对正则方程组的求解，通过逐次超松弛迭代法完成了解谱，获得元素的产额，与通常的加权最小二乘法相比，元素含量的求解精度得到了提高。（2）建立基于直接解调法的可控源元素测井数据处理流程，通过蒙特卡罗数值模拟获得了可控源地层元素测井的非弹标准谱和俘获标准谱。（3）通过在中国石油测井计量站的模型井群试验表明，使用直接解调法的可控源地层元素测井联合解谱方法可以有效地提高能谱解析的精度。可控源地层元素测井仪在吉林油田现场试验结果表明元素含量和矿物组分含量与区域地质特征相吻合，可准确反映地层岩性变化规律。参考文献：
30、1 RADTKE R J,LORENTE M,ADOLPH B,et al.A new capture and inelastic spectroscopy tool takes geochemical logging to the next level C/SPWLA 53rd Annual Logging Symposium,Cartagena,Colombia,2012.2 袁超，周灿灿，薄亚杰，等.直接确定地层有机碳含量测井方法的数值模拟J.原子能科学技术，2015，49（8）：1497-1503.3 邹才能，翟光明，张光亚，等.全球常规-非常规油气形成分布、资源潜力及趋势预测J.石油
31、勘探与开发，2015，42（1）：13-25.4 FENG Z,LI X T,WU H L,et al.Multimineral optimization processing method based on elemental capture spectroscopy logging J.Applied Geophysics,2014,11(1):41-49,116.5 贾承造，郑民，张永峰.中国非常规油气资源与勘探开发前景J.石油勘探与开发，2012，39（2）：129-136.6 QUINTERO L F,GUO W J,GALES R.Characterization of pulse
32、d responses to Monitor CCUS projects C/SPWLA 63rd Annual Logging Symposium,Norway,2022.7 李惕碚，吴枚.高能天文中成像和解谱的直接方法J.天体物理学报，1993，13（3）：215-224.8 张家宇，王焕玉，张承模，等.直接解调方法在嫦娥一号卫星X射线谱仪地面验证实验数据分析中的应用J.光谱学与光谱分析，2008，28（10）：2450-2454.9 艾宪芸，魏义祥，肖无云.直接解调法在碲锌镉探测器谱分析中的应用J.清华大学学报（自然科学版），2006，46（6）：821-824.10 WU W S,NI
33、U W,LUO L.Quantitative analysis of neutron-capture gamma-ray energy spectra using direct demodulation J.Geophysics,2014,79(2):91-98.11 SHEN Z J,ZHOU J F.Accelerated direct demodulation imaging method J.High Energy Physics and Nuclear Physics,2007,31(11):1016-1021.12 沈宗俊.直接解调方法的研究与应用D.北京：清华大学，2007.13
34、 GALFORD J,QUIREIN J,SHANNON S,et al.Field test results of a new neutron induced gamma ray spectroscopy geochemical logging tool C/SPE Annual Technical Conference and Exhibition,New Orleans,USA,2009.14 HERRON M M,GRAU J A,HERRON S L,et al.Total organic carbon and formation evaluation with wireline l
35、ogs in the green river oil shale C/SPE Annual Technical Conference and Exhibition,Denver,USA,2011.15 GONZALEZ J,LEWIS R,HEMINGWAY J,et al.Determination of formation organic carbon content using a new neutron-induced gamma ray spectroscopy service that directly measures carbon C/SPWLA 54th Annual Log
36、ging Symposium,New Orleans,USA,2013.16 GRAU J A.Combining inelastic and capture gamma ray spectroscopy for determining formation elemental weight fractions:EP2791711 P.2015-09-20.17 HERRON S,GRAU J,ROSCOE B,et al.Methods and apparatus using both nuclear capture and inelastic spectroscopy in deriving elemental concentrations:US7366615 P.2008-04-29.（修改回稿日期：2022-09-20 编辑 林丽丽）
