CN103291279B - 一种气测值微异常信号的优化方法 - Google Patents

一种气测值微异常信号的优化方法 Download PDF

Info

Publication number
CN103291279B
CN103291279B CN201210054262.7A CN201210054262A CN103291279B CN 103291279 B CN103291279 B CN 103291279B CN 201210054262 A CN201210054262 A CN 201210054262A CN 103291279 B CN103291279 B CN 103291279B
Authority
CN
China
Prior art keywords
value
gas logging
gas
well depth
logging value
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
CN201210054262.7A
Other languages
English (en)
Other versions
CN103291279A (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.)
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering
Original Assignee
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering
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 China Petroleum and Chemical Corp, Sinopec Research Institute of Petroleum Engineering filed Critical China Petroleum and Chemical Corp
Priority to CN201210054262.7A priority Critical patent/CN103291279B/zh
Publication of CN103291279A publication Critical patent/CN103291279A/zh
Application granted granted Critical
Publication of CN103291279B publication Critical patent/CN103291279B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明属于石油勘探开发领域,尤其涉及一种对油气井的气测值进行处理优化的方法。方法包括,对待录井的各米探测得到气测值,形成原始气测曲线A;当井深为i米时,其气测值为Ai,取i米上下临近的井深的气测值组[Ai-l,Ai+l],且l取值范围为2-10米,后对气测值组[Ai-l,Ai+l]求平均值Bi;并根据平均值和气测值得到该井深点气测值Ai的方差σi;比较Ai与(Bii)的大小,当Ai≥(Bii),该井深为i点的气测值为高异常信号,输出该点信息;对井深各米依次进行上述优化过程,直至结束,输出所有高异常信号点信息。本发明解决了气测值微异常信号的识别和评价问题,有利于地下油气显示的发现和评价。

Description

