自然伽马能谱测井谱解析方法研究

第32卷2010年第5期10月

西南石油大学学报（自然科学版）

Journal of Southwest Petroleum University （Science ＆Technology Edition ）

Vol．32Oct．

No．52010

文章编号：1674－5086（2010）05－0075－04

自然伽马能谱测井谱解析方法研究

123

刘尊年，任爱阁，贾瑞皋

（1．青岛理工大学理学院，山东青岛266033；2．青岛滨海学院基础部，山东青岛266555；

3．中国石油大学物理科学与技术学院，山东东营257062）

*

摘要：自然伽马能谱测井是重要的核测井方法之一，其地面信息处理系统的关键是实测谱解析技术。低分辨率是

自然伽马能谱测井仪的固有弱点，且测井过程中信噪比较低，造成测井置信度不高。从自然伽马能谱测井原理出发，通过中海油实验数据处理后获得的放射性核素的标准计数率谱，对自然伽马能谱测井4种传统数据处理方法的优劣进行了对比，初步探讨了将灵敏度高、空间分辨率好的直接解调方法移植到自然伽马能谱测井谱解析中，为自然伽马能谱测井谱解析方法的选择提供了依据。

关键词：自然伽马能谱测井；谱解析；直接解调方法；响应函数；物理约束中图分类号：TE19

文献标识码：A

DOI ：10．3863/j．issn．1674－5086．2010．05．012

214232

系中Bi 发射的1．76Mev 伽马射线；表征Th 选用Th 系中208Tl 发射的2．62Mev 伽马射线；表征K 只40［14］

能用K 唯一的1．46Mev 伽马射线。1．2自然伽马能谱测井原理

自然伽马能谱测井是核探测技术在地球物理勘探中的应用，其地面信息处理系统的关键是实测谱解析技术，剥谱法、逆矩阵法、最小二乘逆矩阵法、加

［1－3］

权最小二乘逆矩阵法是4种常用的传统谱解析方法。因为自然伽马能谱测井仪探测器的固有分辨率较低，测井过程中统计涨落比较严重，信噪比较低，而直接解调方法是一种先进的数据处理技术，实践证明，该方法用于处理统计性较差的测量数据时

［4－5］

自然伽马能谱测井可以同时做到强度测量和能

谱分析，根据能量分析确定射线由哪一种核素发射，依据强度确定该种核素及其相应元素的含量。对自然伽马能谱信息利用最为充分的是多功能补偿自然伽马能谱测井仪，其测井过程如图1所示

。

，其灵敏度高于传统谱解析方法。

1自然伽马能谱测井原理

Th 、K 自然伽马能谱测井能直接测量地层中U 、

等元素的含量，以此来判断地层性质，定量计算可用于确定泥质含量、鉴别高渗透率、碳酸盐岩裂缝层的划分和粘土成分的确定；定性分析可用于对泥质地层生油岩特性的识别和沉积环境的研究1．1自然伽马能谱

［1，6－10］

。

图1

能谱仪探测过程

Fig．1

Sound process of Natural Gamma Ray Spectra log

自然伽马射线主要是由地层岩石中的U 放射

23240［11］

系、Th 放射系和K 产生的。原始伽马射线能谱均呈现离散的线谱形态，而伽马射线触发伽马谱仪探测器后却得到连续的谱线，通常把这种连续谱称为伽马实测谱或仪器谱

［12－13］

238

2传统谱解析方法

。

238

在自然伽马能谱测井中，表征

：U 通常选用U

剥谱法、逆矩阵法、最小二乘逆矩阵法和加权最

小二乘逆矩阵法是自然伽马能谱测井谱解析的4种

作者简介：刘尊年（1974－），男（汉族），山东郓城人，讲师，硕士，、。

76

主要方法。

2．1自然伽马标准谱

西南石油大学学报（自然科学版）2010年

先找出一种容易识别的核素，将其谱形求出，并从混

然后从剩余谱中再找出第二种核素并合谱中扣除，

直到求出所有核素为止。在利用中进行同样处理，

海油实验数据进行剥谱计算时，先将谱峰划分成3个能窗，即：98～115道为能窗1；120～140道为能窗2；177～198道为能窗3。各能窗将K ，U 和Th 的特1．76Mev 及2．62Mev 恰当地包征能量峰1．46Mev ，

