CN103345002A - 一种测井资料获取砂岩卤水层钾离子含量的方法 - Google Patents

一种测井资料获取砂岩卤水层钾离子含量的方法 Download PDF

Info

Publication number
CN103345002A
CN103345002A CN2013102867435A CN201310286743A CN103345002A CN 103345002 A CN103345002 A CN 103345002A CN 2013102867435 A CN2013102867435 A CN 2013102867435A CN 201310286743 A CN201310286743 A CN 201310286743A CN 103345002 A CN103345002 A CN 103345002A
Authority
CN
China
Prior art keywords
sandstone
bittern
layer
content
potassium
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
Application number
CN2013102867435A
Other languages
English (en)
Other versions
CN103345002B (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.)
Logging Co of Sinopec Jianghan Petroleum Engineering Co Ltd
Exploration and Development Research Institute of Sinopec Jianghan Oilfield Co
Original Assignee
Logging Co of Sinopec Jianghan Petroleum Engineering Co Ltd
Exploration and Development Research Institute of Sinopec Jianghan Oilfield Co
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 Logging Co of Sinopec Jianghan Petroleum Engineering Co Ltd, Exploration and Development Research Institute of Sinopec Jianghan Oilfield Co filed Critical Logging Co of Sinopec Jianghan Petroleum Engineering Co Ltd
Priority to CN201310286743.5A priority Critical patent/CN103345002B/zh
Publication of CN103345002A publication Critical patent/CN103345002A/zh
Application granted granted Critical
Publication of CN103345002B publication Critical patent/CN103345002B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及到地下卤水勘探开发中卤水品位的评价,具体是利用自然伽马测井获取砂岩卤水层钾离子含量的方法,包括以下步骤:1)在待解释地层进行地球物理测井,包括岩性类测井、物性类测井、自然伽马能谱测井;2)划分砂岩层和纯泥岩,求取砂岩层的泥质含量Vsh;3)计算砂岩层孔隙度φ,划分出砂岩卤水层和干砂层;4)利用自然伽马能谱测井资料取得砂岩卤水层总钾含量K、砂岩卤水层附近纯泥岩的钾含量Ksh和砂岩卤水层附近干砂层的钾含量Kma;5)求取砂岩卤水层的钾离子质量含量。该发明是一种全新的求取深井卤水中的钾离子含量的方法,具有较高的精度,满足钾离子含量的测井解释需求。

Description