一种气测值微异常信号的优化方法
技术领域
本发明属于石油勘探开发领域,尤其涉及一种对油气井的气测值进行处理优化的方法。
背景技术
在石油勘探综合录井中,借助于综合录井设备,可连续测量钻井液内烃类气体的含量,用来发现、解释、评价油气显示。一般来说,气测值相对平稳的井段,没有油气显示,如果出现了波峰-气测值高出基线,则意味着出现了油气显示,如图1所示,图中1011-1014米井段出现了波峰,则该井段很可能是油气层。目前,一般用肉眼对曲线进行分析、识别。
但在实际录井工作中,受各种因素影响,当地下出现油气显示时,气测异常并不明显,如图2中,只出现了几处微小的波峰-气测值稍稍高出基线,几乎无法进行人工识别,给发现油气显示带来了很大困难。在现有技术中,一般采用两种方法:
①用肉眼进行仔细分辨,用该方法时,增加了工作量,且人为因素很大,往往漏掉油气显示。
②简单地将所有气测值乘以一个倍数,这种方法将所有的气测值都放大或缩小了,某微小波峰与上下比较时,仍然没有明显的异常,基本起不到效果。
理论上,当钻头钻遇油气显示时,气测值会出现一个明显的波峰,但是,在各种因素干扰下,该波峰往往并不明显。综上所述,我们需要一种方法,利用科学的手段、将微小的波峰识别出来,并利用定量化的方法,并将波峰信号进行放大,形成处理后气测曲线,有利于油气显示的发现及评价。目前未有相关技术内容的披露。
现有技术中,谢元军在《随钻气测录井全烃值标准化方法探讨》一文中,(录井工程,2007,18(4):13-6)提出了对气测值进行处理。该文分析了钻时、钻井液流量等参数对气测值的影响,用反比例函数法、指数函数法、反三角函数法对气测值进行了整体处理,但没有对异常值进行单独处理,没有解决微小波峰的识别问题。
发明内容
本发明为了识别出这些微小波峰,并将微小波峰信号放大,解决现有技术中的问题,研发了一种气测值微异常信号的优化方法,本发明解决了气测值微异常信号的识别和评价问题,有利于地下油气显示的发现和评价。
本发明所采用的方法包括,
一种气测值微异常信号的优化方法,所述优化方法包括,对待录井的各米探测得到气测值,形成原始气测曲线A;当井深为i米时,其气测值为Ai,取i米上下临近的井深的气测值组[Ai-l,Ai+l],且l取值范围为2-10米;后对气测值组[Ai-l,Ai+l]求平均值Bi;并根据平均值和气测值得到该井深点气测值Ai的方差σi;比较Ai与(Bii)的大小,当Ai≥(Bii),该井深为i点的气测值为高异常信号,输出该点信息;对井深各米依次进行上述优化过程,直至结束,输出所有高异常信号点信息。
在具体的实践中,除了迅速输出高异常信号外,对于输出具有显著特征的新气测值曲线也具有突出的实质性效果,因此本方法还包括输出新气测曲线步骤;步骤包括,根据气测值Ai,平均值Bi和方差σi;比较Ai与(Bii)的大小,当Ai≤(Bii),则将Ai乘以小于1的系数k1,得到低异常点值记为D1;比较Ai与(Bii)的大小,当Ai≥(Bii),则将Ai乘以大于1的系数k2,得到高异常点值记为Ei;当不满足Ai≤(Bii)条件,也不满足Ai≥(Bii)条件时,则将Ai乘以系数1,得到正常点值Fi;对待录井各个点均优化完后,形成低异常点D、高异常点E和正常点F三个数据组,连接D,E,F组形成新气测曲线G;输出优化后的新气测曲线。
应用中,优化步骤包括,
(1)得到待录井各个物理参数,包括井深;
(2)探测得到待录井的气测值曲线:每米井深测量一个气测值,记为Ai,所有的Ai值形成原始气测曲线A;
(3)取井深i处,若i≥井深,则跳转至步骤(9),否则进入步骤(4);
(4)提取Ai处上下临近处气测值:在井深i,提取Ai,并提取井深i上下各l米,共计2l+1米个气测值数,形成[Ai-l,Ai+l];l取值范围为2-10米
(5)求取平均值:对[Ai-l,Ai+l]求取平均值为Bi
(6)求取方差值:对[Ai-l,Ai+l]求取方差值为σi
其中, σ i = 1 l Σ i = 1 l ( A i - B i ) 2
比较步骤:
(71)标记低异常点值:比较Ai与(Bii)的大小,若Ai≤(Bii),则将Ai乘以小于1的系数k1,得到低异常值记为Di
(72)标记高异常点值:比较Ai与(Bii)的大小,若Ai≥(Bii),则将Ai乘以大于1的系数k2,得到高异常值记为Ei
(73)标记正常点值:即不满足Ai≤(Bii)条件,也不满足Ai≥(Bii)条件的,则将Ai乘以系数1,得到正常值记为Fi
(7)井深i+1,判断i+1是否大于井深数值;是则进入步骤(9),否则转入步骤(4);
(8)形成新的气测线:低异常点D、高异常点E和正常点F三个数据组,连接D,E,F组形成新气测曲线G;
(9)输出优化后的新气测曲线G。
为了实现更优的效果,所述步骤(4)中的l取值范围为3-5米;所述步骤(71)中的k1取值范围为0.95-0.98;所述步骤(72)中的k2取值范围为1.10-1.20。
本发明在输出了高异常点进行判断后,又对新气测值曲线进行了刻画,直观而清晰的实现了发明效果,同时,本发明为了验证所述优化步骤的正确性,本优化方法还包括校验步骤:
校验过程为:判断数值,若则Gi为高异常点,应与输出新气测值曲线G一致,若不一致,则重新优化步骤。
本发明成功地将气测值中的微小的波峰-气测值稍稍高出基线识别出来,并将该微小波峰信号放大,可有效解决实际生产中的实际问题。
①有利于发现弱油气显示;
②目前,国内外气测录井仪器型号多样,同时受到不同的钻井条件、钻井液性能、检测设备等差异的影响,即使对同一区块油气性质、地层压力和物性完全相同的油气层,气测结果也会有很大不同。某些本来很好的油气显示,受外界条件影响,也会造成气测值异常值很低,只能出现微小的波峰,即通常所说的弱油气显示,造成油气层发现困难。
③有利于评价弱油气显示
识别出弱油气显示后,根据处理后气测值的大小,可更加准确地评价该显示,有效提高油气层解释符合率,可大大提高录井资料的使用率及勘探开发总体效益。
④有利于录井资料的横向对比
受各种钻井参数、钻井液性能、检测设备等影响,气测值低的油气显示并不一定必气测值高的油气显示差,不同井之间的资料无法对比。经过本发明处理后的气测值,能更客观地反应地层含油气情况,实现不同井之间的资料对比,为下步施工提供科学的依据。
附图说明
图1为某井气测值图-明显波峰示意图;
图2为某井气测值图-不明显波峰示意图;
图3为本发明的流程图;
图4为本发明实施例的气测值效果图。
结合背景技术和具体实施方式加以说明
具体实施方式
本发明具体优化步骤包括,
(1)得到待录井各个物理参数,包括井深;
(2)探测得到待录井的气测值曲线:每米井深测量一个气测值,记为Ai,所有的Ai值形成原始气测曲线A;
(3)取井深i处,若i≥井深,则跳转至步骤(9),否则进入步骤(4);
(4)提取Ai处上下临近处气测值:在井深i,提取Ai,并提取井深i上下各l米,共计2l+1米个气测值数,形成[Ai-l,Ai+l];l取值范围为2-10米
(5)求取平均值:对[Ai-l,Ai+l]求取平均值为Bi
(6)求取方差值:对[Ai-l,Ai+l]求取方差值为σi
其中, σ i = 1 l Σ i = 1 l ( A i - B i ) 2
比较步骤:
(71)标记低异常点值:比较Ai与(Bii)的大小,若Ai≤(Bii),则将Ai乘以小于1的系数k1,得到低异常值记为Di
(72)标记高异常点值:比较Ai与(Bii)的大小,若Ai≥(Bii),则将Ai乘以大于1的系数k2,得到高异常值记为Ei
(73)标记正常点值:即不满足Ai≤(Bii)条件,也不满足Ai≥(Bii)条件的,则将Ai乘以系数1,得到正常值记为Fi
(7)井深i+1,判断i+1是否大于井深数值;是则进入步骤(9),否则转入步骤(4);
(8)形成新的气测线:低异常点D、高异常点E和正常点F三个数据组,连接D,E,F组形成新气测曲线G;
(9)输出优化后的新气测曲线G。
为了实现更优的效果,所述步骤(4)中的l取值范围为3-5米;所述步骤(71)中的k1取值范围为0.95-0.98;所述步骤(72)中的k2取值范围为1.10-1.20。
本发明在输出了高异常点进行判断后,又对新气测值曲线进行了刻画,直观而清晰的实现了发明效果,同时,本发明为了验证所述优化步骤的正确性,本优化方法还包括校验步骤:
校验过程为:判断数值,若则Gi为高异常点,应与输出新气测值曲线G一致,若不一致,则重新优化步骤。
实施例:
本试验在某油田M2、W3、L8三口井进行了试验,M2井试验结果如表1、图4,由表1、图4可看出,处理前的气测数据及气测曲线,相对比较平稳,在1014米、1022米出现了极其微弱的波峰,经过本发明科学的手段、定量化的挑选和放大方法,将1014米、1022米的气测值自动挑选出来,并进行了放大(l取3)。
在1014米, G i A i k 2 = 4.8532 4.4120 * 1.10 = 1.0
在1022米, G i A i k 2 = 5.1185 4.6532 * 1.10 = 1.0
表1气测值处理数据表
放大后的气测值,可轻易地识别出来。后经测试证实,1014米是油层,日产油12.6立方米;1022米是油水同层,日产油5.8立方米、产水8.2立方米。