［15］

括在中间。2．2．2

逆矩阵法

逆矩阵法在划分道区时允许各特征峰计数有重叠，不受道区数和道区宽度的限制，能够更充分地利U 和Th 用谱中的计数。利用逆矩阵法解谱时，在K 、特征能量峰所在谱段连续取3个道区，即：95～115道为能窗1；116～144道为能窗2；145～198道为能窗3。2．2．3

最小二乘逆矩阵法

剥谱法和逆矩阵法都是根据核素个数选择特征峰道区，而最小二乘逆矩阵法选用的道区数可多于待求核素数，且道区不一定是特征峰道区，这样就可

U 和充分利用谱中所有比较重要的计数峰确定K ，Th 的含量。显然，这种解谱方法的统计精度更高。2．2．4

加权最小二乘逆矩阵法

假定了各道区皆为等精度在上述3种方法中，

标准谱是在刻度井中测得的仪器谱，对实测谱

Th 和K 含量与谱数据的关进行解析的依据是U 、

Ra 、Th 、K 的系。表1给出了各层位放射性元素U 、K 均以10－2为单位，含量（如不特殊说明，而U 和Th

－6

均以10为单位）。

表1

Ra 、Th 、K 含量表燕郊刻度井各层位放射性元素U 、

Table 1

U /

10－60．400．790．9722．102．4012．30

The contents sheet of U 、Ra 、Th 、K from

China Offshore Oil Logging Ra /10－12

1．72．32．675．18．339．440．9

（Ra /U）×100%

－－－101．4103．295．6100．1

Th /10－61．1263．501．451．215．4134．8036．00

K /10－20．230．225．600．251．055．175．11

层围岩钍层钾层铀层低混层高混层

高混薄层12．20

在刻度井中，用自然伽马能谱测井仪测量只含U 、Th 或K 一种放射性元素的模拟地层，可得到该种元素的标准仪器谱。解谱的依据是由标准谱确定的

Th 和K 在各道或谱段中标准条件下单位含量的U 、

的计数率。实际测井中，井口尺寸与8．5in．刻度井

相近，因此用8．5in．井所得的数据描述自然伽马标Th 和K 三种核素及高混合层、准谱是合适的。U 、低混合层的5种模块的标准谱如图2所示

。

观测。实际上，低能道区的计数率比高能道区高得多，鉴于此，常采用加权最小二乘法来提高非等精度在测量最优值的求取过程中要观测值的处理精度，

使每个道区计数率权重误差平方和趋于最小。这种方法的缺点是：（1）道计数率低，统计精度不高；（2）元素的标准谱难以测定。2．3

四种传统谱解析方法对比

表3和表4显示：误差最大的是剥谱法；误差相

近、效果较好的是逆矩阵法和加权最小二乘逆矩阵法；最小二乘逆矩阵法优于前者、劣于后者。

3直接解调方法

图2

Fig．2

五种模块的标准谱

20世纪90年代初，李惕碚和吴枚提出了直接

［4－5］

，解调方法目的是为了从统计信息少和信噪比低的信号数据中得到复杂天体的高分辨率空间图

像。相对于传统的反演方法，这种方法对观测信息的利用更为充分，其精度是传统方法不可能达到的。－Standard Spectra of the five Modules

2．22．传统谱解析方法的理论及验证：

第5期刘尊年，等：自然伽马能谱测井谱解析方法研究表3

17#仪器的测量值对高混合层，

U

相对误差/%10．532．234．222．39

测量值/10－611．105011．972912．983812．5666

相对误差/%

9．722．665．562．17

34．210833．238132．184932．7411

Th

测量值/10－6

相对误差/%

1．694．497．515．92

77

Table 3

核素

17#剥谱法逆矩阵法最小二乘逆矩阵法加权最小二乘逆矩阵法

The measure values of the 17#SGR from high mixer K

测量值/10－2

4．62575．05464．95165．0463

表4

Table 4

核素17#剥谱法逆矩阵法最小二乘逆矩阵法加权最小二乘逆矩阵法

K

测量值/10－2

1．07951．18021．09961．1416

17#仪器的测量值对低混合层，

U

相对误差/%