一种测井资料获取砂岩卤水层钾离子含量的方法
技术领域
本发明属于地质勘探领域,涉及到地下卤水勘探开发中卤水品位的评价,特别是利用自然伽马能谱测井资料获取砂岩卤水层钾离子含量方法。
背景技术
在地下卤水勘探开发中,通过钻井钻开地层后,利用测井获得的资料解释卤水层,评价卤水的品位(矿化度和微量元素含量的高低),是下步开展试水的重要工作。目前利用测井资料划分砂岩类型的水层已经十分成熟,包括利用自然电位(SP)、自然伽马(GR)、井径(CAL)、光电吸收截面指数(PE)等岩性测井资料划分砂岩层,然后利用声波时差(AC)、密度(DEN)、中子(CNL)等物性测井资料划分水层和干层。
卤水的利用基本上是以利用其中微量稀散元素(B、Li、K、Br等)为主,因此对其中微量元素的计算,是对深层砂岩卤水品位评价的重要内容。目前未见地球物理测井资料求取卤水层中的微量元素包括钾离子含量的计算方法。
发明内容
本发明为填补在卤水勘探开发中利用测井资料获取砂岩卤水层中钾离子含量方法的空白,提供了一种利用测井资料获取砂岩卤水层钾离子含量的方法,以获得更多在在开展试水工作前对深井砂岩卤水层品位的评价数据。
本发明技术方案如下:
一种测井资料获取砂岩卤水层钾离子含量的方法,包括以下步骤:
1)在待解释地层进行地球物理测井,包括岩性类测井、物性类测井、自然伽马能谱测井;
2)利用岩性类测井资料,划分砂岩层和纯泥岩,求取砂岩层的泥质含量Vsh
3)利用物性类测井资料,计算砂岩层孔隙度φ,划分出砂岩卤水层和干砂层;
建立体积模型:砂岩卤水层体积为泥质、砂岩骨架和砂岩孔隙中卤水体积之和,且三者体积为一个单位,以此得出砂岩体积百分含量1-Vsh-φ;
4)利用自然伽马能谱测井取得砂岩卤水层总钾含量K、砂岩卤水层附近纯泥岩的钾含量Ksh和砂岩卤水层附近干砂层的钾含量Kma
在体积模型下建立方程:K=(1-φ-Vsh)Kma+φKw+VshKsh
5)利用Vsh、φ、K、Ksh、Kma,求取砂岩卤水层的钾离子质量含量Kk
Kk=[K-(1-φ-Vsh)Kma-VshKsh]
进一步求取砂岩卤水层体积含量,还包括以下步骤:
通过试验建立地区卤水矿化度与卤水密度的一般函数关系;
ρ=0.0006C+1.0006,其中ρ为卤水的密度,g/cm3,C为卤水矿化度,g/L;
利用测井资料求取卤水矿化度C,进而计算卤水密度ρ;
再利用步骤5)所得钾离子质量含量求取钾离子体积百分含量Kv:
计算方程为:Kv=Kk.(0.0006C+1.0006)×106,Kv单位mg/L,C单位g/L。
上述步骤4)中选择砂岩卤水层附近纯泥岩的钾离子含量代替体积模型中泥质的钾离子含量Ksh,选择砂岩卤水层附近干砂层的钾离子含量代替砂岩骨架的钾离子含量Kma
地球物理测井应该包括以下几种类型。岩性类测井包括自然伽马(GR)、井径(CAL)、光电吸收截面指数(PE)等测井;物性类测井包括声波(AC)、密度(DEN)、中子孔隙度(CNL)等测井;求取钾离子含量重要测井自然伽马能谱测井;以上测井数据经过测量后,进行测井资料处理。
自然伽马能谱测井可以直接求取地层中铀、钍、钾的含量。自然伽马能谱测井的理论基础地层中存在的放射性核素,主要是天然放射性核素,这些核素又分放射系和非放射系的天然放射性核素。放射系为钍系、铀系和锕铀系,但锕铀系的头一个核素U235在自然界中的丰度很低,其放射性贡献甚微,不予考虑。非放射系的天然放射性核素主要是Rb87和K40,但是Rb87无伽马辐射。所以在研究地层中的自然伽马能谱主要是U238、Th232放射系和K40放射的伽马射线能谱。
铀系和钍系在放射性平衡状态下系内核素的原子核数的比例关系是确定的,因此不同能量伽马的相对强度也是确定的,因此我们可以分别在这两个系中选出某种核素的特征核素伽马射线的能量来分别识别铀和钍。这种被选定的某种核素称为特征核素,它发射的伽马射的能量称为特征能量,在自然伽马能谱测井中,通常选用铀系中的Bi214发射的1.76MeV的伽马射线来识别铀,选用钍系中的Tl208发射的2.62MeV的伽马射线来识别钍,用1.46MeV的伽马射线来识别钾。当我们把伽马射线按我们所选定的特征能量分别计数,那么这就叫测谱。因而将测得的自然伽马能谱转换成地层的铀、钍、钾的含量,并计录在磁带上或以连续测井曲线的形式输出,这就是自然伽马能谱测井。
在此数据基础上,将地层中的砂岩卤水层抽象成一个体积模型,即整个体积是由泥质、砂岩骨架和砂岩孔隙中卤水组成,三者体积为一个单位,其中泥质的体积Vsh在步骤2)中可以求出,卤水的体积φ在步骤3)中可以求出,因此砂岩的体积为1-Vsh-φ(见附图1)。
砂岩卤水层中所测得的钾含量是砂岩孔隙中的卤水层、泥质含量、砂岩骨架所含钾含量的综合响应。因此其测井相应方程为:
K=(1-φ-Vsh)Kma+φKw+VshKsh………………………………………………(1)
式中:K—自然伽马能谱测井中钾的测井值,无量纲;Kma—砂岩骨架钾含量,无量纲;Vsh—地层的泥质含量,无量纲;φ—砂岩孔隙度,无量纲;Kw—卤水中单位体积钾离子的含量,无量纲;Ksh—泥质的钾含量(由于各地区砂岩骨架的钾含量不同,不是一个固定值,选取解释的选择砂岩卤水层附近纯泥岩的钾离子含量代替体积模型中泥质的钾离子含量Ksh,选择砂岩卤水层附近干砂层的钾离子含量代替砂岩骨架的钾离子含量Kma)。
设Kk=φKw……………………………………………………………………(2)
式中:KK—孔隙中卤水钾离子含量,无量纲。
联合(1)和(2)式可以求出:
Kk=[K-(1-φ-Vsh)Kma-VshKsh]…………………………………………………(3)
式中,Vsh、φ、K、Kma、Ksh均已知,以此求出卤水层的钾离子质量百分含量(无量纲)。
特别的,在得到质量含量的基础上,还能进一步求取体积含量。
设卤水的质量为m,密度为ρ,则钾离子的体积含量:
Kv=(Kk.m)/(m/ρ)=Kk.ρ……………………………………(4)
卤水的密度与矿化度密切相关,通过试验制作图板建立卤水矿化度和卤水密度的一般函数关系函数后,利用卤水矿化度,可以求取卤水的密度:
ρ=0.0006C+1.0006…………………………………………………………(5)
式中:ρ—卤水的密度,g/cm3;C—卤水矿化度,g/L。
联合(4)和(5)式后求得:
Kv=Kk.(0.0006C+1.0006)…………………………………………………(6)
式中:Kv—卤水钾离子体积含量,g/cm3;Kk—卤水钾离子质量百分含量,无量纲;C—卤水矿化度,g/L。
按照钾离子体积含量的单位习惯是mg/L,将(6)式变换单位后得到:
Kv=Kk.(0.0006C+1.0006)×106…………………………………………(7)
式中:Kv—卤水钾离子体积含量,mg/L;Kk—卤水钾离子质量百分含量,无量纲;C—卤水矿化度,g/L。
在制作卤水矿化度和卤水的密度一般函数关系函数基础上,利用技术手段测得卤水矿化度C后,再用自然伽马能谱求得的卤水质量含量Kk,两者代入(7)中求取卤水的体积含量。
本发明的有益效果为:
(1)填补了在卤水勘探开发中利用测井资料评价砂岩卤水中钾离子含量方法的空白。
(2)该方法为单井快速评价砂岩卤水层中钾离子提供了较好的方法,为下步开展卤水井优选高品位的富钾卤水层进行试水提供了条件。
附图说明
图1:砂岩卤水层体积模型图;
图2:资深1井自然伽马能谱测井计算卤水层钾离子含量图;
图3:江汉盆地卤水矿化度与卤水密度关系图。
具体实施方式
为了更好地理解本发明,下面结合实施例进一步阐明本发明的内容,但本发明的内容不仅仅局限于下面的实施例。
资深1井3909.4-3911.8m和3916.0-3922.4m为一砂岩卤水层(见附图2),两层声波时差为213~227μs/m,密度为2.52~2.45g/cm3。选取密度计算孔隙度,通过计算泥质含量Vsh平均为0.15,孔隙度平均为9.6%和9.8%(见下表),通过电阻率计算矿化度为253.22g/L;然后利用本地区卤水矿化度和密度的函数关系式(见附图3)计算卤水层密度,求得卤水层密度为1.153g/cm3;砂岩的钾骨架值Kma=1.05%,泥岩的钾测量值Ksh=3.2%,运用以上公式计算得到两层的钾离子体积含量分别为6.75g/L和6.83g/L,两者平均含量为6.79g/L。试水资料实验室测定的钾离子含量为7.47g/L。通过对比试水资料,钾离子含量的绝对误差为0.68g/L,相对误差为9.13%,具有一定的精度,满足钾离子含量的测井解释需求。
Figure BDA00003486901500041