Claims (5)

1.一种气测值微异常信号的优化方法,其特征在于,所述优化方法包括,对待录井的各米探测得到气测值,形成原始气测曲线A;当井深为i米时,其气测值为Ai,取i米上下临近的井深气测值组[Ai-l,Ai+l],且l取值范围为2-10米;后对气测值组[Ai-l,Ai+l]求平均值Bi;并根据平均值和气测值得到该井深点气测值Ai的方差σi;比较Ai与(Bii)的大小,当Ai≥(Bii),该井深为i点的气测值为高异常信号,输出该点信息;对井深各米依次进行上述优化过程,直至结束,输出所有高异常信号点信息。
2.根据权利要求1所述的一种气测值微异常信号的优化方法,其特征在于,所述方法还包括输出新气测曲线步骤:根据气测值Ai,平均值Bi和方差σi;比较Ai与(Bii)的大小,当Ai≤(Bii),则将Ai乘以小于1的系数k1,得到低异常点值,记为Di;比较Ai与(Bii)的大小,当Ai≥(Bii),则将Ai乘以大于1的系数k2,得到高异常点值记为Ei;当不满足Ai≤(Bii)条件,也不满足Ai≥(Bii)条件时,则将Ai乘以系数1,得到正常点值Fi;对待录井各个点均优化完后,形成低异常点D、高异常点E和正常点F三个数据组,连接D,E,F组形成新气测曲线G;输出优化后的新气测值曲线。
3.根据权利要求1或2所述的一种气测值微异常信号的优化方法,其特征在于,优化步骤,
(1)得到待录井各个物理参数,包括井深;
(2)探测得到待录井的气测值曲线:每米井深测量一个气测值,记为Ai,所有的Ai值形成原始气测曲线A;
(3)取井深i处,若i≥井深,则跳转至步骤(9),否则进入步骤(4);
(4)提取Ai处上下临近处气测值:在井深i,提取Ai,并提取井深i上下各l米,共计2l+1米个气测值数,形成[Ai-l,Ai+l];l取值范围为2-10米;
(5)求取平均值:对[Ai-l,Ai+l]求取平均值为Bi
(6)求取方差值:对[Ai-l,Ai+l]求取方差值为σi
其中, σ i = 1 l Σ i = 1 l ( A i - B i ) 2
比较步骤:
(71)标记低异常点值:比较Ai与(Bii)的大小,若Ai≤(Bii),则将Ai乘以小于1的系数k1,得到低异常值记为Di
(72)标记高异常点值:比较Ai与(Bii)的大小,若Ai≥(Bii),则将Ai乘以大于1的系数k2,得到高异常值记为Ei
(73)标记正常点值:即不满足Ai≤(Bii)条件,也不满足Ai≥(Bii)条件的,则将Ai乘以系数1,得到正常值记为Fi
(7)井深i+1,判断i+1是否大于井深数值;是则进入步骤(9),否则转入步骤(4);
(8)形成新的气测线:低异常点D、高异常点E和正常点F三个数据组,连接D,E,F组形成新气测值曲线G;
(9)输出优化后的新气测值曲线G。
4.根据权利要求3所述的一种气测值微异常信号的优化方法,其特征在于,所述步骤(4)中的l取值范围为3-5米;
所述步骤(71)中的k1取值范围为0.95-0.98;
所述步骤(72)中的k2取值范围为1.10-1.20。
5.根据权利要求3所述的一种气测值微异常信号的优化方法,其特征在于,所述优化方法还包括校验步骤:
校验过程为:判断数值,若则Gi为高异常点,应与输出新气测值曲线G一致,若不一致,则重新优化步骤。
CN201210054262.7A 2012-03-02 2012-03-02 一种气测值微异常信号的优化方法 Active CN103291279B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210054262.7A CN103291279B (zh) 2012-03-02 2012-03-02 一种气测值微异常信号的优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210054262.7A CN103291279B (zh) 2012-03-02 2012-03-02 一种气测值微异常信号的优化方法