2．8112．404．728．72

测量值/10－6

2．58732．86073．22562．9841

相对误差/%

7．8019．2034．4024．34

5．91335．47174．97065．3735

Th

测量值/10－6

相对误差/%

9．301．148．120．67

The measure values of the 17#SGR from low mixer

［18］

获得了理想的效果。例如，用直接解调方法重新分析欧洲X 射线天文卫星EXOSAT 的扫描数据，得［19］到以前无法找到的一些银河系X 射线源，这些源已被列入SIMBAD 天体表；直接解调方法应用于光

可以认为不变。在研究刻度井实验过程中，情况下，

由于测量环境一定，胜利油田测井公司17#SL6329能谱仪器的响应函数可视为定值，根据得到的标准计数率谱，可以求出仪器的响应函数，如图3所示。

学相干层析OCT 中，可使图像的纵向分辨率提高一

［20］

个数量级。3．1

直接解调方法的原理

直接解调方法的本质是把合理的非线性约束引入到病态观测方程的求解过程，其核心是：首先建立观测方程，选择适当的物理上下限约束，用迭代方法求解此方程。观测过程属于辐射探测过程，可用如下方程描述

N

P （k ，i ）f （i ）∑i =1式中

=d （k ），（k =1，2，…，M ）（1）

图3

Tool 05017对K 、U 和Th 元素的响应

d （k ）—观测数据；

P （k ，i ）—仪器的响应函数（调制函数）；f （i ）—观测天区辐射强度。3．2

自然伽马能谱测井中的直接解调方法直接解调方法的数学本质与自然伽马能谱测井谱解析一样，都是解病态方程，并且二者具有相同的物理、能谱探测基础等。辐射探测过程可以看成是源和探测系统作卷积，得到像函数

P *f =d

（2）

，Fig．3Response function of K 、U and Th from Tool 05017

如果引入一定约束条件并采用逐次迭代法直接

解谱，既可以抑制解的振荡，又会使系统的分辨率大大提高。

4结语

（1）得到了刻度井中放射性核素的标准计数率谱，自然伽马能谱测井4种传统数据处理方法中效。

78西南石油大学学报（自然科学版）2010年

数，初步探讨了将直接解调方法移植到自然伽马能

为谱解析在方法选择上谱测井谱解析中的可行性，提供了可靠的依据。

．岩性油气藏，2007，19（2）：108－能谱测井方法［J ］111．

［10］王谦．自然伽马能谱测井在塔河奥陶系地层划分中

J ］．石油物探，2004（9）：504－507．的应用［

［11］黄凯，徐群洲，杨晓海，等．地震波能量的衰减及其影

J ］．新疆石油地质，1997，9（3）：212－216．响因素［

［12］庞巨丰，迟云鹏，钟振伟．现代核测井技术与仪器

［M ］．北京：石油工业出版社，1998．

［13］黄隆基，J ］．地胡庆东．自然伽马能谱测井薄层响应［

1997，40（2）：272－278．球物理学报，

［14］Grau J A ，Schweitzer J S．Prompt γ?Ray spectral analysis

of well data obtained with NaI and 14MeV neutrons ［J ］．Nuel Geophys ，1987，1（2）：157－165．

［15］任爱阁．自然伽马能谱测井谱解析方法研究［D ］．山

2007．东东营：中国石油大学，

［16］李惕碚，吴枚，陆柱国，等．调制成像γ射线望远镜

［J ］．天文学报，1994，35（1）：105－108．

［17］张澍，李惕碚，吴枚．Compton γ射线望远镜的直接解

J ］．天体物理学报，1997，17（3）：263－270．调成像［

［18］刘力．直接解调成像方法在核医学图像重建中的应

J ］．核技术，1998，21（3）：188－192．用［

［19］卢方军，李惕碚，吴枚，等．ME/EXOSAT银道扫描数据的

［J ］．天文学报，1997，38（1）：56－65．直接解调成像

［20］孙非，薛平，高湔松，等．光学相关层析成像的图像重

J ］．光学学报，2000，20（8）：1043－1046．建［

（编辑：杜增利）

参考文献：

［1］黄隆基．放射性测井原理［M ］．北京：石油工业出版

1985．社，

