CN103485768A - 声波测井曲线的构建方法 - Google Patents
声波测井曲线的构建方法 Download PDFInfo
- Publication number
- CN103485768A CN103485768A CN201210196084.1A CN201210196084A CN103485768A CN 103485768 A CN103485768 A CN 103485768A CN 201210196084 A CN201210196084 A CN 201210196084A CN 103485768 A CN103485768 A CN 103485768A
- Authority
- CN
- China
- Prior art keywords
- curve
- acoustic logging
- lld
- lls
- resistivity
- 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
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
技术领域
本发明涉及石油测井技术领域,特别涉及一种声波测井曲线的构建方法。
背景技术
在石油勘探和测井技术领域,随着测井技术和测井仪器的发展,测井资料的处理与解释可以更真实地反映地下储层情况,但在一些需要重新解释的老井中,存在测井系列较老,测井曲线严重缺失的情况,甚至有的井只有井径和电位或电阻率曲线,这给测井数字处理和解释工作带来了很大的困难,基于此,开展数据挖掘和测井曲线构建研究是非常必要的。
所谓测井曲线构建即是利用测井项目中已有的测井曲线和所要构建曲线之间的关系来得到测井曲线的方法。储层特征曲线构建是基于不同地质体的属性在不同物理场中都有不同反映,同一地质体的属性在不同物理场中的显示都有相关性,同一地质体的不同属性在不同物理场中都有所侧重。储层特征曲线构建就是根据储层预测目标,综合有利于储层预测的相关信息,得到一条更能够突出储层分辨率的特征曲线。储层特征曲线构建方法同时具有一定的局限性,必须根据区域地质背景、测井响应特征及对地层岩性变化规律的指示程度,优选构建方法,得到较为满意的曲线构建效果。
在实现本发明的过程中,发明人发现现有技术至少存在以下问题:在一些需要重新解释评价的老井中,存在测井系列较老、测井曲线严重缺失的情况,甚至有的井只有井径、自然电位和电阻率曲线,这给测井数字处理和解释工作带来了很大的困难,为此开展数据挖掘和测井曲线构建研究具有迫切需要。
发明内容
为了解决现有技术的问题,本发明实施例提供了一种声波测井曲线的构建方法。所述技术方案如下:
提供了一种声波测井曲线的构建方法,所述方法包括:
步骤一,对待构建声波测井曲线所在区域的地质背景资料进行了解、对地层电阻率曲线特征进行分析,获取待构建声波测井曲线所在区域的地质背景数据和地层电阻率曲线特征数据;
步骤二,根据所述地质背景数据和所述地层电阻率曲线特征数据,确定所述待构建声波测井曲线所在井的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP,则所述待构建声波测井曲线所在井的声波测井曲线为:
其中,a0、a1为待定系数,b1、b2为待定拟合指数,a0、a1、b1、b2利用最小二乘法拟合确定。
这样,当确定了待定系数a0、a1,待定拟合指数b1、b2后,只要将该区域的未测井的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP代入声波测井曲线 中即可得到该未测井的声波测井曲线。
进一步地,针对某一油田待测井,利用最小二乘法拟合的所述待定系数:a0=259.6373,a1=-0.0002;所述待定拟合指数:b1=1.2152,b2=0.9045。即该待测井的声波测井曲线为:AC=259.6373-0.0002×DEP1.2152×(LLD+LLS)0.9045×LLD/LLS。
进一步地,当所述待构建声波测井曲线所在井为重新解释评估的老井时,若所述老井缺少已有的测井曲线,则利用虚拟井预测曲线拟合构建声波测井曲线。
另一方面,本发明还提供给了一种声波测井曲线的构建装置,所述装置包括:
数据获取模块,用于获取待构建声波测井曲线所在井区域的地质背景数据和地层电阻率曲线特征数据;
确定模块,用于根据所述地质背景数据和所述地层电阻率曲线特征数据,确定所述待构建声波测井曲线所在井的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP;
构建模块,用于根据所述深侧向电阻率曲线LLD、所述浅侧向电阻率曲线LLS和所述测井深度曲线DEP获取所述待构建声波测井曲线所在井的声波测井曲线。
进一步地,构建模块具体用于确定所述待构建声波测井曲线所在井的声波测井曲线为:
其中,a0、a1为待定系数,b1、b2为待定拟合指数,a0、a1、b1、b2利用最小二乘法拟合确定。
本发明实施例提供的技术方案带来的有益效果是:
通过对待构建声波测井曲线所在井区域的地质背景资料进行了解、对地层电阻率曲线特征进行分析,构建了该区域的老井或者待测井的声波测井曲线,解决了由于缺少声波测井曲线而无法计算储层孔隙度的问题,能够及时服务或满足实际生产的需要,实现了油气储层物性参数的定量计算,具有较好的应用效果。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例一提供的声波测井曲线构建方法的构建声波时差曲线与实测声波时差曲线对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
实施例一
本实施例提供了一种声波测井曲线的构建方法,所述方法包括如下步骤:
步骤一,对待构建声波测井曲线所在井区域的地质背景资料进行了解、对地层电阻率曲线特征进行分析,获取待构建声波测井曲线所在区域的地质背景数据和地层电阻率曲线特征数据;
步骤二,根据地质背景数据和地层电阻率曲线特征数据,确定待构建声波测井曲线所在井的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP,则待构建声波测井曲线所在井的声波测井曲线为:
其中,a0、a1为待定系数,b1、b2为待定拟合指数,a0、a1、b1、b2利用最小二乘法拟合确定。
这样,当确定了待定系数a0、a1,待定拟合指数b1、b2后,只要将该区域的未测井的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP代入声波测井曲线 中,即可得到该未测井的声波测井曲线。
参见图1,选取某油田已测有声波时差曲线资料的井,在构建声波测井曲线时,利用该油田的已有声波时差曲线与所构建的声波测井曲线进行对比和效果检验,具体的步骤是:
对待构建声波测井曲线所在井的所在区域的地质背景资料进行了解、对地层电阻率曲线特征进行分析,获取待构建声波测井曲线所在井的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP,分析地层电阻率曲线LLD和LLS与声波测井曲线是否很好地反映储层特征信息。
利用所得到的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP构建声波测井曲线:
利用最小二乘法拟合确定待定系数:a0=259.6373,a1=-0.0002;待定拟合指数:b1=1.2152,b2=0.9045。将待定系数和待定拟合指数代入所构建的声波测井曲线中,即得到该油田的声波测井曲线为:
AC=259.6373-0.0002×DEP1.2152×(LLD+LLS)0.9045×LLD/LLS,因此,对于该油田内的某一未测量声波测井曲线的井,只要将该井的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP代入该发明的声波测井曲线计算公式,即可得到该井的声波测井曲线。
实施例二
本实施例提供了一种声波测井曲线的构建装置,该装置包括:
数据获取模块,用于获取待构建声波测井曲线所在区域的地质背景数据和地层电阻率曲线特征数据;
确定模块,用于根据所述地质背景数据和所述地层电阻率曲线特征数据,确定所述待构建声波测井曲线所在井的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP;
构建模块,用于根据所述深侧向电阻率曲线LLD、所述浅侧向电阻率曲线LLS和所述测井深度曲线DEP获取所述待构建声波测井曲线所在井的声波测井曲线。
优选地,构建模块具体用于确定所述待构建声波测井曲线所在井的声波测井曲线为:
其中,a0、a1为待定系数,b1、b2为待定拟合指数,a0、a1、b1、b2利用最小二乘法拟合确定。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种声波测井曲线的构建方法,其特征在于,所述方法包括:
获取待构建声波测井曲线所在区域的地质背景数据和地层电阻率曲线特征数据;
根据所述地质背景数据和所述地层电阻率曲线特征数据,确定所述待构建声波测井曲线所在区域的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP,则所述待构建声波测井曲线所在区域的声波测井曲线为:
式(1)中,a0、a1为待定系数,b1、b2为待定拟合指数,a0、a1、b1、b2利用最小二乘法拟合确定。
2.根据权利要求1所述的方法,其特征在于,将未测声波曲线所在井的深侧向电阻率曲线LLD、浅侧向电阻率曲线LLS和测井深度曲线DEP代入式(1)中,即可计算得到未测井的声波测井曲线AC。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210196084.1A CN103485768B (zh) | 2012-06-13 | 2012-06-13 | 声波测井曲线的构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210196084.1A CN103485768B (zh) | 2012-06-13 | 2012-06-13 | 声波测井曲线的构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103485768A true CN103485768A (zh) | 2014-01-01 |
CN103485768B CN103485768B (zh) | 2016-03-09 |
Family
ID=49826313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210196084.1A Active CN103485768B (zh) | 2012-06-13 | 2012-06-13 | 声波测井曲线的构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103485768B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103850679A (zh) * | 2014-04-01 | 2014-06-11 | 北京师范大学 | 一种利用多种测井曲线对声波时差曲线进行重构的方法 |
CN104977616A (zh) * | 2014-04-03 | 2015-10-14 | 中国石油化工股份有限公司 | 针对风化粘土层的定量识别方法 |
CN105092813A (zh) * | 2015-08-21 | 2015-11-25 | 中国石油天然气集团公司 | 一种确定煤层含气量的方法 |
CN105158803A (zh) * | 2015-07-27 | 2015-12-16 | 中国石油天然气集团公司 | 一种选取优势储层的方法 |
CN105569647A (zh) * | 2014-11-06 | 2016-05-11 | 中国石油化工股份有限公司 | 一种储层特征值的生成方法 |
CN107290800A (zh) * | 2016-03-30 | 2017-10-24 | 中国石油化工股份有限公司 | 实用的钻前测井曲线预测方法 |
CN108756867A (zh) * | 2018-05-11 | 2018-11-06 | 中国地质调查局油气资源调查中心 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
CN109116440A (zh) * | 2017-06-23 | 2019-01-01 | 中国石油化工股份有限公司 | 致密灰岩储层裂缝识别方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA995751A (en) * | 1971-10-13 | 1976-08-24 | Alberto Khatchikian | Well logging data processing methods |
CN201747364U (zh) * | 2010-06-01 | 2011-02-16 | 中国石油天然气集团公司 | 一种声波时差曲线重构设备 |
CN102353985A (zh) * | 2011-06-01 | 2012-02-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于非下采样Contourlet变换的拟声波曲线构建方法 |
-
2012
- 2012-06-13 CN CN201210196084.1A patent/CN103485768B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA995751A (en) * | 1971-10-13 | 1976-08-24 | Alberto Khatchikian | Well logging data processing methods |
CN201747364U (zh) * | 2010-06-01 | 2011-02-16 | 中国石油天然气集团公司 | 一种声波时差曲线重构设备 |
CN102353985A (zh) * | 2011-06-01 | 2012-02-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于非下采样Contourlet变换的拟声波曲线构建方法 |
Non-Patent Citations (1)
Title |
---|
宋维琪等: "测井声波时差反演重构技术研究及应用", 《地震地质》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103850679A (zh) * | 2014-04-01 | 2014-06-11 | 北京师范大学 | 一种利用多种测井曲线对声波时差曲线进行重构的方法 |
CN103850679B (zh) * | 2014-04-01 | 2019-12-17 | 北京师范大学 | 一种利用多种测井曲线对声波时差曲线进行重构的方法 |
CN104977616A (zh) * | 2014-04-03 | 2015-10-14 | 中国石油化工股份有限公司 | 针对风化粘土层的定量识别方法 |
CN105569647A (zh) * | 2014-11-06 | 2016-05-11 | 中国石油化工股份有限公司 | 一种储层特征值的生成方法 |
CN105158803A (zh) * | 2015-07-27 | 2015-12-16 | 中国石油天然气集团公司 | 一种选取优势储层的方法 |
CN105158803B (zh) * | 2015-07-27 | 2017-10-17 | 中国石油天然气集团公司 | 一种选取优势储层的方法 |
CN105092813A (zh) * | 2015-08-21 | 2015-11-25 | 中国石油天然气集团公司 | 一种确定煤层含气量的方法 |
CN107290800A (zh) * | 2016-03-30 | 2017-10-24 | 中国石油化工股份有限公司 | 实用的钻前测井曲线预测方法 |
CN109116440A (zh) * | 2017-06-23 | 2019-01-01 | 中国石油化工股份有限公司 | 致密灰岩储层裂缝识别方法 |
CN108756867A (zh) * | 2018-05-11 | 2018-11-06 | 中国地质调查局油气资源调查中心 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
CN108756867B (zh) * | 2018-05-11 | 2021-11-19 | 中国地质调查局油气资源调查中心 | 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103485768B (zh) | 2016-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103485768B (zh) | 声波测井曲线的构建方法 | |
Willis et al. | Spatial orientation and distribution of reservoir fractures from scattered seismic energy | |
Shemeta et al. | It's a matter of size: Magnitude and moment estimates for microseismic data | |
CN103149589B (zh) | 火成岩油气勘探方法及装置 | |
Close et al. | Integrated workflows for shale gas and case study results for the Horn River Basin, British Columbia, Canada | |
CN104345346B (zh) | 一种获取裂缝宽度的方法 | |
Langston | Local earthquake wave propagation through Mississippi embayment sediments, part I: Body-wave phases and local site responses | |
CN104345345A (zh) | 一种页岩储层总有机碳toc含量预测方法 | |
EP2420863A2 (en) | Determining a Position of a Geological Layer Relative to a Wavelet Response in Seismic Data | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
CN102650701A (zh) | 一种基于偶极子波进行储层波阻抗预测的方法 | |
Baird et al. | Frequency-dependent seismic anisotropy due to fractures: Fluid flow versus scattering | |
US6327538B1 (en) | Method and apparatus for evaluating stoneley waves, and for determining formation parameters in response thereto | |
CN105277982A (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
CN101285381B (zh) | 一种泄漏模式波反演软地层横波速度的方法 | |
Hornby et al. | Comparison of fracture apertures computed from electrical borehole scans and reflected Stoneley waves: an integrated interpretation | |
CN104155701A (zh) | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 | |
CN102819038B (zh) | 一种碳酸盐岩内油水识别的方法及系统 | |
CN104316965A (zh) | 一种裂缝方位和强度的预测方法及系统 | |
CN103821495A (zh) | 测井方法 | |
Woolery et al. | The central United States seismic observatory: Site characterization, instrumentation, and recordings | |
CN101315428A (zh) | 基于横波速度的高孔隙碎屑岩沉积地层的高分辨率勘探方法 | |
Cai et al. | Intelligent calculation method of relative sonic attenuation and its application to fracture evaluation in tight sandstone reservoir | |
Somasundaram et al. | Seismic attribute analysis for fracture detection and porosity prediction: A case study from tight volcanic reservoirs, Barmer Basin, India | |
CN104375166B (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 |