CN101363916B - 一种高精度折射静校正数据反演方法 - Google Patents
一种高精度折射静校正数据反演方法 Download PDFInfo
- Publication number
- CN101363916B CN101363916B CN2007101201328A CN200710120132A CN101363916B CN 101363916 B CN101363916 B CN 101363916B CN 2007101201328 A CN2007101201328 A CN 2007101201328A CN 200710120132 A CN200710120132 A CN 200710120132A CN 101363916 B CN101363916 B CN 101363916B
- Authority
- CN
- China
- Prior art keywords
- curve
- refraction
- depth
- data
- static correction
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及地球物理勘探方法,是一种高精度折射静校正数据反演方法,实现步骤是:地震采集工区的生产炮初至波资料;用以下步骤计算折射波延迟时;微测井控制点资料与折射延迟时信息取得工区的低降速层“时-深曲线”。微测井控制点提供目标层深度(HD),折射成果提供延迟时tD,综合tD(HD)曲线,转换t(HD)曲线;得到可靠的表层结构模型及高精度低降速带静校正成果剖面图。本发明用精确的低降速层时深曲线对延迟时间约束反演,消除了由于速度不准造成的深度误差,应用条件宽,适应性强,最终提高成果精度。
Description
所属领域
本发明涉及地球物理勘探方法,是用于地震资料处理中的一种高精度折射静校正数据反演方法。
背景技术
地震勘探是用于寻找油气田的主要物探方法,它用炸药等方法在地面产生人工地震,由专用仪器在一系列观测点上采集地震信号,经过室内资料处理、解释,提取有效信息,调查与油气有关的地质结构情况。
为了得到高信噪比的地震成果资料,目前地震技术普遍采用“多次覆盖”方法,在采集中对测线上每个测点的地下情况,进行共中心点多次观测,每次观测炮点—检波点距离不同。资料处理时通过抽道集——静校正——速度分析——动校正——水平叠加、偏移迭加等流程,得到高质量成果剖面。
上述方法从几何地震学原理出发,假定地面水平,近地表介质均匀。但实际情况常与理论不一致,如地形起伏,低降速带厚度、速度变化大。这种“不一致”,会严重影响地震资料处理效果。静校正即对以上近地表因素校正,改善处理效果。在地表条件复杂的低信噪比地区,静校正工作是决定地震成果质量的首要环节,直接影响速度分析和动校正的精度,最终影响成果资料的信噪比、分辨率,地质构造的精度。
基准面校正是普遍采用的静校正方法。由表层调查取得近地表参数,建立表层模型;通过低降速层时间剥离、高速层时间充填,将炮点、检波点位置校正到地表附近的水平基准面上。建立精确的近地表模型,是提高静校正精度和效果的关键。
“大折射—时深曲线法”,是近年来沙漠区用于表层结构调查及高精度静校正的有效方法。该方法在沙漠区的成果精度与微测井资料相当。采集用直测线、大仪器、大道距、点激发、点接收,炮点与检波点重合的严格互换观测系统;延迟时计算用互换法公式,如:
tD(G)—G点延迟时
h—法线深度
V0—盖层平均速度
VC—折射层界面速度
目的层折射互换时。
由微测井取得“时—深曲线”,用于时—深转换,最终得到可靠的表层模型和高精度静校正反演。
但是大折射—时深曲线法没有利用生产炮的折射初至波,需要单独沿测线采集大折射资料,采集工作量大,费用较高。
发明内容
本发明目的是提供一种省去“大折射-时深曲线法”野外单独采集浅层折射资料,节约了费用的高精度折射静校正数据反演方法。
为实现本发明目的,提供如下技术方案:
1)地震采集工区的生产炮初至波资料;
所述的初至波资料是炸药震源或可控震源资料。
所述的初至波资料二维或三维资料。
所述的初至波资料包括炮记录的初至波时间、炮点、接收点坐标、高程及井深、组合参数。
2)用以下步骤计算折射波延迟时(图2)。
(1)对各地震道进行时间校正;
采用以下A-D公式,逐炮将炮点、检波点不在一条直线分布的二维生产炮资料,校正为直测线、点激发、点接收,炮点和接收点重合的资料;
A.归位校正值Δts
式中:ΔS—归位前后的炮检距离差,Vc—高速层速度,
B.偏移炮点的高程校正Δth
式中:Δh—归位前后炮点高程差,V1—降速层速度
C、组合长度校正ΔtL
式中:ds—炮点组合长度,dr—检波点组合长度,
D、井深校正Δtj
式中:Δhj—井深,V0—低速层速度
(2)根据时差为常数,初至波的高速层速度和波形特征稳定,无明显的干涉确定每炮目标折射层的可靠追踪段:
(3)延长每炮目的层的追踪段长度,由相邻炮目的层初至时距曲线平行移动得到综合时距曲线;
所述的移动的时间值由两炮的追逐时差曲线取得。
(4)按“大折射—时深曲线法”直测线的互换法计算延迟时;
(5)对于三维资料,抽取每条接收测线两侧的近炮,按前述(1)—(4)步骤计算延迟时;
3)微测井控制点资料与折射延迟时信息取得工区的低降速层“时-深曲线”。微测井控制点提供目标层深度(HD),折射成果提供延迟时tD,综合tD(HD)曲线,由tD(HD)及以下公式转换t(HD)曲线;
tD、HD——tD(HD)曲线上任意一点的延迟时间、深度值,t—对应的垂直时间;
所述的控制点资料为:微测井点位置——线号、点号、坐标,垂直时距曲线(t-H),低降速层底界位置(t,HD)。
所述的时-深曲线包括延迟时—深度曲线(tD-HD)和垂直时—深度曲线(t-HD)。
所述的t(HD)曲线也可由多口微测井资料综合得出,再转换tD(HD)曲线。
4)得到可靠的表层结构模型及高精度低降速带静校正成果剖面图。
本发明与“大折射—时深曲线方法”的原始数据不同,使用生产炮的折射初至波,不需要单独沿测线采集折射资料,节约了费用。
本发明首先对生产炮目的层折射初至波作相应的各种校正,消除炮点、检波点位置不规则、偏移炮点高程与检波点不一致及地形、生产因素的影响,使其变为与“大折射-时深曲线法”相当的理想条件,提高了静校正成果的精度。
本发明将单炮目的层追踪段延长,使方法的应用条件放宽,即只要求在归位后的直线上相邻炮记录目的层追踪段有重迭。方法的适应性强,扩大了应用范围。
与扩展广义互换法(EGRM)比较,本法明拟互换(SRM)的延迟时间计算中,增加了偏移炮点的地形校正和施工参数校正,减少随机和系统误差,提高了计算精度。
本发明用精确的低降速层时深曲线对延迟时间约束反演,消除了由于速度不准造成的深度误差,最终提高成果精度。
附图说明
图1EGRM算法示意图;
图2本法明算法示意图。
具体实施方式
本发明实施例1。表层调查采用生产炮抽线计算13条,线距1.5km;布设大折射线8条、线距3km,其中6条线与三维生产炮抽线重合、用于成果对比,验证用于生产炮的计算精度。区内布设7口微测井,由微测井提供表层速度及控制点深度资料。本法明6条重合线的大折射、生产炮初至成果对比,二者的延迟时和低降速层底界面基本重合。延迟时和低降速层底界面高程误差统计:延迟时误差的算术平均值—0.32ms,误差绝对平均值1.77ms;底界面高程误差的算术平均值0.45m,误差绝对平均值1.09m。三维静校正最终成果—低降速层底界面高程及炮、检静校正量数据用于资料处理,单炮及剖面信噪比明显改善。
本发明实施例2。表层调查对二维生产炮初至计算折射成果剖面23条。布设微测井73口。折射成果剖面低降速层底界面高程与微测井资料对比,沙漠区高程误差<2m,边缘区高程误差2~5m。静校正最终成果数据用于资料处理,单炮及剖面信噪比明显改善。
广义互换法(EGRM)与拟互换法(SRM),都是互换法RM的变种。EGRM与SRM法,可用于弯曲测线、炮点与检波点不重合的排列形式,见图1、图2。
图1中▲:归位前炮点位置;△:归位后炮点位置,延迟时综合表达式为:
公式(2-1),相当于将互换排列的炮点(O1、O2)校正到靠近检波点(A、C)—互换点,排列上的解释点(B)校正到互换点直线AC上,在校正后的直线排列上用互换法公式计算延迟时。
归位校正后的直线排列资料与“大折射”有以下区别:由于未作偏移炮点的高程校正,可能互换点时间误差较大;未作组合长度等施工参数校正,造成折射初至时间误差;综合效果,折射延迟时成果有系统和随机误差。EGRM方法的应用条件:测线条带范围内,目的层界面水平(三种方法相同);单炮目的层折射必须有较长的可靠追踪段,保证互换点连接及互换排列目的层有重叠段(这点限制了它的应用范围)。
图2中●:归位前的检波点位置;○:归位后的检波点位置。检波点归位后,其高程不变;炮点归位后,其高程与检波点高程相同。
图2归位直测线的互换排列A1B1C1上,解释点B1的延迟时tD(B1)计算公式为(3-1)式。
包括校正、计算过程的综合表达式为式(3-2)
转换后的一般表达式为(3-3)
公式(3-2)、(3-3)式中右端第一项为互换项,第二项是归位剩余项,第三项是归位炮点高程校正,第四项是施工参数校正,包括组合、井深、仪器延迟时等校正。
具体实现步骤:
本发明对炮点和检波点位置不在一条直线上的二维资料,通过校正将其变为直测线,近似于‘互换法’理想条件中的炮点和接收点位置重合,点激发和点接收;在归位直测线上,确定目的层的追踪段,并延长单炮目的层追踪距离,按互换法公式计算折射成果;拟互换法延迟时计算的综合表达式为:
用微测井资料取得的低降速层时深曲线对折射延迟时间约束反演,得到可靠的表层模型及高精度静校正数据和剖面。
本发明具体技术方案如下:
1)地震采集工区的生产炮初至波资料;
所述的初至波资料是炸药震源或可控震源资料。
所述的初至波资料二维或三维资料。
所述的初至波资料包括炮记录的初至波时间、炮点、接收点坐标、高程及井深、组合参数。
2)用以下步骤计算折射波延迟时,
(1)对各地震道进行时间校正;
采用以下A-D公式,逐炮将炮点、检波点不在一条直线分布的二维生产炮资料,校正为直测线、点激发、点接收,炮点和接收点重合的资料;
A.归位校正值Δts
式中:ΔS一归位前后的炮检距离差,Vc—高速层速度,
B.偏移炮点的高程校正Δth
式中:Δh—归位前后炮点高程差,V1—降速层速度
C、组合长度校正ΔtL
式中:ds—炮点组合长度,dr—检波点组合长度,
D、井深校正Δtj
式中:Δhj—井深,V0—低速层速度
(2)根据时差为常数,初至波的高速层速度和波形特征稳定,无明显的干涉确定每炮目标折射层的可靠追踪段:
(3)延长每炮目的层的追踪段长度,由相邻炮目的层初至时距曲线平行移动得到综合时距曲线;
所述的移动的时间值由两炮的追逐时差曲线取得。
(4)按“大折射—时深曲线法”直测线的互换法计算延迟时;
(5)对于三维资料,抽取每条接收测线两侧的近炮,按前述(1)—(4)步骤计算延迟时;
3)微测井控制点资料与折射延迟时信息取得工区的低降速层“时-深曲线”。微测井控制点提供目标层深度(HD),折射成果提供延迟时tD,综合tD(HD)曲线,由tD(HD)及以下公式转换t(HD)曲线;
tD、HD——tD(HD)曲线上任意一点的延迟时间、深度值,t—对应的垂直时间;
所述的控制点资料为:微测井点位置——线号、点号、坐标,垂直时距曲线(t-H),低降速层底界位置(t,HD)。
所述的时-深曲线包括延迟时—深度曲线(tD-HD)和垂直时—深度曲线(t-HD)。
所述的t(HD)曲线也可由多口微测井资料综合得出,再转换tD(HD)曲线。
4)得到可靠的表层结构模型及高精度低降速带静校正成果剖面图。
Claims (8)
1.一种高精度折射静校正数据反演方法,其特征在于采用以下步骤:
1)地震采集工区的生产炮初至波资料;
2)用以下步骤计算折射波延迟时;
(1)对各地震道进行时间校正;
采用以下A-D公式,逐炮将炮点、检波点不在一条直线分布的二维生产炮资料,校正为直测线、点激发、点接收,炮点和接收点重合的资料;
A.归位校正值Δts,
式中:ΔS-归位前后的炮检距离差,Vc-高速层速度,
B.偏移炮点的高程校正Δth,
式中:Δh-归位前后炮点高程差,V1-降速层速度,
C、组合长度校正ΔtL,
式中:ds-炮点组合长度,dr-检波点组合长度,
D、井深校正Δtj,
式中:Δhj-井深,V0-低速层速度;
(2)根据时差为常数,初至波的高速层速度和波形特征稳定,无明显的干涉确定每炮目的层折射的可靠追踪段:
(3)延长每炮目的层的追踪段长度,由相邻炮目的层初至时距曲线平行移动得到综合时距曲线;
(4)按“大折射-时深曲线法”直测线的互换法计算延迟时;
(5)对于三维资料,抽取每条接收测线两侧的近炮,按前述(1)-(4)步骤计算延迟时;
3)微测井控制点资料与折射延迟时信息取得工区的低降速层“时-深曲线”;微测井控制点提供目的层深度HD,折射成果提供延迟时tD,综合延迟时-深度曲线tD(HD),由延迟时-深度曲线tD(HD)及以下公式转换垂直时-深度曲线t(HD);
tD、HD——延迟时-深度曲线tD(HD)曲线上任意一点的延迟时间、深度值,t-对应的垂直时间;
4)得到可靠的表层结构模型及高精度低降速带静校正成果剖面图。
2.根据权利要求1所述的一种高精度折射静校正数据反演方法,其特征在于,步骤1)所述的初至波资料是炸药震源或可控震源资料。
3.根据权利要求1所述的一种高精度折射静校正数据反演方法,其特征在于,步骤1)所述的初至波资料是二维或三维资料。
4.根据权利要求1所述的一种高精度折射静校正数据反演方法,其特征在于,步骤1)所述的初至波资料包括炮记录的初至波时间、炮点、接收点坐标、高程及井深、组合参数。
5.根据权利要求1所述的一种高精度折射静校正数据反演方法,其特征在于,步骤2)所述移动的时间值由两炮的追逐时差曲线取得。
6.根据权利要求1所述的一种高精度折射静校正数据反演方法,其特征在于,步骤3)所述的控制点资料为:微测井点位置——线号、点号、坐标,垂直时距曲线(t-H),低降速层底界位置(t,HD)。
7.根据权利要求1所述的一种高精度折射静校正数据反演方法,其特征在于,步骤3)所述的时-深曲线包括延迟时-深度曲线tD(HD)和垂直时-深度曲线t(HD)。
8.根据权利要求1所述的一种高精度折射静校正数据反演方法,其特征在于,步骤3)所述的垂直时-深度曲线t(HD)也可由多口微测井资料综合得出,再转换延迟时-深度曲线tD(HD)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101201328A CN101363916B (zh) | 2007-08-10 | 2007-08-10 | 一种高精度折射静校正数据反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101201328A CN101363916B (zh) | 2007-08-10 | 2007-08-10 | 一种高精度折射静校正数据反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101363916A CN101363916A (zh) | 2009-02-11 |
CN101363916B true CN101363916B (zh) | 2012-05-30 |
Family
ID=40390388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101201328A Expired - Fee Related CN101363916B (zh) | 2007-08-10 | 2007-08-10 | 一种高精度折射静校正数据反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101363916B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101533103B (zh) * | 2009-04-13 | 2012-03-21 | 中国石油天然气集团公司 | 地震波速度场构建中井约束横向可变h-v曲线构建方法 |
CN101980054B (zh) * | 2010-10-22 | 2012-07-18 | 中国石油化工股份有限公司 | 一种在高密度地震静校正处理中建立近地表速度模型的方法 |
CN102890288B (zh) * | 2011-07-20 | 2015-05-13 | 中国石油天然气集团公司 | 一种地震波层速度反演方法 |
CN102890289B (zh) * | 2011-07-20 | 2015-05-27 | 中国石油天然气集团公司 | 一种高精度地震波层速度反演方法 |
CN102998699B (zh) * | 2011-09-08 | 2015-09-23 | 中国石油天然气集团公司 | 一种含有套管谐波的垂直地震剖面资料层速度反演方法 |
CN104375176B (zh) * | 2013-08-15 | 2017-05-10 | 中国石油天然气集团公司 | 一种利用初至时间误差确定同相叠加最高截至频率的方法 |
CN103616720B (zh) * | 2013-10-28 | 2016-04-06 | 中国石油天然气集团公司 | 一种折射速度获取方法及装置 |
CN103605163B (zh) * | 2013-10-28 | 2016-05-11 | 中国石油天然气集团公司 | 一种静校正量定量评价方法及装置 |
CN104570091B (zh) * | 2014-12-15 | 2018-02-02 | 中国石油天然气集团公司 | 一种获得初至波射线的方法和装置 |
CN104656134A (zh) * | 2015-03-11 | 2015-05-27 | 中国电建集团北京勘测设计研究院有限公司 | 一种地震波初至读时与时距曲线动态关联的方法和系统 |
CN107656314B (zh) * | 2017-10-20 | 2020-03-10 | 中国石油天然气集团公司 | 道集记录处理方法、装置及计算机存储介质 |
CN109507730B (zh) * | 2018-12-29 | 2020-05-01 | 辽宁工程技术大学 | 一种改进的确定折射面法线深度的方法 |
CN113376697B (zh) * | 2020-03-10 | 2024-04-30 | 中国石油天然气股份有限公司 | 基于三维近地表模型的静校正方法及装置 |
CN111505711A (zh) * | 2020-06-05 | 2020-08-07 | 中国科学技术大学 | 3d折射波波场偏移成像方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6501703B2 (en) * | 2001-01-21 | 2002-12-31 | Xixiang Zhou | Static correction method for exploration seismic data using first arrivals of seismic waves |
CN1138158C (zh) * | 2001-01-19 | 2004-02-11 | 中国石油天然气股份有限公司 | 一种地震多域迭代静校正方法 |
CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
-
2007
- 2007-08-10 CN CN2007101201328A patent/CN101363916B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1138158C (zh) * | 2001-01-19 | 2004-02-11 | 中国石油天然气股份有限公司 | 一种地震多域迭代静校正方法 |
US6501703B2 (en) * | 2001-01-21 | 2002-12-31 | Xixiang Zhou | Static correction method for exploration seismic data using first arrivals of seismic waves |
CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
Non-Patent Citations (3)
Title |
---|
杨海申等.复杂区三维折射静校正技术与应用效果.《石油地球物理勘探》.2005,第40卷(第2期),219-225页. * |
闫杰等.用于折射波静校正计算的拟互换法.《中国地球物理第二十一届年会论文集》.2005,第69页. * |
闫杰等.表层复杂区一种高精度折射静校正方法.《中国地球物理.2003——中国地球物理学会第十九届年会论文集》.2003,第36页. * |
Also Published As
Publication number | Publication date |
---|---|
CN101363916A (zh) | 2009-02-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101363916B (zh) | 一种高精度折射静校正数据反演方法 | |
CN101980054B (zh) | 一种在高密度地震静校正处理中建立近地表速度模型的方法 | |
CN100487488C (zh) | 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法 | |
CN105022031B (zh) | 一种区域岩体微震震源的分层速度定位方法 | |
CN102213769A (zh) | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 | |
CN100549731C (zh) | 检波点二次定位方法 | |
KR20200014387A (ko) | 지하 구조물의 검출 | |
CN106597533A (zh) | 一种用于山前带地震资料处理的深度域速度建模方法 | |
CN104570102B (zh) | 一种近地表速度模型与中深层速度模型的融合方法 | |
US4695984A (en) | Method for establishing a surface consistent correction for the effects of the low velocity layer in seismic data processing | |
CN105510880A (zh) | 一种基于双差法的微地震震源定位方法 | |
CN104656141A (zh) | 共姿态道集地表一致性振幅补偿方法及装置 | |
CN105093320A (zh) | 针对高速结晶盐壳覆盖区层析静校正初至拾取方法 | |
CN101630016A (zh) | 一种提高垂直地震剖面成像质量的方法 | |
CN101986172B (zh) | 一种根据实际钻井轨迹校正vsp下行波的处理方法 | |
CN102053263A (zh) | 调查表层结构的方法 | |
CN101349763B (zh) | 用于起伏地表地震勘探数据处理的Kirchhoff叠前时间偏移方法 | |
Juhlin et al. | Reflection seismic studies over the end-glacial Burträsk fault, Skellefteå, Sweden | |
Shynar et al. | Methodology of assessment and prediction of critical condition of natural-technical systems | |
CN103149588A (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN105093318A (zh) | 一种自适应波动方程波场延拓静校正方法 | |
Xu et al. | Optimal design of microseismic monitoring networking and error analysis of seismic source location for rock slope | |
CN105866839A (zh) | 基于共姿态道集的静校正方法及装置 | |
Juhlin et al. | Imaging of groundwater resources in glacial deposits using high-resolution reflection seismics, Sweden | |
CN102087368B (zh) | 一种近于零震-检距的反射法测取波速的方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120530 Termination date: 20210810 |
|
CF01 | Termination of patent right due to non-payment of annual fee |