［2］黄隆基，首祥云，王瑞平．自然伽马能谱测井原理及应

M ］．北京：石油工业出版社，1995：176－195．用［

［3］黄隆基．核测井原理［M ］．山东东营：石油工业出版

2000．社，

［4］李惕碚，吴枚．高能天文学中成像和解谱的直接方法

［J ］．天体物理学报，1993，17（3）：263－270．

［5］李惕碚．能不能用低分辨率仪器实现高分辨率观测

［J ］．大自然探索，1999，18（67）．

［6］彭曙光，潘和平，王霁霞，等．自然伽马能谱井在储层

．工程地球物理学报，2009，6（4）：评价中的应用［J ］480－484．

［7］赵发展，王尚明，尉中良，等．自然伽马能谱测井确定

J ］．国外测井技术，粘土矿物类型及分布的应用研究［2004，19（2）：47－50．

［8］孙建孟，李召成．应用自然伽马能谱测井确定粘土矿

J ］．石油大学学报：自然科学版，1999，物类型和含量［23（4）：29－32．

［9］王祝文，刘菁华，黄茜．确定粘土矿物含量的自然伽马

·7·

Journal of Southwest Petroleum University （Science ＆Technology Edition ）No．5

Zhanjiang Company ，CNOOC ，Zhanjiang Guangdong 524057，China ；3．Research Institute of Exploration and Develop-Chuanqing Drilling ＆Exploration Company ，CNPC ，Chengdu Sichuan 610051，China ）JOURNAL OF SOUTH-ment ，

VOL．32，NO．5，70－74，2010（ISSN WEST PETROLEUM UNIVERSITY （SCIENCE ＆TECHNOLOGY EDITION ），1674-5086，in Chinese ）

Abstract ：The reservoir of Z1oil group in W Oilfield has low resistance ，high-porosity and low permeability．The geological formation of low resistivity reservoir has many reasons．In view of the developing features of Z oil group ，the clay particles distribution ，clay content and the influence of them on the low-resistivity are studied．clay minerals types ，

It can be found that the clay types in low resistivity layer segment of Z1oil group is mainly the mixed layer of illite and montomorillonite and the content is relatively high．The clay particles are evenly distributed in a thin film and in flocculating form with high average clay content．Additional conductivity caused by clay minerals and high clay content are the main factors affecting the low resistivity of Z1oil group．

Key words ：low-resistivity reservoir ；clay minerals ；additional conductivity ；clay content ；genesis of low resistivity

STUDY ON SPECTRAL ANALYSIS OF NATURAL GAMMA RAY SPECTRA LOG

LIU Zun-nian 1，REN Ai-ge 2，JIA Rui-gao 3（1．School of Science ，Qingdao Technological University ，Qingdao Shan-China ；2．Department of Basic Courses ，Qingdao Bin Hai University ，Qingdao Shandong 266555，Chi-dong 266033，

na ；3．College of Physics Science and Technology ，China University of Petroleum ，Dongying Shandong 257062，Chi-VOL．32，na ）JOURNAL OF SOUTHWEST PETROLEUM UNIVERSITY （SCIENCE ＆TECHNOLOGY EDITION ），NO．5，75－78，2010（ISSN 1674－5086，in Chinese ）

Abstract ：Natural Gamma Ray Spectra log （SGR ）is an important nuclear logging method ，and real spectral analysis is a key technique of the surface processing system in SGR．Due to the inherent low resolution of Natural Gamma there is a low signal-to-noise ratio and low confidence level during the logging process．Based on Ray Spectra log ，

the principle of Natural Gamma Ray Spectra log ，the normal count rate spectra are obtained from the test pit ，and the the authors perform preliminary research on applying the four traditional algorithm methods are tested．Besides ，

direct demodulation method which has a high sensitivity and spatial resolution in gamma-ray log．This research of-fers a credible basis for algorithm selecting．

Key words ：Natural Gamma Ray spectrometry log ；spectrum unscrambling ；direct demodulation method ；response function ；physical constrains

GENETIC MECHANISM OF LOW RESISTIVITY RESERVOIR IN BANQIAO OIL FIELD

YE Xing-shu ，WU Guo-hai （GreatWall Drilling Company ，CNPC ，Chaoyang Beijing 100101，China ）JOURNAL OF ）5

