CN1595202A - 一种氯能谱测井方法 - Google Patents
一种氯能谱测井方法 Download PDFInfo
- Publication number
- CN1595202A CN1595202A CN 03156760 CN03156760A CN1595202A CN 1595202 A CN1595202 A CN 1595202A CN 03156760 CN03156760 CN 03156760 CN 03156760 A CN03156760 A CN 03156760A CN 1595202 A CN1595202 A CN 1595202A
- Authority
- CN
- China
- Prior art keywords
- neutron
- chlorine
- formula
- gamma
- porosity
- 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
本发明涉及在石油开采的盐水地层油藏注水高含水期测定井下含水饱和度的一种氯能谱测井方法。采用热中子探测器和伽马探测器同时测量中子源子发射中子与地层作用后产生的热中子和弹性散射伽马信息,通过测量中子与中子伽马相对于饱和淡水或油地层的差异,求解出产层的视孔隙度和含水饱和度。因而可以克服地层中泥质和骨架的影响,双向差异测量方法也提高了对地层水矿化度变化测量响应的灵敏度。由于使用的仪器结构简单,稳定可靠,测井成本低。
Description
所属领域
本发明涉及在石油开采的盐水地层油藏注水高含水期测定井下含水饱和度的一种氯能谱测井方法。
背景技术
在油田注水开发后期,水驱采油使得含水率不断升高,可达到80%以上。目前,在矿化度较高地区,(地层水氯离子Cl-含量>4×104ppm),主要采用两种测井方法:一是中子寿命测井方法;二是氯能谱测井方法。前者测量值是地层中总的热中子吸收截面,在扣除泥质和骨架的收截面值后,以差值部分来对应孔隙流体的热中子吸收截面,通过孔隙流体的吸收截面值最终来确定油层含水饱和度。这种方法由于受到泥质和骨架不确定性的影响,定量解释精度不高。另外,采用的人造中子发生管结构复杂,造价高,维修困难。一只中子管价值十几万,寿命仅有一百多小时,而每次测井就要使用2小时,仪器寿命短,测井成本高。
目前采用的氯能谱测井方法是测量单一的伽马能段计数,用单一的伽马能段计数的变化反映地层中Cl-矿化度,以Cl-矿化度的来确定地层的含水饱和度。这种Cl-含量的变化所产生伽马计数的差异,对于定量解释含水饱和度,特别是高含水期测定地层含水饱和度来说,其相对灵敏度较差。
发明内容
本发明目的是提供一种能在盐水地层油藏注水开发高含水期确定含水饱和度且仪器测井成本低、测量精度高的氯能谱测井方法。
本发明采用的技术方法是:
采用热中子探测器和伽马探测器同时测量中子源子发射中子与地层作用后产生的热中子和弹性散射伽马信息,通过测量中子与中子伽马相对于饱和淡水或油地层的差异,求解出产层的视孔隙度和含水饱和度。
根据已知地层矿化度、孔隙度、岩性建立刻度模型井,在模拟条件建立视孔隙度关系函数,采用如下公式计算视孔隙度:Φ=aG3+bG2+cG+d,式中:Φ是视孔隙度;a、b、c、d是拟和公式的系数,是使用的热中子和中子伽马测量仪器结构、精度的参数;G=k1Nz+K2Nr,k1,k2是加权系数,Nz、Nr分别表示热中子和伽马氯能段测量的计数。
根据已知地层矿化度、孔隙度、岩性建立刻度模型井,在模拟条件建立含水饱和度关系函数,求含水饱和度,步骤是:
建立归一函数关系,对应地层,每一个中子测量值都有相应的氯测量值与其对应,在淡水地层中对应不同的空隙度,氯与中子的归一关系式;N’z=K1NCl+K2,式中:K1,K2是拟合系数,N’z,NCl分别是测量的热中子和伽马氯能段计数值;
将测井所得氯计数值Nr代入上式得出归一的相应中子计数值N’z;
求氯曲线与中子曲线的双向归一差异值Dt=N’z-Nz;
确定相对差Δ:
在对应地层孔隙度为Φ时,含水饱和度Sw为100%时氯与中子的归一差D,其关系D=KPΦ,式中:K是实验关系刻度系数,P是地层水矿化度;
求出相对差Δ=1-Dt/D;
为了对不同的矿化度值给出一个统一的求解含水饱和度的公式,先对前面求出的Δ值作一个矿化度校正,关系式:δ=Δ+RP+M,式中:R,M是实验关系系数,是使用的测量仪器结构、精度的参数;
含水饱和度Sw=S3δ3+S2δ2+S1δ+S,式中S3,S2,S1S是实验拟合关系系数,是使用的测量仪器结构、精度的参数。
视孔隙度Φ公式可以写成对数形式。Sw的三次多项式公式变为分段的直线函数形式。
本发明的优点是:由热中子探测器和伽马探测器进行组合测量,利用在正源距条件下,矿化度对热中子和中子伽马测量产生的相反影响,即随着矿化度的增加热中子数减少而中子俘获伽马数增加;矿化度减少时则相反。因而利用热中子与俘获伽马加权和确定空隙度就可以克服矿化度的影响。建立了用中子与伽马双差异的变化求解含水饱和度的测井方法。因为所确定的中子与伽马的差异主要是对应于孔隙流体矿化度的变化,因而可以克服地层中泥质和骨架的影响,双向差异测量方法也提高了对地层水矿化度变化测量响应的灵敏度。由于使用的仪器结构简单,稳定可靠,测井成本低。
附图说明:
附图1是本发明测量时仪器探测器的相对位置示意图。
图中:1是仪器外壳,2是传输电路,3是伽马探测器,4是热中子探测器,5是钨镍合金屏蔽体,6是同位素Am-Bi中子源6。
发明实施例
本发明采用热中子探测器和伽马探测器同时测量中子源子发射中子与地层作用后产生的热中子和弹性散射热中子俘获伽马信息。通过在已建立的已知空隙度、矿化度和岩性的系列模型井中刻度得到在套管条件下定量求解地层视空隙度和含水饱和度的数学解释模型。视空隙度是通过中子与伽马氯能段计数的加权和与空隙度的响应关系来确定的。在没有扩径且固井质量良好的情况下与裸眼井的空隙度值相等。含水饱和度是利用矿化度对热中子和中子俘获伽马的双重影响,即随着矿化度的增加热中子数减少而中子俘获伽马数增加;矿化度减少时则相反。通过双重影响所带来的双向差异的变化,确定地层的含水饱和度。
本发明应用双向差异方法确定地层的含水饱和度,是与地层的含氢指数及井眼状况密切相关的。裸眼井的孔隙度资料有两种情况会不可靠:一是当生产井开发时间较长,地层孔隙可能会发生变化;二是有些年代较长的老井,测井资料不完整。在套管井条件下,氯能谱测井可以给出一条现有状态下的视孔隙度指示曲线。当固井状况良好时,这个视孔隙度值应当与裸眼井测量的地层孔隙度相等,可以应用裸眼井的孔隙度资料;而当固井状况不好或是有明显扩径时,这种视孔隙度与裸眼井的孔隙度资料就会产生偏差,而当这两者孔隙度相差较大时,这时就需要根据两者间的差异,对计算方法进行经验校正。
本发明采用热中子探测器和伽马探测器组合并同时测量中子源子发射中子与地层作用后产生的热中子和弹性散射伽马信息,通过测量中子与中子伽马相对于饱和淡水或油地层的差异,求解出产层的含水饱和度。
根据已知地层矿化度、孔隙度、岩性建立刻度模型井,在模拟条件建立视孔隙度、含水饱和度关系函数,
计算视孔隙度
由于矿化度对于热中子和中子伽马测量结果的影响是相反的,因而将热中子与中子伽马的加权和与空隙度建立关系可以基本不受矿化度的影响。依据实验关系可建立如下的确定的视孔隙度的公式:
Φ=aG3+bG2+cG+d……(1)
式中:
Φ是视孔隙度;a、b、c是拟和公式的系数,是使用的热中子和中子伽马测量仪器结构、精度的参数;
G=k1Nz+K2Nr,
k1,k2是加权系数,Nz、Nr分别表示热中子和伽马氯能段测量的计数;
求含水饱和度
1)确定归一双向差值
a)建立归一函数关系
对应着某一地层,每一个中子测量值都有相应的氯测量值与其对应,在淡水地层中对应不同的空隙度,我们可以建立如下的中子和伽马的响应关系
N’z=K1NCl+K2……(2)
式中:K1,K2是拟合系数,N’z,NCl分别是测量的热中子和伽马氯能段计数值;(2)式就是氯与中子的归一关系式。
b)测井所得氯计数值Nr代入公式(2)得出归一的相应中子计数值N’z;
c)氯曲线与中子曲线的双向归一差异值Dt=N’z-Nz;……(3)
2)确定相对差Δ
a)由实验关系可以确定在对应地层空隙度为Φ时,含水饱和度Sw为100%时氯与中子的归一差D,其关系如下式下:
D=KPΦ ……(4)
式中:K是实验关系刻度系数,P是地层水矿化度,Φ地层的矿化度。
b)求出相对差Δ=1-Dt/D
3)求含水饱和度Sw
为了对不同的矿化度值给出一个统一的求解含水饱和度的公式,先对前面求出的Δ值作一个矿化度校正,其关系式如下:
δ=Δ+RP+M ……(5)
式中:R,M是实验关系系数,是使用的热中子和中子伽马测量仪器结构、精度的参数。
根据实验关系,得到含水饱和度得相应函数,其公式如下:
Sw=S3δ3+S2δ2+S1δ+S …… (6)
式中S3,S2,S1S是实验拟合关系系数,是使用的热中子和中子伽马探测器器结构、探测效率的参数。
本发明采用的数学公式,其形式不是唯一的,通常类似的公式可以有多种,例如将Sw的三次多项式变为分段的直线函数,空隙度公式写成对数形式等。
本发明采用热中子探测器和伽马探测器是常规的市售仪器,只要将其排列组合,所使用的记录仪器和方法也是与上述仪器有关的公知手段。
本发明可以克服测井时泥质和骨架的影响的主要机理由以下计算证明:
若以Xh、Xw、Xma和Xsh,分别表示油气、水、骨架和泥质的某一物理参数,则该物理参数的测井值Xlog是上述物理参数与响应体积Vi乘积之和。即:
Xlog=φXwSk+φXwSw+VshXsh+(1-φ-Vsh)Xma
式中等号后第一项为油气对测量响应的贡献;第二项为含水对测井值的影响;第三项为泥质对测井值的影响;第四项为骨架测井值的影响。而采用差值法就有:
ΔXlog=Δ(φXwSk+φXwSw)
就是说,如果在淡水地层(在中子测井上等效100%的纯油层)测量得到一条基线,那么地层由于注水,使测井响应值产生的差异,在理论上就只与孔隙流体中的含水饱和度及注入水的矿化度相关。
Claims (5)
1、一种氯能谱测井方法,其特征在于:采用热中子探测器和伽马探测器同时测量中子源子发射中子与地层作用后产生的热中子和弹性散射伽马信息,通过测量中子与中子伽马相对于饱和淡水或油地层的差异,求解出产层的视孔隙度和含水饱和度。
2、如权利要求1所述的一种氯能谱测井方法,其特征在于根据已知地层矿化度、孔隙度、岩性建立刻度模型井,在模拟条件建立视孔隙度关系函数,采用如下公式计算视孔隙度:Φ=aG3+bG2+cG+d,式中:Φ是视孔隙度;a、b、c、d是拟和公式的系数,是使用的热中子和中子伽马测量仪器结构、精度的参数;G=k1Nz+K2Nr,k1,k2是加权系数,Nz、Nr分别表示热中子和伽马氯能段测量的计数。
3、如权利要求1所述的一种氯能谱测井方法,其特征在于根据已知地层矿化度、孔隙度、岩性建立刻度模型井,在模拟条件建立含水饱和度关系函数,求含水饱和度,步骤是:
建立归一函数关系,对应地层,每一个中子测量值都有相应的氯测量值与其对应,在淡水地层中对应不同的空隙度,氯与中子的归一关系式;N’z=K1NCl+K2,式中:K1,K2是拟合系数,N’z,NCl分别是测量的热中子和伽马氯能段计数值;
将测井所得氯计数值Nr代入上式得出归一的相应中子计数值N’z;
求氯曲线与中子曲线的双向归一差异值Dt=N’z-Nz;
确定相对差△:
在对应地层孔隙度为φ时,含水饱和度Sw为100%时氯与中子的归一差D,其关系D=KPΦ,式中:K是实验关系刻度系数,P是地层水矿化度;
求出相对差△=1-Dt/D;
为了对不同的矿化度值给出一个统一的求解含水饱和度的公式,先对前面求出的△值作一个矿化度校正,关系式:δ=△+RP+M,式中:R,M是实验关系系数,是使用的测量仪器结构、精度的参数;
含水饱和度Sw=S3δ3+S2δ2+S1δ+S,式中S3,S2,S1S是实验拟合关系系数,是使用的测量仪器结构、精度的参数。
4、如权利要求1、2所述的一种氯能谱测井方法,其特征在于可以将视孔隙度Φ公式写成对数形式。
5、如权利要求1、2所述的一种氯能谱测井方法,其特征在于可以将Sw的三次多项式变为分段的直线函数形式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB031567606A CN100492055C (zh) | 2003-09-10 | 2003-09-10 | 一种氯能谱测井方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB031567606A CN100492055C (zh) | 2003-09-10 | 2003-09-10 | 一种氯能谱测井方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1595202A true CN1595202A (zh) | 2005-03-16 |
CN100492055C CN100492055C (zh) | 2009-05-27 |
Family
ID=34660065
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB031567606A Expired - Fee Related CN100492055C (zh) | 2003-09-10 | 2003-09-10 | 一种氯能谱测井方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100492055C (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102330552A (zh) * | 2010-07-13 | 2012-01-25 | 普拉德研究及开发股份有限公司 | 对中子伽马密度测量的校正 |
CN102518431A (zh) * | 2011-12-26 | 2012-06-27 | 中国石油大学(华东) | 基于可控中子源的随钻多参数测井方法 |
CN102084271B (zh) * | 2008-06-25 | 2013-11-06 | 普拉德研究及开发股份有限公司 | 由核谱学确定绝对元素浓度 |
CN108333638A (zh) * | 2018-04-18 | 2018-07-27 | 四川久天测控技术有限公司 | 一种小直径碳氢比饱和度井下仪及测井仪系统 |
CN112302622A (zh) * | 2020-10-29 | 2021-02-02 | 大庆油田有限责任公司 | 一种脉冲中子全谱测井综合解释剩余油饱和度的方法 |
CN113622904A (zh) * | 2020-04-21 | 2021-11-09 | 中国石油天然气股份有限公司 | 用元素俘获谱测井氯离子含量识别测井流体类型的方法 |
RU2771438C1 (ru) * | 2021-09-20 | 2022-05-04 | Общество с ограниченной ответственностью "Институт нефтегазовых технологий "ГеоСпектр" (ООО "ИНТ "ГеоСпектр") | Способ определения лития в рапонасыщенных интервалах геологических разрезов скважин газоконденсатных месторождений по данным мультиметодного многозондового нейтронного каротажа |
CN116066069A (zh) * | 2021-02-11 | 2023-05-05 | 中国石油化工股份有限公司 | 用脉冲中子源和双功能探测器进行地层孔隙度测量的方法 |
RU2815325C1 (ru) * | 2023-09-19 | 2024-03-13 | Общество с ограниченной ответственностью "Институт нефтегазовых технологий "ГеоСпектр" (ООО "ИНТ "ГеоСпектр") | Способ оценки нефтенасыщенности коллекторов в обсаженных нефтегазовых и нефтегазоконденсатных скважинах с высокой минерализацией пластовых вод методом мультиметодного многозондового нейтронного каротажа - ММНК |
-
2003
- 2003-09-10 CN CNB031567606A patent/CN100492055C/zh not_active Expired - Fee Related
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102084271B (zh) * | 2008-06-25 | 2013-11-06 | 普拉德研究及开发股份有限公司 | 由核谱学确定绝对元素浓度 |
CN102330552A (zh) * | 2010-07-13 | 2012-01-25 | 普拉德研究及开发股份有限公司 | 对中子伽马密度测量的校正 |
CN102330552B (zh) * | 2010-07-13 | 2016-06-08 | 普拉德研究及开发股份有限公司 | 对中子伽马密度测量的校正 |
CN102518431A (zh) * | 2011-12-26 | 2012-06-27 | 中国石油大学(华东) | 基于可控中子源的随钻多参数测井方法 |
CN102518431B (zh) * | 2011-12-26 | 2015-04-22 | 中国石油大学(华东) | 基于可控中子源的随钻多参数测井方法 |
CN108333638A (zh) * | 2018-04-18 | 2018-07-27 | 四川久天测控技术有限公司 | 一种小直径碳氢比饱和度井下仪及测井仪系统 |
CN113622904A (zh) * | 2020-04-21 | 2021-11-09 | 中国石油天然气股份有限公司 | 用元素俘获谱测井氯离子含量识别测井流体类型的方法 |
CN113622904B (zh) * | 2020-04-21 | 2024-06-25 | 中国石油天然气股份有限公司 | 用元素俘获谱测井氯离子含量识别测井流体类型的方法 |
CN112302622A (zh) * | 2020-10-29 | 2021-02-02 | 大庆油田有限责任公司 | 一种脉冲中子全谱测井综合解释剩余油饱和度的方法 |
CN116066069A (zh) * | 2021-02-11 | 2023-05-05 | 中国石油化工股份有限公司 | 用脉冲中子源和双功能探测器进行地层孔隙度测量的方法 |
RU2771438C1 (ru) * | 2021-09-20 | 2022-05-04 | Общество с ограниченной ответственностью "Институт нефтегазовых технологий "ГеоСпектр" (ООО "ИНТ "ГеоСпектр") | Способ определения лития в рапонасыщенных интервалах геологических разрезов скважин газоконденсатных месторождений по данным мультиметодного многозондового нейтронного каротажа |
RU2815325C1 (ru) * | 2023-09-19 | 2024-03-13 | Общество с ограниченной ответственностью "Институт нефтегазовых технологий "ГеоСпектр" (ООО "ИНТ "ГеоСпектр") | Способ оценки нефтенасыщенности коллекторов в обсаженных нефтегазовых и нефтегазоконденсатных скважинах с высокой минерализацией пластовых вод методом мультиметодного многозондового нейтронного каротажа - ММНК |
Also Published As
Publication number | Publication date |
---|---|
CN100492055C (zh) | 2009-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1158541C (zh) | 测定地层密度的方法及仪器 | |
US8311744B2 (en) | Use of chemically and geologically consistent constraints for obtaining elemental chemistry downhole | |
US4585939A (en) | Multi-function natural gamma ray logging system | |
CN1209500A (zh) | 用于测量钻井特征和地层特性的方法和器具 | |
EP2831646A1 (en) | Neutron porosity based on one or more gamma ray detectors and a pulsed neutron source | |
CN101906963A (zh) | 利用c/o和地层宏观俘获截面交会技术确定饱和度方法 | |
CN103088803B (zh) | 鉴别岩溶塌陷致塌因素的方法 | |
GB2347497A (en) | Calibrating readings of a downhole tool | |
CN106250619B (zh) | 一种确定地层矿物含量的方法和装置 | |
CN110192125A (zh) | 根据井眼测量的地层水盐度 | |
CN1595202A (zh) | 一种氯能谱测井方法 | |
CN1619100A (zh) | 确定井孔周围结构孔隙度的中子测量方法 | |
CN103744121B (zh) | 碳氢比地层流体饱和度测井方法 | |
CN108661632B (zh) | 一种高精度碳氧比测井方法 | |
CN112523741B (zh) | 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法 | |
Schweitzer et al. | Elemental Concentrations From Gamma Ray Spectroscopy Logs | |
Scott et al. | Response of a multidetector pulsed neutron porosity tool | |
CN1151385C (zh) | 确定岩层中油饱和度的方法 | |
Liu et al. | Methods for evaluating elemental concentration and gas saturation by a three-detector pulsed-neutron well-logging tool | |
Odom et al. | A New 1. 625" Diameter Pulsed Neutron Capture And Inelastic/Capture Spectral Combination System Provides Answers In Complex Reservoirs | |
CN112523742B (zh) | 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法 | |
CN106958443B (zh) | 一种碳氧比测井钻井液侵入影响因素的校正方法及系统 | |
CN115614022B (zh) | 铀裂变瞬发中子测井的井眼水层原位修正方法 | |
Gilchrist Jr et al. | Initial Field Applications of a New 1.7-Inch Pulsed Neutron Instrument | |
CN112593923B (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 | ||
C56 | Change in the name or address of the patentee | ||
CP02 | Change in the address of a patent holder |
Address after: 710077, No. 85, No. two, No. 50, Kam Industrial Zone, Xi'an hi tech Development Zone, Shaanxi, China Patentee after: Technique Center of China Petroleum Logging Co., Ltd. Address before: Fengcheng Weiyang District four road 710021 Shaanxi city of Xi'an Province Patentee before: Technique Center of China Petroleum Logging Co., Ltd. |
|
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090527 Termination date: 20140910 |
|
EXPY | Termination of patent right or utility model |