CN111045108B - 随钻横波计算方法及系统 - Google Patents
随钻横波计算方法及系统 Download PDFInfo
- Publication number
- CN111045108B CN111045108B CN201811184200.1A CN201811184200A CN111045108B CN 111045108 B CN111045108 B CN 111045108B CN 201811184200 A CN201811184200 A CN 201811184200A CN 111045108 B CN111045108 B CN 111045108B
- Authority
- CN
- China
- Prior art keywords
- interval data
- wave
- drilling
- data
- equation
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
公开了一种随钻横波计算方法及系统。该方法可以包括:根据实时随钻录井和测井数据,获得电阻率曲线与纵波声波时差;根据电阻率曲线判断岩石岩性,获得不同岩性的间隔数据Δmi;根据间隔数据的纵波声波时差,获得间隔数据对应的纵波速度;根据间隔数据的岩性,选择转换公式;根据间隔数据对应的纵波速度与转换公式,计算间隔数据对应的横波速度;根据所有间隔数据对应的横波速度,计算最终的横波速度。本发明通过将实时测量获得的声波时差转为实时的横波数据,随着钻遇不同的岩层,采用不同的经验公式将声波时差数据实时的转为横波数据,实时准确,从而为随钻叠前地震反演提供数据基础。
Description
技术领域
本发明涉及石油化工产业地球物理勘探领域,更具体地,涉及一种随钻横波计算方法及系统。
背景技术
传统的基于随钻测量和随钻测井的地质导向技术,在随钻过程中能够得到实时的测井测量数据包括:电阻率曲线、声波时差等,但是在进行随钻叠前地震反演过程中,需要用到横波数据,而横波在随钻过程中并不会通过实时的测井测量获取,因此需要通过岩石物理模型来进行求解获得,目前在随钻测井过程中还没有相关的横波预测技术。因此,有必要开发一种随钻横波计算方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种随钻横波计算方法及系统,其能够通过将实时测量获得的声波时差转为实时的横波数据,随着钻遇不同的岩层,采用不同的经验公式将声波时差数据实时的转为横波数据,实时准确,从而为随钻叠前地震反演提供数据基础。
根据本发明的一方面,提出了一种随钻横波计算方法。所述方法可以包括:根据实时随钻录井和测井数据,获得电阻率曲线与纵波声波时差;根据所述电阻率曲线判断岩石岩性,获得不同岩性的间隔数据Δmi;根据所述间隔数据的纵波声波时差,获得所述间隔数据对应的纵波速度;根据所述间隔数据的岩性,选择转换公式;根据所述间隔数据对应的纵波速度与所述转换公式,计算所述间隔数据对应的横波速度;根据所有间隔数据对应的横波速度,计算最终的横波速度。
优选地,所述转换公式包括科斯塔那(Castagna)一次方程、Castagna二次方程和克里夫(Krief)方程;当所述间隔数据的岩性为砂岩时,选择Castagna一次方程或Krief方程;当所述间隔数据的岩性为石灰岩时,选择Castagna二次方程。
优选地,所述Castagna一次方程为:
优选地,所述Castagna二次方程为:
优选地,所述Krief方程为:
根据本发明的另一方面,提出了一种随钻横波计算系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:根据实时随钻录井和测井数据,获得电阻率曲线与纵波声波时差;根据所述电阻率曲线判断岩石岩性,获得不同岩性的间隔数据Δmi;根据所述间隔数据的纵波声波时差,获得所述间隔数据对应的纵波速度;根据所述间隔数据的岩性,选择转换公式;根据所述间隔数据对应的纵波速度与所述转换公式,计算所述间隔数据对应的横波速度;根据所有间隔数据对应的横波速度,计算最终的横波速度。
优选地,所述转换公式包括Castagna一次方程、Castagna二次方程和Krief方程;当所述间隔数据的岩性为砂岩时,选择Castagna一次方程或Krief方程;当所述间隔数据的岩性为石灰岩时,选择Castagna二次方程。
优选地,所述Castagna一次方程为:
优选地,所述Castagna二次方程为:
优选地,所述Krief方程为:
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的随钻横波计算方法的步骤的流程图。
图2示出了根据本发明的一个实施例的计算不同间隔数据的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的随钻横波计算方法的步骤的流程图。
在该实施例中,根据本发明的随钻横波计算方法可以包括:步骤101,根据实时随钻录井和测井数据,获得电阻率曲线与纵波声波时差;步骤102,根据电阻率曲线判断岩石岩性,获得不同岩性的间隔数据Δmi;步骤103,根据间隔数据的纵波声波时差,获得间隔数据对应的纵波速度;步骤104,根据间隔数据的岩性,选择转换公式;步骤105,根据间隔数据对应的纵波速度与转换公式,计算间隔数据对应的横波速度;步骤106,根据所有间隔数据对应的横波速度,计算最终的横波速度。
在一个示例中,转换公式包括Castagna一次方程、Castagna二次方程和Krief方程;当间隔数据的岩性为砂岩时,选择Castagna一次方程或Krief方程;当间隔数据的岩性为石灰岩时,选择Castagna二次方程。
在一个示例中,Castagna一次方程为:
在一个示例中,Castagna二次方程为:
在一个示例中,Krief方程为:
具体地,根据本发明的随钻横波计算方法可以包括:
根据实时随钻录井和测井数据,获得电阻率曲线与纵波声波时差;根据电阻率曲线判断岩石岩性,获得不同岩性的间隔数据Δmi。
根据间隔数据的纵波声波时差,通过公式(4)获得间隔数据对应的纵波速度:
根据间隔数据的岩性,选择转换公式,其中,转换公式包括Castagna一次方程、Castagna二次方程和Krief方程,其中,Castagna一次方程为公式(1),Castagna二次方程为公式(2),Krief方程为公式(3)。
根据间隔数据对应的纵波速度与转换公式,计算间隔数据对应的横波速度,当间隔数据的岩性为砂岩时,选择Castagna一次方程或Krief方程;当间隔数据的岩性为石灰岩时,选择Castagna二次方程;当间隔数据的岩性为页岩时,选择Krief方程;当间隔数据的岩性为白云岩时,选择Castagna一次方程,其中,不同岩性的计算参数取值为:砂岩:A=0.804、B=0.856;白云岩:A=0.583、B=-0.078;页岩:A=0.770、B=-0.867;石灰岩:A=-0.055、B=1.017、C=-1.03。
根据所有间隔数据对应的横波速度,通过公式(5)计算最终的横波速度:
其中,Vs为最终的横波速度。
本方法通过将实时测量获得的声波时差转为实时的横波数据,随着钻遇不同的岩层,采用不同的经验公式将声波时差数据实时的转为横波数据,实时准确,从而为随钻叠前地震反演提供数据基础。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
图2示出了根据本发明的一个实施例的计算不同间隔数据的示意图。
根据本发明的随钻横波计算方法可以包括:
根据实时随钻录井和测井数据,获得电阻率曲线与纵波声波时差;根据电阻率曲线判断岩石岩性,获得不同岩性的间隔数据Δmi。根据间隔数据的纵波声波时差,通过公式(4)获得间隔数据对应的纵波速度。
根据间隔数据的岩性,选择转换公式,其中,转换公式包括Castagna一次方程、Castagna二次方程和Krief方程,其中,Castagna一次方程为公式(1),Castagna二次方程为公式(2),Krief方程为公式(3)。
根据间隔数据对应的纵波速度与转换公式,计算间隔数据对应的横波速度,如图2所示,当间隔数据的岩性为砂岩时,选择Castagna一次方程;当间隔数据的岩性为石灰岩时,选择Castagna二次方程;当间隔数据的岩性为页岩时,选择Krief方程;当间隔数据的岩性为白云岩时,选择Castagna一次方程,其中,不同岩性的计算参数取值为:砂岩:A=0.804、B=0.856;白云岩:A=0.583、B=-0.078;页岩:A=0.770、B=-0.867;石灰岩:A=-0.055、B=1.017、C=-1.03。
根据所有间隔数据对应的横波速度,通过公式(5)计算最终的横波速度,如图2所示。
综上所述,本发明通过将实时测量获得的声波时差转为实时的横波数据,随着钻遇不同的岩层,采用不同的经验公式将声波时差数据实时的转为横波数据,实时准确,从而为随钻叠前地震反演提供数据基础。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例,提供了一种随钻横波计算系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:根据实时随钻录井和测井数据,获得电阻率曲线与纵波声波时差;根据电阻率曲线判断岩石岩性,获得不同岩性的间隔数据Δmi;根据间隔数据的纵波声波时差,获得间隔数据对应的纵波速度;根据间隔数据的岩性,选择转换公式;根据间隔数据对应的纵波速度与转换公式,计算间隔数据对应的横波速度;根据所有间隔数据对应的横波速度,计算最终的横波速度。
在一个示例中,转换公式包括Castagna一次方程、Castagna二次方程和Krief方程;当间隔数据的岩性为砂岩时,选择Castagna一次方程或Krief方程;当间隔数据的岩性为石灰岩时,选择Castagna二次方程。
在一个示例中,Castagna一次方程为:
在一个示例中,Castagna二次方程为:
在一个示例中,Krief方程为:
本系统通过将实时测量获得的声波时差转为实时的横波数据,随着钻遇不同的岩层,采用不同的经验公式将声波时差数据实时的转为横波数据,实时准确,从而为随钻叠前地震反演提供数据基础。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (6)
1.一种随钻横波计算方法,其特征在于,包括:
根据实时随钻录井和测井数据,获得电阻率曲线与纵波声波时差;
根据所述电阻率曲线判断岩石岩性,获得不同岩性的间隔数据Δmi;
根据所述间隔数据的纵波声波时差,获得所述间隔数据对应的纵波速度;
根据所述间隔数据的岩性,选择转换公式;
根据所述间隔数据对应的纵波速度与所述转换公式,计算所述间隔数据对应的横波速度;
根据所有间隔数据对应的横波速度,计算最终的横波速度;
当所述间隔数据的岩性为砂岩时,选择克里夫方程为转换方程;
当所述间隔数据的岩性为石灰岩时,选择科斯塔那二次方程。
4.一种随钻横波计算系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
根据实时随钻录井和测井数据,获得电阻率曲线与纵波声波时差;
根据所述电阻率曲线判断岩石岩性,获得不同岩性的间隔数据Δmi;
根据所述间隔数据的纵波声波时差,获得所述间隔数据对应的纵波速度;
根据所述间隔数据的岩性,选择转换公式;
根据所述间隔数据对应的纵波速度与所述转换公式,计算所述间隔数据对应的横波速度;
根据所有间隔数据对应的横波速度,计算最终的横波速度;
当所述间隔数据的岩性为砂岩时,选择克里夫方程为转换方程;
当所述间隔数据的岩性为石灰岩时,选择科斯塔那二次方程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811184200.1A CN111045108B (zh) | 2018-10-11 | 2018-10-11 | 随钻横波计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811184200.1A CN111045108B (zh) | 2018-10-11 | 2018-10-11 | 随钻横波计算方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111045108A CN111045108A (zh) | 2020-04-21 |
CN111045108B true CN111045108B (zh) | 2022-12-13 |
Family
ID=70229204
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811184200.1A Active CN111045108B (zh) | 2018-10-11 | 2018-10-11 | 随钻横波计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111045108B (zh) |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8964503B2 (en) * | 2009-04-28 | 2015-02-24 | Baker Hughes Incorporated | Petrophysics-guided processing of LWD acoustic data |
CN102359924B (zh) * | 2011-09-19 | 2013-03-06 | 中国地质大学(北京) | 一种基于多波地震数据的煤岩强度的检测方法 |
CN103775067B (zh) * | 2012-10-20 | 2016-08-31 | 中国石油化工股份有限公司 | 随钻方位声波测井装置 |
CN104181585A (zh) * | 2014-06-27 | 2014-12-03 | 中国石油化工股份有限公司 | 一种地球物理勘探横波估算方法和系统 |
CN105525917B (zh) * | 2014-09-28 | 2019-03-29 | 中国石油化工股份有限公司 | 识别储层的方法以及识别储层中流体类型的方法 |
CN106443770A (zh) * | 2015-08-07 | 2017-02-22 | 中国石油化工股份有限公司 | 一种页岩气地质甜点的预测方法 |
CN106680878B (zh) * | 2016-12-29 | 2019-03-22 | 中国石油大学(华东) | 一种基于改进Biot系数分析的横波速度估算方法 |
CN106896410B (zh) * | 2017-03-09 | 2019-08-23 | 成都理工大学 | 利用声波测井资料解释岩石的变形模量和脆性指数的方法 |
CN107339099B (zh) * | 2017-07-19 | 2020-06-09 | 中国石油天然气集团公司 | 一种确定储层岩性的方法和装置 |
CN107907910B (zh) * | 2017-11-16 | 2019-12-13 | 山东胜软科技股份有限公司 | 一种不同岩性油藏横波测井确定方法 |
CN108104805B (zh) * | 2017-12-13 | 2021-02-05 | 西南石油大学 | 一种利用测井和岩芯资料获取横波时差剖面的方法 |
-
2018
- 2018-10-11 CN CN201811184200.1A patent/CN111045108B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111045108A (zh) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11353615B2 (en) | Method of and apparatus for carrying out acoustic well logging | |
US10429530B2 (en) | Deghosting with adaptive operators | |
US20170299745A1 (en) | Prestack egs migration method for seismic wave multi-component data | |
EP1660914B1 (en) | Method for attenuating water layer multiples in a seismic survey | |
US20070265785A1 (en) | Interpolation and Extrapolation Method for Seismic Recordings | |
US10670750B2 (en) | Multistage full wavefield inversion process that generates a multiple free data set | |
WO2010139998A2 (en) | An improved process for characterising the evolution of an oil or gas reservoir over time | |
US11579323B2 (en) | Noise attenuation | |
Chen | Borehole wave propagation in isotropic and anisotropic media: Three-dimensional finite difference approach | |
Watremez et al. | Traveltime tomography of a dense wide-angle profile across Orphan Basin | |
Munoz et al. | Automatically tying well logs to seismic data | |
US10459102B2 (en) | Signal detection in semblance methods | |
US20160320507A1 (en) | Time lapse seismic data processing | |
Khoshnavaz et al. | Velocity-independent estimation of kinematic attributes in vertical transverse isotropy media using local slopes and predictive painting | |
Cichostępski et al. | Simultaneous inversion of shallow seismic data for imaging of sulfurized carbonates | |
CN111045108B (zh) | 随钻横波计算方法及系统 | |
CN107918153B (zh) | 一种地震信号相干性高精度检测方法 | |
CN105408772A (zh) | 多分量微震数据的噪声衰减 | |
CN113552624B (zh) | 孔隙度预测方法及装置 | |
US20160047929A1 (en) | Methods and systems that combine wavefields associated with generalized source activation times and near-continuously recorded seismic data | |
US9488745B2 (en) | Method for estimating the water speed of an acoustic node | |
AU2015201786B2 (en) | Methods and systems to separate wavefields using pressure wavefield data | |
Vesnaver et al. | Tomographic joint inversion of direct arrivals, primaries and multiples for monochannel marine surveys | |
US10338253B2 (en) | Method of suppressing spectral artefacts of wavefield decomposition caused by imperfect extrapolation | |
US9829593B2 (en) | Determination of an impulse response at a subsurface image level |
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 |