第32卷2010年第5期10月

西南石油大学学报（自然科学版）

Journal of Southwest Petroleum University （Science ＆Technology Edition ）

Vol．32Oct．

No．52010

文章编号：1674－5086（2010）05－0075－04

自然伽马能谱测井谱解析方法研究

123

刘尊年，任爱阁，贾瑞皋

（1．青岛理工大学理学院，山东青岛266033；2．青岛滨海学院基础部，山东青岛266555；

3．中国石油大学物理科学与技术学院，山东东营257062）

*

摘要：自然伽马能谱测井是重要的核测井方法之一，其地面信息处理系统的关键是实测谱解析技术。低分辨率是

自然伽马能谱测井仪的固有弱点，且测井过程中信噪比较低，造成测井置信度不高。从自然伽马能谱测井原理出发，通过中海油实验数据处理后获得的放射性核素的标准计数率谱，对自然伽马能谱测井4种传统数据处理方法的优劣进行了对比，初步探讨了将灵敏度高、空间分辨率好的直接解调方法移植到自然伽马能谱测井谱解析中，为自然伽马能谱测井谱解析方法的选择提供了依据。

关键词：自然伽马能谱测井；谱解析；直接解调方法；响应函数；物理约束中图分类号：TE19

文献标识码：A

DOI ：10．3863/j．issn．1674－5086．2010．05．012

214232

系中Bi 发射的1．76Mev 伽马射线；表征Th 选用Th 系中208Tl 发射的2．62Mev 伽马射线；表征K 只40［14］

能用K 唯一的1．46Mev 伽马射线。1．2自然伽马能谱测井原理

自然伽马能谱测井是核探测技术在地球物理勘探中的应用，其地面信息处理系统的关键是实测谱解析技术，剥谱法、逆矩阵法、最小二乘逆矩阵法、加

［1－3］

权最小二乘逆矩阵法是4种常用的传统谱解析方法。因为自然伽马能谱测井仪探测器的固有分辨率较低，测井过程中统计涨落比较严重，信噪比较低，而直接解调方法是一种先进的数据处理技术，实践证明，该方法用于处理统计性较差的测量数据时

［4－5］

自然伽马能谱测井可以同时做到强度测量和能

谱分析，根据能量分析确定射线由哪一种核素发射，依据强度确定该种核素及其相应元素的含量。对自然伽马能谱信息利用最为充分的是多功能补偿自然伽马能谱测井仪，其测井过程如图1所示

。

，其灵敏度高于传统谱解析方法。

1自然伽马能谱测井原理

Th 、K 自然伽马能谱测井能直接测量地层中U 、

等元素的含量，以此来判断地层性质，定量计算可用于确定泥质含量、鉴别高渗透率、碳酸盐岩裂缝层的划分和粘土成分的确定；定性分析可用于对泥质地层生油岩特性的识别和沉积环境的研究1．1自然伽马能谱

［1，6－10］

。

图1

能谱仪探测过程

Fig．1

Sound process of Natural Gamma Ray Spectra log

自然伽马射线主要是由地层岩石中的U 放射

23240［11］

系、Th 放射系和K 产生的。原始伽马射线能谱均呈现离散的线谱形态，而伽马射线触发伽马谱仪探测器后却得到连续的谱线，通常把这种连续谱称为伽马实测谱或仪器谱

［12－13］

238

2传统谱解析方法

。

238

在自然伽马能谱测井中，表征

：U 通常选用U

剥谱法、逆矩阵法、最小二乘逆矩阵法和加权最

小二乘逆矩阵法是自然伽马能谱测井谱解析的4种

作者简介：刘尊年（1974－），男（汉族），山东郓城人，讲师，硕士，、。

76

主要方法。

2．1自然伽马标准谱

西南石油大学学报（自然科学版）2010年

先找出一种容易识别的核素，将其谱形求出，并从混

然后从剩余谱中再找出第二种核素并合谱中扣除，

直到求出所有核素为止。在利用中进行同样处理，

海油实验数据进行剥谱计算时，先将谱峰划分成3个能窗，即：98～115道为能窗1；120～140道为能窗2；177～198道为能窗3。各能窗将K ，U 和Th 的特1．76Mev 及2．62Mev 恰当地包征能量峰1．46Mev ，