Claims (3)

1.一种测井资料获取砂岩卤水层钾离子含量的方法,包括以下步骤:
1)在待解释地层进行地球物理测井,包括岩性类测井、物性类测井、自然伽马能谱测井;
2)利用岩性类测井资料,划分砂岩层和纯泥岩,求取砂岩层的泥质含量Vsh
3)利用物性类测井资料,计算砂岩层孔隙度φ,划分出砂岩卤水层和干砂层;
建立体积模型:砂岩卤水层体积为泥质、砂岩骨架和砂岩孔隙中卤水体积之和,且三者体积为一个单位,以此得出砂岩体积1-Vsh-φ;
4)利用自然伽马测井取得砂岩卤水层总钾含量K、砂岩卤水层附近纯泥岩的钾含量Ksh和砂岩卤水层附近干砂层的钾含量Kma
在体积模型下建立方程:K=(1-φ-Vsh)Kma+φKw+VshKsh
5)利用Vsh、φ、K、Ksh、Kma,求取砂岩卤水层的钾离子质量含量Kk
方程为:Kk=[K-(1-φ-Vsh)Kma-VshKsh]。
2.如权利要求1所述一种测井资料获取砂岩卤水层钾离子含量的方法,其特征在于还包括以下步骤:
通过试验建立地区卤水矿化度与卤水密度的一般函数关系:
ρ=0.0006C+1.0006,其中ρ为卤水的密度,g/cm3,C为卤水矿化度,g/L;
利用测井资料求取卤水矿化度C,进而计算卤水密度ρ;
再利用步骤5)所得钾离子质量含量Kk计算得到钾离子体积含量Kv:;
计算方程为:Kv=Kk.(0.0006C+1.0006)×106,Kv单位mg/L,C单位g/L。
3.如权利要求1或2所述一种测井资料获取砂岩卤水层钾离子含量的方法,其特征在于步骤4)中选择砂岩卤水层附近纯泥岩的钾离子含量代替体积模型中泥质的钾离子含量Ksh,选择砂岩卤水层附近干砂层的钾离子含量代替砂岩骨架的钾离子含量Kma
CN201310286743.5A 2013-07-09 2013-07-09 一种测井资料获取砂岩卤水层钾离子含量的方法 Active CN103345002B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310286743.5A CN103345002B (zh) 2013-07-09 2013-07-09 一种测井资料获取砂岩卤水层钾离子含量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310286743.5A CN103345002B (zh) 2013-07-09 2013-07-09 一种测井资料获取砂岩卤水层钾离子含量的方法

