CN112649866A - 基于面元均化的拟三维处理方法及系统 - Google Patents

基于面元均化的拟三维处理方法及系统 Download PDF

Info

Publication number
CN112649866A
CN112649866A CN201910969329.1A CN201910969329A CN112649866A CN 112649866 A CN112649866 A CN 112649866A CN 201910969329 A CN201910969329 A CN 201910969329A CN 112649866 A CN112649866 A CN 112649866A
Authority
CN
China
Prior art keywords
dimensional
determining
pseudo
surface element
observation system
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.)
Pending
Application number
CN201910969329.1A
Other languages
English (en)
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 Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201910969329.1A priority Critical patent/CN112649866A/zh
Publication of CN112649866A publication Critical patent/CN112649866A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/53Statics correction, e.g. weathering layer or transformation to a datum
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time

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

公开了一种基于面元均化的拟三维处理方法及系统。该方法可以包括:根据目标工区的二维测线资料,确定检波点道距,进而确定测线覆盖次数;确定三维网格面元,建立拟三维观测系统;根据三维网格面元,确定覆盖次数不均匀的测线;针对覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统。本发明通过面元均化的拟三维方法,解决了面元不一致的问题,按照拟三维的方式处理二维测线,统一进行静校正计算和速度分析,提高了处理效率,减少了二维测线的闭合差,处理剖面信噪比和保真度有所提高,为解释人员提供拟三维资料,使三维解释技术可以应用到二维。

Description