［15］

括在中间。2．2．2

逆矩阵法

逆矩阵法在划分道区时允许各特征峰计数有重叠，不受道区数和道区宽度的限制，能够更充分地利U 和Th 用谱中的计数。利用逆矩阵法解谱时，在K 、特征能量峰所在谱段连续取3个道区，即：95～115道为能窗1；116～144道为能窗2；145～198道为能窗3。2．2．3

最小二乘逆矩阵法

剥谱法和逆矩阵法都是根据核素个数选择特征峰道区，而最小二乘逆矩阵法选用的道区数可多于待求核素数，且道区不一定是特征峰道区，这样就可

U 和充分利用谱中所有比较重要的计数峰确定K ，Th 的含量。显然，这种解谱方法的统计精度更高。2．2．4

加权最小二乘逆矩阵法

假定了各道区皆为等精度在上述3种方法中，

标准谱是在刻度井中测得的仪器谱，对实测谱

Th 和K 含量与谱数据的关进行解析的依据是U 、

Ra 、Th 、K 的系。表1给出了各层位放射性元素U 、K 均以10－2为单位，含量（如不特殊说明，而U 和Th

－6

均以10为单位）。

表1

Ra 、Th 、K 含量表燕郊刻度井各层位放射性元素U 、

Table 1

U /

10－60．400．790．9722．102．4012．30

The contents sheet of U 、Ra 、Th 、K from

China Offshore Oil Logging Ra /10－12

1．72．32．675．18．339．440．9

（Ra /U）×100%

－－－101．4103．295．6100．1

Th /10－61．1263．501．451．215．4134．8036．00

K /10－20．230．225．600．251．055．175．11

层围岩钍层钾层铀层低混层高混层

高混薄层12．20

在刻度井中，用自然伽马能谱测井仪测量只含U 、Th 或K 一种放射性元素的模拟地层，可得到该种元素的标准仪器谱。解谱的依据是由标准谱确定的

Th 和K 在各道或谱段中标准条件下单位含量的U 、

的计数率。实际测井中，井口尺寸与8．5in．刻度井

相近，因此用8．5in．井所得的数据描述自然伽马标Th 和K 三种核素及高混合层、准谱是合适的。U 、低混合层的5种模块的标准谱如图2所示

。

观测。实际上，低能道区的计数率比高能道区高得多，鉴于此，常采用加权最小二乘法来提高非等精度在测量最优值的求取过程中要观测值的处理精度，

使每个道区计数率权重误差平方和趋于最小。这种方法的缺点是：（1）道计数率低，统计精度不高；（2）元素的标准谱难以测定。2．3

四种传统谱解析方法对比

表3和表4显示：误差最大的是剥谱法；误差相

近、效果较好的是逆矩阵法和加权最小二乘逆矩阵法；最小二乘逆矩阵法优于前者、劣于后者。

3直接解调方法

图2

Fig．2

五种模块的标准谱

20世纪90年代初，李惕碚和吴枚提出了直接

［4－5］

，解调方法目的是为了从统计信息少和信噪比低的信号数据中得到复杂天体的高分辨率空间图

像。相对于传统的反演方法，这种方法对观测信息的利用更为充分，其精度是传统方法不可能达到的。－Standard Spectra of the five Modules

2．22．传统谱解析方法的理论及验证：

第5期刘尊年，等：自然伽马能谱测井谱解析方法研究表3

17#仪器的测量值对高混合层，

U

相对误差/%10．532．234．222．39

测量值/10－611．105011．972912．983812．5666

相对误差/%

9．722．665．562．17

34．210833．238132．184932．7411

Th

测量值/10－6

相对误差/%

1．694．497．515．92

77

Table 3

核素

17#剥谱法逆矩阵法最小二乘逆矩阵法加权最小二乘逆矩阵法

The measure values of the 17#SGR from high mixer K

测量值/10－2

4．62575．05464．95165．0463

表4

Table 4

核素17#剥谱法逆矩阵法最小二乘逆矩阵法加权最小二乘逆矩阵法

K

测量值/10－2

1．07951．18021．09961．1416

17#仪器的测量值对低混合层，

U

相对误差/%

2．8112．404．728．72