Publications (2)

Publication Number Publication Date
CN103291279A CN103291279A (zh) 2013-09-11
CN103291279B true CN103291279B (zh) 2016-01-13

Family

ID=49092742

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210054262.7A Active CN103291279B (zh) 2012-03-02 2012-03-02 一种气测值微异常信号的优化方法

Country Status (1)

Country Link
CN (1) CN103291279B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105466927B (zh) * 2014-07-02 2020-04-10 深圳迈瑞生物医疗电子股份有限公司 一种比浊法异常反应曲线的识别、修正及报警方法
CN108898117A (zh) * 2018-06-30 2018-11-27 防灾科技学院 一种滑动阈值的自适应随机信号异常提取方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3821472B2 (ja) * 2002-01-29 2006-09-13 松下電器産業株式会社 異常波形検出回路および情報再生装置
CN101369152A (zh) * 2008-10-17 2009-02-18 中国安全生产科学研究院 油气开采集输作业安全监控预警及安全管理系统及其方法
CN101696627A (zh) * 2009-10-30 2010-04-21 湖北工业大学 一种钻井过程事故监测方法
CN101976494B (zh) * 2010-08-26 2013-06-26 中国石油集团川庆钻探工程有限公司 录井实时数据采集及同步预警方法
CN102312671B (zh) * 2011-09-16 2014-02-12 中国石油化工股份有限公司 能快速解释评价储层流体性质的方法

