CN110888167B - 一种二维测线地震数据闭合差的校正方法及处理终端 - Google Patents
一种二维测线地震数据闭合差的校正方法及处理终端 Download PDFInfo
- Publication number
- CN110888167B CN110888167B CN201911014861.4A CN201911014861A CN110888167B CN 110888167 B CN110888167 B CN 110888167B CN 201911014861 A CN201911014861 A CN 201911014861A CN 110888167 B CN110888167 B CN 110888167B
- Authority
- CN
- China
- Prior art keywords
- measuring
- lines
- measuring lines
- line
- group
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000012545 processing Methods 0.000 title claims abstract description 10
- 238000010276 construction Methods 0.000 claims abstract description 35
- 238000012163 sequencing technique Methods 0.000 claims abstract description 8
- 238000005259 measurement Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 6
- 230000007547 defect Effects 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
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
本发明涉及一种二维测线地震数据闭合差的校正方法及处理终端，所述方法包括如下步骤：步骤1：获得二维测线地震数据，并对二维测线按施工特征进行排序，相同施工特征的二维测线的排序相同；步骤2：根据排序结果，对所有测线进行分组，以使得同一组的二维测线不相交，不同组之间的二维测线相交，从而得到若干组测线；步骤3：计算当前组测线与下一组测线之间的校正量，根据校正量对当前组测线与下一组测线之间进行校正，校正完后将两组测线进行合并为同一组测线，并作为当前组测线继续与下一组测线之间进行校正，依次进行，完成校正。本发明消除闭合差的效率和精度高，能够低成本第自动进行校正，特别适合大量地震数据的闭合差校正。
Description
技术领域
本发明涉及地震数据闭合差校正技术领域，具体涉及一种二维测线地震数据闭合差的校正方法及处理终端。
背景技术
在地震解释工作中，需要在由测线测量的地震剖面上将层位线和断层刻画出来，也即需要在每个地震剖面上画线，画线的结果为包括三轴方向的xyz数据，其中xy数据为测线所在平面的位置信息，z数据为对应平面上纵向值(时间域的时间值或深度域的深度值)。将不同地震剖面同一层位的解释数据经过网格化处理后，可制作多种具有生产指导的工业化图件，如构造图、油藏分布图等。
地震剖面数据分为二维地震数据和三维地震数据，对于二维测线的地震数据，其采集为单线采集，二维测线在测线交点处存在差异，该差异为点的闭合差，也即不同测线上的数据之间存在连续性问题，闭合差的存在会影响后续地震数据进行网格化结果，因此对闭合差进行校正，以消除闭合差。
目前闭合差的校正方法，主要包括两方面的不足：1)常用的闭合差校正方法中，通过计算出校正量，然后进行闭合差校正，对于简单闭合差的情况(例如二维测线只有横纵方向测线)，因为只需要进行一次校正，校正结果误差在允许的误差范围内，但对于复杂闭合差的情况，需要进行多次校正，虽然每一次的校正是在允许的误差范围内，但累加多次的校正后，往往最终的校正结果却超出允许的误差范围，无法满足要求；2)校正过程的测线分组往往依靠人工逐次筛选进行分组，分组简单，效率非常低下。
发明内容
针对现有技术的不足，本发明的目的之一提供一种二维测线地震数据闭合差的校正方法，其能够解决闭合差校正效率低的问题；
本发明的目的之二提供一种处理终端，其能够解决闭合差校正效率低的问题；
实现本发明的目的之一的技术方案为：一种二维测线地震数据闭合差的校正方法，包括如下步骤：
步骤1：获得二维测线地震数据，并对二维测线按施工特征进行排序，相同施工特征的二维测线的排序相同，施工特征是指进行二维测线施工作业时的施工特征；
步骤2：根据排序结果，对所有测线进行分组，以使得同一组的二维测线不相交，不同组之间的二维测线相交，从而得到若干组测线；
步骤3：从第一组测线开始，计算当前组测线与下一组测线之间的校正量，根据校正量对当前组测线与下一组测线之间进行校正，校正完后将两组测线进行合并为同一组测线，合并后的同一组测线作为当前组测线继续与下一组测线之间进行校正，依次进行，直至当前组测线与最后一组测线完成校正，从而消除闭合差。
进一步地，所述施工特征为施工时间、施工时的测线角度、施工批次中的任一个。
进一步地，所述步骤2中，对所有测线进行分组的具体过程包括如下步骤：
步骤21：从所有测线中选取排序最靠前的二维测线作为第一根测线，从第一根测线开始，选取第i根测线作为未分组测线，i＝1,2,3,……，n，n为二维测线总数，比较第i+1根测线和第i根测线，若第i+1根测线和第i根测线不相交，则将第i+1根测线标记为与第i根测线为同一组测线，否则，将第i+1根测线标记为与第i根测线为相异组测线，依次进行，遍历所有测线，找出所有与第i根测线不相交的测线，从而得到与第i根测线均不相交的所有测线并作为第一组测线，剩下的测线均为与第i根测线均相交的相异组测线；
步骤22：从与第i根测线相交的相异组测线中选取排序最靠前的二维测线作为第一根测线，重复执行步骤21,得到第二组测线，依次进行，直至将所有测线均分组完成，分别得到第一组测线、第二组测线、……、第m组测线。
进一步地，所述根据校正量对当前组测线与下一组测线之间进行校正，其具体实现过程包括如下步骤：根据当前组测线与下一组测线的相交关系得到交点处的校正量，将所述校正量线性插值到每根测线的所有数据点上，然后将每根测线的数据点加上对应的校正量，完成校正。
实现本发明目的之二的技术方案为：一种处理终端，其包括，
存储器，用于存储程序指令；
处理器，用于运行所述程序指令，以执行所述二维测线地震数据闭合差的校正方法中的步骤。
本发明的有益效果为：本发明消除闭合差的效率和精度高，能够低成本的自动进行校正，特别适合大量地震数据的闭合差校正。
附图说明
图1为本发明的流程示意图；
图2为未对测线进行分组的原始二维测线示意图；
图3为对图2的二维测线进行分组后的二维测线示意图；
图4为对各组测线进行校正的流程示意图；
图5为未消除闭合差直接成图的地震剖面示意图；
图6为消除闭合差的地震剖面示意图；
图7为本发明一种处理终端的结构示意图。
具体实施方式
下面，结合附图以及具体实施方案，对本发明做进一步描述。
如图1至图6所示，一种二维测线地震数据闭合差的校正方法，包括如下步骤：
步骤1：获得二维测线地震数据，并对二维测线按施工特征进行排序，相同施工特征的二维测线的排序相同，施工特征是指进行二维测线施工作业时的施工特征，包括施工时间、施工时的测线角度、施工批次等。例如，一批二维测线地震数据中，包括十条二维测线，分别记为测线g0-g9，其中测线g0、g2、g6、g7、g8为同一批次的测线，施工批次名记为0，g1、g5和g9为同一批次的测线，施工批次名记为1，剩下的测线为同一批次的测线，施工批次名记为2。10条二维测线按施工批次得到的排序如表一所示，同样的，若按施工的测线角度则得到表二：
| 二维测线名 | 施工批次 |
| g0 | 0 |
| g2 | 0 |
| g6 | 0 |
| g7 | 0 |
| g8 | 0 |
| g1 | 1 |
| g5 | 1 |
| g9 | 1 |
| g3 | 2 |
| g4 | 2 |
表一
| 二维测线名 | 施工的二维测线角度 |
| g1 | 45.1 |
| g3 | 45.6 |
| g8 | 45.6 |
| g0 | 45.7 |
| g6 | 45.8 |
| g5 | 135.6 |
| g2 | 135.7 |
| g4 | 135.8 |
| g7 | 136.7 |
| g9 | 179.1 |
表二
步骤2：根据步骤1的排序结果，对所有测线进行分组，以使得同一组的二维测线不相交，也即同一组的二维测线之间不存在闭合差，只有组与组之间的二维测线之间存在闭合差，从而得到若干组测线。分组的具体过程包括如下步骤：
步骤21：从所有测线中选取排序最靠前的二维测线作为第一根测线，从第一根测线开始，选取第i根测线作为未分组测线，i＝1,2,3,……，n，n为二维测线总数，比较第i+1根测线和第i根测线，若第i+1根测线和第i根测线不相交，则将第i+1根测线标记为与第i根测线为同一组测线，否则，将第i+1根测线标记为与第i根测线为相异组测线，依次进行，遍历所有测线，找出所有与第i根测线不相交的测线，从而得到与第i根测线均不相交的所有测线并作为第一组测线，剩下的测线均为与第i根测线均相交的相异组测线；
步骤22：从与第i根测线相交的相异组测线中选取排序最靠前的二维测线作为第一根测线，重复执行步骤21,得到第二组测线，直至将所有测线均分组完成，分别得到第一组测线、第二组测线、……、第m组测线，其中，同一组测线内的各根测线不相交，组与组之间的测线相交。
例如，经过步骤1进行排序后有5根测线，各根测线排序均不相同，按排序顺序的各根测线分别记为0号测线、1号测线、2号测线、3号测线、4号测线，从排序最靠前的0号测线开始，比较1号测线与0号测线是否相交，若1号测线与0号测线不相交，则1号测线与0号测线为同一组测线，否则，1号测线与0号测线为相异组测线，依次进行，分别比较2号线与0号线、3号线与0号线、4号线与0号线，假设经此轮比较后，只有1号线和3号线与0号为同一组，与0号测线相交的相异组测线包括2号测线和4号测线，即第一组测线为包括0号测线、1号测线和3号测线，
然后，从相异组测线中的2号测线和4号测线，选取排序靠前的2号测线开始，比较4号测线和2号测线是否相交，若不相交，则2号测线和4号测线为同一组测线，得到第二组测线，否则，2号测线和4号测线分别作为第二组测线和第三组测线。
如图2和图3所示，图2为未对测线进行分组的原始二维测线示意图，图3为对图2的二维测线进行分组后的二维测线示意图。
步骤3：如图4所示，从第一组测线开始，计算当前组测线与下一组测线之间的校正量，根据校正量对当前组测线与下一组测线之间进行校正，校正完后将两组测线进行合并为同一组测线，合并后的同一组测线作为当前组测线继续与下一组测线之间进行校正，依次进行，直至当前组测线与最后一组测线完成校正，从而消除闭合差。其中，当前组测线与下一组测线之间的校正量采用线性插值方法生成，得到每根测线上的所有数据点的校正量，也即根据当前组测线与下一组测线的相交关系得到交点处的校正量，将所述校正量线性插值到每根测线的所有数据点上，然后将每根测线的数据点加上对应的校正量，完成校正。
以一个示例性的实例进一步解释本步骤的处理过程，假如第一组测线有3条测线，第二组测线有3条测线，第一组测线和第二组测线的校正过程如下：
步骤31：第二组测线中的每一根测线都会与第一组测线中的测线相交，只有交点处有闭合差，也即只有交点处有闭合差的校正量，测线上的其余点只需要使用反距离加权或者线性插值等插值方法，计算出对应的校正量，每个点加上对应的校正量，完成一根测线上所有点的校正，进而完成这根测线与第一组测线的闭合差；
步骤32：重复步骤31，将第二组测线中剩余的测线完成闭合差校正，使得第二组中的任何一根测线与第一组中的任何一根测线均不存在闭合差，完成第一组测线和第二组测线的校正。
步骤33：将完成校正的第一组测线和第二组测线合并为一组测线，并作为当前测线组，重复步骤31-33，完成所有测线组的校正。
如图5和图6所示，图5为未消除闭合差直接成图的地震剖面示意图，图6为按本发明对与图5相同的原始二维测线地震数据消除闭合差的地震剖面示意图。从图6可以看出，本发明能够有效进行闭合差校正，消除闭合差，消除闭合差的效率和精度高。
如图7所示，本发明还涉及二维测线地震数据闭合差的校正方法的实体实现处理终端100，其包括，
存储器101，用于存储程序指令；
处理器102，用于运行所述程序指令，以执行所述地震构造图的建立方法中的步骤。
本说明书所公开的实施例只是对本发明单方面特征的一个例证，本发明的保护范围不限于此实施例，其他任何功能等效的实施例均落入本发明的保护范围内。对于本领域的技术人员来说，可根据以上描述的技术方案以及构思，做出其它各种相应的改变以及变形，而所有的这些改变以及变形都应该属于本发明权利要求的保护范围之内。
Claims (4)
1.一种二维测线地震数据闭合差的校正方法，其特征在于，包括如下步骤：
步骤1：获得二维测线地震数据，并对二维测线按施工特征进行排序，相同施工特征的二维测线的排序相同，施工特征是指进行二维测线施工作业时的施工特征；
步骤2：根据排序结果，对所有测线进行分组，以使得同一组的二维测线不相交，不同组之间的二维测线相交，从而得到若干组测线；
步骤3：从第一组测线开始，计算当前组测线与下一组测线之间的校正量，根据校正量对当前组测线与下一组测线之间进行校正，校正完后将两组测线进行合并为同一组测线，合并后的同一组测线作为当前组测线继续与下一组测线之间进行校正，依次进行，直至当前组测线与最后一组测线完成校正，从而消除闭合差，
所述步骤2中，对所有测线进行分组的具体过程包括如下步骤：
步骤21：从所有测线中选取排序最靠前的二维测线作为第一根测线，从第一根测线开始，选取第i根测线作为未分组测线，i＝1,2,3,……，n，n为二维测线总数，比较第q根测线和第i根测线，q＝1,2,3,…,n且q≠i，若第q根测线和第i根测线不相交，则将第q根测线标记为与第i根测线为同一组测线，否则，将第q根测线标记为与第i根测线为相异组测线，依次进行，遍历所有测线，找出所有与第i根测线不相交的测线，从而得到与第i根测线均不相交的所有测线并作为第一组测线，剩下的测线均为与第i根测线均相交的相异组测线；
步骤22：从与第i根测线相交的相异组测线中选取排序最靠前的二维测线作为第一根测线，重复执行步骤21,得到第二组测线，依次进行，直至将所有测线均分组完成，分别得到第一组测线、第二组测线、……、第m组测线。
2.根据权利要求1所述的二维测线地震数据闭合差的校正方法，其特征在于，所述施工特征为施工时间、施工时的测线角度、施工批次中的任一个。
3.根据权利要求1所述的二维测线地震数据闭合差的校正方法，其特征在于，所述根据校正量对当前组测线与下一组测线之间进行校正，其具体实现过程包括如下步骤：根据当前组测线与下一组测线的相交关系得到交点处的校正量，将所述校正量线性插值到每根测线的所有数据点上，然后将每根测线的数据点加上对应的校正量，完成校正。
4.一种处理终端，其特征在于，其包括，
存储器，用于存储程序指令；
处理器，用于运行所述程序指令，以执行如权利要求1-3任一项所述的二维测线地震数据闭合差的校正方法中的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911014861.4A CN110888167B (zh) | 2019-10-23 | 2019-10-23 | 一种二维测线地震数据闭合差的校正方法及处理终端 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911014861.4A CN110888167B (zh) | 2019-10-23 | 2019-10-23 | 一种二维测线地震数据闭合差的校正方法及处理终端 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110888167A CN110888167A (zh) | 2020-03-17 |
| CN110888167B true CN110888167B (zh) | 2021-07-30 |
Family
ID=69746437
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911014861.4A Active CN110888167B (zh) | 2019-10-23 | 2019-10-23 | 一种二维测线地震数据闭合差的校正方法及处理终端 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110888167B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113514880B (zh) * | 2020-04-10 | 2024-02-20 | 中国石油化工股份有限公司 | 一种二维测线交点处叠加速度构建方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101750628A (zh) * | 2008-12-11 | 2010-06-23 | 中国石油天然气股份有限公司 | 二维叠加速度及均方根速度场闭合差校正方法 |
| CN102323617A (zh) * | 2011-06-13 | 2012-01-18 | 中国石油化工股份有限公司 | 一种复杂地表的二维地震资料连片处理方法 |
| CN102353993A (zh) * | 2011-06-20 | 2012-02-15 | 中国石油集团川庆钻探工程有限公司 | 二维工区速度分析时速度场闭合差的实时检查方法 |
| CN103869356A (zh) * | 2012-12-10 | 2014-06-18 | 中国石油天然气集团公司 | 一种地震反射界面闭合方法 |
| CN108828669A (zh) * | 2018-04-02 | 2018-11-16 | 中国石油天然气集团有限公司 | 一种二维相交测线静校正处理方法、装置及系统 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6757615B2 (en) * | 2001-09-14 | 2004-06-29 | Exxonmobil Upstream Research Company | Method for automated horizon transfer and alignment through the application of time-shift volumes |
| US8594939B2 (en) * | 2009-01-16 | 2013-11-26 | Fairfield Industries, Incorporated | Method for calibrating seismic imaging velocities |
-
2019
- 2019-10-23 CN CN201911014861.4A patent/CN110888167B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101750628A (zh) * | 2008-12-11 | 2010-06-23 | 中国石油天然气股份有限公司 | 二维叠加速度及均方根速度场闭合差校正方法 |
| CN102323617A (zh) * | 2011-06-13 | 2012-01-18 | 中国石油化工股份有限公司 | 一种复杂地表的二维地震资料连片处理方法 |
| CN102353993A (zh) * | 2011-06-20 | 2012-02-15 | 中国石油集团川庆钻探工程有限公司 | 二维工区速度分析时速度场闭合差的实时检查方法 |
| CN103869356A (zh) * | 2012-12-10 | 2014-06-18 | 中国石油天然气集团公司 | 一种地震反射界面闭合方法 |
| CN108828669A (zh) * | 2018-04-02 | 2018-11-16 | 中国石油天然气集团有限公司 | 一种二维相交测线静校正处理方法、装置及系统 |
Non-Patent Citations (7)
| Title |
|---|
| Correcting amplitude, time, and phase mis-ties in seismic data;Thomas N. Bishop et al.;《GEOPHYSICS》;19940630;第59卷(第6期);第946-953页 * |
| Michael D. Harper et al..Seismic mis-tie resolution technique.《GEOPHYSICS》.1991,第56卷(第11期), * |
| Seismic mis-tie resolution technique;Michael D. Harper et al.;《GEOPHYSICS》;19911130;第56卷(第11期);第1829-1830页 * |
| 二维地震数据闭合方法新探;沈克非等;《石油地球物理勘探》;20081231;第43卷;第155-157页 * |
| 交互闭合差校正方法;苏永昌;《石油地球物理勘探》;19941231;第29卷(第6期);第793-798页 * |
| 地震数据闭合差精细校正三步法;肖春燕;《1996年中国地球物理学会第十二届学术年会论文集》;19960930;第48页 * |
| 基于二维地震数据的宏观三维速度建模—以南海西北部海域为例;赵斌等;《海洋地质前沿》;20190731;第35卷(第7期);第18-24页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110888167A (zh) | 2020-03-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105787281B (zh) | 一种海浪有效波高场的融合反演方法和装置 | |
| CN109710893B (zh) | 一种用于修正矿山边坡变形监测异常数据的时空插值方法 | |
| CN110333543B (zh) | 基于反射系数分析的低阻体解释及成像方法与系统 | |
| CN104374317A (zh) | 一种基于激光跟踪仪多点测量技术的机床误差标定方法 | |
| CN110888167B (zh) | 一种二维测线地震数据闭合差的校正方法及处理终端 | |
| CN107831516B (zh) | 融合gnss和地面监测网获取大坝实时高精度位移的方法 | |
| CN109782355B (zh) | Obs检波点漂移的检测方法及装置 | |
| Manning et al. | 14C record and wiggle-match placement for the Anatolian (Gordion area) juniper tree-ring chronology~ 1729 to 751 cal BC, and typical Aegean/Anatolian (growing season related) regional 14C offset assessment | |
| CN108802860B (zh) | 数据修正方法、数据修正装置 | |
| CN110376644B (zh) | 叠前地震数据拼接方法及装置 | |
| CN116029036B (zh) | 用于运营普速铁路的平面坐标测量数据搭接方法和系统 | |
| CN110261820B (zh) | 一种多测量站的时差定位方法及装置 | |
| Adamczak et al. | An analysis of strategies of measurement of 3D rotary elements | |
| CN110567489B (zh) | 角度交会测量系统的动态误差的获取方法及系统 | |
| CN108387194A (zh) | 一种井下闭合导线坐标测量方法 | |
| CN105140164B (zh) | 一种用于承载芯片的定位平台的旋转中心标定方法 | |
| CN113673759A (zh) | 一种水文资料实时整编方法及终端 | |
| CN109116391B (zh) | 一种基于改进正交分解的区域划分方法 | |
| CN107451992A (zh) | 一种sar图像变化检测的方法与装置 | |
| CN111624655A (zh) | 初至波剩余静校正量的确定方法及装置 | |
| CN110968656B (zh) | 地形图处理方法及处理装置 | |
| CN110096674B (zh) | 一种基于参数依赖的仪器参数校正方法 | |
| CN111629321A (zh) | 基于离群点剥离的中值卡尔曼滤波的定位方法 | |
| CN109344545A (zh) | 建筑设计软件中构件尺寸的精确度调整方法、装置 | |
| CN116819573B (zh) | 一种削弱gps相关lod系统误差的长期模型构建方法 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |