CN103291279B - 一种气测值微异常信号的优化方法 - Google Patents
一种气测值微异常信号的优化方法 Download PDFInfo
- 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
Links
Abstract
本发明属于石油勘探开发领域,尤其涉及一种对油气井的气测值进行处理优化的方法。方法包括,对待录井的各米探测得到气测值,形成原始气测曲线A;当井深为i米时,其气测值为Ai,取i米上下临近的井深的气测值组[Ai-l,Ai+l],且l取值范围为2-10米,后对气测值组[Ai-l,Ai+l]求平均值Bi;并根据平均值和气测值得到该井深点气测值Ai的方差σi;比较Ai与(Bi+σi)的大小,当Ai≥(Bi+σi),该井深为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与(Bi+σi)的大小,当Ai≥(Bi+σi),该井深为i点的气测值为高异常信号,输出该点信息;对井深各米依次进行上述优化过程,直至结束,输出所有高异常信号点信息。
在具体的实践中,除了迅速输出高异常信号外,对于输出具有显著特征的新气测值曲线也具有突出的实质性效果,因此本方法还包括输出新气测曲线步骤;步骤包括,根据气测值Ai,平均值Bi和方差σi;比较Ai与(Bi-σi)的大小,当Ai≤(Bi-σi),则将Ai乘以小于1的系数k1,得到低异常点值记为D1;比较Ai与(Bi+σi)的大小,当Ai≥(Bi+σi),则将Ai乘以大于1的系数k2,得到高异常点值记为Ei;当不满足Ai≤(Bi-σi)条件,也不满足Ai≥(Bi+σi)条件时,则将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;
其中,
比较步骤:
(71)标记低异常点值:比较Ai与(Bi-σi)的大小,若Ai≤(Bi-σi),则将Ai乘以小于1的系数k1,得到低异常值记为Di;
(72)标记高异常点值:比较Ai与(Bi+σi)的大小,若Ai≥(Bi+σi),则将Ai乘以大于1的系数k2,得到高异常值记为Ei;
(73)标记正常点值:即不满足Ai≤(Bi-σi)条件,也不满足Ai≥(Bi+σi)条件的,则将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;
其中,
比较步骤:
(71)标记低异常点值:比较Ai与(Bi-σi)的大小,若Ai≤(Bi-σi),则将Ai乘以小于1的系数k1,得到低异常值记为Di;
(72)标记高异常点值:比较Ai与(Bi+σi)的大小,若Ai≥(Bi+σi),则将Ai乘以大于1的系数k2,得到高异常值记为Ei;
(73)标记正常点值:即不满足Ai≤(Bi-σi)条件,也不满足Ai≥(Bi+σi)条件的,则将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米,
在1022米,
表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与(Bi+σi)的大小,当Ai≥(Bi+σi),该井深为i点的气测值为高异常信号,输出该点信息;对井深各米依次进行上述优化过程,直至结束,输出所有高异常信号点信息。
2.根据权利要求1所述的一种气测值微异常信号的优化方法,其特征在于,所述方法还包括输出新气测曲线步骤:根据气测值Ai,平均值Bi和方差σi;比较Ai与(Bi-σi)的大小,当Ai≤(Bi-σi),则将Ai乘以小于1的系数k1,得到低异常点值,记为Di;比较Ai与(Bi+σi)的大小,当Ai≥(Bi+σi),则将Ai乘以大于1的系数k2,得到高异常点值记为Ei;当不满足Ai≤(Bi-σi)条件,也不满足Ai≥(Bi+σi)条件时,则将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;
其中,
比较步骤:
(71)标记低异常点值:比较Ai与(Bi-σi)的大小,若Ai≤(Bi-σi),则将Ai乘以小于1的系数k1,得到低异常值记为Di;
(72)标记高异常点值:比较Ai与(Bi+σi)的大小,若Ai≥(Bi+σi),则将Ai乘以大于1的系数k2,得到高异常值记为Ei;
(73)标记正常点值:即不满足Ai≤(Bi-σi)条件,也不满足Ai≥(Bi+σi)条件的,则将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一致,若不一致,则重新优化步骤。
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)
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)
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 | 中国石油化工股份有限公司 | 能快速解释评价储层流体性质的方法 |
-
2012
- 2012-03-02 CN CN201210054262.7A patent/CN103291279B/zh active Active
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 |