CN102944570A - 一种基于x射线元素荧光的岩石样品岩性定名的分析方法 - Google Patents
一种基于x射线元素荧光的岩石样品岩性定名的分析方法 Download PDFInfo
- Publication number
- CN102944570A CN102944570A CN2012105038237A CN201210503823A CN102944570A CN 102944570 A CN102944570 A CN 102944570A CN 2012105038237 A CN2012105038237 A CN 2012105038237A CN 201210503823 A CN201210503823 A CN 201210503823A CN 102944570 A CN102944570 A CN 102944570A
- Authority
- CN
- China
- Prior art keywords
- lithology
- sample
- unknown
- template
- rock sample
- 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.)
- Granted
Links
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种基于X射线元素荧光的岩石样品岩性定名的分析方法。本方法利用X射线荧光光谱仪对未知岩石样品岩性进行分析定名。首先确定各种岩性在分析时使用的元素种类,然后将此地区测得的已知岩性样品的元素含量数据存入模板数据库。软件根据使用的元素种类及模板数据库自动识别元素正态分布的期望值与标准差。根据样品元素含量及期望值与标准差,计算出该样品这种岩性的概率值。最终取概率值最高的岩性名为此岩石样品的岩性名。本方法在对油田现场录井复杂地层的高维元素分析中,通过对已知岩性的统计、分类、分析后的模板建立,提高了岩性定名的准确性,实施本方法不需要对模板建立人员的苛刻要求,便于在油田现场录井作业中推广和应用。
Description
技术领域
本发明涉及油田现场对岩石样品进行岩性定名方法,特别涉及一种基于X射线元素荧光的岩石样品岩性定名的分析方法。
背景技术
油田现场对岩石样品进行岩性定名需要进行岩石样品元素含量的检测与记录,并通过相应的解释分析方法对岩性定名进行判断。目前,各油田录井公司采用的元素分析方法包括曲线交汇法,图版分析法等。这些方法,使用元素种类少,维度利用率低,不适合高维的元素分析;而且模板建立人员需要一定的理论基础及实践经验,模板难以准确建立。因此,目前的分析方法对复杂地层的高维元素分析并不适用。
发明内容
鉴于以上现有技术存在的问题,本发明提供一种基于X射线元素荧光的岩石样品岩性定名的分析方法。本方法利用X射线荧光光谱仪和基于统计规律对未知岩石样品岩性进行分析定名。在同一地区,同种岩性的岩石各元素含量具有正态分布的规律。本方法首先确定各种岩性在分析时使用的元素种类,然后将在此地区测得的已知岩性岩石样品的元素含量数据存入模板数据库。软件根据使用的元素种类及模板数据库自动识别元素正态分布的期望值与标准差。软件分析时,根据岩石样品元素含量及模板期望值与标准差,计算出该岩石样品这种岩性的概率值。最终,取概率值最高的岩性名为此岩石样品的岩性名。
本发明采取的技术方案是:.一种基于X射线元素荧光的岩石样品岩性定名的分析方法,其特征在于,利用X射线荧光光谱仪对未知岩石样品岩性进行分析定名,其分析步骤如下:
一、建立模板数据库:首先人工选定各岩性使用的分析元素种类,对同一地区已知岩性的岩石样品进行数据采集并汇总存入模板数据库;
二、建立模板:根据选定的分析元素的含量通过模板数据库获得已知岩性样品各岩性各元素含量分布的期望值与标准差,并保存到模板文件;
三、检测并读取未知岩性样品元素含量数据;
四、根据已知岩石样品各岩性各元素含量分布的期望值与标准差以及未知岩性样品元素数据,使用正态分布模型计算未知岩性样品各岩性各元素的概率值;
五、对未知岩性样品各岩性各元素概率值取几何平均值;
六、比较未知岩性样品各岩性各元素概率值的几何平均值,其中最高值的岩性即为未知岩性的岩性名。
本发明所产生的有益效果是:本方法在对油田现场录井复杂地层的高维元素分析中,通过对已知岩性的统计、分类、分析后的模板建立,提高了岩性定名的准确性,实施本方法不需要对模板建立人员的苛刻要求,便于在油田现场录井作业中推广和应用。
附图说明
图1是本发明的软件流程图。
具体实施方式
以下结合实施例和附图对本发明作进一步说明:
实施例:
一、建立模板数据库:首先汇总同一地区的安山岩、煌斑岩、玄武质泥岩的各元素含量并存入模板数据库;安山岩、煌斑岩、玄武质泥岩的各元素含量数据见表1、表2和表3(由于篇幅所限,本表仅列出8种元素及5条数据)。
表1安山岩数据
表2煌斑岩数据
表3玄武质泥岩数据
二、建立模板:通过观察,优选出Al、Si、K、Ca、Ti、Fe这6种主要元素作为分析元素。根据选定的分析元素通过模板数据库获得安山岩、煌斑岩、玄武质泥岩各元素的期望值与标准差,见表4、表5、表6,并保存到模板文件。
表4安山岩模板
岩性 | Al期望 | Al标准差 | Si期望 | Si标准差 | K期望 | K标准差 |
安山岩 | 15.15 | 0.55 | 55.00 | 1.67 | 7.05 | 0.58 |
岩性 | Ca期望 | Ca标准差 | Ti期望 | Ti标准差 | Fe期望 | Fe标准差 |
安山岩 | 4.00 | 1.00 | 1.25 | 0.15 | 12.80 | 1.13 |
表5煌斑岩模板
岩性 | Al期望 | Al标准差 | Si期望 | Si标准差 | K期望 | K标准差 |
煌斑岩 | 7.10 | 1.20 | 36.00 | 1.33 | 14.90 | 0.63 |
岩性 | Ca期望 | Ca标准差 | Ti期望 | Ti标准差 | Fe期望 | Fe标准差 |
煌斑岩 | 7.00 | 0.83 | 3.40 | 0.27 | 18.15 | 0.58 |
表6玄武质泥岩模板
岩性 | Al期望 | Al标准差 | Si期望 | Si标准差 | K期望 | K标准差 |
玄武质泥岩 | 2.10 | 0.70 | 26.50 | 1.83 | 5.30 | 0.77 |
岩性 | Ca期望 | Ca标准差 | Ti期望 | Ti标准差 | Fe期望 | Fe标准差 |
玄武质泥岩 | 23.50 | 3.83 | 5.55 | 0.58 | 29.35 | 2.68 |
三、检测并读取未知岩性样品元素含量数据:检测未知岩性样品元素含量数据见表7。
表7未知岩性样品数据
四、根据已知岩石样品各岩性各元素含量分布的期望值与标准差以及未知岩性样品元素含量数据,使用正态分布模型计算未知岩性样品各岩性各元素的概率值:
式中:ƒμ , σ(y)为未知岩性样品此岩性此元素的概率值;
y为未知岩性样品此元素的含量值;
μ为未知岩性此元素的期望值;
σ为未知岩性此元素的标准差。
经过计算得出:
安山岩Al=1.64*10-57
安山岩Si=2.65*10-26
安山岩K=8.63*10-24
安山岩Ca=1.67*10-6
安山岩Ti=5.94*10-82
安山岩Fe=2.29*10-6
煌斑岩Al=0.78
煌斑岩Si=0.79
煌斑岩K=1.48*10-2
煌斑岩Ca=3.50*10-2
煌斑岩Ti=1.89*10-2
煌斑岩Fe=0.77
玄武质泥岩Al=2.22*10-8
玄武质泥岩Si=9.92*10-8
玄武质泥岩K=5.53*10-23
玄武质泥岩Ca=9.13*10-4
玄武质泥岩Ti=5.64*10-2
玄武质泥岩Fe=3.16*10-4。
五、对未知岩性样品各岩性各元素概率值取几何平均值:
安山岩总=(安山岩Al*安山岩Si*安山岩K*安山岩Ca*安山岩Ti*安山岩Fe)1/6=9.74*10-34。
煌斑岩总=(煌斑岩Al*煌斑岩Si*煌斑岩K*煌斑岩Ca*煌斑岩Ti*煌斑岩Fe)1/6=0.13。
玄武质泥岩总=(玄武质泥岩Al*玄武质泥岩Si*玄武质泥岩K*玄武质泥岩Ca*玄武质泥岩Ti*玄武质泥岩Fe)1/6=3.54*10-8。
六、岩性定名:比较未知岩性样品各岩性各元素概率值的几何平均值,其中最高值的岩性即为未知岩性的岩性名。
因为:煌斑岩总>玄武质泥岩总>安山岩总;
所以:未知样品的岩性定名为煌斑岩。
参照图1,本发明软件流程包括以下步骤:
1)数据分析:根据岩石样品的元素含量进行岩性分析,并保存入库。
2)数据回放:可对已经保存入库的元素数据进行查看和修改,并进行保存。
3)显示设置:设置软件界面显示数据表格及饼图显示的元素种类及次序。
4)生成报告:根据数据库生成岩石分析解释报告。
5)模板建立:确定各岩性的使用元素,并导入模板数据库。对生成的模板期望值与标准差进行修改及保存。
6)曲线显示:确定绘制曲线的元素种类,含量范围和颜色。读取分析数据库,绘制元素曲线,进行查看。
7)生成长图:确定长图的通道数,每个通道的元素数,元素种类,含量范围和颜色,图头,图底等,进行长图绘制。
Claims (1)
1.一种基于X射线元素荧光的岩石样品岩性定名的分析方法,其特征在于,利用X射线荧光光谱仪对未知岩石样品岩性进行分析定名,其分析步骤如下:
一、建立模板数据库:首先人工选定各岩性使用的分析元素种类,对同一地区已知岩性的岩石样品进行数据采集并汇总存入模板数据库;
二、建立模板:根据选定的分析元素的含量通过模板数据库获得已知岩性样品各岩性各元素含量分布的期望值与标准差,并保存到模板文件;
三、检测并读取未知岩性样品元素含量数据;
四、根据已知岩石样品各岩性各元素含量分布的期望值与标准差以及未知岩性样品元素数据,使用正态分布模型计算未知岩性样品各岩性各元素的概率值;
五、对未知岩性样品各岩性各元素概率值取几何平均值;
六、比较未知岩性样品各岩性各元素概率值的几何平均值,其中最高值的岩性即为未知岩性的岩性名。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210503823.7A CN102944570B (zh) | 2012-12-02 | 2012-12-02 | 一种基于x射线元素荧光的岩石样品岩性定名的分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210503823.7A CN102944570B (zh) | 2012-12-02 | 2012-12-02 | 一种基于x射线元素荧光的岩石样品岩性定名的分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102944570A true CN102944570A (zh) | 2013-02-27 |
CN102944570B CN102944570B (zh) | 2014-11-26 |
Family
ID=47727532
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210503823.7A Active CN102944570B (zh) | 2012-12-02 | 2012-12-02 | 一种基于x射线元素荧光的岩石样品岩性定名的分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102944570B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104612675A (zh) * | 2015-02-02 | 2015-05-13 | 中石化西南石油工程有限公司地质录井分公司 | 一种碳酸盐岩地层随钻岩性快速识别方法 |
CN105044140A (zh) * | 2015-06-30 | 2015-11-11 | 北京奥意尔工程技术有限公司 | 基于x射线荧光元素录井技术的盐底卡取方法 |
CN107367520A (zh) * | 2017-06-27 | 2017-11-21 | 中国石油天然气股份有限公司 | 一种基于xrf识别细粒沉积岩岩性的方法 |
CN107489417A (zh) * | 2017-07-25 | 2017-12-19 | 中国海洋石油总公司 | 基于xrf录井的岩矿组成反演法砂泥岩岩性识别方法 |
CN108240999A (zh) * | 2016-12-26 | 2018-07-03 | 中国石油天然气股份有限公司 | 一种碎屑岩岩性的识别方法 |
CN109490266A (zh) * | 2018-11-15 | 2019-03-19 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种岩样无损取样方法 |
CN110286141A (zh) * | 2019-07-15 | 2019-09-27 | 中国石油集团渤海钻探工程有限公司 | 一种基于逻辑回归的自动岩性定名方法 |
CN110689830A (zh) * | 2019-10-10 | 2020-01-14 | 辽宁工程技术大学 | 一种地质岩样标本采集展示系统 |
CN111860983A (zh) * | 2020-07-06 | 2020-10-30 | 中国煤炭地质总局勘查研究总院 | 煤中伴生金属资源量估算方法 |
US20220291154A1 (en) * | 2021-03-10 | 2022-09-15 | Eiffage Gc Infra Lineaires | Method for identifying the lithological layer of an excavated material for the purpose of its valorisation |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4510573A (en) * | 1981-05-06 | 1985-04-09 | Unc Nuclear Industries, Inc. | Method for X-ray fluorescence spectroscopy |
CN1441244A (zh) * | 2003-03-01 | 2003-09-10 | 中国电子科技集团公司第二十二研究所 | 岩屑描述仪及岩屑岩性和含油丰度的定量分析方法 |
CN101082277A (zh) * | 2007-07-05 | 2007-12-05 | 北京奥能瑞科石油技术有限公司重庆分公司 | 石油钻井地质x射线荧光岩屑录井方法 |
CN101354362A (zh) * | 2008-07-21 | 2009-01-28 | 中国石化集团华北石油局 | 石油钻井中x射线荧光泥质含量的分析方法 |
US20090132372A1 (en) * | 2007-11-20 | 2009-05-21 | Nhn Corporation | Method for calculating predicted charge amount of advertisement for each keyword and system for executing the method |
WO2010129280A2 (en) * | 2009-04-28 | 2010-11-11 | Baker Hughes Incorporated | Petrophysics-guided processing of lwd acoustic data |
CN101936929A (zh) * | 2010-07-27 | 2011-01-05 | 中国石化集团华北石油局 | X射线荧光元素分析装置及其录井仪 |
WO2012149003A2 (en) * | 2011-04-26 | 2012-11-01 | Baker Hughes Incorporated | Apparatus and method for estimating formation lithology using x-ray flourescence |
-
2012
- 2012-12-02 CN CN201210503823.7A patent/CN102944570B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4510573A (en) * | 1981-05-06 | 1985-04-09 | Unc Nuclear Industries, Inc. | Method for X-ray fluorescence spectroscopy |
CN1441244A (zh) * | 2003-03-01 | 2003-09-10 | 中国电子科技集团公司第二十二研究所 | 岩屑描述仪及岩屑岩性和含油丰度的定量分析方法 |
CN101082277A (zh) * | 2007-07-05 | 2007-12-05 | 北京奥能瑞科石油技术有限公司重庆分公司 | 石油钻井地质x射线荧光岩屑录井方法 |
US20090132372A1 (en) * | 2007-11-20 | 2009-05-21 | Nhn Corporation | Method for calculating predicted charge amount of advertisement for each keyword and system for executing the method |
CN101354362A (zh) * | 2008-07-21 | 2009-01-28 | 中国石化集团华北石油局 | 石油钻井中x射线荧光泥质含量的分析方法 |
WO2010129280A2 (en) * | 2009-04-28 | 2010-11-11 | Baker Hughes Incorporated | Petrophysics-guided processing of lwd acoustic data |
CN101936929A (zh) * | 2010-07-27 | 2011-01-05 | 中国石化集团华北石油局 | X射线荧光元素分析装置及其录井仪 |
WO2012149003A2 (en) * | 2011-04-26 | 2012-11-01 | Baker Hughes Incorporated | Apparatus and method for estimating formation lithology using x-ray flourescence |
Non-Patent Citations (2)
Title |
---|
朱根庆等: "X射线荧光录井资料基本解释方法", 《录井工程》 * |
李一超等: "X射线荧光分析在岩屑录井中的应用", 《岩石矿物学杂志》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104612675B (zh) * | 2015-02-02 | 2017-07-14 | 中石化西南石油工程有限公司地质录井分公司 | 一种碳酸盐岩地层随钻岩性快速识别方法 |
CN104612675A (zh) * | 2015-02-02 | 2015-05-13 | 中石化西南石油工程有限公司地质录井分公司 | 一种碳酸盐岩地层随钻岩性快速识别方法 |
CN105044140B (zh) * | 2015-06-30 | 2018-03-16 | 北京奥意尔工程技术有限公司 | 基于x射线荧光元素录井技术的盐底卡取方法 |
CN105044140A (zh) * | 2015-06-30 | 2015-11-11 | 北京奥意尔工程技术有限公司 | 基于x射线荧光元素录井技术的盐底卡取方法 |
CN108240999A (zh) * | 2016-12-26 | 2018-07-03 | 中国石油天然气股份有限公司 | 一种碎屑岩岩性的识别方法 |
CN107367520A (zh) * | 2017-06-27 | 2017-11-21 | 中国石油天然气股份有限公司 | 一种基于xrf识别细粒沉积岩岩性的方法 |
CN107489417A (zh) * | 2017-07-25 | 2017-12-19 | 中国海洋石油总公司 | 基于xrf录井的岩矿组成反演法砂泥岩岩性识别方法 |
CN109490266A (zh) * | 2018-11-15 | 2019-03-19 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种岩样无损取样方法 |
CN109490266B (zh) * | 2018-11-15 | 2021-08-27 | 中国石油化工股份有限公司 | 一种岩样无损取样方法 |
CN110286141A (zh) * | 2019-07-15 | 2019-09-27 | 中国石油集团渤海钻探工程有限公司 | 一种基于逻辑回归的自动岩性定名方法 |
CN110689830A (zh) * | 2019-10-10 | 2020-01-14 | 辽宁工程技术大学 | 一种地质岩样标本采集展示系统 |
CN111860983A (zh) * | 2020-07-06 | 2020-10-30 | 中国煤炭地质总局勘查研究总院 | 煤中伴生金属资源量估算方法 |
US20220291154A1 (en) * | 2021-03-10 | 2022-09-15 | Eiffage Gc Infra Lineaires | Method for identifying the lithological layer of an excavated material for the purpose of its valorisation |
US11879858B2 (en) * | 2021-03-10 | 2024-01-23 | Eiffage Gc Infra Lineaires | Method for identifying the lithological layer of an excavated material for the purpose of its valorisation |
Also Published As
Publication number | Publication date |
---|---|
CN102944570B (zh) | 2014-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102944570B (zh) | 一种基于x射线元素荧光的岩石样品岩性定名的分析方法 | |
CN102759748B (zh) | 基于目标分析的交汇分析地震储层预测方法 | |
CN103913781B (zh) | 一种预测铀矿体埋深的组合探测方法 | |
CN103775075B (zh) | 一种全井段岩性识别方法 | |
CN114813903B (zh) | 一种基于石榴石微区化学成分判别矿种的方法 | |
CN106295113A (zh) | 一种复杂油气储层渗透率的定量评价方法 | |
CN103678778B (zh) | 一种放射性物化探信息综合方法 | |
CN104806232A (zh) | 一种确定碳酸盐岩储层孔隙度下限的方法 | |
CN106682377A (zh) | 一种钻孔岩芯蚀变矿物的高光谱提取与定量估算方法 | |
CN115267928B (zh) | 一种用于随钻元素测井的能谱智能处理方法 | |
CN107505344A (zh) | 利用“最小二乘积”法的岩性解释方法 | |
CN103336305A (zh) | 一种基于灰色理论划分致密砂岩储层岩石物理相的方法 | |
CN105114067A (zh) | 一种岩性测井相的方法 | |
WO2019205725A1 (zh) | 一种基于电成像测井分形维数的火山岩岩性识别方法 | |
CN104143115A (zh) | 地质雷达技术实现土壤含水性分类识别的技术方法 | |
CN106777707B (zh) | 一种利用改进的蜘蛛网图进行测井岩性定量识别的方法 | |
CN114609675A (zh) | 基于高频旋回对碳酸盐岩地层沉积微地貌的定量恢复方法 | |
CN104459790B (zh) | 含油气性盆地有效储层的分析方法和装置 | |
CN110208874B (zh) | 一种致密砂岩油藏有效储层识别方法 | |
CN109283577B (zh) | 一种地震层位标定方法 | |
CN106468171A (zh) | 多信息耦合的递进式砂砾岩岩相判别方法 | |
CN105590018A (zh) | 砂泥岩薄互层油藏的油水层判别方法 | |
CN106872669A (zh) | 一种花岗岩潜山裂缝储层电性分类方法 | |
CN107991705B (zh) | 基于二维统计特征的测井曲线校正方法和装置 | |
CN116027454A (zh) | 一种砂岩厚度等值线图的成图方法及智能绘制装置 |
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 | ||
CP02 | Change in the address of a patent holder |
Address after: 300450 1st floor, no.6, Haitai Huake 10th Road, Huayuan Industrial Zone, Binhai New Area, Tianjin Patentee after: TIANJIN LAND & OCEAN PETROLEUM EQUIPMENT SYSTEMS ENGINEERING Co.,Ltd. Address before: 300192 No. 16 West research road, Tianjin, Nankai District Patentee before: TIANJIN LAND & OCEAN PETROLEUM EQUIPMENT SYSTEMS ENGINEERING Co.,Ltd. |
|
CP02 | Change in the address of a patent holder |