测量值/10－6

2．58732．86073．22562．9841

相对误差/%

7．8019．2034．4024．34

5．91335．47174．97065．3735

Th

测量值/10－6

相对误差/%

9．301．148．120．67

The measure values of the 17#SGR from low mixer

［18］

获得了理想的效果。例如，用直接解调方法重新分析欧洲X 射线天文卫星EXOSAT 的扫描数据，得［19］到以前无法找到的一些银河系X 射线源，这些源已被列入SIMBAD 天体表；直接解调方法应用于光

可以认为不变。在研究刻度井实验过程中，情况下，

由于测量环境一定，胜利油田测井公司17#SL6329能谱仪器的响应函数可视为定值，根据得到的标准计数率谱，可以求出仪器的响应函数，如图3所示。

学相干层析OCT 中，可使图像的纵向分辨率提高一

［20］

个数量级。3．1

直接解调方法的原理

直接解调方法的本质是把合理的非线性约束引入到病态观测方程的求解过程，其核心是：首先建立观测方程，选择适当的物理上下限约束，用迭代方法求解此方程。观测过程属于辐射探测过程，可用如下方程描述

N

P （k ，i ）f （i ）∑i =1式中

=d （k ），（k =1，2，…，M ）（1）

图3

Tool 05017对K 、U 和Th 元素的响应

d （k ）—观测数据；

P （k ，i ）—仪器的响应函数（调制函数）；f （i ）—观测天区辐射强度。3．2

自然伽马能谱测井中的直接解调方法直接解调方法的数学本质与自然伽马能谱测井谱解析一样，都是解病态方程，并且二者具有相同的物理、能谱探测基础等。辐射探测过程可以看成是源和探测系统作卷积，得到像函数

P *f =d

（2）

，Fig．3Response function of K 、U and Th from Tool 05017

如果引入一定约束条件并采用逐次迭代法直接

解谱，既可以抑制解的振荡，又会使系统的分辨率大大提高。

4结语

（1）得到了刻度井中放射性核素的标准计数率谱，自然伽马能谱测井4种传统数据处理方法中效。

78西南石油大学学报（自然科学版）2010年

数，初步探讨了将直接解调方法移植到自然伽马能

为谱解析在方法选择上谱测井谱解析中的可行性，提供了可靠的依据。

．岩性油气藏，2007，19（2）：108－能谱测井方法［J ］111．

［10］王谦．自然伽马能谱测井在塔河奥陶系地层划分中

J ］．石油物探，2004（9）：504－507．的应用［

［11］黄凯，徐群洲，杨晓海，等．地震波能量的衰减及其影

J ］．新疆石油地质，1997，9（3）：212－216．响因素［

［12］庞巨丰，迟云鹏，钟振伟．现代核测井技术与仪器

［M ］．北京：石油工业出版社，1998．

［13］黄隆基，J ］．地胡庆东．自然伽马能谱测井薄层响应［

1997，40（2）：272－278．球物理学报，

［14］Grau J A ，Schweitzer J S．Prompt γ?Ray spectral analysis

of well data obtained with NaI and 14MeV neutrons ［J ］．Nuel Geophys ，1987，1（2）：157－165．

［15］任爱阁．自然伽马能谱测井谱解析方法研究［D ］．山

2007．东东营：中国石油大学，

［16］李惕碚，吴枚，陆柱国，等．调制成像γ射线望远镜

［J ］．天文学报，1994，35（1）：105－108．

［17］张澍，李惕碚，吴枚．Compton γ射线望远镜的直接解

J ］．天体物理学报，1997，17（3）：263－270．调成像［

［18］刘力．直接解调成像方法在核医学图像重建中的应

J ］．核技术，1998，21（3）：188－192．用［

［19］卢方军，李惕碚，吴枚，等．ME/EXOSAT银道扫描数据的

［J ］．天文学报，1997，38（1）：56－65．直接解调成像

［20］孙非，薛平，高湔松，等．光学相关层析成像的图像重

J ］．光学学报，2000，20（8）：1043－1046．建［

（编辑：杜增利）

参考文献：

［1］黄隆基．放射性测井原理［M ］．北京：石油工业出版

1985．社，

［2］黄隆基，首祥云，王瑞平．自然伽马能谱测井原理及应