Publications (2)

Publication Number Publication Date
CN103345002A true CN103345002A (zh) 2013-10-09
CN103345002B CN103345002B (zh) 2016-06-01

Family

ID=49279815

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310286743.5A Active CN103345002B (zh) 2013-07-09 2013-07-09 一种测井资料获取砂岩卤水层钾离子含量的方法

Country Status (1)

Country Link
CN (1) CN103345002B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103852475A (zh) * 2014-03-07 2014-06-11 中国科学院青海盐湖研究所 一种基于伽马射线的多道测钾仪
CN104612673A (zh) * 2015-01-15 2015-05-13 中国地质科学院矿产资源研究所 一种寻找深层富钾卤水矿床的矿物-化学方法及其应用
CN107045154A (zh) * 2017-02-08 2017-08-15 中国海洋石油总公司 一种水平井环境中的识别地层产状的方法和装置
CN111045114A (zh) * 2019-12-31 2020-04-21 核工业北京地质研究院 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法
CN112878998A (zh) * 2021-01-28 2021-06-01 成都理工大学 地下卤水型钾矿和锂矿的储卤层预测和资源量评价方法
CN116931062A (zh) * 2023-07-21 2023-10-24 成都理工大学 基于敏感参数分析的优质富锂钾卤水储层识别方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4436997A (en) * 1981-08-21 1984-03-13 Mobil Oil Corporation Method for the determination of clay and mica concentrations in subsurface sandstone formations through radioactive logging
US5065016A (en) * 1990-06-11 1991-11-12 Texaco Inc. Radioactive well logging to determine vertical brine flow
US20100193676A1 (en) * 2007-08-14 2010-08-05 Halliburton Energy Services, Inc. Determining formation characteristics
WO2011001269A2 (en) * 2009-07-02 2011-01-06 Sclumberger Technology B.V. System and method for drilling using drilling fluids
CN102012526A (zh) * 2010-09-09 2011-04-13 四川德阳西德电器有限公司 采用电阻率资料进行储层流体类型判别的方法
CN102518428A (zh) * 2011-11-11 2012-06-27 中国石油天然气股份有限公司 油水层识别方法及装置
CN102576091A (zh) * 2009-11-11 2012-07-11 雪佛龙美国公司 分析和变换地球物理和岩石物理数据的系统和方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4436997A (en) * 1981-08-21 1984-03-13 Mobil Oil Corporation Method for the determination of clay and mica concentrations in subsurface sandstone formations through radioactive logging
US5065016A (en) * 1990-06-11 1991-11-12 Texaco Inc. Radioactive well logging to determine vertical brine flow
US20100193676A1 (en) * 2007-08-14 2010-08-05 Halliburton Energy Services, Inc. Determining formation characteristics
WO2011001269A2 (en) * 2009-07-02 2011-01-06 Sclumberger Technology B.V. System and method for drilling using drilling fluids
CN102576091A (zh) * 2009-11-11 2012-07-11 雪佛龙美国公司 分析和变换地球物理和岩石物理数据的系统和方法
CN102012526A (zh) * 2010-09-09 2011-04-13 四川德阳西德电器有限公司 采用电阻率资料进行储层流体类型判别的方法
CN102518428A (zh) * 2011-11-11 2012-06-27 中国石油天然气股份有限公司 油水层识别方法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
李英华: "油田水地化指标研究的新认识", 《中国海上油气(地质)》, vol. 12, no. 01, 15 February 1998 (1998-02-15) *
胡素华等: "成都盐盆富钾卤水储层测井解释方法初探", 《中国井矿盐》, vol. 39, no. 06, 25 November 2008 (2008-11-25) *
陈建立: "地球物理测井在安棚天然碱矿勘查中的应用", 《地质找矿论丛》, vol. 25, no. 3, 15 September 2010 (2010-09-15), pages 252 - 259 *
陈科贵等: "运用测井资料进行K2O视含量计算方法研究", 《石油天然气学报》, vol. 34, no. 3, 15 March 2012 (2012-03-15), pages 96 - 100 *
黄华等: "潜江凹陷潜江组砂岩卤水矿床测井识别方法研究", 《化工矿产地质》, vol. 35, no. 2, 15 June 2013 (2013-06-15), pages 65 - 71 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103852475A (zh) * 2014-03-07 2014-06-11 中国科学院青海盐湖研究所 一种基于伽马射线的多道测钾仪
CN103852475B (zh) * 2014-03-07 2017-01-11 中国科学院青海盐湖研究所 一种基于伽马射线的多道测钾仪
CN104612673A (zh) * 2015-01-15 2015-05-13 中国地质科学院矿产资源研究所 一种寻找深层富钾卤水矿床的矿物-化学方法及其应用
CN107045154A (zh) * 2017-02-08 2017-08-15 中国海洋石油总公司 一种水平井环境中的识别地层产状的方法和装置
CN111045114A (zh) * 2019-12-31 2020-04-21 核工业北京地质研究院 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法
CN111045114B (zh) * 2019-12-31 2021-11-30 核工业北京地质研究院 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法
CN112878998A (zh) * 2021-01-28 2021-06-01 成都理工大学 地下卤水型钾矿和锂矿的储卤层预测和资源量评价方法
CN116931062A (zh) * 2023-07-21 2023-10-24 成都理工大学 基于敏感参数分析的优质富锂钾卤水储层识别方法
CN116931062B (zh) * 2023-07-21 2023-12-22 成都理工大学 基于敏感参数分析的优质富锂钾卤水储层识别方法

