CN107255834A - 一种基于地震约束的声波测井曲线校正方法 - Google Patents
一种基于地震约束的声波测井曲线校正方法 Download PDFInfo
- Publication number
- CN107255834A CN107255834A CN201710357237.9A CN201710357237A CN107255834A CN 107255834 A CN107255834 A CN 107255834A CN 201710357237 A CN201710357237 A CN 201710357237A CN 107255834 A CN107255834 A CN 107255834A
- Authority
- CN
- China
- Prior art keywords
- data
- well
- sound wave
- wavelet
- curve
- 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.)
- Pending
Links
- 239000002131 composite material Substances 0.000 claims abstract description 29
- 238000004088 simulation Methods 0.000 claims description 2
- 238000012937 correction Methods 0.000 abstract description 21
- 238000012545 processing Methods 0.000 abstract description 4
- 238000000034 method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于地震约束的声波测井曲线校正方法,涉及油气田勘探开发的测井资料处理领域,利用地震记录校正测井曲线的思路或流程,将声波测井资料转换成高分辨率的合成记录,将此合成记录与对应位置的地震记录进行匹配,根据匹配情况优选匹配程度高的井段,将声波曲线段与地震记录段建立关系,再将该关系作用于所有井的地震数据,获得校正后的声波数据。本发明从全新的角度来进行声波曲线的校正,在较大程度上克服该问题,得到反映横向上相对变化的声波曲线。
Description
技术领域
本发明涉及油气田勘探开发的测井资料处理领域,确切第说涉及一种基于地震约束的声波测井曲线校正方法。
背景技术
在一个区域的测井资料中,不同井的同一参数曲线一般都存在非地质因素引起的差别,有些井间差别较大。造成这些差别的原因很多,如井眼垮塌及泥浆污染、裸眼测井与套管测井方式不同、测井仪器类型不一以及设备精度误差、测量时间变化等的不同,而且不同参数受到的影响因素不同。这样造成测井解释出现偏差,甚至错误。此外,还造成测井资料不能有效应用于地震资料。特别是,当反映时深关系的声波测井曲线出现较大误差时,与地震剖面就建不起正确对应关系,使得各条测井曲线都不能与地震资料正确对应,以致地震资料难以正确应用。为此,需要对测井曲线进行归一化处理或进行校正,校正后的曲线应基本克服非地质因素引起的误差。
测井资料归一化或校正处理的方法很多,一般是采取多井统计、回归方式,但这些方式受测井资料本身质量和井口数的影响,特别是在横向上出现局部地质情况变化时,难以判断测井曲线是否是正确的,可能造成统计失误而得不到正确的校正结果。通常的校正还会采取减少条件或分类的方式来校正,如选择同一种仪器所测量的、同为裸眼测量的、同一批次测量的等来提高校正精度。进一步,还可以在分类校正后再统一校正。
发明内容
本发明的目的在于克服现有测井资料的校正方法存在横向上变化不易把握、校正结果可靠性受影响的技术问题,从全新的角度来进行声波曲线的校正,在较大程度上克服该问题,得到反映横向上相对变化的声波曲线。
本发明是通过采用下述技术方案实现的:
一种基于地震约束的声波测井曲线校正方法,其特征在于:利用地震记录校正测井曲线的思路或流程,将声波(可含密度)测井资料转换成高分辨率的合成记录,将此合成记录与对应位置的地震记录进行匹配,根据匹配情况优选匹配程度高(相关系数高)的井段,将声波曲线段与地震记录段建立关系,再将该关系作用于所有井的地震数据,获得校正后的声波数据。
具体步骤如下:
a、输入声波测井资料,将声波数据转换成速度数据;
b、将速度与密度相乘,若没有密度曲线,则密度取值1,得到波阻抗,再计算反射系数;
c、给定子波及子波主频;
d、将反射系数与子波进行褶积运算,即得到合成记录;
e、将研究区域内的井的合成记录与井旁地震数据进行相关运算;
f、选择相关系数高,即合成记录与地震数据匹配好的井;
g、再在相关系数高的井中,再进一步优选匹配最优的声波曲线段;
h、将所选的声波曲线段与对应的地震数据,依次输入到第i步中;
i、建立声波曲线与地震数据的数学关系,用函数表示为y=f(x),其中,x为输入参数-地震数据,y为输出参数-声波曲线;
j、将建立的关系应用于所有井的井旁地震数据,得到新的声波曲线。
所述步骤a中,输入声波测井资料包含密度测井曲线。
所述步骤i中的函数包含各种可以实现本发明目的的数学物理函数。
合成记录的子波形态可以用固定的,也可以用任意的,包括雷克子波、测井与地震联合提取的子波、直接从地震数据上提取的子波、模拟的子波和假设子波。
子波主频根据地震数据主频给定。可以任意选择算法,只要能建立声波测井曲线与地震数据的关系都可以应用,如,统计法、拟合法、随机反演法等。
与现有技术相比,本发明所达到的有益效果如下:
1、本发明可充分利用地震数据具有横向连续性及反映地质横向相对变化的特点,能够对声波测井资料横向不连续、不稳定的缺陷进行有效弥补。
2、本发明为一种基于地震约束的声波测井曲线校正方法,采用地震资料约束来校正声波测井资料,与仅用测井资料校正的方式形成差异:
要应用声波测井数据(可含密度测井数据),其中声波数据是必须的,有密度曲线时应使用。先用声波与密度测井资料制作合成记录。至少2口井有声波曲线,至少1口井的合成记录与地震资料匹配较好(由应用研究人员确定匹配效果)。必须有地震数据,以地震剖面的形式呈现,地震剖面由地震记录或地震道构成。建立关系时要根据合成记录与地震记录匹配情况来定位。定位后,要选择最优的井段来建立地震记录与声波之间的关系。
附图说明
下面将结合说明书附图和具体实施方式对本发明作进一步的详细说明,其中:
图1为本发明流程示意图;
图2为合成记录与地震剖面匹配对比与优选示意图;
图3为声波测井曲线校正前后合成记录与地震剖面匹配对比示意图
(图中,上部为校正前,下部为校正后,各自中间部分的短图为合成记录,两侧为地震剖面。合成记录与地震剖面的一致性上差下好,相关系数分别为0.265和0.914);
图4为原始声波曲线(左)和校正后的声波曲线(右)对比示意图(最大变化在中部和上部)。
具体实施方式
本发明利用地震资料作约束,以一种新方式进行测井曲线的校正。
地震资料与测井资料的显著不同在于,前者记录地层空间(三维)信息,后者记录井筒(线)信息。一般情况下,对于一个区域而言,地震资料的采集、处理参数在横向上一般是保持不变的,得到的地震资料在横向上就会保持其相对关系的一致和稳定,将其与测井结合,可校正测井资料非地质原因引起的横向变化,即在地震资料控制下,对测井资料进行校正。
因测井声波与地震纵波形成机制相近,可以将声波曲线转换成合成记录,与地震资料进行匹配。优选匹配较好的井段,建立声波测井曲线与地震数据的数学、物理关系,将该关系作用于匹配不好的井段的地震资料,得到校正后的声波曲线。
本发明可以用于不同参数类型曲线的校正,其中声波曲线可以转换为合成记录与地震资料对比,确定该曲线是否正确,然后确定是否用本发明方法进行校正,但其它曲线因事先不能与地震资料直观对比做出判断,故本发明不涉及除声波以外的曲线。
实施例1
根据图1,可以完成声波测井曲线校正。
输入声波测井资料(也可以包含密度测井曲线),将声波数据转换成速度数据;
将速度与密度相乘(若没有密度曲线,则密度取值1),得到波阻抗,再计算反射系数;
给定子波及子波主频;
将反射系数与子波进行褶积运算,即得到合成记录;
将研究区域内的井的合成记录与井旁地震数据进行相关运算;
选择相关系数高,即合成记录与地震数据匹配好的井;
再在相关系数高的井中,再进一步优选匹配最优的声波曲线段;
将所选的声波曲线段与对应的地震数据,依次输入到第9步中;
建立声波曲线与地震数据的数学关系,用函数表示为y=f(x)。其中,x为输入参数(地震数据),y为输出参数(声波曲线)。这里的函数包含各种可以实现本发明目的的数学物理函数。
将建立的关系应用于所有井的井旁地震数据,得到新的声波曲线。
实施例2
具体实施过程如下:
资料前提:在一个勘探开发的地区,或一个地质构造,已经有2口以上的钻井,这些钻井已经做过测井,测井资料中必须有声波、密度(最好有密度,没有也可以)。在此地区(或地质构造),已经做过地震勘探采集,采集的地震资料经过处理形成了地震剖面数据。钻井与地震资料在空间上位置关系是明确的。
提取每口井的声波、密度曲线,该曲线从上到下对应不同深度地层,即一个声波、密度数据点对应一个深度。每个点的声波是声波在地层中传播时间的倒数,可转换为速度,各点的速度与密度相乘为对应地层的波阻抗,求上一个地层与下一个地层的波阻抗差与波阻抗和之比,得到反射系数。
给定一个类似于震源的子波,子波可以用很多方法计算,最常用的是雷克子波,具有一个单峰和弱小的旁瓣,其谱为白化谱,分辨能力最强。
将反射系数序列与子波进行褶积运算,即得到合成记录,褶积运算方法与现有方法相同。
将合成记录与井旁地震记录匹配,确定位置,同时选出最优段。图2中,井名13-48:702,其合成记录(中部短图)与地震记录相关系数为0.528,上部、下部匹配较差,优选中间段(方框中),相关系数达到0.675。其它井类似。
将优选井段的声波曲线与对应位置的地震记录建立数学关系。再将关系应用于其它井(包括本井),即得到各井校正后的声波曲线。
图3中,井名13-21:732,其校正前后的声波曲线制作的合成记录与地震剖面对比。可见校正后匹配效果明显变好,相关系数从0.265达到0.914。校正后,合成记录波形(含振幅、相位)与地震记录基本一致。
图4中,井名13-21:732,为校正前后的声波曲线对比,声波曲线转换为了速度曲线,横轴为速度,单位m/s,纵轴为深度,单位m。曲线的总体面貌变化不大,但细节变化多,主要在上部和中部变化明显。与图3是对应的。
Claims (6)
1.一种基于地震约束的声波测井曲线校正方法,其特征在于:利用地震记录校正测井曲线的思路或流程,将声波测井资料转换成高分辨率的合成记录,将此合成记录与对应位置的地震记录进行匹配,根据匹配情况优选匹配程度高的井段,将声波曲线段与地震记录段建立关系,再将该关系作用于所有井的地震数据,获得校正后的声波数据。
2.根据权利要求1所述的一种基于地震约束的声波测井曲线校正方法,其特征在于:具体步骤如下:
a、输入声波测井资料,将声波数据转换成速度数据;
b、将速度与密度相乘,若没有密度曲线,则密度取值1,得到波阻抗,再计算反射系数;
c、给定子波及子波主频;
d、将反射系数与子波进行褶积运算,即得到合成记录;
e、将研究区域内的井的合成记录与井旁地震数据进行相关运算;
f、选择相关系数高,即合成记录与地震数据匹配好的井;
g、再在相关系数高的井中,再进一步优选匹配最优的声波曲线段;
h、将所选的声波曲线段与对应的地震数据,依次输入到第i步中;
i、建立声波曲线与地震数据的数学关系,用函数表示为y=f(x),其中,x为输入参数-地震数据,y为输出参数-声波曲线;
j、将建立的关系应用于所有井的井旁地震数据,得到新的声波曲线。
3.根据权利要求2所述的一种基于地震约束的声波测井曲线校正方法,其特征在于:所述步骤a中,输入声波测井资料包含密度测井曲线。
4.根据权利要求2所述的一种基于地震约束的声波测井曲线校正方法,其特征在于:所述步骤i中的函数包含各种可以实现本发明目的的数学物理函数。
5.根据权利要求2所述的一种基于地震约束的声波测井曲线校正方法,其特征在于:合成记录的子波形态包括雷克子波、测井与地震联合提取的子波、直接从地震数据上提取的子波、模拟的子波和假设子波。
6.根据权利要求2所述的一种基于地震约束的声波测井曲线校正方法,其特征在于:子波主频根据地震数据主频给定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710357237.9A CN107255834A (zh) | 2017-05-19 | 2017-05-19 | 一种基于地震约束的声波测井曲线校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710357237.9A CN107255834A (zh) | 2017-05-19 | 2017-05-19 | 一种基于地震约束的声波测井曲线校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107255834A true CN107255834A (zh) | 2017-10-17 |
Family
ID=60027642
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710357237.9A Pending CN107255834A (zh) | 2017-05-19 | 2017-05-19 | 一种基于地震约束的声波测井曲线校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107255834A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001815A (zh) * | 2018-08-27 | 2018-12-14 | 中国石油化工股份有限公司 | 基于地震和井径约束的声波曲线校正方法 |
CN109946739A (zh) * | 2019-03-15 | 2019-06-28 | 成都理工大学 | 一种基于压缩感知理论的地震剖面增强方法 |
CN110133717A (zh) * | 2019-04-15 | 2019-08-16 | 长江大学 | 确定区域地震波相位的方法及设备 |
CN111208558A (zh) * | 2018-11-22 | 2020-05-29 | 中国石油天然气股份有限公司 | 超深低幅度三维地质构造的建立方法及装置 |
CN111650647A (zh) * | 2020-05-26 | 2020-09-11 | 中国石油天然气集团有限公司 | 一种基于地震数据约束的声波测井曲线重构方法 |
CN111650646A (zh) * | 2020-05-26 | 2020-09-11 | 中国石油天然气集团有限公司 | 一种用于地震处理中的声波曲线校正方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002071101A1 (en) * | 2001-03-01 | 2002-09-12 | Baker Hughes Incorporated | An inhomogeneous background based focusing method for multi-array induction measurements in a deviated well |
CN103389513A (zh) * | 2013-08-02 | 2013-11-13 | 中国石油集团川庆钻探工程有限公司 | 应用声波测井资料约束反演提高地震资料分辨率的方法 |
WO2014084752A1 (en) * | 2012-11-30 | 2014-06-05 | Schlumberger Holdings Limited | A method for processing acoustic waveforms |
CN105301657A (zh) * | 2015-10-29 | 2016-02-03 | 中国石油天然气股份有限公司 | 一种基于岩石物理意义的曲线校正方法 |
CN106405645A (zh) * | 2016-08-30 | 2017-02-15 | 英得赛斯科技(北京)有限公司 | 一种基于资料品质分析的信噪比可控的地震拓频处理方法 |
-
2017
- 2017-05-19 CN CN201710357237.9A patent/CN107255834A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002071101A1 (en) * | 2001-03-01 | 2002-09-12 | Baker Hughes Incorporated | An inhomogeneous background based focusing method for multi-array induction measurements in a deviated well |
WO2014084752A1 (en) * | 2012-11-30 | 2014-06-05 | Schlumberger Holdings Limited | A method for processing acoustic waveforms |
CN103389513A (zh) * | 2013-08-02 | 2013-11-13 | 中国石油集团川庆钻探工程有限公司 | 应用声波测井资料约束反演提高地震资料分辨率的方法 |
CN105301657A (zh) * | 2015-10-29 | 2016-02-03 | 中国石油天然气股份有限公司 | 一种基于岩石物理意义的曲线校正方法 |
CN106405645A (zh) * | 2016-08-30 | 2017-02-15 | 英得赛斯科技(北京)有限公司 | 一种基于资料品质分析的信噪比可控的地震拓频处理方法 |
Non-Patent Citations (1)
Title |
---|
张绍辉 等: "最优相关分析技术在精细层位标定中的应用", 《特种油气藏》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001815A (zh) * | 2018-08-27 | 2018-12-14 | 中国石油化工股份有限公司 | 基于地震和井径约束的声波曲线校正方法 |
CN109001815B (zh) * | 2018-08-27 | 2019-12-31 | 中国石油化工股份有限公司 | 基于地震和井径约束的声波曲线校正方法 |
CN111208558A (zh) * | 2018-11-22 | 2020-05-29 | 中国石油天然气股份有限公司 | 超深低幅度三维地质构造的建立方法及装置 |
CN109946739A (zh) * | 2019-03-15 | 2019-06-28 | 成都理工大学 | 一种基于压缩感知理论的地震剖面增强方法 |
CN110133717A (zh) * | 2019-04-15 | 2019-08-16 | 长江大学 | 确定区域地震波相位的方法及设备 |
CN111650647A (zh) * | 2020-05-26 | 2020-09-11 | 中国石油天然气集团有限公司 | 一种基于地震数据约束的声波测井曲线重构方法 |
CN111650646A (zh) * | 2020-05-26 | 2020-09-11 | 中国石油天然气集团有限公司 | 一种用于地震处理中的声波曲线校正方法 |
CN111650647B (zh) * | 2020-05-26 | 2023-06-13 | 中国石油天然气集团有限公司 | 一种基于地震数据约束的声波测井曲线重构方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107255834A (zh) | 一种基于地震约束的声波测井曲线校正方法 | |
US10802171B2 (en) | High resolution seismic data derived from pre-stack inversion and machine learning | |
US7840394B2 (en) | Method for generating a 3D earth model | |
US11230922B2 (en) | Fracture interpretation with resistivity and sonic logs in biaxial anisotropic formations | |
US20070168133A1 (en) | Computer-based method for while-drilling modeling and visualization of layered subterranean earth formations | |
US20210089897A1 (en) | High-resolution earth modeling using artificial intelligence | |
CN107219553B (zh) | 基于gr分频反演的暗河充填预测方法 | |
AU2019340410B2 (en) | Reservoir characterization utilizing inversion of resampled seismic data | |
CA3010894C (en) | Methods of selecting an earth model from a plurality of earth models | |
DE112016005567T5 (de) | Erzeugen eines Erdmodells anhand räumlicher Korrelationen von äquivalenten Erdmodellen | |
CN104090303B (zh) | 一种地震反演方法和装置 | |
US11703608B2 (en) | Reservoir characterization using machine-learning techniques | |
Munoz et al. | Automatically tying well logs to seismic data | |
US20070064529A1 (en) | Method for obtaining a seismic wavelet using seismic data and well logs in a deviated well | |
Leahy et al. | Uncertainty in subsurface interpretation: a new workflow | |
CN111650646A (zh) | 一种用于地震处理中的声波曲线校正方法 | |
CN106199715A (zh) | 基于模型的强反射分离方法和装置 | |
CN107991705A (zh) | 基于二维统计特征的测井曲线校正方法和装置 | |
US11914089B2 (en) | Method and device for determining sonic slowness | |
US11693141B2 (en) | Methods and systems for processing borehole dispersive waves with a physics-based machine learning analysis | |
Donald et al. | Digital Transformation of Borehole Sonic Services in the Oil and Gas Industry | |
CN107783192B (zh) | 碳酸盐岩古河道充填物性质地震预测方法及系统 | |
US20240201404A1 (en) | Methods and systems for analyzing seismic amplitude | |
US20240168190A1 (en) | Well Log Data Conditioning using Quantified Uncertainties in Rock Physics and Seismic Inversion Results | |
CN107561580A (zh) | 基于多层多井建立初始地质模型的方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171017 |