CN108756867B - 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 - Google Patents
基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 Download PDFInfo
- Publication number
- CN108756867B CN108756867B CN201810447124.2A CN201810447124A CN108756867B CN 108756867 B CN108756867 B CN 108756867B CN 201810447124 A CN201810447124 A CN 201810447124A CN 108756867 B CN108756867 B CN 108756867B
- Authority
- CN
- China
- Prior art keywords
- curve
- frequency
- acoustic
- resistivity
- low
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims abstract description 5
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 5
- 238000013507 mapping Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 239000011435 rock Substances 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 150000002430 hydrocarbons Chemical class 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 239000004215 Carbon black (E152) Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 238000003556 assay Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于声波测井曲线和电阻率测井曲线进行压裂选层的方法,通过声波测井获得声波曲线,通过电阻率测井获得电阻率曲线,包括如下步骤:步骤1:对声波曲线和电阻率曲线分别进行频谱分解,获得声波低频分量、声波高频分量、电阻率低频分量、电阻率高频分量;步骤2:将低频声波曲线和高频电阻率曲线合成广义声波曲线;步骤3:将步骤2合成的广义声波曲线与声波曲线作比,生成指示曲线;步骤4:根据指示曲线选取压裂层段。
Description
技术领域
本发明涉及一种基于声波测井曲线和电阻率测井曲线进行压裂选层的方法。
背景技术
储层改造是页岩油增产的重要措施,油井的压裂选层是储层改造中的重要技术手段。常规压裂选层主要有两类方法:
第一类方法,成熟探区,统计已有井压裂效果与压裂效果影响参数的规律,对目的井优选压裂层段。成熟探区,选择压裂影响参数,一般有目前产量、原始和目前的油藏压力、井的位置、射孔信息、完井信息等。利用数理统计的方法,找出压后产能与影响参数的规律,找出最主要的影响因素,指导目的井的压裂选层。此种方法存在的不足之处在于,需要压裂后的参数,对初探区或者压裂井较少的探区不适用。优选出的压裂效果影响参数一般为组段,厚度量级一般为几十米,而压裂施工段级一般为几米,精度不够。
第二类方法,对初探区缺少压裂井资料的情况下,优选与产能相关的烃源岩优势参数作为压裂选层依据。初探区,对目的层段烃源岩取心,进行实验室有机地化参数分析,分析与产量密切相关的烃源岩参数:有机碳含量(TOC)、生烃量(S1+S2)、孔隙度、渗透率、脆性矿物含量等,优选较好的层段进行压裂。此种方法存在的不足之处在于钻井取心,分析化验成本高。
发明内容
本发明目的在于提供一种基于声波测井曲线和电阻率测井曲线进行压裂选层的方法,对成熟探区和初探区均适用,且检测方便,成本低。
实现本发明目的的技术方案:
一种基于声波测井曲线和电阻率测井曲线进行压裂选层的方法,通过声波测井获得声波曲线,通过电阻率测井获得电阻率曲线,其特征在于,包括如下步骤:
步骤1:对声波曲线和电阻率曲线分别进行频谱分解,获得声波低频分量、声波高频分量、电阻率低频分量、电阻率高频分量;
步骤2:将低频声波曲线和高频电阻率曲线合成广义声波曲线;
步骤3:将步骤2合成的广义声波曲线与声波曲线作比,生成指示曲线;
步骤4:根据指示曲线选取压裂层段。
进一步地,步骤1中,将声波曲线和电阻率曲线的深度域信号作为时间序列通过傅里叶变换得到频率域信号,利用滤波器获得曲线的低频分量、高频分量。
进一步地,步骤1中,通过均值方差法对电阻率高频分量和声波高频分量进行标准化校正。
进一步地,步骤2中,分别将频率域的声波低频分量和校正的电阻率高频分量通过傅里叶逆变换将频率映射到时域,在时域内将声波低频分量与电阻率高频分量相加得到广义声波曲线。
进一步地,步骤4中,指示曲线值大于1.15时,对应油层;指示曲线值大于0.97小于1.15时,对应油水层。
进一步地,声波曲线采用声速测井曲线
本发明具有的有益效果:
本发明对声波曲线和电阻率曲线分别进行频谱分解,获得声波低频分量、声波高频分量、电阻率低频分量、电阻率高频分量;将低频声波曲线和高频电阻率曲线合成广义声波曲线;将广义声波曲线与声波曲线作比,生成指示曲线;根据指示曲线选取压裂层段。本发明基于常规的声波曲线和电阻率测井曲线,数据容易获取,不需要成熟探区的压裂数据,不需要对初探区进行钻井取心,分析化验,本发明可广泛适用于成熟探区和初探区,在保证检测精度的前提下,检测更加方便,有效降低检测成本。本发明对电阻率高频分量和声波高频分量进行标准化校正,再生成指示曲线,能够进一步保证检测的准确率。
附图说明
图1是本发明方法流程图;
图2是本发明方法用matlab编程后运行结果图;
图3是本发明所述的声波曲线和电阻率曲线的高低频分量、广义声波曲线、指示曲线与试油层的对应关系图;
图4是本发明指示曲线与烃源岩分析化验烃源岩(TOC)、生烃量(S1、S2)、孔隙度、渗透率、脆性指数等参数的对应关系图;
图5是本发明指示曲线与全井试油层段的对应关系图;
图6是本发明指示曲线与试油层段的对应关系局部放大图;
图7是本发明指示曲线与试油层段的对应关系图。
具体实施方式
如图1所示,本发明通过声波测井获得声波曲线,通过电阻率测井获得电阻率曲线,包括如下步骤:
步骤1:对声波曲线和电阻率曲线分别进行频谱分解,获得声波低频分量、声波高频分量、电阻率低频分量、电阻率高频分量;
步骤2:将低频声波曲线和高频电阻率曲线合成广义声波曲线;
步骤3:将步骤2合成的广义声波曲线与声波曲线作比,生成指示曲线,即指示曲线=广义声波曲线/声波曲线;
其中,声波(时差)曲线(DT)和电阻率曲线(RT)均通过常规的测井方法获得。声波测井是利用声波在不同岩石的中传播时,速度、幅度及频率的变化等声学特性不相同来研究钻井的地质剖面,判断固井质量的一种测井方法。将一个受控声波振源放入井中,声源发出的声波引起周围质点的振动,在地层中产生体波即纵波和横波,在井壁一钻井液界面上产生诱导的界面波即伪瑞利波和斯通莱波。这些波作为地层信息的载体,被井下接收器接收,送至地面的记录下来,就是声波测井。声波在井内地层中传播由于地层岩石成分、结构、孔隙中流体成分的变化其波的速度、幅度甚至频率都会发生变化。只记录声波速度变化的称为声速测井(AC),而记录声幅度变化的则称为声幅测井。实施时,选用声速测井曲线(AC),纵坐标参数us/m,即微秒每米。电阻率测井是在钻孔中采用布置在不同部位的供电电极和测量电极来测定岩石(包括其中的流体)电阻率的方法,纵坐标参数:Ω·m或ohmm欧姆米。
步骤1中,声波曲线和电阻率曲线的深度域信号通过傅里叶变换得到频率域信号,利用滤波器获得曲线的低频分量、高频分量;通过均值方差法对电阻率高频分量和声波高频分量进行标准化校正。
步骤2中,分别将频率域的声波低频分量和校正的电阻率高频分量通过傅里叶逆变换将频率映射到时域,在时域内将声波低频分量与电阻率高频分量相加得到广义声波曲线。
步骤1和步骤2中涉及的公式如下:
傅里叶变换:
傅里叶逆变换:
式中:F(w)为频率域函数,f(t)为深度域函数,w为角速度。
滤波器的频谱函数H(f)表达式如下:
式中:fL为滤波器的截止频率;f为频率。
步骤4中,如图7所示,指示曲线值大于1.15时,对应油层;指示曲线值大于0.97小于1.15时,对应油水层。本发明检测精度达到0.125m,满足压裂施工要求。
Claims (4)
1.一种基于声波测井曲线和电阻率测井曲线进行压裂选层的方法,通过声波测井获得声波曲线,通过电阻率测井获得电阻率曲线,其特征在于,包括如下步骤:
步骤1:对声波曲线和电阻率曲线分别进行频谱分解,获得声波低频分量、声波高频分量、电阻率低频分量、电阻率高频分量;
通过均值方差法对电阻率高频分量和声波高频分量进行标准化校正;
步骤2:将低频声波曲线和高频电阻率曲线合成广义声波曲线;
步骤3:将步骤2合成的广义声波曲线与声波曲线作比,生成指示曲线;
步骤4:根据指示曲线选取压裂层段;
步骤4中,指示曲线值大于1.15时,对应油层;指示曲线值大于0.97小于1.15时,对应油水层。
2.根据权利要求1所述的方法,其特征在于:步骤1中,声波曲线和电阻率曲线的深度域信号作为时间序列通过傅里叶变换得到频率域信号,利用滤波器获得曲线的低频分量、高频分量。
3.根据权利要求1所述的方法,其特征在于:步骤2中,分别将频率域的声波低频分量和校正的电阻率高频分量通过傅里叶逆变换将频率映射到时域,在时域内将声波低频分量与电阻率高频分量相加得到广义声波曲线。
4.根据权利要求1至3任何一项所述的方法,其特征在于:声波曲线采用声速测井曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810447124.2A CN108756867B (zh) | 2018-05-11 | 2018-05-11 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810447124.2A CN108756867B (zh) | 2018-05-11 | 2018-05-11 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108756867A CN108756867A (zh) | 2018-11-06 |
CN108756867B true CN108756867B (zh) | 2021-11-19 |
Family
ID=64009731
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810447124.2A Expired - Fee Related CN108756867B (zh) | 2018-05-11 | 2018-05-11 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108756867B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110346416B (zh) * | 2019-07-17 | 2022-04-15 | 北京金海能达科技有限公司 | 基于声波和电阻率的特征参数曲线识别火山岩岩性的方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101915091A (zh) * | 2010-06-22 | 2010-12-15 | 中国石油天然气股份有限公司 | 一种测井频率匹配薄层校正方法与设备 |
CN102353985A (zh) * | 2011-06-01 | 2012-02-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于非下采样Contourlet变换的拟声波曲线构建方法 |
CN103389513A (zh) * | 2013-08-02 | 2013-11-13 | 中国石油集团川庆钻探工程有限公司 | 应用声波测井资料约束反演提高地震资料分辨率的方法 |
CN103485768A (zh) * | 2012-06-13 | 2014-01-01 | 中国石油天然气集团公司 | 声波测井曲线的构建方法 |
CN103850679A (zh) * | 2014-04-01 | 2014-06-11 | 北京师范大学 | 一种利用多种测井曲线对声波时差曲线进行重构的方法 |
CN104047597A (zh) * | 2014-01-26 | 2014-09-17 | 中国石油大学(华东) | 含油气泥页岩地层测井曲线标准化方法 |
US8848484B2 (en) * | 2010-12-08 | 2014-09-30 | Schlumberger Technology Corporation | Filtering acoustic waveforms in downhole environments |
CN105089652A (zh) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | 一种拟声波曲线重构与稀疏脉冲联合反演方法 |
CN105223608A (zh) * | 2015-08-07 | 2016-01-06 | 中国石油大学(华东) | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 |
CN106405642A (zh) * | 2016-07-19 | 2017-02-15 | 西安石油大学 | 基于脱压实声波速度的地震反演储层预测方法 |
CN107065011A (zh) * | 2017-06-22 | 2017-08-18 | 东北石油大学 | 一种应用于陆相盆地储层反演的曲线频率融合方法 |
-
2018
- 2018-05-11 CN CN201810447124.2A patent/CN108756867B/zh not_active Expired - Fee Related
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101915091A (zh) * | 2010-06-22 | 2010-12-15 | 中国石油天然气股份有限公司 | 一种测井频率匹配薄层校正方法与设备 |
US8848484B2 (en) * | 2010-12-08 | 2014-09-30 | Schlumberger Technology Corporation | Filtering acoustic waveforms in downhole environments |
CN102353985A (zh) * | 2011-06-01 | 2012-02-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于非下采样Contourlet变换的拟声波曲线构建方法 |
CN103485768A (zh) * | 2012-06-13 | 2014-01-01 | 中国石油天然气集团公司 | 声波测井曲线的构建方法 |
CN103389513A (zh) * | 2013-08-02 | 2013-11-13 | 中国石油集团川庆钻探工程有限公司 | 应用声波测井资料约束反演提高地震资料分辨率的方法 |
CN103389513B (zh) * | 2013-08-02 | 2016-05-04 | 中国石油集团川庆钻探工程有限公司 | 应用声波测井资料约束反演提高地震资料分辨率的方法 |
CN104047597A (zh) * | 2014-01-26 | 2014-09-17 | 中国石油大学(华东) | 含油气泥页岩地层测井曲线标准化方法 |
CN103850679A (zh) * | 2014-04-01 | 2014-06-11 | 北京师范大学 | 一种利用多种测井曲线对声波时差曲线进行重构的方法 |
CN105089652A (zh) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | 一种拟声波曲线重构与稀疏脉冲联合反演方法 |
CN105223608A (zh) * | 2015-08-07 | 2016-01-06 | 中国石油大学(华东) | 一种含煤强屏蔽缝洞型储层的地震预测与描述方法 |
CN106405642A (zh) * | 2016-07-19 | 2017-02-15 | 西安石油大学 | 基于脱压实声波速度的地震反演储层预测方法 |
CN107065011A (zh) * | 2017-06-22 | 2017-08-18 | 东北石油大学 | 一种应用于陆相盆地储层反演的曲线频率融合方法 |
Non-Patent Citations (7)
Title |
---|
《Application of well logs integration and wavelet transform to improve fracture zones detection in metamorphic rocks》;Huaijie Yang,etc;《Journal of Petroleum Science and Engineering》;20170721;全文 * |
《Identifying Reservoir Fluids by Wavelet Transform of Well Logs》;Wenzheng Yue, Guo Tao, Zhengwu Liu;《Society of Petroleum Engineers (SPE)》;20041020;全文 * |
《东营北带深层砂砾岩测井层序对比与岩性识别研究》;王鑫;《中国优秀硕士学位论文全文数据库基础科学辑》;20110615(第6期);第14-22页 * |
《声波曲线重构反演技术在霍多莫尔地区储层预测中的应用》;朱国军;《物探化探计算技术》;20170531;第39卷(第3期);全文 * |
《拟声波曲线重构反演在柳屯储层预测中的应用》;李素贞,商建瓴;《复杂油气藏》;20111231;第4卷(第4期);全文 * |
《拟声波重构技术在巴彦呼舒凹陷储层预测中的应用》;卢伟;《断块油气田》;20160531;第23卷(第3期);第310-313页 * |
《英台断陷龙深2井区营城组二段储层预测》;胡立,邓宏文,邵明礼;《地下水》;20161130;第38卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108756867A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6681185B1 (en) | Method of seismic signal processing | |
WO2012092131A2 (en) | Stress in formations from azimuthal variation in acoustic and other properties | |
CN111060986B (zh) | 一种地层压力预测方法及岩性油藏评价方法 | |
CN100455768C (zh) | 确定扩径层段地层密度的方法 | |
CN113006779B (zh) | 一种砾岩地层岩石可钻性级值的评价方法和校正方法 | |
CN108756867B (zh) | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 | |
WO2016209822A1 (en) | Predicting hydraulic fracture propagation | |
CN112560246A (zh) | 一种目标井散点地层压力系数的预测方法 | |
CN112180464B (zh) | 一种储层物性的识别方法 | |
CN110967742A (zh) | 一种孔隙度反演方法及系统 | |
CN108732621B (zh) | 一种基于ffc-电阻率的随钻精细时深预测方法 | |
US11035966B2 (en) | Method for estimating reservoir formation quality using seismic velocities | |
CN113720745A (zh) | 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法 | |
CN114086938A (zh) | 一种非均质性砂岩储层的含气饱和度预测方法 | |
CN104699975A (zh) | 从声电效应井下探测器测量数据中提取参数的方法 | |
CN112709566B (zh) | 一种碳酸盐岩气藏产层的预测方法 | |
CN108107482A (zh) | 一种少井条件下烃源岩多参数预测方法及系统 | |
CN113945970B (zh) | 一种致密砂岩储层预测方法 | |
RU2225020C1 (ru) | Способ геофизической разведки для определения нефтепродуктивности трещинных глинистых коллекторов в межскважинном пространстве | |
Donald et al. | Wellsite Full Waveform Sonic Interpretation | |
CN107622450B (zh) | 基于测井资料快速判别水平井具备经济产量的方法 | |
CN116856921A (zh) | 一种潜山裂缝性储层的产能评价方法、系统及存储介质 | |
Carpenter | Seismic Work Flow Characterizes Deep, Thin Formation Layers Offshore South China Sea | |
Kerimova | NEW APPROACH FOR CORRECTION OF BOUNDARIES OF HORIZONS AND LAYER SERIES BY THE WELL DATA | |
Jugastreanu et al. | CONSIDERATION ON APPRAISAL OF THE RESERVOIR PRESSURE BY GEOPHYSICAL METHODS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211119 |