基于面元均化的拟三维处理方法及系统
技术领域
本发明涉及地球物理勘探地震资料处理领域,更具体地,涉及一种基于面元均化的拟三维处理方法及系统。
背景技术
同一工区不同年份采集的多条二维线资料,由于观测系统的差异,会造成面元大小不一致的问题,进行连片统一面元处理,会出现资料覆盖次数不均匀的现象,对振幅保真等带来不利影响。
并且,在进行同一工区多条二维线连片处理时,由于静校正量、速度分析等差异,往往会造成严重的交点闭合差,有时可达几十毫秒,处理人员需要对每条线进行闭合平差处理,这样往往会造成构造失真问题,影响解释人员的构造解释等。
对于面元不一致的二维测线资料,每条测线分别采用的不同面元处理,不能同时解决工区资料面元不一致、覆盖次数不均匀和交点闭合差等问题,不利于三维处理、解释技术作用的发挥,不利于二维资料潜力的挖掘。因此,有必要开发一种基于面元均化的拟三维处理方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种基于面元均化的拟三维处理方法及系统,其能够通过面元均化的拟三维方法,解决了面元不一致的问题,按照拟三维的方式处理二维测线,统一进行静校正计算和速度分析,提高了处理效率,减少了二维测线的闭合差,处理剖面信噪比和保真度有所提高,为解释人员提供拟三维资料,使三维解释技术可以应用到二维。
根据本发明的一方面,提出了一种基于面元均化的拟三维处理方法。所述方法可以包括:根据目标工区的二维测线资料,确定检波点道距,进而确定测线覆盖次数;确定三维网格面元,建立拟三维观测系统;根据三维网格面元,确定覆盖次数不均匀的测线;针对所述覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统。
优选地,还包括:针对目标工区的所有二维线分别进行观测系统定义,进而确定三维网格面元。
优选地,针对所述覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统包括:确定面元搜索范围;在所述覆盖次数不均匀的测线的相邻面元中确定面元缺失的炮检距组,进而进行均化处理,获得搜索均化后的拟三维观测系统。
优选地,还包括:针对所述均化后的拟三维观测系统,计算静校正量与速度。
优选地,所述三维网格面元为5-30m。
根据本发明的另一方面,提出了一种基于面元均化的拟三维处理系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:根据目标工区的二维测线资料,确定检波点道距,进而确定测线覆盖次数;确定三维网格面元,建立拟三维观测系统;根据三维网格面元,确定覆盖次数不均匀的测线;针对所述覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统。
优选地,还包括:针对目标工区的所有二维线分别进行观测系统定义,进而确定三维网格面元。
优选地,针对所述覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统包括:确定面元搜索范围;在所述覆盖次数不均匀的测线的相邻面元中确定面元缺失的炮检距组,进而进行均化处理,获得搜索均化后的拟三维观测系统。
优选地,还包括:针对所述均化后的拟三维观测系统,计算静校正量与速度。
优选地,所述三维网格面元为5-30m。
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的基于面元均化的拟三维处理方法的步骤的流程图。
图2示出了根据本发明的一个实施例的20米道距的测线定义为15米*15米面元的覆盖次数的分布示意图。
图3示出了根据图2的面元均化后的覆盖次数分布示意图。
图4示出了根据本发明的一个实施例的均化后的拟三维观测系统的检波点层析静校正量的示意图。
图5示出了根据本发明的一个实施例的拟三维观测系统的炮点层析静校正量的示意图。
图6示出了根据本发明的一个实施例的相交测线交点闭合差质量控制的示意图。
图7示出了根据本发明的一个实施例的二维处理后的叠加剖面的示意图。
图8示出了根据本发明的一个实施例的拟三维处理后的叠加剖面的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的基于面元均化的拟三维处理方法的步骤的流程图。
在该实施例中,根据本发明的基于面元均化的拟三维处理方法可以包括:步骤101,根据目标工区的二维测线资料,确定检波点道距,进而确定测线覆盖次数;步骤102,确定三维网格面元,建立拟三维观测系统;步骤103,根据三维网格面元,确定覆盖次数不均匀的测线;步骤104,针对覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统。
在一个示例中,还包括:针对目标工区的所有二维线分别进行观测系统定义,进而确定三维网格面元。
在一个示例中,针对覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统包括:确定面元搜索范围;在覆盖次数不均匀的测线的相邻面元中确定面元缺失的炮检距组,进而进行均化处理,获得搜索均化后的拟三维观测系统。
在一个示例中,还包括:针对均化后的拟三维观测系统,计算静校正量与速度。
在一个示例中,三维网格面元为5-30m。
具体地,根据本发明的基于面元均化的拟三维处理方法可以包括:
根据目标工区的二维测线资料,确定检波点道距,进而确定测线覆盖次数;针对目标工区的所有二维线分别进行观测系统定义,确定三维网格面元为5-30m,建立拟三维观测系统。根据三维网格面元,确定覆盖次数不均匀的测线,拟三维网格定义所有测线的面元必须一致,如检波点道距分别为20m、30m的二维测线资料,若统一定义拟三维网格面元为15m,则道距为20m的测线覆盖次数分布不均匀,覆盖次数会在90次和180次左右来回跳动,如图2所示。
针对覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统。面元均化是指在地震采集中,为适应扩大面元范围,从相邻的范围内选出面元内所缺少的反射点道,使不同面元的共中心点基本相同的处理方法。以当前面元为中心定义一个长方形的大面元为面元搜索范围;在一个圆范围内的相邻面元中确定面元缺失的炮检距组,进而进行均化处理,获得搜索均化后的拟三维观测系统,其中圆的半径是炮检距组的线性函数,其范围限定在面元搜索范围内。
针对均化后的拟三维观测系统,按照层析静校正方法来计算静校正量,由于每条线的速度是闭合的,可以按照三维方式进行速度分析,计算速度。
本方法通过面元均化的拟三维方法,解决了面元不一致的问题,按照拟三维的方式处理二维测线,统一进行静校正计算和速度分析,提高了处理效率,减少了二维测线的闭合差,处理剖面信噪比和保真度有所提高,为解释人员提供拟三维资料,使三维解释技术可以应用到二维。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
图2示出了根据本发明的一个实施例的20米道距的测线定义为15米*15米面元的覆盖次数的分布示意图。
图3示出了根据图2的面元均化后的覆盖次数分布示意图。
根据目标工区的二维测线资料,确定检波点道距,进而确定测线覆盖次数;针对目标工区的所有二维线分别进行观测系统定义,确定三维网格面元为15m,建立拟三维观测系统。根据三维网格面元,确定覆盖次数不均匀的测线,拟三维网格定义所有测线的面元必须一致,检波点道距为20m的二维测线资料的测线覆盖次数分布不均匀,覆盖次数会在90次和180次左右来回跳动,如图2所示;针对覆盖次数不均匀的测线进行面元均化,均化后的测线覆盖次数如图3所示,获得均化后的拟三维观测系统,针对均化后的拟三维观测系统,计算静校正量与速度。
图4示出了根据本发明的一个实施例的均化后的拟三维观测系统的检波点层析静校正量的示意图。
图5示出了根据本发明的一个实施例的拟三维观测系统的炮点层析静校正量的示意图。
图4与图5中所有二维线都是地表一致性的,每条测线同一点处(交点处)的静校正量都是相同的,构造更加符合实际地下的情况,更好合理,同时不用单独进行每条线的闭合差计算,提高了处理效率。
图6示出了根据本发明的一个实施例的相交测线交点闭合差质量控制的示意图。因为静校正量、速度等的差异,每两条二维线的交点,交点处可能会有时差,表现为交点处的同相轴不连续,有上下位移。图6中黑线即为两条二维线交点处,可以看到,黑线左右两边没有时差,就是说明交点闭合差几乎为0,闭合比较好,与实际情况相符。
图7示出了根据本发明的一个实施例的二维处理后的叠加剖面的示意图。
图8示出了根据本发明的一个实施例的拟三维处理后的叠加剖面的示意图。
图8与图7对比,进行了面元均化处理后的叠加剖面的信噪比更高,构造与图7的相同,说明采用拟三维网格计算后提高了信噪比,构造符合实际地下的情况。
综上所述,本发明通过面元均化的拟三维方法,解决了面元不一致的问题,按照拟三维的方式处理二维测线,统一进行静校正计算和速度分析,提高了处理效率,减少了二维测线的闭合差,处理剖面信噪比和保真度有所提高,为解释人员提供拟三维资料,使三维解释技术可以应用到二维。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例,提供了一种基于面元均化的拟三维处理系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:根据目标工区的二维测线资料,确定检波点道距,进而确定测线覆盖次数;确定三维网格面元,建立拟三维观测系统;根据三维网格面元,确定覆盖次数不均匀的测线;针对覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统。
在一个示例中,还包括:针对目标工区的所有二维线分别进行观测系统定义,进而确定三维网格面元。
在一个示例中,针对覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统包括:确定面元搜索范围;在覆盖次数不均匀的测线的相邻面元中确定面元缺失的炮检距组,进而进行均化处理,获得搜索均化后的拟三维观测系统。
在一个示例中,还包括:针对均化后的拟三维观测系统,计算静校正量与速度。
在一个示例中,三维网格面元为5-30m。
本系统通过面元均化的拟三维方法,解决了面元不一致的问题,按照拟三维的方式处理二维测线,统一进行静校正计算和速度分析,提高了处理效率,减少了二维测线的闭合差,处理剖面信噪比和保真度有所提高,为解释人员提供拟三维资料,使三维解释技术可以应用到二维。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (10)

