CN104375181B - 一种快速生成三维共检波点初至道集的方法 - Google Patents

一种快速生成三维共检波点初至道集的方法 Download PDF

Info

Publication number
CN104375181B
CN104375181B CN201410682365.7A CN201410682365A CN104375181B CN 104375181 B CN104375181 B CN 104375181B CN 201410682365 A CN201410682365 A CN 201410682365A CN 104375181 B CN104375181 B CN 104375181B
Authority
CN
China
Prior art keywords
file
common receiver
certain
dimensional
arrival
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
CN201410682365.7A
Other languages
English (en)
Other versions
CN104375181A (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 National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201410682365.7A priority Critical patent/CN104375181B/zh
Publication of CN104375181A publication Critical patent/CN104375181A/zh
Application granted granted Critical
Publication of CN104375181B publication Critical patent/CN104375181B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提出了一种快速生成三维共检波点初至道集的方法，包括：建立实验地区的SPS文件；获取每一炮点的一定偏移距范围内的初至时间，生成初至数据文件；选取要生成三维共检波点初至道集的检波点作为共检波点；建立三维检波点线的索引指示文件及三维炮点线的索引指示文件；建立每一炮点关于检波点线及检波点的排列关系的索引关系文件；设定线序号搜索范围，在炮点文件中选取并生成某一共检波点对应的炮点线序号搜索范围，选取并生成某一共检波点对应的炮点序号范围；分别判断每一炮点的排列关系中是否包含某一共检波点，如果包含，计算生成某一共检波点的初至数据道集；重复前述两步计算每一共检波点的初至数据道集，生成三维共检波点初至道集。

Description

一种快速生成三维共检波点初至道集的方法
技术领域
本发明涉及地球物理勘探方法，尤指一种快速生成三维共检波点初至道集的方法。
背景技术
共检波点道集是地震资料处理常用的一个数据体，野外数据是共炮点道集，因此共检波点道集需要根据SPS文件(包含炮点文件、检波点文件及关系文件)自主生成。最简单的方法就是从炮点中逐一寻找，提取相关信息，然后生成共检波点道集。当然这样做其计算量也是最大的，尤其是三维，其计算量更是十分惊人。在一般基于初至计算静校正的数据处理中，也会用到共检波点的相关信息。假设某三维，共有检波点M个，炮点N个，每炮有L道接收，则生成1个共检波点道集其要搜索及判断的次数大致为：N×L，则全部共检波点道集搜索及判断的次数为：M×N×L，再加上要提取的每个地震道的相关信息(有些信息需要临时计算，如通过炮点和检波点的坐标信息计算炮检距)，其数据量十分惊人，所生成共检波点道集的时间很难满足实际生产资料的时效要求。
在目前并未有比较快捷的生成共检波点道集的方法，因此依照上述分析，生成三维共检波点道集的计算量十分巨大，在普通微机上所需时间相当长，实际生产中很难满足地震资料处理的时效要求。
发明内容
针对现有技术中生成共检波点道集的方法计算量大，效率低等问题，本发明提出了一种快速生成三维共检波点初至道集的方法，能在普通微机上高速生成共检波点初至道集文件，缩短数据生成时间，从而能满足实际生产的需要。
为达到上述目的，本发明提出了1、一种快速生成三维共检波点初至道集的方法，包括：步骤1，根据野外生产资料，获取并建立实验地区的SPS文件，其中包括：检波点文件、炮点文件及关系文件；步骤2，根据所述SPS文件，获取每一炮点的一定偏移距范围内的初至时间，生成初至数据文件；步骤3，在所述实验地区的检波点中，选取要生成所述三维共检波点初至道集的检波点作为共检波点；步骤4，根据所述检波点文件及炮点文件，分别建立三维检波点线的索引指示文件及三维炮点线的索引指示文件，并根据三维检波点线的索引指示文件，计算生成所述共检波点在检波点文件中的位置数据；步骤5，根据所述关系文件，建立每一炮点关于检波点线及检波点的排列关系的索引关系文件；步骤6，设定线序号搜索范围，根据所述线序号搜索范围在所述炮点文件中选取并生成某一共检波点对应的炮点线序号搜索范围，再根据所述炮点线序号搜索范围及所述三维炮点线的索引指示文件，在所述炮点文件中选取并生成所述某一共检波点对应的炮点序号范围；步骤7，根据所述某一共检波点对应的炮点序号范围以及每一炮点的索引关系文件，分别判断每一炮点的排列关系中是否包含所述某一共检波点，如果包含，根据所述共检波点在检波点文件中的位置数据、索引关系文件及初至数据文件，提取相关数据并计算生成所述某一共检波点的初至数据道集；步骤8，重复上述步骤6、7，计算每一共检波点的初至数据道集，生成三维共检波点初至道集。
本发明提出的快速生成三维共检波点初至道集的方法计算速度快，生成共检波点道集时间短，能提高三维剩余静校正的计算效率，缩短地震资料处理周期，具有广泛的适用性和通用性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1为本发明一实施例的快速生成三维共检波点初至道集的方法流程图。
图2为一具体实施例对某一地区三维项目生成三维共检波点初至道集采用常规方法和本发明方法用时的对比图。
具体实施方式
以下配合图式及本发明的较佳实施例，进一步阐述本发明为达成预定发明目的所采取的技术手段。
图1为本发明一实施例的快速生成三维共检波点初至道集的方法流程图。如图1所示，该方法包括：
步骤1，根据野外生产资料，获取并建立实验地区的SPS文件，其中包括：检波点文件、炮点文件及关系文件；
在步骤1中，检波点文件按检波点桩号升序排序、炮点文件按炮点桩号升序排序。
步骤2，根据SPS文件，获取每一炮点的一定偏移距范围内的初至时间，生成初至数据文件；
步骤3，在实验地区的检波点中，选取要生成三维共检波点初至道集的检波点作为共检波点；
步骤4，根据检波点文件及炮点文件，分别建立三维检波点线的索引指示文件及三维炮点线的索引指示文件，并根据三维检波点线的索引指示文件，计算生成共检波点在检波点文件中的位置数据；
步骤5，根据关系文件，建立每一炮点关于检波点线及检波点的排列关系的索引关系文件；
步骤6，设定线序号搜索范围，根据线序号搜索范围在炮点文件中选取并生成某一共检波点对应的炮点线序号搜索范围，再根据炮点线序号搜索范围及三维炮点线的索引指示文件，在炮点文件中选取并生成某一共检波点对应的炮点序号范围；
步骤7，根据某一共检波点对应的炮点序号范围以及每一炮点的索引关系文件，分别判断每一炮点的排列关系中是否包含某一共检波点，如果包含，根据共检波点在检波点文件中的位置数据、索引关系文件及初至数据文件，提取相关数据并计算生成某一共检波点的初至数据道集；
步骤8，重复上述步骤6、7，计算每一共检波点的初至数据道集，生成三维共检波点初至道集。
在本实施例中的步骤4中，根据三维检波点线的索引指示文件，计算生成共检波点在检波点文件中的位置数据还包括：
根据三维检波点线的索引指示文件，计算某一共检波点在检波点文件中的位置，并提取某一共检波点的坐标，计算某一共检波点在检波点文件中的位置的公式如下：
Rloci,j＝Rloci,1+Rstai,1-Rstai,j；
其中，Rloci,j为某一共检波点在检波点文件中的位置，i为检波点线号，j为检波点桩号；
Rloci,1为检波点线第i线第1个检波点在三维检波点文件中的位置，其值在建立三维检波点线的索引指示文件时获得；
Rstai,1为检波点线第i线第1个检波点桩号；
Rstai,j为某一共检波点桩号。
在本实施例中，步骤7还包括：
如果某一炮点的排列关系中包含某一共检波点，根据炮点文件提取某一炮点的坐标，根据三维检波点线的索引指示文件，计算某一共检波点在检波点文件中的位置，并提取某一共检波点的坐标，计算炮检距；
根据索引关系文件，计算某一共检波点相对某一炮点的排列位置，根据排列位置从初至数据文件中提取某一共检波点及某一炮点对应地震道的初至时间数据；
根据炮检距及初至时间数据，生成某一共检波点的初至数据道集。
其中，步骤7是利用以下公式计算某一共检波点在某一炮点内的排列位置：
SR＝SRnumi,1+SRstai,1-Rstai,j；
其中，SR为某一共检波点相对该炮点的排列位置；
SRnumi,1表示在索引关系文件中，检波点线第i线第1个检波点相对某一炮点的排列位置，其值在建立索引关系文件时获得；
SRstai,1表示在索引关系文件中，检波点线第i线第1个检波点相对某一炮点的检波点桩号，其值在建立索引关系文件时获得；
Rstai,j为某一共检波点桩号。
在本实施例中，步骤8中生成的三维共检波点初至道集可以供计算炮点剩余静校正应用。
本发明的目的是针对在计算三维初至波剩余静校正时，需要事先生成共检波点初至道集数据供计算炮点剩余静校正用。由于三维数据一般数据量较大，按常规方法生成共检波点道集的计算量特别巨大，采用普通计算机计算费时长，很难满足地震资料处理的时效要求；对此，本发明提出了一种快速生成三维共检波点初至道集的方法，能在普通微机上高速生成共检波点初至道集文件，缩短数据生成时间，从而能满足实际生产的需要。
为了对上述快速生成三维共检波点初至道集的方法进行更为清楚的解释，下面结合一个具体的实施例来进行说明，然而值得注意的是该实施例仅是为了更好地说明本发明，并不构成对本发明不当的限定。
在某一地区的三维项目中，利用本发明的方法快速生成三维共检波点初至道集，以供计算炮点剩余静校正应用。
首先，结合步骤1，根据野外生产资料，获取并建立该地区的SPS文件，其中包括：检波点文件、炮点文件及关系文件；检波点文件按检波点桩号升序排序、炮点文件按炮点桩号升序排序。
结合步骤2，根据SPS文件，获取每一炮点的偏移距在100米至2800米范围内的地震道的初至时间，生成初至数据文件；
结合步骤3，以每条检波点线的起始桩号为起点，检波点增量10(也可以是1至20内的任一数)为间隔，选取并确定该三维的共检波点初至道集的共检波点；
结合步骤4，根据检波点文件及炮点文件，分别建立三维检波点线的索引指示文件及三维炮点线的索引指示文件,并根据三维检波点线的索引指示文件，计算生成每一共检波点在检波点文件中的位置数据
结合步骤5，根据关系文件，建立每一炮点关于检波点线及检波点的排列关系的索引关系文件；
结合步骤6，根据实验观测系统的相关参数，设定线序号搜索范围为-10至10(也可以根据实际情况改变，如：-4，8；-6，6等)，根据线序号搜索范围，确定某一个共检波点在炮点文件中要搜索的炮点线号范围在该检波点线号的±10之内，也就可以在炮点文件中确定某一个检波点要对应搜索的炮点范围。
结合步骤7，根据对应的炮点搜索范围以及每一炮点的索引关系文件，分别判断每一炮点的排列关系中是否包含某一共检波点，如果包含，根据共检波点在检波点文件中的位置数据、索引关系文件及初至数据文件，计算生成某一共检波点的初至数据道集；
结合步骤8，重复上述步骤6、7，计算每一共检波点的初至数据道集，生成三维共检波点初至道集；生成的三维共检波点初至道集可供计算炮点剩余静校正应用。
本发明的实现过程主要为：先分别对检波点文件和炮点文件建立检波点线和炮点线的索引文件，对排列关系文件建立索引指示文件，然后对每个炮点根据其排列关系与检波点索引文件，建立该炮点的检波点索引信息文件；最后根据要生成的共检波点逐一抽取相关的炮点信息生成共检波点的初至数据道集，供后续计算炮点剩余静校正应用。
图2为一具体实施例对某一地区三维项目生成三维共检波点初至道集采用常规方法和本发明方法用时的对比图。如图2所示，二者在同一微机上运行。从对比图上可以看出，采用常规方法需要264小时，而采用本发明的方法只需要40分钟(0.67小时)，因此采用本方法后，运行速度提高了396倍。从而很好满足了三维处理的时间要求。
本发明提出的快速生成三维共检波点初至道集的方法计算速度快，生成共检波点道集时间短，能提高三维剩余静校正的计算效率，缩短地震资料处理周期，具有广泛的适用性和通用性。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种快速生成三维共检波点初至道集的方法，其特征在于，包括：
步骤1，根据野外生产资料，获取并建立实验地区的SPS文件，其中包括：检波点文件、炮点文件及关系文件；
步骤2，根据所述SPS文件，获取每一炮点的一定偏移距范围内的初至时间，生成初至数据文件；
步骤3，在所述实验地区的检波点中，选取要生成所述三维共检波点初至道集的检波点作为共检波点；
步骤4，根据所述检波点文件及炮点文件，分别建立三维检波点线的索引指示文件及三维炮点线的索引指示文件，并根据所述三维检波点线的索引指示文件，计算生成所述共检波点在检波点文件中的位置数据；
步骤5，根据所述关系文件，建立每一炮点关于检波点线及检波点的排列关系的索引关系文件；
步骤6，设定线序号搜索范围，根据所述线序号搜索范围在所述炮点文件中选取并生成某一共检波点对应的炮点线序号搜索范围，再根据所述炮点线序号搜索范围及所述三维炮点线的索引指示文件，在所述炮点文件中选取并生成所述某一共检波点对应的炮点序号范围；
步骤7，根据所述某一共检波点对应的炮点序号范围以及每一炮点的索引关系文件，分别判断每一炮点的排列关系中是否包含所述某一共检波点，如果包含，根据所述共检波点在检波点文件中的位置数据、索引关系文件及初至数据文件，提取相关数据并计算生成所述某一共检波点的初至数据道集；
步骤8，重复所述步骤6、步骤7，计算每一共检波点的初至数据道集，生成三维共检波点初至道集。
2.根据权利要求1所述的快速生成三维共检波点初至道集的方法，其特征在于，在所述步骤1中，检波点文件按检波点桩号升序排序、炮点文件按炮点桩号升序排序。
3.根据权利要求2所述的快速生成三维共检波点初至道集的方法，其特征在于，在所述步骤4中，根据三维检波点线的索引指示文件，计算生成所述共检波点在检波点文件中的位置数据还包括：
根据所述三维检波点线的索引指示文件，计算某一共检波点在检波点文件中的位置，并提取所述某一共检波点的坐标，计算某一共检波点在检波点文件中的位置的公式如下：
Rloci,j＝Rloci,1+Rstai,1-Rstai,j；
其中，Rloci,j为所述某一共检波点在所述检波点文件中的位置，i为检波点线号，j为检波点桩号；
Rloci,1为检波点线第i线第1个检波点在三维检波点文件中的位置，其值在建立三维检波点线的索引指示文件时获得；
Rstai,1为检波点线第i线第1个检波点桩号；
Rstai,j为所述某一共检波点桩号。
4.根据权利要求3所述的快速生成三维共检波点初至道集的方法，其特征在于，所述步骤7还包括：
如果某一炮点的排列关系中包含所述某一共检波点，根据炮点文件提取所述某一炮点的坐标，根据所述三维检波点线的索引指示文件，计算某一共检波点在检波点文件中的位置，并提取所述某一共检波点的坐标，计算炮检距；
根据索引关系文件，计算所述某一共检波点相对所述某一炮点的排列位置，根据所述排列位置从所述初至数据文件中提取所述某一共检波点及某一炮点对应地震道的初至时间数据；
根据所述炮检距及初至时间数据，生成所述某一共检波点的初至数据道集。
5.根据权利要求3所述的快速生成三维共检波点初至道集的方法，其特征在于，所述步骤7还包括：
利用以下公式计算所述某一共检波点在所述某一炮点内的排列位置：
SR＝SRnumi,1+SRstai,1-Rstai,j；
其中，SR为所述某一共检波点相对该炮点的排列位置；
SRnumi,1表示在所述索引关系文件中，检波点线第i线第1个检波点相对所述某一炮点的排列位置，其值在建立索引关系文件时获得；
SRstai,1表示在所述索引关系文件中，检波点线第i线第1个检波点相对所述某一炮点的检波点桩号，其值在建立索引关系文件时获得；
Rstai,j为所述某一共检波点桩号。
CN201410682365.7A 2014-11-24 2014-11-24 一种快速生成三维共检波点初至道集的方法 Active CN104375181B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410682365.7A CN104375181B (zh) 2014-11-24 2014-11-24 一种快速生成三维共检波点初至道集的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410682365.7A CN104375181B (zh) 2014-11-24 2014-11-24 一种快速生成三维共检波点初至道集的方法

Publications (2)

Publication Number Publication Date
CN104375181A CN104375181A (zh) 2015-02-25
CN104375181B true CN104375181B (zh) 2017-03-08

Family

ID=52554209

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410682365.7A Active CN104375181B (zh) 2014-11-24 2014-11-24 一种快速生成三维共检波点初至道集的方法

Country Status (1)

Country Link
CN (1) CN104375181B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107677237B (zh) * 2017-09-27 2020-05-15 中国石油集团东方地球物理勘探有限责任公司 一种三维观测系统方位角的计算方法
CN108107471B (zh) * 2017-11-17 2019-10-11 中国石油天然气集团公司 一种分方位初至数据体的获取方法及装置
CN112241023B (zh) * 2019-07-16 2023-04-25 中国石油天然气集团有限公司 陆上三维炮点二次定位及定量偏移监控方法和装置
CN112764092B (zh) * 2019-10-21 2023-11-28 中国石油天然气集团有限公司 检波点位置定位方法及系统
CN115586575B (zh) * 2022-06-14 2024-01-30 中国地质大学（北京） 一种3d被动源拟炮集计算方法、系统、设备及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4596005A (en) * 1983-04-20 1986-06-17 Chevron Research Company Method of seismic collection utilizing multicomponent processing receivers and processing resultant conventional and converted P- or S-wave data
CN1365008A (zh) * 2001-01-19 2002-08-21 中国石油天然气股份有限公司 一种地震多域迭代静校正方法
CN101520516A (zh) * 2008-02-25 2009-09-02 中国石油集团东方地球物理勘探有限责任公司 一种三维地震记录的叠前似二维变换的方法
CN101630015A (zh) * 2008-07-16 2010-01-20 中国石油天然气集团公司 一种提高初至波拾取精度和效率的方法
CN102707318A (zh) * 2012-06-18 2012-10-03 彭苏萍 虚共炮点射线追踪动态静校正方法
CN103837895A (zh) * 2014-03-10 2014-06-04 中国石油集团川庆钻探工程有限公司地球物理勘探公司 拟合初至波获取短波长静校正量的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO170441C (no) * 1985-05-06 1992-10-14 Western Atlas Int Inc Fremgangsmaate for undertrykkelse av koherent stoey i marine seismikk-data

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4596005A (en) * 1983-04-20 1986-06-17 Chevron Research Company Method of seismic collection utilizing multicomponent processing receivers and processing resultant conventional and converted P- or S-wave data
CN1365008A (zh) * 2001-01-19 2002-08-21 中国石油天然气股份有限公司 一种地震多域迭代静校正方法
CN101520516A (zh) * 2008-02-25 2009-09-02 中国石油集团东方地球物理勘探有限责任公司 一种三维地震记录的叠前似二维变换的方法
CN101630015A (zh) * 2008-07-16 2010-01-20 中国石油天然气集团公司 一种提高初至波拾取精度和效率的方法
CN102707318A (zh) * 2012-06-18 2012-10-03 彭苏萍 虚共炮点射线追踪动态静校正方法
CN103837895A (zh) * 2014-03-10 2014-06-04 中国石油集团川庆钻探工程有限公司地球物理勘探公司 拟合初至波获取短波长静校正量的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"SPS文件恢复与地震资料重处理";马强 等;《石油物探》;20100525;第49卷(第3期);第306-311页 *
"交互迭代静校正方法";王彦春 等;《石油物探》;19980630;第37卷(第2期);第63-70页 *
"共炮（检）点剩余静校正方法";潘树林 等;《石油地球物理勘探》;20110215;第46卷(第1期);第83-88页 *

Also Published As

Publication number Publication date
CN104375181A (zh) 2015-02-25

Similar Documents

Publication Publication Date Title
CN104375181B (zh) 一种快速生成三维共检波点初至道集的方法
CN104483705B (zh) 一种三维剩余静校正方法
CN107290722B (zh) 微震源的定位方法和装置
CN104076392A (zh) 基于网格搜索和牛顿迭代的微震震源定位联合反演方法
CN103837895B (zh) 拟合初至波获取短波长静校正量的方法
CN102565857B (zh) 自动剩余动校正方法
CN104422963B (zh) 一种变偏移距vsp资料时差校正方法
CN106772577A (zh) 基于微地震数据和spsa优化算法的震源反演方法
CN104122587B (zh) 一种基于炮检距域的异常初至识别方法以及系统
CN104155691A (zh) 转换波各向异性速度分析方法及装置
CN103439738A (zh) 地震勘探单炮记录异常道识别方法
CN103713318A (zh) 对三维叠后地震数据进行拼接的方法
CN104182651B (zh) 用于三分量检波器接收的微地震事件方位角自动质控方法
CN105093319A (zh) 基于三维地震数据的地面微地震静校正方法
CN105319587A (zh) 一种基于初至子波提取的可控震源整形滤波方法
CN104391324A (zh) 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术
CN105116447A (zh) 一种基于曲率异常条带的地质河道方向判别方法
GB2526709A (en) Migration velocity analysis method for vertical seismic profile data
CN111273352A (zh) 地质构造智能检测方法、装置与电子设备
CN103592691B (zh) 一种反向剩余磁性磁异常识别方法及装置
CN105277981A (zh) 基于波场延拓补偿的非一致性时移地震面元匹配方法
CN109581506A (zh) 一种数据驱动的无拉伸动校正方法
CN112419493B (zh) 页岩储层三维属性模型建立方法及装置
CN105445788A (zh) 一种基于模型和全局寻优的速度谱自动解释方法
CN105353420B (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