CN101393270B - 一种连续介质地区的表层高精度静校正方法 - Google Patents
一种连续介质地区的表层高精度静校正方法 Download PDFInfo
- Publication number
- CN101393270B CN101393270B CN2007101220808A CN200710122080A CN101393270B CN 101393270 B CN101393270 B CN 101393270B CN 2007101220808 A CN2007101220808 A CN 2007101220808A CN 200710122080 A CN200710122080 A CN 200710122080A CN 101393270 B CN101393270 B CN 101393270B
- Authority
- CN
- China
- Prior art keywords
- static correction
- surface layer
- high precision
- depth
- formula
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及地球物理勘探技术,是一种地震勘探数据的表层连续介质地区的表层高精度静校正方法。步骤是利用微测井方法采集得到高速顶界面高程和低速带厚度;对区域每一微测井的时深关系进行二次曲线拟合得到每一口微测井点的近地表初始速度;按平面位置的变化进行插值,得到整个区域空间上连续变化的时深关系,得到时深关系的三维空间数据体;计算低速带静校正量以及其他静校正量;用通常的方法将静校正量进行地震资料处理得到高分辨反射波成像。本发明可以避开难以分层近地表的人为分层误差,与沙丘曲线相比具有更广的适用性和通用性,精度更高。
Description
技术领域
本发明涉及地球物理勘探方法,是一种地震勘探数据的表层连续介质地区的表层高精度静校正方法。
背景技术
在地震勘探中,静校正主要是解决地震波在近地表地层旅行时间变化影响反射波叠加成像的问题。由于低速带的存在地震反射波出射地表时接近垂直地表方向。近地表简单的地区近地表结构呈层状,静校正的一般算法是按地震波在层状介质的垂直旅行时间计算的。在近地表呈连续变化的地区表层结构分层困难,常规静校正算法不适用,一般做法是利用统计规律统计出某一地区的近地表的时深关系,利用时深关系得到低速带校正量,如适用沙漠区的沙丘曲线静校正方法。
该方法是解决大沙漠区中、长波长静校正问题的较好方法。沙漠腹地物性成因简单,应用一条沙丘曲线就可以得到较好的效果;而沙漠边缘因物性成因比较复杂,近地表在垂向上虽然呈连续变化,但在横向上变化很大,不再具有统计规律,单一时深曲线不在适用。原理就是利用微测井或折射原理作沙丘曲线调查在一定的误差范围内统计生产沙漠区的时深关系曲线,然后根据表层调查得到近地表高速顶即表层厚度,利用沙丘曲线求取表层静校正量。该方法在表层比较稳定的沙漠区或相对比较小的地区适用,不管是分段拟合还是分块都没有突破单一时深曲线的思路,不适宜具有相似表层特点的巨厚砾石区、风成沙漠区。巨厚砾石区、风成沙漠区的表层在单点或较小的区域是符合上述规律的,但直接应用以上方法建立时深曲线(沙漠区叫沙丘曲线、砾石区可称为砾石曲线)不再适用,否则将产生较大的静校正误差。
发明的内容
本发明目的是在不稳定的表层连续介质地区地震勘探时,通过高精度的静校正提高地震资料的反射波成像效果。
本发明通过以下具体题步骤实现:
1)利用微测井方法采集得到高速顶界面高程(Hg),低速带厚度(h);
步骤1)所述的低速带厚度(h)为地表高程减高速顶界面高程。
2)利用最小二乘原理对区域每一口微测井的时深关系进行二次曲线拟合;
3)按照公式 (1)与拟合的二次曲线得到每一口微测井点的近地表初始速度V0,速度随深度的变化率β,其中h为微测井实测测量点的深度(m),T为实测测量点到井口的单程旅行时间(s);
4)将V0和β按按平面位置的变化进行插值,得到整个区域空间上连续变化的时深关系,根据公式(1)得到时深关系的三维空间数据体;
步骤4)所述的插值采用有理Lagrange插值曲面方法。
5)由h,V0和β,用以下公式得到低速带静校正量tw以及其他静校正量;
式中:td:静校正量(s)
tw:低速带静校正量
Hd:静校正计算的基准面(m)
Hg:高速顶界面高程(m)
Vd:基准面替换速度(m/s)
6)将静校正结果与地震资料处理结合,在同一基准面采集、接收,通过处理叠加,得到高分辨反射波成像。
本发明比较常规静校正方法相比可以避开难以分层近地表的人为分层误差,与沙丘曲线相比具有更广的适用性和通用性,精度更高。
本发明图2为发明应用目标区所有微测井时深散点叠合图。从图中可以看出砾石区和沙漠区时深关系是有明显差别的;同一类型地表不同位置的时深关系也存在一定的差别。利用每一口微测井拟合一条时深曲线,然后空间内插到每一个物理点。空变的时深曲线比单一时深曲线更符合该地区的表层情况。地震资料剖面效果如图3,为目标区LD三维工区现场处理剖面对比,常规时深曲线法大部分成像效果好,大沙漠区成像精度略差;沙漠戈壁交界部位构造形态异常,可能存在长波长静校正问题,空变时深曲线静校正方法的应用,沙漠戈壁构造形态自然,成像效果相当,大沙漠区成像局部变好。
附图说明
图1是微测井时深关系二次曲线拟合图;
图2是工区微测井时深散点叠合图;
图3是三维常规时深曲线法(左)与本发明方法(右)剖面对比。
具体实施方式
针对不稳定表层连续介质地区静校正存在的问题,本发明依据近地表横向连续变化,将时深曲线进行按空间不同位置而变化而形成一种不稳定的表层连续介质地区的表层高精度静校正方法。
本发明实现过程是表层连续变化,横向连续变化,对横向求极限,在某一点上有纵向连续变化,该点可以得到时深关系曲线,这个点的数据得来依靠微测井(单一沙丘曲线可以来源于微测井,也可以利用折射原理作沙丘曲线调查在一定的误差范围内统计),当有数量足够多的微测井就可以得到分布在不同位置的时深曲线,然后对这些离散的曲线形成平面过渡。最后形成空间的立体时深曲线关系,利用空变的时深关系曲线求取低速带校正量。
本发明通过以下具体题步骤实现:
1)利用通常的微测井方法采集得到高速项界面高程(Hg),低速带厚度(h,地表高程减高速顶界面高程);
2)对每一口微测井的时深关系进行二次曲线拟合(图1);
由于表层呈连续变化,微测井时深关系连续变化时深关系为二次曲线关系,利用最小二乘方法进行二次曲线拟合;
在地区微测井中,单程旅行时间与激发点深度存在二次变化关系,符合连续介质定义,通用连续介质定义 当n=2,有 介质速度随深度非线性变化。
建立直角坐标,纵坐标为微测井的激发点深度,横坐标为该点激发到地表接收的垂直旅行时间,将微测井实测数据展布在该坐标系中,利用最小二乘原理,拟合二次曲线,深度是时间的二次函数,且深度为零时时间也为零。
3)按照公式
对比2)中拟合的二次曲线提取每一口微测井点的近地表初始速度V0和速度随深度的变化率β,
公式中h:为微测井实测测量点的深度(m),
T:为实测测量点到井口的单程旅行时间(s);
4)将V0和β按“有理Lagrange插值曲面”原理按平面位置的变化进行插值,得到整个区域空间上连续变化的时深关系,根据公式(1)可知即得到时深关系的三维空间数据体;
5)由1)得到h,4)得到V0和β,取得低速带静校正量tw;
6)计算静校正量
公式中td:静校正量(s)
tw:低速带静校正量
Hd:静校正计算的基准面(m)
Hg:高速顶界面高程(m)
Vd:基准面替换速度(m/s)
7)通过1)-4)解决不稳定的表层连续介质地区的表层问题,由5)-6)得到高精度的静校正;将静校正应用到地震资料处理中,达到在同一基准面采集、接收,通过处理叠加,反射波成像效果得到提高。
Claims (3)
1.一种连续介质地区的表层高精度静校正方法,其特征在于:具体通过以下步骤实现:
1)利用微测井方法采集得到高速顶界面高程Hg,低速带厚度;
2)利用最小二乘原理对区域每一口微测井的时深关系进行二次曲线拟合;
每一口微测井点的近地表初始速度V0,速度随深度的变化率β,其中h为微测井实测测量点的深度(m),T为实测测量点到井口的单程旅行时间(s);
4)将V0和β按平面位置的变化进行插值,得到整个区域空间上连续变化的时深关系,根据公式(1)得到时深关系的三维空间数据体;
5)由h,V0和β,用以下公式得到低速带静校正量tw以及其他静校正量;
式中:td:静校正量(s)
tw:低速带静校正量
Hd:静校正计算的基准面(m)
Hg:高速顶界面高程(m)
Vd:基准面替换速度(m/s)
6)用通常的方法将静校正量进行地震资料处理,在同一基准面对地震资料采集、接收,通过对地震资料处理叠加,得到高分辨反射波成像。
2.根据权利要求1所述的连续介质地区的表层高精度静校正方法,其特征在于步骤1)所述的低速带厚度为地表高程减高速顶界面高程。
3.根据权利要求1所述的连续介质地区的表层高精度静校正方法,其特征在于步骤4)所述的插值采用有理Lagrange插值曲面方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101220808A CN101393270B (zh) | 2007-09-21 | 2007-09-21 | 一种连续介质地区的表层高精度静校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101220808A CN101393270B (zh) | 2007-09-21 | 2007-09-21 | 一种连续介质地区的表层高精度静校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101393270A CN101393270A (zh) | 2009-03-25 |
CN101393270B true CN101393270B (zh) | 2011-04-20 |
Family
ID=40493654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101220808A Active CN101393270B (zh) | 2007-09-21 | 2007-09-21 | 一种连续介质地区的表层高精度静校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101393270B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101819279B (zh) * | 2010-03-26 | 2012-05-23 | 中国石油集团西北地质研究所 | 基于变差函数拟合重构的高精度三维静校正方法 |
CN102236104B (zh) * | 2010-05-06 | 2014-12-24 | 中国石油天然气集团公司 | 三维密度变速地震反射空校方法 |
CN102003172B (zh) * | 2010-10-11 | 2013-04-03 | 青岛杰瑞自动化有限公司 | 测井系统中的测井深度校正计算方法 |
CN101980054B (zh) * | 2010-10-22 | 2012-07-18 | 中国石油化工股份有限公司 | 一种在高密度地震静校正处理中建立近地表速度模型的方法 |
CN102798893B (zh) * | 2011-05-25 | 2015-05-20 | 中国石油天然气集团公司 | 三维叠合层析反演静校正方法 |
CN103728666B (zh) * | 2012-10-16 | 2017-08-18 | 中国石油化工股份有限公司 | 一种确定静校正中间参考面替换速度的方法 |
CN103954997B (zh) * | 2014-04-18 | 2016-09-14 | 中国石油化工股份有限公司 | 一种地震勘探资料近地表静校正处理方法 |
CN104133249B (zh) * | 2014-07-31 | 2016-08-31 | 中国石油天然气集团公司 | 一种微测井资料与声波测井资料联合解释的方法及装置 |
CN104156624B (zh) * | 2014-08-27 | 2017-08-01 | 中国电子科技集团公司第二十二研究所 | 测井设备的时深数据的处理方法和装置 |
CN105572741B (zh) * | 2014-10-08 | 2019-02-01 | 中国石油化工股份有限公司 | 一种计算3d高频静校正量的方法 |
CN104570122B (zh) * | 2015-01-14 | 2017-06-23 | 中国石油集团东方地球物理勘探有限责任公司 | 基于反射波的地震数据静校正方法及装置 |
CN105277978B (zh) * | 2015-10-15 | 2018-01-05 | 中国石油天然气集团公司 | 一种确定近地表速度模型的方法及装置 |
CN109581503B (zh) * | 2017-09-28 | 2020-08-25 | 中国石油化工股份有限公司 | 一种解决沙漠探区层析静校正边界问题的方法及系统 |
CN111323818B (zh) * | 2018-12-13 | 2022-06-03 | 中国石油天然气股份有限公司 | 陆上地震资料的静校正方式的筛选方法及装置 |
CN109884707B (zh) * | 2019-03-20 | 2021-09-07 | 中国石油化工股份有限公司 | 近地表分层时深曲线静校正方法 |
RU2722861C1 (ru) * | 2019-07-09 | 2020-06-04 | Общество с ограниченной ответственностью "СИГМА-ГЕО" | Способ расчета статических поправок |
CN112630838B (zh) * | 2019-09-24 | 2024-06-18 | 中国石油化工股份有限公司 | 一种时间域浮动面转换深度域计算方法及计算系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5663928A (en) * | 1994-10-19 | 1997-09-02 | Elf Aquitaine Production | Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction |
FR2801387A1 (fr) * | 1999-11-23 | 2001-05-25 | Geophysique Cie Gle | Perfectionnements aux procedes de traitement sismique mettant en oeuvre une determination des corrections dynamiques |
CN1308240A (zh) * | 2001-01-21 | 2001-08-15 | 周熙襄 | 用初至折射波对地震勘探资料进行短波长静校正处理的方法 |
CN1365008A (zh) * | 2001-01-19 | 2002-08-21 | 中国石油天然气股份有限公司 | 一种地震多域迭代静校正方法 |
CN1595201A (zh) * | 2004-06-25 | 2005-03-16 | 大庆油田有限责任公司 | 一种提高地震分辨率的方法 |
CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
CN1948999A (zh) * | 2005-10-10 | 2007-04-18 | 中国石油天然气股份有限公司 | 近似层替换静校正方法 |
-
2007
- 2007-09-21 CN CN2007101220808A patent/CN101393270B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5663928A (en) * | 1994-10-19 | 1997-09-02 | Elf Aquitaine Production | Method for analysing and processing seismic reflection data for the determination of a high resolution spatial velocity field for hyperbolicity correction |
FR2801387A1 (fr) * | 1999-11-23 | 2001-05-25 | Geophysique Cie Gle | Perfectionnements aux procedes de traitement sismique mettant en oeuvre une determination des corrections dynamiques |
CN1365008A (zh) * | 2001-01-19 | 2002-08-21 | 中国石油天然气股份有限公司 | 一种地震多域迭代静校正方法 |
CN1308240A (zh) * | 2001-01-21 | 2001-08-15 | 周熙襄 | 用初至折射波对地震勘探资料进行短波长静校正处理的方法 |
CN1595201A (zh) * | 2004-06-25 | 2005-03-16 | 大庆油田有限责任公司 | 一种提高地震分辨率的方法 |
CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
CN1948999A (zh) * | 2005-10-10 | 2007-04-18 | 中国石油天然气股份有限公司 | 近似层替换静校正方法 |
Non-Patent Citations (8)
Title |
---|
吴永富.石油勘探地震资料处理质量控制.《甘肃科技》.2007,(第2期),全文. |
安树杰等.中国西部YKB复杂山地山前带三维地震采集实例.《CPS/SEG2004国际地球物理会议论文集》.2004,全文. |
张玉铭 |
张玉铭;吴永富.石油勘探地震资料处理质量控制.《甘肃科技》.2007,(第2期),全文. * |
陈学强 |
陈耿毅.撒哈拉沙漠腹地静校正方法研究及应用.《中国博士学位论文全文数据库基础科学辑》.2007,(第2期),全文. * |
黄显国 |
黄显国;陈学强;安树杰等.中国西部YKB复杂山地山前带三维地震采集实例.《CPS/SEG2004国际地球物理会议论文集》.2004,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN101393270A (zh) | 2009-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101393270B (zh) | 一种连续介质地区的表层高精度静校正方法 | |
CN104091064B (zh) | 基于优化解空间搜索法的PS‑DInSAR地表形变测量参数估计方法 | |
CN103869368B (zh) | 一种无表层调查资料约束的大炮初至综合建模静校正方法 | |
CN101776768B (zh) | 一种各向异性速度分析和动校正方法 | |
CN101533103B (zh) | 地震波速度场构建中井约束横向可变h-v曲线构建方法 | |
CN101750628B (zh) | 二维叠加速度及均方根速度场闭合差校正方法 | |
CN103292792B (zh) | 一种适用海底探测与假地形处理的实测svp重构方法 | |
CN103399346B (zh) | 一种井震联合初始波阻抗建模方法 | |
CN103472483A (zh) | 基于真地表或浮动基准面的速度建模方法 | |
CN102590862B (zh) | 补偿吸收衰减的叠前时间偏移方法 | |
CN102590860B (zh) | 基于地震波初至信息的反射波建模方法 | |
CN104330827B (zh) | 一种地表模型静校正的处理方法 | |
CN105093307B (zh) | 下古生界斜地层真厚度求取方法 | |
CN106646601A (zh) | 多信息联合约束的浅中深层三维q体建立方法 | |
CN104090301B (zh) | 一种求取三维高频静校正量的方法 | |
CN103630934A (zh) | 一种确定转换波检波点大的横波静校正量的方法 | |
CN104268848A (zh) | 一种海洋内波波速监测的方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN102043167B (zh) | 一种曲面地表地震资料速度分析的方法 | |
CN104570103A (zh) | 一种低信噪比地震资料的井约束速度谱拾取方法 | |
CN101452082A (zh) | 一种分形地震波初至拾取的方法 | |
CN104678439A (zh) | 一种多波连续介质速度分析方法 | |
CN104155690B (zh) | 基于椭球展开的三维地震数据叠加速度求取方法 | |
CN107870374B (zh) | 一种地层剥蚀演化恢复方法 | |
CN114839678A (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 |