1.一种基于面元均化的拟三维处理方法,其特征在于,包括:
根据目标工区的二维测线资料,确定检波点道距,进而确定测线覆盖次数;
确定三维网格面元,建立拟三维观测系统;
根据三维网格面元,确定覆盖次数不均匀的测线;
针对所述覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统。
2.根据权利要求1所述的基于面元均化的拟三维处理方法,其中,还包括:
针对目标工区的所有二维线分别进行观测系统定义,进而确定三维网格面元。
3.根据权利要求1所述的基于面元均化的拟三维处理方法,其中,针对所述覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统包括:
确定面元搜索范围;
在所述覆盖次数不均匀的测线的相邻面元中确定面元缺失的炮检距组,进而进行均化处理,获得搜索均化后的拟三维观测系统。
4.根据权利要求1所述的基于面元均化的拟三维处理方法,其中,还包括:
针对所述均化后的拟三维观测系统,计算静校正量与速度。
5.根据权利要求1所述的基于面元均化的拟三维处理方法,其中,所述三维网格面元为5-30m。
6.一种基于面元均化的拟三维处理系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
根据目标工区的二维测线资料,确定检波点道距,进而确定测线覆盖次数;
确定三维网格面元,建立拟三维观测系统;
根据三维网格面元,确定覆盖次数不均匀的测线;
针对所述覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统。
7.根据权利要求6所述的基于面元均化的拟三维处理系统,其中,还包括:
针对目标工区的所有二维线分别进行观测系统定义,进而确定三维网格面元。
8.根据权利要求6所述的基于面元均化的拟三维处理系统,其中,针对所述覆盖次数不均匀的测线进行面元均化,获得均化后的拟三维观测系统包括:
确定面元搜索范围;
在所述覆盖次数不均匀的测线的相邻面元中确定面元缺失的炮检距组,进而进行均化处理,获得搜索均化后的拟三维观测系统。
9.根据权利要求6所述的基于面元均化的拟三维处理系统,其中,还包括:
针对所述均化后的拟三维观测系统,计算静校正量与速度。
10.根据权利要求6所述的基于面元均化的拟三维处理系统,其中,所述三维网格面元为5-30m。
CN201910969329.1A 2019-10-12 2019-10-12 基于面元均化的拟三维处理方法及系统 Pending CN112649866A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910969329.1A CN112649866A (zh) 2019-10-12 2019-10-12 基于面元均化的拟三维处理方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910969329.1A CN112649866A (zh) 2019-10-12 2019-10-12 基于面元均化的拟三维处理方法及系统