Also Published As

Publication number Publication date
CN103345002B (zh) 2016-06-01

Similar Documents

Publication Publication Date Title
CN104636588B (zh) 计算烃源岩中有机碳含量的方法及装置
CN101930082B (zh) 采用电阻率资料进行储层流体类型判别的方法
CN103345002A (zh) 一种测井资料获取砂岩卤水层钾离子含量的方法
CN102012526A (zh) 采用电阻率资料进行储层流体类型判别的方法
KR101415199B1 (ko) 전기비저항 및 밀도 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 슬로우니스값, 영률, 포아송비 및 취성 추정방법
CN109521487B (zh) 一种利用元素伽马能谱测井识别气层的方法
Nansheng et al. Geothermal regime and Jurassic source rock maturity of the Junggar Basin, northwest China
CN105487135B (zh) 沉积盆地型地热资源类型划分方法
CN105467465B (zh) 利用孔隙度差值计算粘土含量的方法
CA2867583C (en) Fracking method for fracking intervals of a horizontal drilling zone in a sweet spot range based on measurements of resistivity and neutron logging data in the horizontal drillingzone
KR101415198B1 (ko) 전기비저항 및 중성자 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 슬로우니스값, 영률, 포아송비 및 취성 추정방법
CN101906963A (zh) 利用c/o和地层宏观俘获截面交会技术确定饱和度方法
CN103343687B (zh) 一种测井资料获取砂岩卤水层等效氯化钠矿化度的方法
McKenna et al. Radiogenic Heat Production in Sendimentary Rocks of the Gulf of Mexico Basin, South Texas
CN102183788B (zh) 地层条件下横波测井曲线的合成方法
CN104297810A (zh) 一种地层元素测井中净非弹性散射伽马能谱的获取方法
CN103765247A (zh) 用于快速正演中子模型的地层表征
KR101415197B1 (ko) 전기비저항 및 중성자 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 파쇄구간 선정방법
Ferronsky Nuclear geophysics: applications in hydrology, hydrogeology, engineering geology, agriculture and environmental science
Manning et al. Direct observation of the depth of active groundwater circulation in an alpine watershed
CN103306663A (zh) 铀矿测井方法
Zhang et al. Monte Carlo simulation for determining gas saturation using three-detector pulsed neutron logging technology in tight gas reservoir and its application
Fu et al. A new method of gas reservoir evaluation based on neutron cross section logging
Wang et al. Exploration process and genesis mechanism of deep geothermal resources in the North Jiangsu Basin, East China: From nothing to something
Evans et al. A sourceless alternative to conventional LWD nuclear logging

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