CN103777236A - 一种可控震源同步扫描波场正演数值模拟方法 - Google Patents
一种可控震源同步扫描波场正演数值模拟方法 Download PDFInfo
- Publication number
- CN103777236A CN103777236A CN201210408060.8A CN201210408060A CN103777236A CN 103777236 A CN103777236 A CN 103777236A CN 201210408060 A CN201210408060 A CN 201210408060A CN 103777236 A CN103777236 A CN 103777236A
- Authority
- CN
- China
- Prior art keywords
- vibroseis
- simulation
- time
- parameter
- synchronous scanning
- 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
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种可控震源同步扫描波场正演数值模拟方法,属于石油勘探领域。本方法包括:(1)输入模拟参数;(2)根据步骤(1)输入的模拟参数计算获得可控震源同步扫描模拟参数,并将这些参数输出到参数文件中;所述可控震源同步扫描模拟参数包括:每一个震源位置、起始与终止扫描时间及其属于哪个同步源、连续记录时间长度和每一个同步元的扫描信号;(3)通过高阶有限差分波场延拓计算得到模拟记录;(4)输出模拟结果。本发明方法能够根据实际模拟的需要,通过输入参数的可控,模拟滑动扫描、独立同步扫描、伪随机扫描、距离拉开同步扫描技术和V1等多种方式的可控震源同步扫描、任意多个可控震源的正演模拟连续记录。
Description
技术领域
本发明属于石油勘探领域,具体涉及一种可控震源同步扫描波场正演数值模拟方法。
背景技术
可控震源同步扫描(又称为同步激发)技术是近年来发展起来的一种高效地震采集技术。由于在采集过程中不考虑各个震源扫描信号之间的相互影响,因此该技术可以大大缩短采集周期,降低采集成本。目前该技术正逐步走向实用,但其资料处理仍然存在着各种方法技术的攻关难题。而可控震源同步扫描正演数值模拟技术的开发对于这种处理技术的攻关具有重要的意义。
目前常规地震波场数值模拟技术已发展到岩性与油藏模拟阶段,基于弹性波方程的炸药震源和可控震源扫描的常规地震波场正演数值模拟技术已极为完善,但未见有关可控震源同步扫描正演模拟技术的文献和报道,目前在进行多震源地震方法技术研究时大多采用由常规炸药震源模拟数据(或常规采集数据)的合成多震源地震数据或实际野外可控震源同步扫描的实际数据,由于这种合成数据不同于真正的多震源数据,而实际资料又过于复杂,从而影响了多震源地震处理技术的进度。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种可控震源同步扫描波场正演数值模拟方法,开发基于弹性波方程的2D、3D可控震源同步扫描连续记录波场正演数值模拟技术,可以模拟滑动扫描、独立同步扫描、伪随机扫描、距离拉开同步扫描技术和V1等多种主要可控震源同步扫描连续地震记录,为这些可控震源同步扫描地震技术的方法研发提供实验数据。
本发明是通过以下技术方案实现的:
一种可控震源同步扫描波场正演数值模拟方法,所述方法包括以下步骤:
(1)输入模拟参数:
所述模拟参数包括常规观测系统参数和可控震源同步扫描参数;所述常规观测系统参数包括炮间距、道间距、最大、最小偏移距和时间采样间隔;所述可控震源同步扫描参数包括同步震源个数,扫描信号长度,最低频率,最高频率、听时间长度,最大随机时差以及各个震源扫描信号的分段数;
(2)根据步骤(1)输入的模拟参数计算获得可控震源同步扫描模拟参数,并将这些参数输出到参数文件中;所述可控震源同步扫描模拟参数包括:每一个震源位置、起始与终止扫描时间及其属于哪个同步源、连续记录时间长度和每一个同步元的扫描信号;
(3)通过高阶有限差分波场延拓计算得到模拟记录:
根据步骤(1)输入的模拟参数和步骤(2)计算得到的可控震源同步扫描模拟参数,利用高阶有限差分形式的波动方程进行正演模拟计算得到模拟结果;
(4)输出模拟结果:
将步骤(3)得到的模拟结果采用野外地震采集解编前的时间切片记录方式进行记录并输出。
所述步骤(2)包括以下步骤:
(A)根据同步源个数将地面炮点区域划分为与同步源个数相同的、形状、大小近似相同的子区域,并对这些子区域进行编号;每个子区域采用一个扫描信号;所述扫描信号根据子区域编号、最低频率,最高频率和总的分段时间自动生成;
(B)按照常规观测系统参数计算每一个子区域中的震源位置,并按顺序对同一个子区域的震源位置进行独立编号;同时利用随机生成函数,根据最大随机时差生成每一炮的随机时差,并与其扫描信号长度和听时间长度相加生成每一炮的完成时间长度;按照震源位置编号将每一炮的完成时间长度进行累加,从而得到同一子区域的各炮的起始与终止扫描时间;完成所有子区域的计算后得到所有炮的起始与终止时间(即所有震源位置扫描信号的起始与终止时间)以及每一个子区域的最长连续记录时间;
(C)比较所有子区域的最长连续记录时间的长度,选择其中最长的作为模拟连续记录时间长度。
所述步骤(3)中,在正演模拟计算时模拟计算区域的边界采用的是高阶有限差分形式最佳匹配层吸收边界条件方程;在正演模拟计算时的循环顺序是直接进行时间层的循环,且循环时间长度为步骤(C)中得到的模拟连续记录时间长度。
与现有技术相比,本发明的有益效果是:本发明方法能够根据实际模拟的需要,通过输入参数的可控,可以模拟滑动扫描、独立同步扫描、伪随机扫描、距离拉开同步扫描技术和V1等多种主要方式的可控震源同步扫描、任意多个可控震源的正演模拟连续记录,从而为多震源地震方法技术的研发提供实验用数据和检测方法技术正确性的数据。
附图说明
图1是本发明可控震源同步扫描波场正演数值模拟方法的步骤框图。
图2是本发明实施例中的两个同步扫描震源分别在两个不同位置扫描弹性波模拟垂直分量连续记录。
图3-1是本发明实施例中第一个震源的扫描信号。
图3-2是本发明实施例中第二个震源的扫描信号。
具体实施方式
下面结合附图对本发明作进一步详细描述:
本发明以波动理论为指导,以弹性波方程为基础(体现在本发明方法的步骤(3)中,即:采用波动方程高阶有限差分正演模拟技术进行模拟计算),利用交错网格高阶有限差分法对弹性波方程进行离散化,边界条件采用最佳匹配层吸收边界条件,开发出可控震源同步扫描高精度正演数值模拟技术,同步源之间可以有任意时差,扫描信号可以是线性的,也可以是非线性的,还可以是分段的不连续信号(如图3-1和图3-2所示的扫描信号)。可以提供多种主要可控震源同步扫描模拟连续地震记录。本发明方法的主要创新点是:利用波动方程高阶有限差分法模拟各种多震源地震野外采集原始地震记录,为多震源地震技术的研发提供实验数据。
本发明可控震源同步扫描波场正演数值模拟方法的步骤如图1所示,包括:
(1).模拟参数输入
输入如下参数:
常规观测系统参数:炮间距,道间距,最大、最小偏移距,时间采样间隔。可控震源同步扫描参数:同步震源个数,扫描信号长度,最低频率,最高频率、听时间长度,最大随机时差以及各个震源扫描信号的分段数和总的分段时间。
(2).可控震源同步扫描模拟参数计算
根据输入参数计算如下参数,并将这些参数输出到参数文件中。
具体计算方法是:
(A)根据同步源个数将地面炮点区域划分为与同步源个数相同的、形状、大小近似相同的子区域(二维情况下,这些子区域都是一个线段;三维情况下,这些子区域都是长方形),并对这些子区域进行编号。每个子区域采用一个扫描信号。扫描信号根据子区域编号、最低频率,最高频率和总的分段时间自动生成。
(B)按照观测系统参数计算每一个子区域中的震源位置,并按顺序对同一个子区域的震源位置进行独立编号。同时利用随机生成函数,根据最大随机时差生成每一炮的随机时差,并与其扫描信号长度和听时间长度相加生成每一炮的完成时间长度;按照震源位置编号将每一炮的完成时间长度进行累加,从而得到同一子区域各炮的起始、终止扫描时间。每一子区域均进行这种计算,从而得到所有震源位置扫描信号的起始、终止时间以及每一个子区域的最长连续记录时间。
(C)对所有子区域的记录时间长度进行对比,从而得到模拟连续记录时间长度。
经过上述三步的计算即得到可控震源同步扫描模拟参数:每一个震源位置、起始与终止扫描时间及其属于哪个同步源,连续记录时间长度,每一个同步源的扫描信号。
(3)高阶有限差分波场延拓计算
地震正演模拟计算都需要使用各种不同形式的边界条件以消除模拟计算边界所产生的反射,因为这种反射在实际地震勘探中是不存在的。这一步主要是根据输入参数和计算得到的可控震源同步扫描模拟参数,利用高阶有限差分形式的波动方程进行模拟正演模拟计算(可参考的文献很多,例如可参考董良国,马在田,曹景中等的“一阶弹性波交错网格高阶差分解法”【J】.地球物理学报,2000,43(3):411~419),在模拟计算区域的边界,采用高阶有限差分形式最佳匹配层吸收边界条件方程(该方程是现有的,参考文献有:Berenger,J.,A perfectly matched layer for the absorption of electromagneticwaves[J],J.Comput.Phys.,1994,114(1),185-200。)来消除边界反射。但与常规单炮模拟相比,其循环顺序有较大变化。常规单炮模拟最外层是炮循环,其次是时间循环,最后是模型空间网格点的循环,而本发明中可控震源同步扫描没有炮循环这一层,直接进行时间层的循环,且循环时间长度由常规单炮模拟的单炮记录长度变为全部模拟结束的连续记录时间长度。
(4)输出模拟结果。
常规单炮模拟形成的是记录,由于多震源同步扫描连续记录时间非常长,受计算机内存的制约,不可能采用一炮一炮的记录形式,因此本发明中的模拟结果是采用野外地震采集解编前的时间切片记录方式进行记录,然后输出。
图2是两个同步源扫描2D弹性波模拟垂直分量连续记录。其中S11和S12所指示为的第一个震源在两个不同位置的起始激发时间,S21和S22所指示为第一个震源在两个不同位置的起始激发时间。每一个震源各自在指定的不同位置各扫描了一次,并连续记录在一个记录中。图3-1是第一个震源的扫描信号,图3-2是第二个震源的扫描信号,两个信号为同步源,图中的纵坐标是扫描信号振幅,横坐标是时间。从图3-1和图3-2可以看出两个震源的扫描信号各不相同,而且每个扫描信号的中间独有间断。实际上这两个扫描信号也可以是相同的,也可以是不间断的,间断次数等都可以由输入参数控制。
可控震源同步扫描技术是近几年来发展起来的地震勘探资料高效野外采集技术,但其资料处理还面临这一系列的技术攻关,本发明可控震源同步扫描波场正演数值模拟方法提供了认识这种复杂地震波场的工具,另一方面为进行可控震源同步扫描资料处理技术研发提供了能计算试验和检验数据的工具。
为此本发明以弹性波方程为基础,以交错网格和旋转交错网格的高阶有限差分算法为工具,结合最佳匹配层吸收边界条件,开发了2D/3D高精度可控震源同步扫描正演数值模拟技术,可以根据实际需要,模拟以滑动扫描、独立同步扫描和伪随机扫描等为主的多种方式的可控震源同步扫描连续记录采集数据
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (3)
1.一种可控震源同步扫描波场正演数值模拟方法,其特征在于:所述方法包括以下步骤:
(1)输入模拟参数:
所述模拟参数包括常规观测系统参数和可控震源同步扫描参数;所述常规观测系统参数包括炮间距、道间距、最大、最小偏移距和时间采样间隔;所述可控震源同步扫描参数包括同步震源个数,扫描信号长度,最低频率,最高频率、听时间长度,最大随机时差以及各个震源扫描信号的分段数;
(2)根据步骤(1)输入的模拟参数计算获得可控震源同步扫描模拟参数,并将这些参数输出到参数文件中;所述可控震源同步扫描模拟参数包括:每一个震源位置、起始与终止扫描时间及其属于哪个同步源、连续记录时间长度和每一个同步元的扫描信号;
(3)通过高阶有限差分波场延拓计算得到模拟记录:
根据步骤(1)输入的模拟参数和步骤(2)计算得到的可控震源同步扫描模拟参数,利用高阶有限差分形式的波动方程进行正演模拟计算得到模拟结果;
(4)输出模拟结果:
将步骤(3)得到的模拟结果采用野外地震采集解编前的时间切片记录方式进行记录并输出。
2.根据权利要求1所述的可控震源同步扫描波场正演数值模拟方法,其特征在于:所述步骤(2)包括以下步骤:
(A)根据同步源个数将地面炮点区域划分为与同步源个数相同的、形状、大小近似相同的子区域,并对这些子区域进行编号;每个子区域采用一个扫描信号;所述扫描信号根据子区域编号、最低频率,最高频率和总的分段时间自动生成;
(B)按照常规观测系统参数计算每一个子区域中的震源位置,并按顺序对同一个子区域的震源位置进行独立编号;同时利用随机生成函数,根据最大随机时差生成每一炮的随机时差,并与其扫描信号长度和听时间长度相加生成每一炮的完成时间长度;按照震源位置编号将每一炮的完成时间长度进行累加,从而得到同一子区域的各炮的起始与终止扫描时间;完成所有子区域的计算后得到所有炮的起始与终止时间以及每一个子区域的最长连续记录时间;
(C)比较所有子区域的最长连续记录时间的长度,选择其中最长的作为模拟连续记录时间长度。
3.根据权利要求1所述的可控震源同步扫描波场正演数值模拟方法,其特征在于:所述步骤(3)中,在正演模拟计算时模拟计算区域的边界采用的是高阶有限差分形式最佳匹配层吸收边界条件方程;在正演模拟计算时的循环顺序是直接进行时间层的循环,且循环时间长度为步骤(C)中得到的模拟连续记录时间长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210408060.8A CN103777236A (zh) | 2012-10-24 | 2012-10-24 | 一种可控震源同步扫描波场正演数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210408060.8A CN103777236A (zh) | 2012-10-24 | 2012-10-24 | 一种可控震源同步扫描波场正演数值模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103777236A true CN103777236A (zh) | 2014-05-07 |
Family
ID=50569696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210408060.8A Pending CN103777236A (zh) | 2012-10-24 | 2012-10-24 | 一种可控震源同步扫描波场正演数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103777236A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105259580A (zh) * | 2015-10-30 | 2016-01-20 | 中国石油大学(北京) | 一种可控震源信号低频拓展方法 |
CN106126316A (zh) * | 2016-06-23 | 2016-11-16 | 哈尔滨工业大学 | 基于随机点配置方法的电磁仿真快速参数扫描方法 |
CN106199691A (zh) * | 2015-05-26 | 2016-12-07 | 中石化石油工程地球物理有限公司胜利分公司 | 基于可控震源滑动扫描方法的并行正演模拟方法 |
CN106199713A (zh) * | 2016-07-07 | 2016-12-07 | 崔宏良 | 一种地震勘探数值模拟的分析方法 |
CN106597523A (zh) * | 2016-12-02 | 2017-04-26 | 中国石油天然气集团公司 | 一种可控震源动态扫描方法及装置 |
CN110568485A (zh) * | 2019-09-06 | 2019-12-13 | 广州海洋地质调查局 | 一种基于神经网络的多道地震连续记录分离方法 |
CN112782754A (zh) * | 2019-11-06 | 2021-05-11 | 中国石油天然气集团有限公司 | 多源同步采集数据模拟方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008110743A2 (en) * | 2007-03-09 | 2008-09-18 | Geco Technology B.V. | Vibroseis acquisition method |
CN101980052A (zh) * | 2010-09-28 | 2011-02-23 | 中国科学院地质与地球物理研究所 | 叠前逆时偏移成像的方法及装置 |
-
2012
- 2012-10-24 CN CN201210408060.8A patent/CN103777236A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008110743A2 (en) * | 2007-03-09 | 2008-09-18 | Geco Technology B.V. | Vibroseis acquisition method |
CN101980052A (zh) * | 2010-09-28 | 2011-02-23 | 中国科学院地质与地球物理研究所 | 叠前逆时偏移成像的方法及装置 |
Non-Patent Citations (6)
Title |
---|
CLAUDIO BAGAINI,ET AL: "Overview of simultaneous Vibroseis acquisition methods", 《SEG/NEW ORLEANS 2006 ANNUAL MEETING》, 31 December 2006 (2006-12-31), pages 70 - 74 * |
佘德平: "多震源地震正演数值模拟技术", 《地球物理学进展》, vol. 27, no. 4, 31 August 2012 (2012-08-31) * |
左莹: "基于高阶交错网格的有限差分地震波场数值模拟", 《中国优秀硕士学位论文全文数据库 基础科学辑》, no. 2, 15 February 2010 (2010-02-15) * |
爱德华 A. 博蒙特 等: "《油气圈闭勘探》", 28 February 2002 * |
王井富等: "可控震源高效采集技术简介及对装备的需求", 《物探装备》, vol. 20, no. 2, 30 April 2010 (2010-04-30), pages 106 - 109 * |
肖忠祥 等: "《数据采集原理》", 28 February 2001 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106199691A (zh) * | 2015-05-26 | 2016-12-07 | 中石化石油工程地球物理有限公司胜利分公司 | 基于可控震源滑动扫描方法的并行正演模拟方法 |
CN105259580A (zh) * | 2015-10-30 | 2016-01-20 | 中国石油大学(北京) | 一种可控震源信号低频拓展方法 |
CN105259580B (zh) * | 2015-10-30 | 2017-12-26 | 中国石油大学(北京) | 一种可控震源信号低频拓展方法 |
CN106126316A (zh) * | 2016-06-23 | 2016-11-16 | 哈尔滨工业大学 | 基于随机点配置方法的电磁仿真快速参数扫描方法 |
CN106126316B (zh) * | 2016-06-23 | 2019-06-11 | 哈尔滨工业大学 | 基于随机点配置方法的电磁仿真快速参数扫描方法 |
CN106199713A (zh) * | 2016-07-07 | 2016-12-07 | 崔宏良 | 一种地震勘探数值模拟的分析方法 |
CN106597523A (zh) * | 2016-12-02 | 2017-04-26 | 中国石油天然气集团公司 | 一种可控震源动态扫描方法及装置 |
CN110568485A (zh) * | 2019-09-06 | 2019-12-13 | 广州海洋地质调查局 | 一种基于神经网络的多道地震连续记录分离方法 |
CN110568485B (zh) * | 2019-09-06 | 2021-02-26 | 广州海洋地质调查局 | 一种基于神经网络的多道地震连续记录分离方法 |
CN112782754A (zh) * | 2019-11-06 | 2021-05-11 | 中国石油天然气集团有限公司 | 多源同步采集数据模拟方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103777236A (zh) | 一种可控震源同步扫描波场正演数值模拟方法 | |
CN103238158B (zh) | 利用互相关目标函数进行的海洋拖缆数据同时源反演 | |
CN106526674A (zh) | 一种三维全波形反演能量加权梯度预处理方法 | |
CN102645670B (zh) | 一种基于叠加响应分析的观测系统优化设计方法 | |
CN102707316B (zh) | 复杂地质目标的2.5次三维地震勘探方法 | |
CN106772579B (zh) | 一种薄煤层中地震叠前反演方法和装置 | |
CN106842295A (zh) | 测井信息约束的波形反演方法 | |
CN108445533A (zh) | 基于obs的长偏移距构建浅水区速度模型的方法 | |
CN106415321A (zh) | 用于储层建模的基于瞬时等时属性的地质体识别 | |
Krahenbuhl et al. | Understanding the applications and limitations of time-lapse gravity for reservoir monitoring | |
CN103592698A (zh) | 一种基于地震属性的观测系统评价方法及装置 | |
CN109212589A (zh) | 一种基于gpu/cpu协同并行的地震采集观测系统设计方法 | |
CN103913768A (zh) | 基于地震波资料对地表中浅层进行建模的方法及装置 | |
CN105277976A (zh) | 基于岩石露头雷达探测的地震正演模拟方法 | |
CN107894618A (zh) | 一种基于模型平滑算法的全波形反演梯度预处理方法 | |
CN103119472B (zh) | 利用同时和顺序源方法进行全波形反演的混合方法 | |
Herrmann et al. | Frugal full-waveform inversion: From theory to a practical algorithm | |
CN104391324A (zh) | 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术 | |
CN102914790B (zh) | 二维观测系统和三维观测系统一次采集的观测系统方法 | |
Wang et al. | Velocity model estimation of karstic fault reservoirs using full-waveform inversion accelerated on graphics processing unit | |
Jiang et al. | Velocity model optimization for surface microseismic monitoring via amplitude stacking | |
CN102096099A (zh) | 折射波射线层析成像方法 | |
Zhao et al. | Re-Exploration Programs for Petroleum-Rich Sags in Rift Basins | |
Zhao et al. | An automatical infill shot method for uniform imaging of target layer | |
Roach et al. | An assessment of the time-lapse seismic repeatability using a permanent array for reservoir monitoring at the Aquistore CO2 storage site, Saskatchewan, Canada |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140507 |