Publications (1)

Publication Number Publication Date
CN112649866A true CN112649866A (zh) 2021-04-13

Family

ID=75343312

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910969329.1A Pending CN112649866A (zh) 2019-10-12 2019-10-12 基于面元均化的拟三维处理方法及系统

Country Status (1)

Country Link
CN (1) CN112649866A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102707316A (zh) * 2012-06-06 2012-10-03 中国石油天然气股份有限公司 复杂地质目标的2.5次三维地震勘探方法
US20190080512A1 (en) * 2017-09-12 2019-03-14 Microsoft Technology Licensing, Llc Three-dimensional graphics image processing
CN110187390A (zh) * 2019-06-17 2019-08-30 中国矿业大学(北京) 一种煤矿巷道平行测线立体地震观测与成像方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102707316A (zh) * 2012-06-06 2012-10-03 中国石油天然气股份有限公司 复杂地质目标的2.5次三维地震勘探方法
US20190080512A1 (en) * 2017-09-12 2019-03-14 Microsoft Technology Licensing, Llc Three-dimensional graphics image processing
CN110187390A (zh) * 2019-06-17 2019-08-30 中国矿业大学(北京) 一种煤矿巷道平行测线立体地震观测与成像方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张振平: "海上三维地震资料处理面元覆盖均化技术", 石油地球物理勘探, vol. 24, no. 6, pages 686 - 688 *
杨文君等: "拟三维层析反演技术在鄂尔多斯盆地黄土山地区研究与应用", 2018油气田勘探开发国际会议(IFEDC2018)论文集, pages 3 *
王西文等: "相对保幅的地震资料连片处理方法研究", 石油物探, vol. 45, no. 2, pages 105 - 110 *

Similar Documents

Publication Publication Date Title
CN110007357B (zh) 一种航空tem和航空mt联合反演方法
CN105093299A (zh) 一种基于炮检距向量片技术优化观测系统的方法及装置
US12007515B2 (en) Optimal survey design
CN104237940A (zh) 一种基于动力学特征的绕射波成像方法及装置
CN105301641B (zh) 一种方位各向异性速度反演的方法及装置
US9366771B2 (en) System and method for migration velocity modeling
US8706462B2 (en) System and method for providing a physical property model
CN103454680A (zh) Walk-away VSP观测系统垂向覆盖次数的计算方法
CN110361792A (zh) 一种地球物理数据融合及成像方法、介质与设备
CN104570088A (zh) 一种在共偏移距共方位角域进行多方位层析速度建模方法
Cui et al. The accurate estimation of GPR migration velocity and comparison of imaging methods
CN104122586A (zh) 一种地震勘探无线节点数据监控方法
CN105093300A (zh) 一种地质体边界识别方法及装置
CN104155703A (zh) 评价三维观测系统的方法和设备
CN112649866A (zh) 基于面元均化的拟三维处理方法及系统
CN109581494B (zh) 叠前偏移方法及装置
CN115754007A (zh) 一种基于声发射技术和层析成像技术的损伤检测方法
CN108957539A (zh) 层析偏移速度分析中的射线追踪方法及装置
KR102205080B1 (ko) 가시성이 향상된 gpr 데이터의 처리방법
CN110058318B (zh) 一种基于测量电极坐标的天然电场的一维方向校正方法
EP2784553A2 (en) Predicting sensitivity to positioning for seismic surveys
Sadala Valente et al. Time-to-depth conversion and velocity estimation by image-wavefront propagation
CN111045090B (zh) 一种磁异常网格快速缝合方法
CN115993641B (zh) 一种提取被动源面波频散曲线的方法
Zhang et al. Research on Crack Detection Technology of Buildings After Earthquake Based on Structured Light

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