CN112882097B - 一种大斜度井和水平井的标定方法 - Google Patents
一种大斜度井和水平井的标定方法 Download PDFInfo
- Publication number
- CN112882097B CN112882097B CN202110062180.6A CN202110062180A CN112882097B CN 112882097 B CN112882097 B CN 112882097B CN 202110062180 A CN202110062180 A CN 202110062180A CN 112882097 B CN112882097 B CN 112882097B
- Authority
- CN
- China
- Prior art keywords
- well
- seismic
- synthetic
- channel
- waveform
- 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
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 239000003086 colorant Substances 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims description 19
- 239000004973 liquid crystal related substance Substances 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 5
- 206010028896 Needle track marks Diseases 0.000 abstract description 3
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000003786 synthesis reaction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 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/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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
本发明公开了一种大斜度井和水平井的标定方法,包括以下步骤:S1、通过相控建模模拟井轨迹的地震道波阻抗曲线;S2、对地震道波阻抗曲线进行拟值井合成记录,得到合成道;S3、合成道覆盖于井旁地震剖面对应的地震道;S4、合成地震记录垂直显示在井轨迹附近,且井轨迹标注测量深度与垂直深度。通过相控建模技术将井轨迹所处地震道的波阻抗都模拟出来转换为直井合成记录,所有合成记录显示都是以0作为中心大于零部分(波峰)填充颜色,小于零部分(波谷)直接线条连接,标记测量深度与垂直深度,解决了斜井段合成记录波形失真,倾角越大,失真度越大,尤其是在水平段,无法得到合成记录的问题。
Description
技术领域
本发明涉及地震资料时深标定技术领域,更具体的说是涉及一种大斜度井和水平井的标定方法。
背景技术
目前,合成地震记录是用声波与密度测井或垂直地震剖面资料经过人工合成转换成的地震记录(地震道)。它是地震模型技术中应用非常广泛的一种,也是层位标定、油藏描述等工作的基础,是把地质模型转化为地震信息的中间媒介。通过合成地震记录,可以描述不同的地质层位在地震时间剖面上在位置,以及分析不同岩性组合或流体性质的波形特征。
但是,针对斜井、水平井合成记录,在斜井段,其处理方式与直井相同,但是,斜井段合成记录波形失真,倾角越大,失真度越大,而且一旦进入水平段,没有水平段上下测井资料,无法进行合成记录。
因此,针对斜井或水平井合成记录采用相控建模技术将井轨迹所处地震道的波阻抗都模拟出来转换为直井合成记录,是本领域技术人员亟需解决的问题。
发明内容
有鉴于此,本发明提供了一种大斜度井和水平井的标定方法,通过相控建模技术将井轨迹所处地震道的波阻抗都模拟出来转换为直井合成记录,所有合成记录显示都是以0作为中心大于零部分(波峰)涂黑,小于零(波谷)直接线条连接,标记测量深度与垂直深度,解决了斜井段合成记录波形失真,倾角越大,失真度越大,尤其是在水平段,无法得到合成记录的问题。
为了实现上述目的,本发明采用如下技术方案:
一种大斜度井和水平井的标定方法,包括以下步骤:
S1、通过相控建模模拟井轨迹的地震道波阻抗曲线;
S2、对地震道波阻抗曲线进行拟值井合成记录,得到合成道;
S3、合成道覆盖于井旁地震剖面对应的地震道;
S4、合成地震记录垂直显示在井轨迹附近,且井轨迹标注测量深度与垂直深度。
优选的,所述步骤1具体包括:
S11、建立目的层段地震-测井数据库;
S12、将目的层段采样点序号P置零,总采样点数为M,M>0;
S13、计算未知属性地震波形与数据库中波形的相关性;
S14、选取最大相关系数对应的侧井值作为当前点预测值;
S15、判断目的层段采样点是否达到总采样点数,若否,则令目的层段采样点序号P=P+1后返回步骤S13,若是,则进行步骤S16;
S6、输出地震道波阻抗曲线。
优选的,所述步骤S2具体包括:
S21、标记测量深度与处置深度;
S22、通过声波测井与密度测井,计算出地层反射系数;
S23、通过地震资料振幅普分析或自相关系数,得到地震子波;
S24、地震子波与反射系数褶积,得到合成地震记录:
其中,为合成地震记录,/>为反射系数,/>为地震子波;
S25、以当前位置井轨迹为中心,大于零部分填充颜色,小于零部分线条连接。
优选的,所述步骤S25中所述的大于零的部分具体指代地震道波形的波峰,所述的小于零的部分具体指导地震道波形的波谷。
经由上述的技术方案可知,与现有技术相比,本发明公开提供了一种大斜度井和水平井的标定方法,通过相控建模技术将井轨迹所处地震道的波阻抗都模拟出来转换为直井合成记录,所有合成记录显示都是以0作为中心大于零部分(波峰)填充颜色,小于零部分(波谷)直接线条连接,标记测量深度与垂直深度,解决了斜井段合成记录波形失真,倾角越大,失真度越大,尤其是在水平段,无法得到合成记录的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1附图为本发明提供的流程结构示意图。
图2附图为本发明提供的步骤1的流程结构示意图。
图3附图为本发明提供的直井合成记录示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种大斜度井和水平井的标定方法,包括以下步骤:
S1、通过相控建模模拟井轨迹的地震道波阻抗曲线;
S2、对地震道波阻抗曲线进行拟值井合成记录,得到合成道;
S3、合成道覆盖于井旁地震剖面对应的地震道;
S4、合成地震记录垂直显示在井轨迹附近,且井轨迹标注测量深度与垂直深度。
为进一步优化上述技术方案,步骤1具体包括:
S11、建立目的层段地震-测井数据库;
S12、将目的层段采样点序号P置零,总采样点数为M,M>0;
S13、计算未知属性地震波形与数据库中波形的相关性;
S14、选取最大相关系数对应的侧井值作为当前点预测值;
S15、判断目的层段采样点是否达到总采样点数,若否,则令目的层段采样点序号P=P+1后返回步骤S13,若是,则进行步骤S16;
S6、输出地震道波阻抗曲线。
为进一步优化上述技术方案,步骤S2具体包括:
S21、标记测量深度与处置深度;
S22、通过声波测井与密度测井,计算出地层反射系数;
S23、通过地震资料振幅普分析或自相关系数,得到地震子波;
S24、地震子波与反射系数褶积,得到合成地震记录:
其中,为合成地震记录,/>为反射系数,/>为地震子波;
S25、以当前位置井轨迹为中心,大于零部分填充颜色,小于零部分线条连接。
为进一步优化上述技术方案,步骤S25中所述的大于零的部分具体指代地震道波形的波峰,所述的小于零的部分具体指导地震道波形的波谷。
步骤S11中建立目的层段地震-测井数据库的具体过程为:
S111、设定时窗长度为T1~T2,设T0为时窗中点,T0为目的层段;
S112、选取已知测井T1~T2时窗长度内的地震波形,以及T0点的测井值,建立地震波形与测井的对应关系;
S113、选取已知测井T1~T2时窗长度内的地震波形代表T0点的测井值,构建目的层段地震-测井数据库。
步骤S13中的计算未知属性地震波形与数据库中波形的相关性具体过程为:
将未知属性地震波形中位于T1~T2时窗长度内的地震波形与已知测井的(T1-Δt)~(T2+Δt)段时窗长度内的地震波形比较,滑动计算相关系数并记录最大相关系数Cmax及其对应位置Tm,其中,Δt为允许漂移量,相关系数的计算公式为:
其中,X为已知测井的(T1-Δt)~(T2+Δt)段时窗长度内的地震波形,Y为未知属性地震波形中位于T1~T2时窗长度内的地震波形,N为T1~T2段采样点。
步骤S14包括:根据最大相关系数所对应的位置Tm,获取测井值,并将测井值作为未知属性地震波形中T0处的测井值。
(1)建立地震与测井的非一一对应数据库,分别选取每一口井一定时窗地震波形,与指定时窗指定位置某位置(中点)的测井值作为一组波形与测井的对应关系(图2),如图2所示,ABC分别是三口已钻井,井名下左侧为地震波形,右侧为测井曲线,设时窗长度为T1~T2,T0是T1~T2的中点,是研究目的层,分别选取ABC三口井T1~T2段地震波形代表T0处测井值,作为大数据分析的样本数据库。
(2)计算未知属性地震波形与所有已知波形的相关性,如图2中所示,D处只有地震波形资料,我们需要预测T0时刻的测井物性。将D位置T1~T2波形分别与A、B、C井(T1-△t)~(T2+△t)(其中△t为预先设定的允许漂移量)段波形滑动计算相关系数且记录最大相关系数Cmax及其对应位置Tm。设X为(T1-△t)~(T2+△t)某段时窗长度与(T1~T2)相同已知地震波形,Y为D位置T1~T2段地震波形,T1~T2段采样点为N,相关系数的计算公式为:
(3)筛选最大相关系数,寻找相关系数最大那口井Tm位置处的测井值为D处T0预测测井值。
(4)向下滑动时窗,滑动间隔为一采样点,重复(1)-(3),直到目的层所有地震采样点都实现。
如图3所示,图3A为水平井轨迹及其波阻抗曲线,图3B为将沿井轨迹波阻抗曲线通过相控建模技术模拟的“直井”波阻抗曲线及其合成记录,图3C为只在井轨迹附近显示的合成记录。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (3)
1.一种大斜度井和水平井的标定方法,其特征在于,包括以下步骤:
S1、通过相控建模模拟井轨迹的地震道波阻抗曲线;
S2、对地震道波阻抗曲线进行拟值井合成记录,得到合成道;
S3、合成道覆盖于井旁地震剖面对应的地震道;
S4、合成地震记录垂直显示在井轨迹附近,且井轨迹标注测量深度与垂直深度;
所述步骤1具体包括:
S11、建立目的层段地震-测井数据库;
S12、将目的层段采样点序号P置零,总采样点数为M,M>0;
S13、计算未知属性地震波形与数据库中波形的相关性;
S14、选取最大相关系数对应的侧井值作为当前点预测值;
S15、判断目的层段采样点是否达到总采样点数,若否,则令目的层段采样点序号P=P+1后返回步骤S13,若是,则进行步骤S16;
S16、输出地震道波阻抗曲线。
2.根据权利要求1所述的一种大斜度井和水平井的标定方法,其特征在于,所述步骤S2具体包括:
S21、标记测量深度与处置深度;
S22、通过声波测井与密度测井,计算出地层反射系数;
S23、通过地震资料振幅普分析或自相关系数,得到地震子波;
S24、地震子波与反射系数褶积,得到合成地震记录:
其中,为合成地震记录,/>为反射系数,/>为地震子波;
S25、以当前位置井轨迹为中心,大于零部分填充颜色,小于零部分线条连接。
3.根据权利要求2所述的一种大斜度井和水平井的标定方法,其特征在于,所述步骤S25中所述的大于零的部分具体指代地震道波形的波峰,所述的小于零的部分具体指导地震道波形的波谷。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110062180.6A CN112882097B (zh) | 2021-01-18 | 2021-01-18 | 一种大斜度井和水平井的标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110062180.6A CN112882097B (zh) | 2021-01-18 | 2021-01-18 | 一种大斜度井和水平井的标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112882097A CN112882097A (zh) | 2021-06-01 |
CN112882097B true CN112882097B (zh) | 2023-07-25 |
Family
ID=76048987
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110062180.6A Active CN112882097B (zh) | 2021-01-18 | 2021-01-18 | 一种大斜度井和水平井的标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112882097B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117523035B (zh) * | 2024-01-05 | 2024-03-26 | 中国石油集团东方地球物理勘探有限责任公司 | 水平井剖面中曲线变密填充方法、装置、设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2029016A (en) * | 1978-07-24 | 1980-03-12 | Schlumberger Ltd | Vertical seismic exploration and profiling technique |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100557464C (zh) * | 2007-01-15 | 2009-11-04 | 中国石油大港油田勘探开发研究院 | 基于叠前波场模拟的地震勘探层位标定方法 |
CN104656142B (zh) * | 2013-11-19 | 2017-06-06 | 中国石油天然气集团公司 | 一种利用垂直地震剖面与测井联合的地震层位标定方法 |
CN111175820B (zh) * | 2020-01-10 | 2022-08-30 | 杨林海 | 一种合成地震记录清晰断点标定确定时深关系的分析方法 |
CN111474578A (zh) * | 2020-04-24 | 2020-07-31 | 成都捷科思石油天然气技术发展有限公司 | 一种利用模拟测井资料进行随钻跟踪的方法 |
CN111897004B (zh) * | 2020-08-31 | 2023-03-24 | 北京瑞码恒杰科技有限公司 | 一种基于大数据分析技术的测井预测方法 |
-
2021
- 2021-01-18 CN CN202110062180.6A patent/CN112882097B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2029016A (en) * | 1978-07-24 | 1980-03-12 | Schlumberger Ltd | Vertical seismic exploration and profiling technique |
Non-Patent Citations (1)
Title |
---|
合成地震记录制作原理与方法;杜萌,等;山东化工;第45卷(第17期);第159-161页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112882097A (zh) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108802812A (zh) | 一种井震融合的地层岩性反演方法 | |
US8531914B2 (en) | Method of imaging a target area of the subsoil from walkaway type data | |
US7525874B2 (en) | Method for converting seismic data from the time domain to the depth domain | |
CN107121699A (zh) | 一种地震相控制下的沉积微相识别方法 | |
US7768870B2 (en) | Method for adjusting a seismic wave velocity model according to information recorded in wells | |
MXPA05013185A (es) | Metodos y sistemas para el procedimiento de ondas acusticas. | |
CN108562938B (zh) | 一种消除频散效应的方法、装置及系统 | |
CN112882097B (zh) | 一种大斜度井和水平井的标定方法 | |
CN106772570A (zh) | 火成岩开启缝地震预测方法 | |
CN111487692A (zh) | 一种盐间页岩油韵律层地震响应特征及储层厚度预测方法 | |
Bell | AAPG Memoir 76, Chapter 18: Velocity Estimation for Pore-Pressure Prediction | |
US4817059A (en) | Borehole logging methods for detection and imaging of formation structural features | |
CN108375785A (zh) | 裂缝带位置校正方法及装置 | |
US4218766A (en) | Method of seismic wave amplitude normalization | |
CN107238859B (zh) | 计算虚拟井合成地震记录的方法及系统 | |
CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
CN109324344A (zh) | 基于纯纵波和拟声波反演的页岩厚度预测方法及系统 | |
CN108897039B (zh) | 基于层位控制的dtw多波匹配方法 | |
CN113589365B (zh) | 基于时频域信息的储层尖灭线描述方法 | |
CN116500679A (zh) | 一种井震标定方法、装置和相关设备 | |
CN110568490B (zh) | 一种高速层顶薄储层的识别方法 | |
CN110673211B (zh) | 一种基于测井与地震数据的品质因子建模方法 | |
CN112505767A (zh) | 深度域层位的标定方法 | |
CN111897004A (zh) | 一种基于大数据分析技术的测井预测方法 | |
CN112198554B (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 |