M ］．北京：石油工业出版社，1995：176－195．用［

［3］黄隆基．核测井原理［M ］．山东东营：石油工业出版

2000．社，

［4］李惕碚，吴枚．高能天文学中成像和解谱的直接方法

［J ］．天体物理学报，1993，17（3）：263－270．

［5］李惕碚．能不能用低分辨率仪器实现高分辨率观测

［J ］．大自然探索，1999，18（67）．

［6］彭曙光，潘和平，王霁霞，等．自然伽马能谱井在储层

．工程地球物理学报，2009，6（4）：评价中的应用［J ］480－484．

［7］赵发展，王尚明，尉中良，等．自然伽马能谱测井确定

J ］．国外测井技术，粘土矿物类型及分布的应用研究［2004，19（2）：47－50．

［8］孙建孟，李召成．应用自然伽马能谱测井确定粘土矿

J ］．石油大学学报：自然科学版，1999，物类型和含量［23（4）：29－32．

［9］王祝文，刘菁华，黄茜．确定粘土矿物含量的自然伽马

·7·

Journal of Southwest Petroleum University （Science ＆Technology Edition ）No．5

Zhanjiang Company ，CNOOC ，Zhanjiang Guangdong 524057，China ；3．Research Institute of Exploration and Develop-Chuanqing Drilling ＆Exploration Company ，CNPC ，Chengdu Sichuan 610051，China ）JOURNAL OF SOUTH-ment ，

VOL．32，NO．5，70－74，2010（ISSN WEST PETROLEUM UNIVERSITY （SCIENCE ＆TECHNOLOGY EDITION ），1674-5086，in Chinese ）

Abstract ：The reservoir of Z1oil group in W Oilfield has low resistance ，high-porosity and low permeability．The geological formation of low resistivity reservoir has many reasons．In view of the developing features of Z oil group ，the clay particles distribution ，clay content and the influence of them on the low-resistivity are studied．clay minerals types ，

It can be found that the clay types in low resistivity layer segment of Z1oil group is mainly the mixed layer of illite and montomorillonite and the content is relatively high．The clay particles are evenly distributed in a thin film and in flocculating form with high average clay content．Additional conductivity caused by clay minerals and high clay content are the main factors affecting the low resistivity of Z1oil group．

Key words ：low-resistivity reservoir ；clay minerals ；additional conductivity ；clay content ；genesis of low resistivity

STUDY ON SPECTRAL ANALYSIS OF NATURAL GAMMA RAY SPECTRA LOG

LIU Zun-nian 1，REN Ai-ge 2，JIA Rui-gao 3（1．School of Science ，Qingdao Technological University ，Qingdao Shan-China ；2．Department of Basic Courses ，Qingdao Bin Hai University ，Qingdao Shandong 266555，Chi-dong 266033，

na ；3．College of Physics Science and Technology ，China University of Petroleum ，Dongying Shandong 257062，Chi-VOL．32，na ）JOURNAL OF SOUTHWEST PETROLEUM UNIVERSITY （SCIENCE ＆TECHNOLOGY EDITION ），NO．5，75－78，2010（ISSN 1674－5086，in Chinese ）

Abstract ：Natural Gamma Ray Spectra log （SGR ）is an important nuclear logging method ，and real spectral analysis is a key technique of the surface processing system in SGR．Due to the inherent low resolution of Natural Gamma there is a low signal-to-noise ratio and low confidence level during the logging process．Based on Ray Spectra log ，

the principle of Natural Gamma Ray Spectra log ，the normal count rate spectra are obtained from the test pit ，and the the authors perform preliminary research on applying the four traditional algorithm methods are tested．Besides ，

direct demodulation method which has a high sensitivity and spatial resolution in gamma-ray log．This research of-fers a credible basis for algorithm selecting．

Key words ：Natural Gamma Ray spectrometry log ；spectrum unscrambling ；direct demodulation method ；response function ；physical constrains

GENETIC MECHANISM OF LOW RESISTIVITY RESERVOIR IN BANQIAO OIL FIELD

YE Xing-shu ，WU Guo-hai （GreatWall Drilling Company ，CNPC ，Chaoyang Beijing 100101，China ）JOURNAL OF ）5