CN104375172B - 火山岩下伏地层的构造形态正确成像的方法 - Google Patents
火山岩下伏地层的构造形态正确成像的方法 Download PDFInfo
- Publication number
- CN104375172B CN104375172B CN201310356114.5A CN201310356114A CN104375172B CN 104375172 B CN104375172 B CN 104375172B CN 201310356114 A CN201310356114 A CN 201310356114A CN 104375172 B CN104375172 B CN 104375172B
- Authority
- CN
- China
- Prior art keywords
- volcanic
- depth
- volcanic rock
- speed
- rock
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
Abstract
本发明是地震数据处理中火山岩下伏地层的构造形态正确成像的方法。利用火山岩顶、底层位,建立叠前深度偏移沿层速度模型,在深度偏移后的目标线上沿真实的地质层位追踪火山岩顶层位深度Ht和火山岩底层位深度Hb;然后在火山岩底层位深度Hb基础上,根据钻井资料、地质资料确定预期的火山岩底界层位深度Heb;利用(Heb‑Ht)与(Hb‑Ht)两者的厚度之比求得火山岩的层速度系数;将火山岩层速度系数与火山岩层速度相乘,得到深度偏移速度模型中火山岩层新速度,将火山岩层的新速度填充到深度速度模型中,进行目标线深度偏移,最终使火山岩下伏地层的构造形态正确成像。本发明可以解决二叠系火山岩发育区下伏地层随火山岩厚度变化起伏问题。
Description
技术领域
本发明属于地震数据处理技术,具体是一种火山岩下伏地层的构造形态正确成像的方法。
背景技术
由于火山岩速度与厚度的空间变化,时间偏移得到的偏移成果不能正确描述火山岩下伏地层的构造形态,通常的解决方案是叠前深度偏移,而对于火山岩速度和厚度的准确成像则是能否解决问题的关键。
目前在深度偏移速度建模中通常应用的方法主要有两种:一是首先在时间域的偏移剖面上进行速度层位的解释;然后在井速度趋势的约束下,将叠前时间偏移的均方根速度转换为层速度,建立初始速度模型;在深度偏移处理道集校平、剩余速度归零准则的约束下,迭代修正火山岩的层速度,得到最终速度场。第二种方法是不考虑速度层位的约束,整个速度场的速度是点点变化的,通常被称为网格层析速度反演技术;该技术首先需要一个初始速度模型(可以有叠前时间偏移的均方根速度转换得到,也可以采用第一种方法中得到的初始或最终速度模型),然后通过网格层析速度反演迭代技术获得最终速度模型。
在实际应用中,前面两种方法都不能很好的校正火山岩下伏地层跟随火山岩厚度、速度变化的问题。第一种速度建模方法人为因素较多,在进行火山岩层位拾取时,如果不能很精确的拾取火山岩顶面和底面的变化(火山岩顶面通常起伏剧烈,有些地段识别困难),仅靠“道集校平、剩余速度归零”准则约束速度建模,是不能准确给出火山岩速度的空间变化的,那么下伏地层的起伏是不能消除的。第二种方法虽然减少了人为因素的影响,但在应用中最终速度模型对初始速度趋势的依赖性很大;即使通过网格层析的多次迭代,也很难改变初始速度模型的趋势。综上所述,对地层沿层分布、厚度不稳定、岩性多变的高速火山岩,即使地震资料处理中采用了叠前深度偏移处理技术,也不能使下伏地层的构造形态正确成像问题依然很难解决。
发明内容
本发明的目的就是提供一种能与钻井吻合的火山岩下伏地层的构造形态正确成像的方法。
本发明通过以下步骤实现:
1)采集地震数据,预处理得到叠前时间偏移剖面,在叠前时间偏移剖面上以较密的间隔追踪火山岩顶、底层位;
所述的顶、底层位要求与火山岩顶、底界面的起伏保持一致。
在火山岩顶、底层位解释的过程中,如附图1中深黑色和浅灰色所示火山岩顶、底层位解释都是不符合要求的,只有白色是符合火山岩顶、底层位解释的要求。
2)利用火山岩顶、底层位,按照常规叠前深度偏移速度的方法建立火山岩层速度的叠前深度偏移沿层速度模型;
3)用步骤2)得到的速度模型进行目标线的Kirchhoff积分法深度偏移;
所述的目标线的密度要求能控制火山岩界面的起伏变化。
4)首先在深度偏移后的目标线上沿真实的地质层位追踪火山岩顶层位深度Ht和火山岩底层位深度Hb;
所述的火山岩底层位深度Hb反映火山岩底界的起伏变化;
然后在火山岩底层位深度Hb基础上,根据钻井资料、地质资料确定预期的火山岩底界层位深度Heb;
所述的确定预期的火山岩底界层位深度Heb是:
(1)在火山岩变厚并且下伏地层形态扭曲的部位重新拾取,其它部位与原层位保持一致(如附图2中(a)所示);
(2)遇到断层时,保持断层断距大小变化(如附图2中(b)所示)。
5)采用以下公式校正火山岩层速度:
Vc=(Heb-Ht)/(Hb-Ht)
Vc:火山岩层速度系数,Ht:火山岩顶层位深度(单位:m);
Hb:火山岩底原层位深度(单位:m),Heb:火山岩底预期层位深度(单位:m);
将火山岩层速度系数与步骤2)中得到的火山岩层速度相乘,得到深度偏移速度模型中火山岩层新速度;
6)将火山岩层的新速度填充到步骤2)中的深度速度模型中,进行目标线深度偏移,当偏移结果满足条件,火山岩速度迭代结束,反之,重新从步骤4)开始进行迭代,其迭代流程见附图3,满足条件时迭代结束,最终使火山岩下伏地层的构造形态正确成像。
所述的偏移结果满足条件是:
a、道集拉平、剩余谱归零;
b、下伏地层构造形态符合地质认识,并且与钻井资料吻合;
c、奥陶系碳酸盐岩缝洞成像精度得到提高。
本发明可以解决二叠系火山岩发育区下伏地层随火山岩厚度变化起伏问题。以往深度偏移成果在火山岩厚度大的区域志留系底界面(埋藏深度5400米)与钻井深度误差大(约90米左右),而火山岩薄的区域志留系底界面与钻井误差小(约50米左右),验证钻井误差达到40米以内。本发明在消除火山岩下伏地层构造形态影响的同时,还使深部奥陶系溶洞反射的归位更加收敛,从而进一步说明了正确性和可行性。
附图说明
图1时间域火山岩顶、底层位解释实例;
图2火山岩底原层位与预期层位解释,(a)在火山岩变厚并且下伏地层形态扭曲的部位重新拾取,其它部位与原层位保持一致,(b)遇到断层时,保持断层断距大小变化的合理性;
图3准确描述火山岩层速度迭代流程图;
图4时间偏移剖面:二叠系火成岩分布不均、厚度、岩性、速度横向变化;
图5所示叠前深度偏移基本流程图;
图6火山岩速度校正前后的平面图及剖面对比;
图7准确描述火山岩层速度前深度偏剖面;
图8准确描述火山岩层速度后深度偏剖面。
具体实施方式
本发明具体实施方式如下:
1)首先经过野外地震采集获得原始地震数据,经过预处理以及kirchhoff积分法偏移得到叠前时间偏移剖面,在叠前时间偏移剖面上以能够控制构造变化的间隔追踪火山岩顶、底层位;
所述的顶、底层位要求与火山岩顶、底界面的起伏保持一致。
本发明试验实例是奥陶系凸起的延伸,区域位置十分有利。全区发育一套二叠系的火成岩,分布不均、厚度、岩性、速度横向变化剧烈(如附图4中所示),主要目的层受上覆二叠系火成岩影响较大,火成岩速度和厚度纵横向上的剧烈变化影响了其下伏地层低幅度构造的准确落实,以石炭系东河砂岩为例,T0图表现的低幅度构造经时深转换后低幅度构造形态发生了变化,说明研究火成岩速度场的变化对下伏地层构造的准确落实非常重要。
本发明在该区叠前时间偏移剖面上追踪火山岩顶、底层位,在火山岩顶、底层位解释的过程中,如附图1中深黑色和浅灰色所示火山岩顶、底层位解释都是不符合要求的,只有白色是符合火山岩顶、底层位解释的要求。
2)按照附图5所示叠前深度偏移基本流程,建立叠前深度偏移沿层速度模型,该速度模型利用火山岩顶、底层位控制火山岩层速度,如附图6中所示是火山岩速度校正前层速度剖面。
3)用步骤2)得到的速度模型进行目标线的Kirchhoff积分法深度偏移;
所述的目标线的密度要求能控制火山岩界面的起伏变化。
4)首先在深度偏移后的目标线上沿真实的地质层位追踪火山岩顶层位深度Ht和火山岩底层位深度Hb;
所述的火山岩底层位深度Hb反映火山岩底界的起伏变化(如附图2中原层位解释所示);
然后在火山岩底层位深度Hb基础上,根据钻井资料、地质资料确定预期的火山岩底界层位深度Heb;
所述的确定预期的火山岩底界层位深度Heb是:
(1)在火山岩变厚并且下伏地层形态扭曲的部位重新拾取,其它部位与原层位保持一致(如附图2中(a)所示);
(2)遇到断层时,保持断层断距大小变化(如附图2中(b)所示)。
5)采用以下公式校正火山岩层速度:
Vc=(Heb-Ht)/(Hb-Ht)
Vc:火山岩层速度系数,Ht:火山岩顶层位深度(单位:m);
Hb:火山岩底原层位深度(单位:m),Heb:火山岩底预期层位深度(单位:m);
将火山岩层速度系数与步骤2)中得到的火山岩层速度相乘,得到深度偏移速度模型中火山岩层新速度;
6)将火山岩层的新速度填充到步骤2)中的深度速度模型中,进行目标线深度偏移,当偏移结果满足条件,火山岩速度迭代结束,反之,重新从步骤4)开始进行迭代,其迭代流程见附图3,满足条件时迭代结束,最终使火山岩下伏地层的构造形态正确成像。
所述的偏移结果满足条件是:
a、道集拉平、剩余谱归零;
b、下伏地层构造形态符合地质认识,并且与钻井资料吻合;
c、奥陶系碳酸盐岩缝洞成像精度得到提高。
在实验区中,按照如前所述步骤4)~6)的火山岩速度校正方法并通过迭代,可以得到如附图6所示的火山岩新速度,将该速度填充到整个深度-速度模型中,利用新的深度-速度模型进行叠前深度偏移。
如图7和8分别为该项技术应用前后的叠前深度偏移剖面,通过对比剖面偏移效果可以看出,应用该项技术之后的深度偏移结果火山岩下伏地层构造形态合理,并且奥陶系碳酸盐岩缝洞的成像效果得到改善。
Claims (5)
1.一种火山岩下伏地层的构造形态正确成像的方法,特点是通过以下步骤实现:
1)采集地震数据,预处理得到叠前时间偏移剖面,在叠前时间偏移剖面上以较密的间隔追踪火山岩顶、底层位;
2)利用火山岩顶、底层位,按照常规叠前深度偏移速度的方法建立火山岩层速度的叠前深度偏移沿层速度模型;
3)用步骤2)得到的速度模型进行目标线的Kirchhoff积分法深度偏移;
4)首先在深度偏移后的目标线上沿真实的地质层位追踪火山岩顶层位深度Ht和火山岩底原层位深度Hb;然后在火山岩底原层位深度Hb基础上,根据钻井资料、地质资料确定预期的火山岩底界层位深度Heb;
5)采用以下公式校正火山岩层速度:
Vc=(Heb-Ht)/(Hb-Ht)
Vc:火山岩层速度系数,Ht:火山岩顶层位深度;
Hb:火山岩底原层位深度,Heb:火山岩底预期层位深度;
将火山岩层速度系数与步骤2)中得到的火山岩层速度相乘,得到深度偏移速度模型中火山岩层新速度;
6)将火山岩层的新速度填充到步骤2)中的叠前深度偏移沿层速度模型中,进行目标线深度偏移,当偏移结果满足条件,火山岩层速度迭代结束,反之,重新从步骤4)开始进行迭代,直到满足条件时迭代结束,最终使火山岩下伏地层的构造形态正确成像,所述的偏移结果满足条件是:a、道集拉平或剩余谱归零;b、下伏地层构造形态符合地质认识,并且与钻井资料吻合;c、奥陶系碳酸盐岩缝洞成像精度得到提高。
2.根据权利要求1的方法,特点是步骤1)所述的顶、底层位与火山岩顶、底界面的起伏保持一致。
3.根据权利要求1的方法,特点是步骤3)所述的目标线的密度要求能控制火山岩界面的起伏变化。
4.根据权利要求1的方法,特点是步骤4)所述的火山岩底原层位深度Hb反映火山岩底界的起伏变化。
5.根据权利要求1的方法,特点是步骤4)所述的确定预期的火山岩底界层位深度Heb是:
(1)在火山岩变厚并且下伏地层形态扭曲的部位重新拾取,其它部位与原层位保持一致;
(2)遇到断层时,保持断层断距大小变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310356114.5A CN104375172B (zh) | 2013-08-15 | 2013-08-15 | 火山岩下伏地层的构造形态正确成像的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310356114.5A CN104375172B (zh) | 2013-08-15 | 2013-08-15 | 火山岩下伏地层的构造形态正确成像的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104375172A CN104375172A (zh) | 2015-02-25 |
CN104375172B true CN104375172B (zh) | 2017-02-15 |
Family
ID=52554200
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310356114.5A Active CN104375172B (zh) | 2013-08-15 | 2013-08-15 | 火山岩下伏地层的构造形态正确成像的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104375172B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106896421B (zh) * | 2017-03-29 | 2019-01-08 | 中国海洋石油总公司 | 基于计算机图形学的喷发相火山岩地质体三维建模方法 |
CN108333627B (zh) * | 2018-01-16 | 2019-01-08 | 成都理工大学 | 火成岩地区断裂真伪的识别方法及装置 |
CN109343114A (zh) * | 2018-11-10 | 2019-02-15 | 北京科胜伟达石油科技股份有限公司 | 复杂岩相火山岩的层拉平的断裂解释方法 |
CN109740212A (zh) * | 2018-12-21 | 2019-05-10 | 中科恒道(北京)科技股份有限公司 | 基于层剥离技术的时间域速度建模方法 |
CN109709605B (zh) * | 2019-02-28 | 2020-09-01 | 中国石油化工股份有限公司 | 火成岩复杂区碳酸盐岩断溶体地震数据成像方法及设备 |
CN112147686B (zh) * | 2019-06-28 | 2022-06-21 | 中国石油化工股份有限公司 | 多期发育火成岩叠前深度偏移成像速度建模方法及系统 |
CN112859167B (zh) * | 2019-11-28 | 2023-08-22 | 中国石油天然气股份有限公司 | 畸变地质体速度场的校正方法及装置 |
CN113589368B (zh) * | 2020-04-30 | 2023-12-01 | 中国石油化工股份有限公司 | 火成岩影响地震反射特征变化及规律的定量分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2093591A1 (en) * | 2008-02-22 | 2009-08-26 | PGS Geophysical AS | Method for Three Dimensional Seismic Travel Time Tomography in Transversely Isotropic Media |
CN101937100A (zh) * | 2010-08-17 | 2011-01-05 | 中国科学院地质与地球物理研究所 | 一种叠前深度偏移方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6864890B2 (en) * | 2002-08-27 | 2005-03-08 | Comoco Phillips Company | Method of building and updating an anisotropic velocity model for depth imaging of seismic data |
-
2013
- 2013-08-15 CN CN201310356114.5A patent/CN104375172B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2093591A1 (en) * | 2008-02-22 | 2009-08-26 | PGS Geophysical AS | Method for Three Dimensional Seismic Travel Time Tomography in Transversely Isotropic Media |
CN101937100A (zh) * | 2010-08-17 | 2011-01-05 | 中国科学院地质与地球物理研究所 | 一种叠前深度偏移方法 |
Non-Patent Citations (5)
Title |
---|
ADDRESSING ANISOTROPY IN PRESTACK DEPTH MIGRATION:A SOUTHERN NORTH SEA CASE STUDY;K. HAWKINS,et al.;《EAGE 63rd Conference & Technical Exhibition》;20010615;第1-4页 * |
VTI各向异性叠前深度偏移技术应用研究——以JZ工区三维地震资料处理为例;赵玲芝等;《石油物探》;20110331;第50卷(第2期);摘要及第202-203页 * |
兴城火山岩地区三维地震资料叠前深度偏移;周振兴等;《勘探地球物理进展》;20061031;第29卷(第5期);第367页 * |
叠前深度偏移技术在深层火山岩成像中的应用;任朝发等;《内蒙古石油化工》;20091231(第14期);摘要及第90-92页 * |
各向异性叠前偏移技术及应用;田彦灿等;《新疆石油地质》;20111031;第32卷(第5期);第533-536页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104375172A (zh) | 2015-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104375172B (zh) | 火山岩下伏地层的构造形态正确成像的方法 | |
CN106597533B (zh) | 一种用于山前带地震资料处理的深度域速度建模方法 | |
CN105317431B (zh) | 一种用于解释和评价水平井测井参数的方法 | |
CN104847264B (zh) | 一种实现地质导向的方法和装置 | |
CN104459782B (zh) | 使用薄层的沿层层析反演地震速度建模方法以及建模单元 | |
CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
CN103077558A (zh) | 碳酸盐岩缝洞型油藏大型溶洞储集体分布模型的建模方法 | |
CN104849762B (zh) | 利用水平井入层点信息校正油藏顶面微构造的方法 | |
CN107132578B (zh) | 一种微地震地面监测速度模型校正算法 | |
CN104793247A (zh) | 近地表速度建模方法及建模装置 | |
CN105093307B (zh) | 下古生界斜地层真厚度求取方法 | |
CN108121009B (zh) | 开发后期复杂断块变速构造成图方法 | |
CN103149588B (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN105093320A (zh) | 针对高速结晶盐壳覆盖区层析静校正初至拾取方法 | |
CN105700013B (zh) | 构建页岩气甜点产能模型的方法 | |
CN109164487A (zh) | 一种基于模型建立平均速度场的方法及精细成图的方法 | |
CN104597494A (zh) | 地震地层体分析方法及装置 | |
CN102943669B (zh) | 适用于水平井轨迹设计的高精度深度域随机模拟反演方法 | |
CN105114063A (zh) | 一种斜井双侧向测井曲线校正方法 | |
CN105549084A (zh) | 一种三维高精度速度建模方法和系统 | |
Han et al. | Characteristics and genesis of the polygonal fault system in southern slope of the Qiongdongnan Basin, South China Sea | |
CN105911592A (zh) | 一种实钻约束的三维地震资料长波长静校正方法 | |
CN105116447A (zh) | 一种基于曲率异常条带的地质河道方向判别方法 | |
US20210396897A1 (en) | Computer implemented method for correcting a reservoir model of a reservoir geological formation based on seismic images | |
Eisner et al. | Effective VTI anisotropy for consistent monitoring of microseismic events |
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 |