CN109001815B - 基于地震和井径约束的声波曲线校正方法 - Google Patents

基于地震和井径约束的声波曲线校正方法 Download PDF

Info

Publication number
CN109001815B
CN109001815B CN201810984266.2A CN201810984266A CN109001815B CN 109001815 B CN109001815 B CN 109001815B CN 201810984266 A CN201810984266 A CN 201810984266A CN 109001815 B CN109001815 B CN 109001815B
Authority
CN
China
Prior art keywords
seismic
curve
acoustic
time difference
calculating
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
Application number
CN201810984266.2A
Other languages
English (en)
Other versions
CN109001815A (zh
Inventor
王侠
石好果
徐佑德
任新成
程长领
刘德志
刘博奥
刘婧
刘振阳
翟晓薇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petrochemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Original Assignee
China Petrochemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petrochemical Corp, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petrochemical Corp
Priority to CN201810984266.2A priority Critical patent/CN109001815B/zh
Publication of CN109001815A publication Critical patent/CN109001815A/zh
Application granted granted Critical
Publication of CN109001815B publication Critical patent/CN109001815B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48—Processing data
    • G01V1/50—Analysing data

Abstract

本发明提供一种基于地震和井径约束的声波曲线校正方法，该基于地震和井径约束的声波曲线校正方法包括：步骤1，进行合成地震记录标定，确定两个包含井壁垮塌段的地震反射轴a和地震反射轴b，并计算地震反射轴a、b间的地震双程时差Δt1；步骤2，计算合成地震记录上与地震反射轴a、b对应的反射轴间的双程时差Δt2；步骤3，分别计算原始声波曲线及校正后的声波曲线与地震反射轴a、b对应的深度段间的值的总和Σac、ΣAC；步骤4，计算校正系数K，进行声波时差曲线的校正。该基于地震和井径约束的声波曲线校正方法能够准确的对井壁垮塌段的声波曲线进行校正，而且操作相对简单、便于实现。

Description

基于地震和井径约束的声波曲线校正方法
技术领域
本发明涉及油田勘探技术领域，特别是涉及到一种基于地震和井径约束的声波曲线校正方法。
背景技术
井壁垮塌(扩径)是导致测井过程中声波时差曲线出现异常的重要原因。目前针对井壁垮塌进行校正的方法主要有两类：
第一类是基于声波测井的基本原理出发,对声波测井仪器经过井径扩径段时声波到达接收器的射线路径进行分析,找出产生声波时差异常的原因，从而进行校正。
第二类是分析统计校正法。这种方法利用选取的井径状况良好的标准井或者井径状况良好的某一层段，统计不同岩性的声波时差，然后对垮塌井段的声波曲线进行校正。
第一类从测井原理和方法出发，校正结果相对准确一些，但是专业性比较强，普通地质人员难以操作和实施。第二类原理简单，容易操作和实现，但是由于涉及大量的统计和分析，工作量较大，此外，如果垮塌井段若出现特殊岩性，则容易出现误差。为此我们发明了一种基于地震和井径约束的声波曲线校正方法，解决了以上技术问题。
发明内容
本发明的目的是提供一种能够准确的对井壁垮塌段的声波曲线进行校正，而且操作相对简单、便于实现的基于地震和井径约束的声波曲线校正方法。
本发明的目的可通过如下技术措施来实现：基于地震和井径约束的声波曲线校正方法，该基于地震和井径约束的声波曲线校正方法包括：步骤1，进行合成地震记录标定，确定两个包含井壁垮塌段的地震反射轴a和地震反射轴b，并计算地震反射轴a、b间的地震双程时差Δt1；步骤2，计算合成地震记录上与地震反射轴a、b对应的反射轴间的双程时差Δt2；步骤3，分别计算原始声波曲线及校正后的声波曲线与地震反射轴a、b对应的深度段间的值的总和Σac、ΣAC；步骤4，计算校正系数K，进行声波时差曲线的校正。
本发明的目的还可通过如下技术措施来实现：
在步骤1中，利用原始声波曲线进行合成地震记录标定，在过井地震剖面上选取两个能够包含井壁垮塌段且地震反射特征容易识别的地质界面所对应的地震反射轴a和地震反射轴b，并计算地震反射轴a、b间的地震双程时差Δt1。
在步骤2中，利用原始声波曲线所做的合成地震记录标定的结果，计算合成地震记录上与步骤1中选取的两个地质界面对应的地震反射轴a’、b’间的双程时差Δt2。
在步骤3中，计算原始声波曲线与步骤1选定的两个地质界面对应的深度段间的声波时差值的总和Σac,然后根据公式1计算校正后与步骤1选定的两个地质界面对应的深度段间的声波时差总和ΣAC，
式中：Δt1-选定的两个地质界面对应的地震反射轴a、b间的地震双程时差，单位为毫秒；
Δt2-合成地震记录上与选定的两个地质界面对应的地震反射轴a’、b’间的双程时差，单位为毫秒；
ac-原始声波时差；
AC-校正后的声波时差；
Σac-原始声波曲线与选定的两个地质界面对应的深度段间的声波值的总和；
ΣAC-与选定的两个地质界面对应的深度段间的校正后的声波时差值的总和。
在步骤4中，根据需要校正井段的井径基线确定D0值，D0值等于井径基线值，D0值的获取为在测井综合柱状图的井径曲线中进行读取。
在步骤4中，确定D0值后，依据公式2计算校正系数K：
式中：D0为井径基线对应的井径值，单位为厘米；
D为井径值，单位为厘米；
K为校正系数，无单位；
ac-原始声波时差；
AC-校正后的声波时差；
Σac-原始声波曲线与选定的两个地质界面对应的深度段间的声波值的总和；
ΣAC-与选定的两个地质界面对应的深度段间的校正后的声波时差值的总和。
在步骤4中，在求得校正系数K的基础上，根据公式3对原始声波时差曲线进行校正：
本发明的基于地震和井径约束的声波曲线校正方法，涉及到由于井径的变化引起的声波曲线异常的校正方法，能够准确的对井壁垮塌段的声波曲线进行校正，从根本上解决了由于井壁垮塌引起的声波时差曲线异常造成合成记录标定不准确的问题。相对于基于声波测井原理出发,对声波测井仪器经过井径扩径段时声波到达接收器的射线路径进行分析和校正的方法而言，该方法的原理简单，同时，该方法是以合成地震记录标定为基础的，因此更易于地震解释人员理解和操作。与统计分析校正法相比较，该方法省去了大量的统计工作，而且校正更加准确和快速。
附图说明
图1为本发明的基于地震和井径约束的声波曲线校正方法的一具体实施例的流程图；
图2为本发明的一具体实施例中董701井综合柱状图；
图3为本发明的一具体实施例中董701原始曲线合成地震记录标定的示意图；
图4为本发明的一具体实施例中董701过井地震剖面的示意图；
图5为本发明的一具体实施例中董701校正后合成地震记录标定的示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合所附图式，作详细说明如下。
如图1所示，图1为本发明的基于地震和井径约束的声波曲线校正方法的一具体实施例的流程图。
在步骤101，首先利用原始声波曲线进行合成地震记录标定。在过井地震剖面上选取两个能够包含井壁垮塌段且地震反射特征容易识别的地质界面所对应的地震反射轴a和地震反射轴b，并计算地震反射轴a、b间的地震双程时差Δt1。流程进入到步骤102。
在步骤102，利用原始声波曲线所做的合成地震记录标定的结果，计算合成地震记录上与步骤101中选取的两个地质界面对应的地震反射轴a’、b’间的双程时差Δt2。流程进入到步骤103。
在步骤103，计算原始声波曲线与步骤1选定的两个地质界面对应的深度段间的声波时差值的总和Σac,然后根据公式1计算校正后与步骤1选定的两个地质界面对应的深度段间的声波时差总和ΣAC。
式中：Δt1选定的两个地质界面对应的地震反射轴a、b间的地震双程时差，单位为毫秒。
Δt2合成地震记录上与选定的两个地质界面对应的地震反射轴a’
、b’间的双程时差，单位为毫秒。。
ac原始声波时差
AC校正后的声波时差
Σac原始声波曲线与选定的两个地质界面对应的深度段间的声波值的总和。
ΣAC与选定的两个地质界面对应的深度段间的校正后的声波时差值的总和。
流程进入到步骤104。
在步骤104，首先根据需要校正井段的井径基线确定D0值，D0值等于井径基线值，D0值的获取可在测井综合柱状图的井径曲线中进行读取。确定D0值后，依据公式2计算校正系数K。
式中：D0为井径基线对应的井径值(单位为厘米)。
D为井径值(单位为厘米)。
K为校正系数，无单位。
在求得校正系数K的基础上，根据公式3对原始声波时差曲线进行校正。
在应用本发明的一具体实施例中，准格尔盆地中部4区块董2井北三维区，区内共有探井8口，测井资料齐全，三维地震资料品质较好，地震反射与地质层位对应关系明确。工区内地震反射标志层为白垩系底强反射及西山窑组煤层强反射。
根据董701井原始声波曲线所做的合成地震记录标定，在对准西山窑组煤层反射时，合成地震记录上白垩系底面位置与区域上白垩系底面反射不对应。此外，齐古组油层(J3q)在地震上为一套强反射(图4)，但合成地震记录却标定在弱反射上(图3)，这与实际情况不吻合。针对以上情况，认为井壁垮塌(图2)引起的声波曲线异常是产生以上问题的主要原因，因此对该井的声波曲线进行校正，包括了以下步骤：
在步骤1中，准格尔盆地中部4区块董2井北三维区董701井进行合成地震记录标定(图3)，并确定地震反射特征明显的白垩系底及西山窑组煤层两个强反射分别为a和b(图4)。a、b时差为Δt1＝480毫秒。
在步骤2中，根据合成地震记录标定的结果，合成地震记录上与白垩系底及西山窑组煤层对应的地震反射为a’、b’，其地震双程时差为Δt2＝504毫秒。(图3)。
在步骤3中，对原始声波时差曲线白垩系底面至西山窑组煤层反射对应深度段3775米-4820米，求和为Σac＝2012289；Δt1＝480，Δt2＝504，根据公式1，可以求得校正后的对应深度段声波时差的和为ΣAC＝1916456.71。
在步骤4中，首先根据井径曲线确定井径的基线值12.4英寸(图2)，换算成厘米制，D0＝31.5厘米；然后对白垩系底面至西山窑组煤层反射对应深度段3775米-4820米深度段每个采样点的井径和井径基线之差与声波时差值相乘并累计求和，结果为7866374.12；最后，根据步骤3可求得校正前、后白垩系底面至西山窑组煤层反射对应深度段3775米-4820米间的声波时差累计求和后的差为95832.29，根据公式2可以求得K＝82.08；把K值带入公式3，即可求得校正后每个采样点声波时差值，也就完成了对声波时差的校正。
用该方法对井曲线进行校正，其基于地震和井径约束的思想体现在两个方面：其一，整条曲线的校正量大小是由地震标志层反射及合成地震记录对应反射间的双程时差决定的(公式1)；其二，曲线上每个采样点的校正量大小与则与井壁垮塌的大小成正比例(公式3)。
从图2中可以看到，在侏罗系顶部，井壁垮塌比较严重，井径曲线与井径基线偏差很大，因此这一垮塌段校正前后的声波曲线差异较大，而白垩系砂岩地层几乎没有垮塌，校正前后二者几乎没有差别。通过校正前后的合成地震记录标定对比(图3与图5)来看，校正后不仅白垩系底及西山窑组煤层与地震反射对应性好，而且校正段的整个合成地震记录与地震剖面对应关系良好，齐古组的油层对应关系也有了非常明显的改善，标定结果既符合区域上的地质认识，对局部砂体的标定也非常准确。校正前，合成地震记录与过井地震的相关度只有32％，校正后达到了71％，可以看出，校正效果是非常明显的。

Claims (1)

1.一种基于地震和井径约束的声波曲线校正方法，其特征在于：包括以下步骤：
步骤1，利用原始声波曲线进行合成地震记录标定，在过井地震剖面上选取两个能够包含井壁垮塌段且地震反射特征容易识别的地质界面所对应的地震反射轴a和地震反射轴b，并计算地震反射轴a、b间的地震双程时差Δt1；
步骤2，利用原始声波曲线所做的合成地震记录标定的结果，计算合成地震记录上与步骤1中选取的两个地质界面对应的地震反射轴a’、b’间的双程时差Δt2；
步骤3，分别计算原始声波曲线和校正后的声波曲线与地震反射轴a、b对应的深度段间的声波时差值的总和Σac、ΣAC；
步骤4，计算校正系数K，进行声波时差曲线的校正；
在步骤3中，计算原始声波曲线与步骤1选定的两个地质界面对应的深度段间的声波时差值的总和Σac,然后根据公式1计算校正后的声波曲线与步骤1选定的两个地质界面对应的深度段间的声波时差值总和ΣAC，
ac-原始声波时差；
AC-校正后的声波时差；
在步骤4中，根据需要校正井段的井径基线确定D0值，D0值为井径基线值，D0值通过在测井综合柱状图的井径曲线中读取；
确定D0值后，依据公式2计算校正系数K：
式中：D为井径值，单位为厘米；
K为校正系数，无单位；
在求得校正系数K的基础上，根据公式3对原始声波时差曲线进行校正：
CN201810984266.2A 2018-08-27 2018-08-27 基于地震和井径约束的声波曲线校正方法 Active CN109001815B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810984266.2A CN109001815B (zh) 2018-08-27 2018-08-27 基于地震和井径约束的声波曲线校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810984266.2A CN109001815B (zh) 2018-08-27 2018-08-27 基于地震和井径约束的声波曲线校正方法

Publications (2)

Publication Number Publication Date
CN109001815A CN109001815A (zh) 2018-12-14
CN109001815B true CN109001815B (zh) 2019-12-31

Family

ID=64593286

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810984266.2A Active CN109001815B (zh) 2018-08-27 2018-08-27 基于地震和井径约束的声波曲线校正方法

Country Status (1)

Country Link
CN (1) CN109001815B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111060986B (zh) * 2019-10-18 2022-01-28 中国石油化工股份有限公司 一种地层压力预测方法及岩性油藏评价方法
CN111650646B (zh) * 2020-05-26 2023-05-12 中国石油天然气集团有限公司 一种用于地震处理中的声波曲线校正方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102454399A (zh) * 2010-10-26 2012-05-16 中国石油化工股份有限公司 测井声波时差信号校正方法
CN103573250A (zh) * 2013-07-22 2014-02-12 中国石油天然气股份有限公司 一种计算水平井井眼到地层上下界面距离的方法
CN107065011A (zh) * 2017-06-22 2017-08-18 东北石油大学 一种应用于陆相盆地储层反演的曲线频率融合方法
CN107092031A (zh) * 2017-06-01 2017-08-25 中国石油天然气集团公司 一种声波测井数据校正方法及装置
CN107255834A (zh) * 2017-05-19 2017-10-17 中国石油集团川庆钻探工程有限公司 一种基于地震约束的声波测井曲线校正方法
CN108037532A (zh) * 2017-12-06 2018-05-15 中国科学院地质与地球物理研究所 用于标定合成地震记录的方法、装置、系统和计算机可读介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102454399A (zh) * 2010-10-26 2012-05-16 中国石油化工股份有限公司 测井声波时差信号校正方法
CN103573250A (zh) * 2013-07-22 2014-02-12 中国石油天然气股份有限公司 一种计算水平井井眼到地层上下界面距离的方法
CN107255834A (zh) * 2017-05-19 2017-10-17 中国石油集团川庆钻探工程有限公司 一种基于地震约束的声波测井曲线校正方法
CN107092031A (zh) * 2017-06-01 2017-08-25 中国石油天然气集团公司 一种声波测井数据校正方法及装置
CN107065011A (zh) * 2017-06-22 2017-08-18 东北石油大学 一种应用于陆相盆地储层反演的曲线频率融合方法
CN108037532A (zh) * 2017-12-06 2018-05-15 中国科学院地质与地球物理研究所 用于标定合成地震记录的方法、装置、系统和计算机可读介质

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
严重垮塌井眼条件下测井资料的校正方法;孙中春 等;《新疆石油地质》;20061031;第27卷(第5期);第559-561页 *
井眼垮塌条件下声波和密度曲线的校正方法;吉伟平 等;《国外测井技术》;20150831(第4期);第32-35页 *
基于井径变化的声波测井时差校正;刘浩杰 等;《石油物探》;20061130;第45卷(第6期);第638-641页 *

Also Published As

Publication number Publication date
CN109001815A (zh) 2018-12-14

Similar Documents

Publication Publication Date Title
CN106842313B (zh) 基于方位叠前地震数据的各向异性参数反演方法
AU2012385250B2 (en) Anisotropy parameter estimation
CN106405651B (zh) 一种基于测井匹配的全波形反演初始速度模型构建方法
Maxwell et al. Anisotropic velocity modeling for microseismic processing: Part 1—Impact of velocity model uncertainty
US20150109885A1 (en) Method for correcting first break arrival time
CN105277982A (zh) 一种泥页岩总有机碳含量地震预测方法
CN109001815B (zh) 基于地震和井径约束的声波曲线校正方法
US4399693A (en) Applications of borehole gravimetric techniques to determine residual oil saturation
WO2020215170A1 (zh) 基于地震岩石物理实验分析的测井与地震速度匹配方法
CN109057784A (zh) 利用岩石切削强度快速确定岩体普通强度参数的方法
CN111650646B (zh) 一种用于地震处理中的声波曲线校正方法
CN113568046B (zh) 裂缝方位约束的多维信息裂缝表征方法
CN111208564A (zh) 一种深度域层位标定方法及装置
CN111308558B (zh) 页岩气水平井纵波时差校正方法
CN106761666A (zh) 一种四探头散射伽马测井及非线性数据反演的方法及装置
CN111089904B (zh) 一种考虑地层特性的室内测量波速围压、频散校正方法
CN106761667B (zh) 一种四探头散射伽马测井及线性数据反演的方法及装置
CN108260359A (zh) 一种地震信号处理方法、装置和系统
CN112925022B (zh) 页岩vti介质各向异性参数的预测方法
CN110850504B (zh) 一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法
CN112666610B (zh) 地层有机碳含量计算方法及系统
Grilli et al. Sonic Log Quality Control: A Mandatory Workflow to Ensure Reliable Results in Different Disciplines
Valero et al. Radial profilings integration for optimal well completion design
Onayemi et al. Rock physics templates modeling for time-lapse seismic analysis of reservoir M1 in Moremi Field, Offshore Niger Delta Basin
CN115903015A (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
GR01 Patent grant
GR01 Patent grant