CN101354362B - 石油钻井中x射线荧光泥质含量的分析方法 - Google Patents

石油钻井中x射线荧光泥质含量的分析方法 Download PDF

Info

Publication number
CN101354362B
CN101354362B CN2008101407236A CN200810140723A CN101354362B CN 101354362 B CN101354362 B CN 101354362B CN 2008101407236 A CN2008101407236 A CN 2008101407236A CN 200810140723 A CN200810140723 A CN 200810140723A CN 101354362 B CN101354362 B CN 101354362B
Authority
CN
China
Prior art keywords
ray fluorescence
shale index
value
analysis
shale
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
Application number
CN2008101407236A
Other languages
English (en)
Other versions
CN101354362A (zh
Inventor
陈英毅
何国贤
周天顺
王飞龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Logging Branch Of Sinopec Jingwei Co ltd
North China Measurement And Control Co Of Sinopec Jingwei Co ltd
Sinopec Oilfield Service Corp
Sinopec North China Petroleum Engineering Corp
Sinopec Jingwei Co Ltd
Original Assignee
Sinopec North China Petroleum Bureau
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sinopec North China Petroleum Bureau filed Critical Sinopec North China Petroleum Bureau
Priority to CN2008101407236A priority Critical patent/CN101354362B/zh
Publication of CN101354362A publication Critical patent/CN101354362A/zh
Application granted granted Critical
Publication of CN101354362B publication Critical patent/CN101354362B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明涉及一种利用X射线荧光分析石油钻井中泥质含量的方法。该方法是采用能量色散型X射线荧光分析仪在石油钻井现场随钻井深度对岩心、岩屑中的化学元素进行连续分析,得出随深度变化的X射线荧光分析数据,采用线性回归方法对所述X射线荧光分析数据进行相关性分析和统计,计算出不同元素之间的相关系数,然后计算某元素与其他11种元素的平均相关系数,从中选出平均相关系数最大的元素;进行泥质含量的计算;绘制出随钻井深度变化的泥质含量曲线。本发明的方法可用于钻井现场,对已钻地层的岩性进行及时、准确地泥质含量分析,并根据泥质含量分析数据进行岩性和储集性能的定性和定量评价,形成基于X射线荧光分析的随钻地质解释评价方法,能正确有效地指导石油天然气勘探开发钻井生产。

Description

石油钻井中X射线荧光泥质含量的分析方法
一、技术领域:
本发明涉及一种石油钻井地质中泥质含量的分析方法,特别涉及一种利用X射线荧光分析石油钻井中泥质含量的方法。
二、背景技术:
在石油钻井地质中,岩石中的“泥质含量”是评价储集层性质的重要参数,及时获得泥质含量参数,可及时评价储集层性质,给油气层保护赢得主动,并且给完井方案的确定和后期勘探开发方案的制定提供可靠数据。
根据GB/T 17412.2-1998标准定义,“泥质”是指沉积岩石中的粘土矿物成分和粒度<0.004mm的碎屑成分,泥质含量>50%的岩石称为泥质岩。泥质含量是评价储集层性质的重要参数,泥质含量越高,储集性能越差,反之,泥质含量越小,储集性越好。大量的实践表明,泥质含量与油气产能有着负相关关系。为了评价储集层性质及计算储量,各油田不同的油气层组,都制定了储集层划分标准,规定了储集层和非储集层的泥质含量界线。如:鄂尔多斯盆地某气田盒2+3气层组,泥质含量<15%为一类储集层,15-30%为二类储集层,30-35%为三类储集层,35-40%为四类储集层,>40%为非储集层。
获得定量的泥质含量参数的传统方法大体分两类:(1)试验室岩石样品化验分析。该方法一般用钻井岩心样品分析,可获得较准确的泥质含量数据。但是,因钻井取心井段少,且不连续,所以获得的泥质含量参数也是不连续的、有限的;另外,因分析周期长,难以满足生产需要。(2)地球物理测井数据计算。该方法是利用地球物理测井(如自然伽马、电阻率等)数据、经验性公式和经验性常数计算获得连续的泥质含量参数,这种方法对系统的解释评价储集层性质很有帮助。但这种方法只能等到完钻(或中途完钻),进行地球物理测井后才能获得泥质含量参数。
传统的钻井地质录井中,没有定量获取泥质含量的方法,只是通过肉眼对岩心、岩屑实物的观察,获得描述性的、定性的泥质含量资料。然而PDC钻头的应用、空气钻井的实施及其他钻井工程工艺技术的发展,造成岩屑细小、甚至呈粉尘状。因此,肉眼观察资料的可靠程度越来越低。
申请号为200710078690.2的专利申请文件公开了一种名称为“石油钻井地质X射线荧光岩屑录井方法”,该文件公开了一种石油钻井地质X射线荧光岩屑录井方法,按钻井深度由浅至深连续取样,通过对样品进行清洗、烘干、研磨、压制成样片;再按钻井深度,用能量色散型X射线荧光分析仪检测样片成份数据;计算机系统按照钻井深度将各岩石成份数据制作成图谱,根据数据及图谱分析获得钻井岩性与层位判断。但该方法主要介绍X射线荧光岩屑录井的工作流程,对经过X射线荧光分析获得的岩石组分(元素)数据的具体应用方法未涉及,更未涉及利用岩心、岩屑X射线荧光分析成果对岩石泥质含量的分析方法。
本发明的X射线荧光分析方法,经中国石油化工股份有限公司石油勘探开发研究院查新表明,利用岩屑X射线荧光进行录井岩性识别和元素分析,除上述公开的专利申请的“石油钻井地质X射线荧光岩屑录井方法”专利外,未见其它公开发表的国内外同类文献资料。
三、发明内容:
本发明要解决的技术问题:克服现有技术的缺点,提出了一种在钻井现场可随着钻井进度实时采集岩心、岩屑样品,实时进行X射线荧光分析和岩石泥质含量的分析方法,并对岩石储集性能进行实时评价。
本发明的技术方案:
本发明的石油钻井地质录井工程中泥质含量的分析方法,其技术核心包括钻井现场的岩心、岩屑X射线荧光分析、化学元素数据统计分析、泥质含量计算及数据处理技术。
本发明的石油钻井中X射线荧光泥质含量的分析方法,包括以下步骤,
(1)采用能量色散型X射线荧光分析仪在石油钻井现场随钻井深度对岩心、岩屑中的化学元素进行连续分析,得出随深度变化的X射线荧光分析数据,所述化学元素为硅、铝、铁、钙、钾、镁、钛、磷、锰、硫、钡和氯;
(2)采用线性回归方法对所述X射线荧光分析数据进行相关性分析和统计;
(3)选出平均相关系数最大的A元素,用该元素的分析数据作为泥质含量计算的基础数据,计算公式为:
Figure S2008101407236D00031
V SH = 2 FC * SH - 1 2 FC - 1 × 100 %
式中:A为A元素X射线荧光分析值(脉冲数),Amax为纯泥岩A元素X射线荧光分析最大值(脉冲数),Amin为纯砂岩A元素X射线荧光分析最小值(脉冲数),FC为泥质含量计算的经验系数,Vsh为泥质含量(%);
(4)进行数据处理,并绘制出泥质含量曲线;
(5)利用泥质含量值和曲线变化进行岩石储集性能的定量评价。
所述的步骤(2)中采用线性回归方法分析X射线荧光数据的步骤为:
(2.1)数据筛选:选择以砂、泥岩为主的X射线荧光分析数据,作为相关性分析的原始基础数据;
(2.2)元素相关性分析与计算:采用一元线性回归方法对不同元素作相关性分析,得出不同元素之间的相关系数,计算公式为:
r = Σ i = 1 n ( x i - x ‾ ) ( y i - y ‾ ) Σ i = 1 n ( x i - x ‾ ) 2 Σ i = 1 n ( y i - y ‾ ) 2
式中:r为相关系数,-1≤r≤1;xi为甲元素测量值(脉冲数);x为甲元素平均值(脉冲数);yi为乙元素测量值(脉冲数);y为乙元素平均值(脉冲数);
(2.3)相关系数统计
先统计出不同元素之间的相关系数,取绝对值,然后计算出某元素与其他元素相关系数的平均值;
(2.4)选择用于泥质含量计算的化学元素
从X射线荧光分析数据中选出与泥质含量有明显正相关关系的元素,然后再从中选出平均相关系数最大的元素。
从大量的X射线荧光分析数据看出,除硅、钙、钡元素之外的其他9种元素含量,都与泥质含量有着较大的正相关关系,从9种元素中挑选出平均相关系数最大的元素,作为计算泥质含量的元素。
所述的明显正相关关系的元素为铝、铁、钾、镁、钛、磷、锰、硫和氯共9种元素。根据鄂尔多斯盆地北部二叠系下石盒子组、山西组和黄骅坳陷东营组元素相关系数统计结果,铁元素一般是平均相关系数最大的元素。地区不同,平均相关系数最大的元素可能不同,如可能是钛或锰等。
数据处理时是将小于0的Vsh值按0计算,将大于100的Vsh值按100计算。
本发明的积极有益效果:
(1)本发明的方法通过钻井岩心、岩屑样品X射线荧光分析,获得岩石样品的X射线荧光图谱,经计算机软件处理获得岩石样品的元素含量,可用于岩性识别。在钻井过程中,利用能量色散型X射线荧光分析仪器,及时对钻井岩心、岩屑样品进行分析,获得岩石样品X射线荧光分析数据,并根据X射线荧光分析结果,及时进行岩性识别和地层判断,并可进行岩石储集性能的定性和定量评价,形成基于X射线荧光分析的随钻地质解释评价方法,能正确有效地指导石油天然气勘探开发钻井生产。
(2)本发明方法利用石油钻井中的岩心、岩屑的化学元素值计算岩石的泥质含量,在钻井现场随着钻井进度实时采集岩心、岩屑样品,实时进行X射线荧光分析和岩石泥质含量分析,实时利用泥质含量数据进行储集层性能评价,因此,该方法可及时地指导油气勘探开发实践,这种实时性是地球物理测井和常规化验分析不可能达到的,具有较大的优越性。
(3)本发明方法利用的岩心、岩屑的化学元素本身就是泥质的重要组成部分,因此求取的岩石泥质含量值比利用地球物理测井信息解释的泥质含量更直接,也更具有代表性和说服力;利用的分析样品主要为岩屑,可按照岩屑取样间距连续取样、连续分析,所获得的泥质含量数据是等间距、系统的分析数据;而常规化验分析化验的只是岩心,得到的泥质含量参数有限且不连续,分析周期较长。本发明对样品只做粉碎、压片和X光照射等物理方法处理,未改变岩石的组分,是对样品的无损分析。因此,同一样品可以进行多次重复性检验,而常规化验分析很难做到重复性检验,地球物理测井也因井筒条件的不断变化难以做到重复性检验。因此,本方法具有系统性上的优势。
(4)本发明利用的分析测量设备和辅助设备都是在地面进行工作的,因此该方法比利用井下分析测量设备的地球物理测井方法的工程技术风险要小得多;分析测量设备和辅助设备的总价格比任何地球物理测井设备的造价都要小得多,设备投资相对较低。因此,本方法具有较好的应用和推广前景。
四、附图说明:
图1:鄂尔多斯盆地某气田某井的X射线荧光泥质含量分析结果与常规录井、电测井及测试结果综合分析对比图。
图2:渤海湾盆地某油田某钻井的X射线荧光泥质含量分析结果与地球物理测井、气测录井结果综合分析图。
五、具体实施方式:
下面结合附图和具体实例对本发明作进一步详细说明。
实施例1:
本发明的石油钻井中X射线荧光泥质含量的分析方法,包括对钻井现场的岩心、岩屑进行X射线荧光分析、化学元素数据统计分析、泥质含量计算及数据处理等过程,本例选取鄂尔多斯盆地某气田的某井X射线荧光分析数据作说明,具体步骤如下:
(一)岩心、岩屑X射线荧光分析
采用能量色散型X射线荧光分析仪,在石油钻井现场对岩心、岩屑样品化学元素进行连续分析。分析时对岩心、岩屑样品先进行粉碎、压片处理,然后进行X射线荧光分析。
1.分析仪器选择
分析仪器采用能量色散型X射线荧光分析仪,现用仪器为四川新先达测控技术有限公司生产的CIT-3000SY型石油X荧光岩屑分析仪。选择依据是仪器的元素分析范围、分析检出限、分析精度、射线能量、检测周期、重复性、稳定性、辐射剂量、中文软件等技术条件的实用性与仪器性价比等方面。
2.分析元素选择
根据地壳中化学元素的克拉克值得知,地壳中主要元素及含量为:氧46.71%、硅27.69%、铝8.07%、铁5.05%、钙3.65%、钠2.75%、钾2.58%、镁2.08%、钛0.62%、氢0.14%、磷0.13%、碳0.094%、锰0.090%、硫0.052%、钡0.050%、氯0.013%,这16种元素占地壳元素总量的99.769%,其他元素仅占0.231%。
元素周期表中1-11号元素的X射线荧光产额低,不易通过能量色散型X射线荧光分析仪测定。因此,X射线荧光分析的元素选择硅、铝、铁、钙、钾、镁、钛、磷、锰、硫、钡和氯,共12种元素。
(二)X射线荧光分析数据统计分析
采用线性回归分析方法,将工区钻井岩心、岩屑的X射线荧光分析数据进行相关性分析,获取最具有代表性的元素种类。
1、数据筛选
选择以砂、泥岩为主的X射线荧光分析数据,然后去掉非砂、泥岩(如煤层、碳酸盐岩、盐岩等)井段数据,筛选后的数据作为相关性分析的原始基础数据。
筛选后的数据如表1所示,表1是鄂尔多斯盆地某气田某井岩屑X射线荧光分析数据表。分析井段为2630-2965米,1点/米,共336点。
分析仪器采用CIT-3000SY型石油X荧光岩屑分析仪,参数设置为:
主板参数:量程2,放大倍数128,测量时间60秒;
控制参数:设置管压120,设置管流20。
2、元素相关性分析
使用一元线性回归方法作不同元素之间的相关性分析,可获得不同元素之间的相关系数。
相关性系数计算公式:
r = Σ i = 1 n ( x i - x ‾ ) ( y i - y ‾ ) Σ i = 1 n ( x i - x ‾ ) 2 Σ i = 1 n ( y i - y ‾ ) 2 · · · ( 1 )
式中:r:样本相关系数,-1≤r≤1;
xi:甲元素测量值(脉冲数);
x:甲元素平均值(脉冲数);
yi:乙元素测量值(脉冲数);
y:乙元素平均值(脉冲数)。
表2为利用表1中数据得出的不同元素之间的相关系数(绝对值)统计表。
3、相关系数统计
先统计不同元素之间的相关系数,然后计算某元素与其他11种元素相关系数(绝对值)的算术平均值。
4、用于泥质含量计算的元素选择
从沉积岩石学理论得知,硅元素为砂质岩的主要成分,硅元素含量代表着砂质含量;钙元素为碳酸盐岩的主要成分,钙元素含量代表着碳酸盐岩含量;钡元素主要矿物为重晶石(BaSO4)和毒重石(BaCO3),在沉积岩中主要存在于碳酸盐岩中。因此,用于泥质含量计算的元素暂不考虑硅、钙、钡这三种元素。
从大量的X射线荧光分析数据可看出,除硅、钙、钡元素之外,其他9种元素含量都与泥质含量有着较大的正相关关系,从这9种元素中挑选出平均相关系数最大的元素,作为泥质含量计算的基础元素。
如表3所示,表3为某元素与其他8种元素的相关系数(绝对值)算术平均值统计表。从表中看出,铁元素的平均相关系数最大,为0.7912。因此,选择铁元素的X射线荧光分析值作为泥质含量计算的基础数据。
(三)X射线荧光泥质含量计算
1、泥质含量计算
选择平均相关系数最大的元素分析数据,作为泥质含量计算的原始基础数据,并用下列公式计算,以铁元素为例:
SH = Fe - Fe min Fe max - Fe min · · · ( 2 )
V SH = 2 FC * SH - 1 2 FC - 1 × 100 % · · · ( 3 )
式中:Fe:铁元素X射线荧光测量值(脉冲数),Femax:纯泥岩铁元素X射线荧光分析最大值(脉冲数),Femin纯砂岩铁元素X射线荧光分析最小值(脉冲数);FC:用铁元素计算泥质含量时的经验系数,鄂尔多斯盆地采用2.0;Vsh:泥质含量(%)。
2、数据处理
纯泥岩元素值和纯砂岩元素值是根据地区研究和该井测量情况选定的常量,当元素测量值大于纯泥岩元素值或小于纯砂岩元素值时,会出现VSH小于0或大于100的情形。若出现此种情形,应首先检查纯泥岩元素值和纯砂岩元素值的取值是否合理,若取值没问题,则直接将小于0的VSH值,作为0处理,将大于100的VSH值,作为100处理。
如表4所示,表4是利用铁元素X射线荧光分析值,采用公式(2)、(3)计算得出的SH、VSH值。其中,公式(2)中纯泥岩铁元素值Femax取600(脉冲数),纯砂岩铁元素值Femin取160(脉冲数),经验系数FC取2.0。
在数据处理中,将VSH值小于0的值,按0值处理;将VSH值大于100的值,按100值处理。
将得到的X射线荧光泥质含量绘制成图,如图1所示,本例的鄂尔多斯盆地某气田某井的X射线荧光泥质含量分析结果与常规录井、电测井及测试结果综合分析对比图。
图1内容包括:1、铁(Fe)元素曲线、按铁元素计算的泥质含量曲线(Fe-SH);2、常规录井的气测全烃曲线、岩性剖面;3、地球物理测井的自然伽马曲线、按自然伽马计算的泥质含量曲线(GR-SH);4、综合解释、射孔位置和测试结果等。
自然伽马测井是成熟的地球物理测井项目,其测量值主要用于泥质含量的解释。从图1中可看出,地球物理测井的自然伽马曲线和X射线荧光分析的铁元素曲线形态相似,变化趋势一致,说明X射线荧光分析的铁元素值也具有解释泥质含量的作用;按自然伽马计算的泥质含量曲线(GR-SH)和按铁元素X射线荧光分析值计算的泥质含量曲线(Fe-SH)形态相似,变化趋势一致。两种方法计算的泥质含量值也相近,地球物理测井2630-2965米井段的泥质含量平均值为52.762%,而该井段的X射线荧光铁元素计算的泥质含量平均值为48.267%,两种方法的相对误差仅为8.519%,二者相对误差较小。
气测全烃值是气测录井工程中利用色谱气测仪器,直接测量钻井液中的天然气烃类气体总量。气测全烃值高,说明地层含气性好。从图1中还可看出,常规录井的气测全烃曲线与铁元素计算的泥质含量曲线(Fe-SH)具有较明显的负相关关系,即泥质含量高,气测全烃值低。这一结果也说明利用X射线荧光元素分析值计算泥质含量的可行性。
利用铁元素计算出泥质含量值较真实地反映了岩石泥质含量情况,这从测试结果也得到很好的证明。图1中看出,射孔井段2733-2737米,泥质含量平均值为6.48%,测试天然气无阻流量为6.095万方/日;射孔井段2875.5-2879米、2903-2906米,泥质含量平均值为9.43%,测试天然气无阻流量为1.521万方/日;射孔井段2949.5-2953.5米,泥质含量平均值为19.06%,测试天然气无阻流量为0.551万方/日;测试证明了泥质含量与天然气产能的负相关关系,也表明利用X射线荧光分析的元素值计算泥质含量的可行性。
实施例2:
如图2所示,图2为渤海湾盆地某油田某钻井根据X射线荧光分析的铁元素计算的泥质含量成果与地球物理测井、气测录井结果综合分析图。从图2中看出,以铁元素计算的泥质含量曲线与地球物理测井的自然伽马曲线具有一定相关性,说明在渤海湾盆地利用铁元素计算泥质含量的可行性;并且以铁元素计算的泥质含量与气测全烃值具有较明显的负相关关系。即泥质含量高,气测全烃值低;反之泥质含量低,则气测全烃值高。因此,从气测录井结果也证明了利用X射线荧光分析的元素值计算泥质含量的可行性。
表1:鄂尔多斯盆地某井钻井岩屑X射线荧光分析数据表
Figure S2008101407236D00111
Figure S2008101407236D00121
续表1
Figure S2008101407236D00131
续表1
Figure S2008101407236D00141
续表1
Figure S2008101407236D00151
续表1
Figure S2008101407236D00161
续表1
  井深〔m〕 Mg Al Si P S Cl K Ca Ba Ti Mn Fe
  2954   3.52   4.61   49.13   2.14   2.17   2.17   27.01   24.59   3.22   14.14   14.73   313.32
  2955   4.08   3.85   46.30   2.40   2.43   2.40   28.01   25.24   3.07   15.37   15.80   337.63
  2956   4.66   6.30   38.13   2.96   6.66   3.50   30.79   19.86   3.08   34.87   15.61   399.20
  2957   4.15   6.79   39.29   2.61   3.47   2.80   36.22   18.43   3.06   38.12   13.71   297.56
  2958   4.42   6.05   42.58   3.03   3.50   2.92   37.91   25.23   3.49   27.80   16.36   409.76
  2959   4.85   5.21   43.93   2.94   3.97   3.45   37.89   26.85   3.70   25.20   20.88   540.80
  2960   4.22   4.50   50.13   2.16   2.51   2.37   31.21   43.02   4.95   20.06   16.01   343.24
  2961   3.72   3.76   49.92   2.58   2.11   2.39   26.58   72.21   8.11   10.41   15.43   295.45
  2962   3.63   4.05   56.42   2.39   2.07   2.55   25.98   56.62   6.12   8.88   13.35   251.99
  2963   3.36   4.15   55.32   2.66   2.29   2.34   25.41   57.76   6.00   9.31   14.00   253.19
  2964   3.65   3.81   53.96   2.38   2.31   2.25   26.18   57.85   6.10   9.68   13.85   255.37
  2965   3.71   3.58   54.45   2.15   2.32   2.23   25.65   57.20   6.31   9.30   13.84   253.32
续表1
表2:硅、铝、铁、钙、钾、镁、钛、磷、锰、硫、钡、氯12种元素相互之间相关系数(绝对值)统计表
  R   Mg   Al   Si   P   S   Cl   K   Ca   Ba   Ti   Mn   Fe
  Mg   1   0.7824   0.6192   0.7208   0.5941   0.8584   0.8103   0.0775   0.1030   0.7951   0.8386   0.8982
  Al   1   0.6951   0.4828   0.5073   0.7082   0.7512   0.2834   0.1404   0.9090   0.6337   0.6888
  Si   1   0.4358   0.5418   0.6526   0.6519   0.2970   0.1507   0.7868   0.5644   0.6833
  P   1   0.5344   0.7964   0.6200   0.0490   0.2311   0.5791   0.6829   0.8040
  S   1   0.6651   0.3803   0.0245   0.1371   0.6012   0.6020   0.6251
  Cl   1   0.7901   0.0447   0.1418   0.7595   0.8049   0.9167
  K   1   0.1158   0.0032   0.7573   0.6253   0.7784
  Ca   1   0.9591   0.2944   0.0173   0.1852
  Ba   1   0.1296   0.2025   0.0141
  Ti   1   0.6700   0.7529
  Mn   1   0.8653
  Fe   1
表3:铝、铁、钾、镁、钛、磷、锰、硫、氯9种元素之间的相关系数算术平均值统计表。
  R   1   2   3   4   5   6   7   8   平均
  Mg   0.7824   0.7208   0.5941   0.8584   0.8103   0.7951   0.8386   0.8982   0.7872
  Al   0.7824   0.4828   0.5073   0.7082   0.7512   0.9090   0.6337   0.6888   0.6829
  P   0.7208   0.4828   0.5344   0.7964   0.6200   0.5791   0.6829   0.8040   0.6526
  S   0.5941   0.5073   0.5344   0.6651   0.3803   0.6012   0.6020   0.6251   0.5637
  Cl   0.8584   0.7082   0.7964   0.6651   0.7901   0.7595   0.8049   0.9167   0.7874
  K   0.8103   0.7512   0.6200   0.3803   0.7901   0.7573   0.6253   0.7784   0.6891
  Ti   0.7951   0.9090   0.5791   0.6012   0.7595   0.7573   0.6700   0.7529   0.7280
  Mn   0.8386   0.6337   0.6829   0.6020   0.8049   0.6253   0.6700   0.8653   0.7153
  Fe   0.8982   0.6888   0.8040   0.6251   0.9167   0.7784   0.7529   0.8653   0.7912
表4:以铁元素计算的泥质含量数据表
Figure S2008101407236D00191
Figure S2008101407236D00201
续表4

Claims (4)

1.一种石油钻井中X射线荧光泥质含量的分析方法,其特征在于:包括以下步骤,
(1)采用能量色散型X射线荧光分析仪在石油钻井现场随钻井深度对岩心或岩屑中的化学元素进行连续分析,得出随深度变化的X射线荧光分析数据,所述化学元素为硅、铝、铁、钙、钾、镁、钛、磷、锰、硫、钡和氯共12种;
(2)采用线性回归方法对所述X射线荧光分析数据进行相关性分析和统计,计算出不同元素之间的相关系数,然后计算某元素与其他11种元素的平均相关系数,从中选出平均相关系数最大的元素A;
(3)用A元素的分析数据作为泥质含量计算的基础数据进行泥质含量的计算,其计算公式为:
Figure FSB00000211636200011
V SH = 2 FC * SH - 1 2 FC - 1 × 100 %
式中:A为A元素X射线荧光分析值,Amax为纯泥岩A元素X射线荧光分析最大值,Amin为纯砂岩A元素X射线荧光分析最小值,FC为泥质含量计算的经验系数,VSH为泥质含量;
(4)进行数据处理,并绘制出随钻井深度变化的泥质含量曲线;
(5)利用泥质含量值和泥质含量曲线变化进行岩石储集性能评价。
2.根据权利要求1所述的石油钻井中X射线荧光泥质含量的分析方法,其特征在于:所述的采用线性回归方法对X射线荧光分析数据进行相关性分析和统计的步骤为:
(1)数据筛选:选择以砂和泥岩为主的X射线荧光分析数据,作为相关性分析的原始基础数据;
(2)元素相关性分析与计算:采用一元线性回归方法对不同元素作相关性分析,得出不同元素之间的相关系数,计算公式为:
r = Σ i = 1 n ( x i - x ‾ ) ( y i - y ‾ ) Σ i = 1 n ( x i - x ‾ ) 2 Σ i = 1 n ( y i - y ‾ ) 2
式中:r为相关系数,-1≤r≤1;xi为甲元素测量值;
Figure FSB00000211636200022
为甲元素平均值;yi为乙元素测量值;
Figure FSB00000211636200023
为乙元素平均值;
(3)相关系数统计
先统计出不同元素之间的相关系数,取绝对值,然后计算出某元素与其他元素相关系数的算术平均值;
(4)选择用于泥质含量计算的化学元素
从X射线荧光分析数据中选出与泥质含量有明显正相关关系的元素,然后再从中选出平均相关系数最大的元素。
3.根据权利要求2所述的石油钻井中X射线荧光泥质含量的分析方法,其特征在于:所述的有明显正相关关系的元素为铝、铁、钾、镁、钛、磷、锰、硫和氯共9种元素,所述的平均相关系数最大的元素为铁元素。
4.根据权利要求1所述的石油钻井中X射线荧光泥质含量的分析方法,其特征在于:所述的数据处理是将小于0的VSH值按0计算,将大于100的VSH值按100计算。
CN2008101407236A 2008-07-21 2008-07-21 石油钻井中x射线荧光泥质含量的分析方法 Active CN101354362B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008101407236A CN101354362B (zh) 2008-07-21 2008-07-21 石油钻井中x射线荧光泥质含量的分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008101407236A CN101354362B (zh) 2008-07-21 2008-07-21 石油钻井中x射线荧光泥质含量的分析方法

Publications (2)

Publication Number Publication Date
CN101354362A CN101354362A (zh) 2009-01-28
CN101354362B true CN101354362B (zh) 2011-01-12

Family

ID=40307273

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008101407236A Active CN101354362B (zh) 2008-07-21 2008-07-21 石油钻井中x射线荧光泥质含量的分析方法

Country Status (1)

Country Link
CN (1) CN101354362B (zh)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102788856B (zh) * 2011-05-17 2015-12-16 中国石油化工股份有限公司 液化石油气中元素硫定性和定量的分析方法
CN102944570B (zh) * 2012-12-02 2014-11-26 天津陆海石油设备系统工程有限责任公司 一种基于x射线元素荧光的岩石样品岩性定名的分析方法
CN103969271B (zh) * 2013-01-30 2016-05-25 中国建材检验认证集团股份有限公司 X射线荧光分析法测定水泥中不同价态硫元素的方法
CN104422706B (zh) * 2013-08-30 2017-03-15 中国石油天然气股份有限公司 一种富含同生泥屑砂岩粘土矿物相对含量的检测方法
CN104089967B (zh) * 2014-07-15 2017-04-19 南京市产品质量监督检验院 一种实体面材产品中铝、钙、硅含量的快速测定方法
CN104612675B (zh) * 2015-02-02 2017-07-14 中石化西南石油工程有限公司地质录井分公司 一种碳酸盐岩地层随钻岩性快速识别方法
CN104730095B (zh) * 2015-02-12 2017-09-19 中国石油大学(华东) 一种基于电成像测井资料的复杂地层粘土含量测量方法
CN105044140B (zh) * 2015-06-30 2018-03-16 北京奥意尔工程技术有限公司 基于x射线荧光元素录井技术的盐底卡取方法
CN106703794B (zh) * 2015-11-12 2021-05-18 中石化石油工程技术服务有限公司 一种碳酸盐岩岩性识别方法
CN106124545A (zh) * 2016-08-31 2016-11-16 中国石油集团川庆钻探工程有限公司 利用岩屑x射线荧光光谱特征进行分层卡层的方法
CN106885815B (zh) * 2017-03-22 2020-01-31 中国石油大学(华东) 一种基于砂岩中元素xrf法判别沉积环境的方法
CN108931545B (zh) * 2017-05-26 2021-07-02 中国石油天然气股份有限公司 矿物种类及含量的确定方法
CN107489417A (zh) * 2017-07-25 2017-12-19 中国海洋石油总公司 基于xrf录井的岩矿组成反演法砂泥岩岩性识别方法
CN109142411A (zh) * 2018-08-07 2019-01-04 巨石集团有限公司 一种玻璃池窑废气治理污泥化学成分检测方法
CN111255435B (zh) * 2020-01-17 2021-09-07 西安石油大学 一种复杂储层泥质含量计算方法
CN111458363A (zh) * 2020-04-21 2020-07-28 中国地质大学(武汉) 一种基于手持式x荧光分析仪快速圈定钛矿床方法
CN114428094A (zh) * 2020-09-25 2022-05-03 中国石油化工股份有限公司 一种岩心矿物组成的分析方法和应用
CN116068000B (zh) * 2021-12-22 2024-05-24 中国石油天然气集团有限公司 一种基于岩屑x射线荧光光谱分析数据拟合伽马曲线的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5778154A (en) * 1993-11-01 1998-07-07 Atlantic Richfield Company Advisor system and method for determining reservoir properties
US6108608A (en) * 1998-12-18 2000-08-22 Exxonmobil Upstream Research Company Method of estimating properties of a multi-component fluid using pseudocomponents
CN101082277A (zh) * 2007-07-05 2007-12-05 北京奥能瑞科石油技术有限公司重庆分公司 石油钻井地质x射线荧光岩屑录井方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5778154A (en) * 1993-11-01 1998-07-07 Atlantic Richfield Company Advisor system and method for determining reservoir properties
US6108608A (en) * 1998-12-18 2000-08-22 Exxonmobil Upstream Research Company Method of estimating properties of a multi-component fluid using pseudocomponents
CN101082277A (zh) * 2007-07-05 2007-12-05 北京奥能瑞科石油技术有限公司重庆分公司 石油钻井地质x射线荧光岩屑录井方法

Also Published As

Publication number Publication date
CN101354362A (zh) 2009-01-28

Similar Documents

Publication Publication Date Title
CN101354362B (zh) 石油钻井中x射线荧光泥质含量的分析方法
CN101344001B (zh) 石油钻井中x射线荧光陆源碎屑岩孔隙度的分析方法
EP3414566B1 (en) Thermal maturity determination of rock formations using mud gas isotope logging
Khaksar et al. Rock strength from core and logs, where we stand and ways to go
Boutt et al. A field study (Massachusetts, USA) of the factors controlling the depth of groundwater flow systems in crystalline fractured-rock terrain
AU2013246056A1 (en) Reservoir and completion quality assessment in unconventional (shale gas) wells without logs or core
CN101082277A (zh) 石油钻井地质x射线荧光岩屑录井方法
CN104514552A (zh) 一种煤层气储层识别与丰度预测的方法
CN103510952A (zh) 碳酸盐岩岩屑岩性组合式鉴定方法
CN107505344A (zh) 利用“最小二乘积”法的岩性解释方法
US20130292111A1 (en) Method of constructing a well log of a quantitative property from sample measurements and log data
US10890066B1 (en) Determination of a rock testability index for formation testing
Jadoon et al. Mineralogical modelling and petrophysical parameters in Permian gas shales from the Roseneath and Murteree formations, Cooper Basin, Australia
CN106484925A (zh) 页岩气水平井压裂选段系统及选段方法
Barson et al. Spectroscopy: the key to rapid, reliable petrophysical answers
Smith Improvements in blast fragmentation using measurement while drilling parameters
CN106907140B (zh) 基于测井资料的水动力封闭型煤层气高产层位确定方法
CN113050168B (zh) 基于阵列声波测井和声波远探测测井资料的裂缝有效性评价方法
FANG et al. Intelligent recognition of tunnel stratum based on advanced drilling tests
AU602656B2 (en) Quantitative determination by elemental logging of subsurface formation properties
Brindle et al. An integrated approach to unconventional resource play reservoir characterization, Thistleton-1 case study, NW England
CN110728074A (zh) 混积细粒岩岩性连续解释的方法及其模型的建模方法
You et al. Stratigraphic identification using real-time drilling data
Nelson et al. Geophysical and geochemical logs from a copper oxide deposit, Santa Cruz project, Casa Grande, Arizona
Zuber et al. Coalbed-methane evaluation techniques-The current state of the art

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: SINOPEC OILFIELD SERVICE CORPORATION

Free format text: FORMER OWNER: CNPC HUABEI PETROLEUM

Effective date: 20150212

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 450042 ZHENGZHOU, HENAN PROVINCE TO: 100029 CHAOYANG, BEIJING

TR01 Transfer of patent right

Effective date of registration: 20150212

Address after: 100029 Beijing City, Chaoyang District Hui Street No. six twelfth floor

Patentee after: SINOPEC OILFIELD SERVICE CORPORATION

Address before: 450042 No. 4 Jia Jia Road, Zhongyuan Town, Zhengzhou City, Henan Province

Patentee before: CNPC Huabei Petroleum

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220211

Address after: 100028 Chaoyang District, Beijing Hui Xin Street 6, Twelfth level.

Patentee after: SINOPEC OILFIELD SERVICE Corp.

Patentee after: Sinopec Jingwei Co.,Ltd.

Patentee after: SINOPEC NORTH CHINA PETROLEUM ENGINEERING Co.,Ltd.

Patentee after: North China Measurement and control company of Sinopec Jingwei Co.,Ltd.

Patentee after: North China logging branch of Sinopec Jingwei Co.,Ltd.

Address before: 100029 Chaoyang District, Beijing Hui Xin Street six, Twelfth level.

Patentee before: SINOPEC OILFIELD SERVICE Corp.