Also Published As

Publication number Publication date
CN103291279A (zh) 2013-09-11

Similar Documents

Publication Publication Date Title
CN103728655B (zh) 一种回采工作面冲击危险性采前预评价方法
Tokhmchi et al. Estimation of the fracture density in fractured zones using petrophysical logs
CN105160414B (zh) 预测全油藏产区类型的方法及装置
Karacan et al. Inference of strata separation and gas emission paths in longwall overburden using continuous wavelet transform of well logs and geostatistical simulation
He et al. Pressure-transient behavior of multisegment horizontal wells with nonuniform production: theory and case study
Clerke et al. Wireline spectral porosity analysis of the Arab Limestone—From Rosetta Stone to CIPHER
CN104899358A (zh) 奥灰岩溶裂隙水网络横向分布的预测方法
CN111502650B (zh) 一种利用气测衍生参数识别凝析气层的方法及其应用
Goodwin et al. Improved water use estimates for drilling and hydrualic fracturing in Northeastern Colorado
CN103291279B (zh) 一种气测值微异常信号的优化方法
Kharrat et al. Rock Typeand Permeability Prediction of a Heterogeneous Carbonate Reservoir Using Artificial Neural Networks Based on Flow Zone Index Approach
CN105064987A (zh) 利用随钻录井q参数进行水层识别的解释评价方法
CN104514553B (zh) 一种岩屑砂岩层中气层的识别方法及其应用
CN103422852B (zh) 一种不同井间气测值转换对比方法
Zhang et al. Study on Pressure Characteristics of Disturbed Wells due to Interwell Fracturing Interference and Its Application in Small Well Spacing Fracturing
Qin et al. Diagnosis of Water-Influx Locations of Horizontal Well Subject to Bottom-Water Drive through Well-Testing Analysis
CN111827994B (zh) 利用气测全烃百分含量和烃组分湿度比解释砂岩油水层的方法
Wei Well test pressure derivatives and the nature of fracture networks
Dastgerdi et al. Estimation of underground interwell connectivity: A data-driven technology
CN108661629A (zh) 一种针对页岩地层的工程甜点定量评价方法
Onur et al. Application of modern well test analysis techniques to pressure transient tests in Kizildere geothermal field, Turkey
CN112855131A (zh) 一种页岩油储层含油率预测方法
Zhu et al. Dipole flow tracer test models for identifying skin properties and ambient groundwater flow effects
Akin Design and analysis of multi-well interference tests
Fang et al. An intelligent production fluctuation monitoring system for giant oilfield development

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