CN112394413A - 三维初至波剩余静校正方法及装置 - Google Patents

三维初至波剩余静校正方法及装置 Download PDF

Info

Publication number
CN112394413A
CN112394413A CN202011202638.5A CN202011202638A CN112394413A CN 112394413 A CN112394413 A CN 112394413A CN 202011202638 A CN202011202638 A CN 202011202638A CN 112394413 A CN112394413 A CN 112394413A
Authority
CN
China
Prior art keywords
correction value
static correction
refraction
shot point
slowness
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
CN202011202638.5A
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 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 CN202011202638.5A priority Critical patent/CN112394413A/zh
Publication of CN112394413A publication Critical patent/CN112394413A/zh
Pending legal-status Critical Current

Links

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. analysis, for interpretation, for correction
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362—Effecting static or dynamic corrections; Stacking
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum

Abstract

本发明提供了一种三维初至波剩余静校正方法及装置，该方法包括：利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。通过建立考虑低速层厚度的线性方程组，将低速层的厚度因素考虑在内，通过利用最小二乘法进行三维折射初至波的旅行时拟合，将不同炮之间初至曲线的变化的因素考虑在内，从而提高了校正的精度；通过对炮点静校正量和检波点静校正量进行去低频处理，从而提高了剩余静校正量的计算精度。

Description

三维初至波剩余静校正方法及装置
技术领域
本发明涉及地震资料处理技术领域，尤其涉及一种三维初至波剩余静校正方法及装置。
背景技术
近些年进行地震资料处理时，发展起来一类静校正方法——初至折射剩余静校正方法，这类方法不需要确定近地表层的厚度和速度，而是用直线或某种曲线对实际初至曲线进行拟合，然后再将实际旅行时曲线与拟合的趋势曲线之间的剩余差时差分解成炮点和检波点静校正量。具体地，如拟合迭代折射静校正、全差分折射静校正、折射初至对齐静校正等。这类方法能较好的改正单炮初至曲线的畸变，使叠加剖面上的同相轴能很好的聚焦，计算速度快、精度高，能较好地解决短波长静校正问题。但是这类方法没有考虑不同炮之间初至曲线的变化和低速层的厚度等因素，因而导致校正的精度不足。
发明内容
本发明实施例提供一种三维初至波剩余静校正方法，用以提高校正精度，该方法包括：
利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，所述线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；
将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；
对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。
本发明实施例还提供一种三维初至波剩余静校正装置，用以提高校正精度，该装置包括：
旅行时拟合模块，用于利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，所述线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；
迭代求解模块，用于将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；
去低频模块，用于对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述三维初至波剩余静校正方法。
本发明实施例也提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述三维初至波剩余静校正方法的计算机程序。
本发明实施例中，通过利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，该线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。通过建立考虑低速层厚度的线性方程组，将低速层的厚度因素考虑在内，通过利用最小二乘法进行三维折射初至波的旅行时拟合，将不同炮之间初至曲线的变化的因素考虑在内，从而提高了校正的精度；为了避免拾取折射初至时间的误差、炮检距范围的限制以及坏道剔除等因素造成炮点和检波点的有效道数小于4时带来的计算误差，通过对炮点静校正量和检波点静校正量进行去低频处理，从而提高了剩余静校正量的计算精度。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中三维初至波剩余静校正方法的示意图。
图2为本发明具体实施例中步骤101的实现方法流程示意图。
图3为本发明具体实施例中二层水平介质折射模型的示意图。
图4为本发明具体实施例中步骤102的实现方法流程示意图。
图5为本发明一具体应用实施中某三维单炮线性动校正初至图。
图6为本发明一具体应用实施中某三维单炮线性动校正叠加剖面示意图。
图7为本发明一具体应用实施中进行三维初至波剩余静校正后的某三维单炮线性动校正初至图。
图8为本发明一具体应用实施中进行三维初至波剩余静校正后的某三维单炮线性动校正叠加剖面示意图。
图9为本发明实施例中三维初至波剩余静校正装置的示意图。
图10为本发明具体实施例中迭代求解模块902的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例提供了一种三维初至波剩余静校正方法，用以提高校正精度，如图1所示，该方法包括：
步骤101：利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；
步骤102：将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；
步骤103：对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。
由图1所示流程可以得知，本发明实施例中，通过利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，该线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。通过建立考虑低速层厚度的线性方程组，将低速层的厚度因素考虑在内，通过利用最小二乘法进行三维折射初至波的旅行时拟合，将不同炮之间初至曲线的变化的因素考虑在内，从而提高了校正的精度；为了避免拾取折射初至时间的误差、炮检距范围的限制以及坏道剔除等因素造成炮点和检波点的有效道数小于4时带来的计算误差，通过对炮点静校正量和检波点静校正量进行去低频处理，从而提高了剩余静校正量的计算精度。
具体实施时，首先利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组。具体实施过程，如图2所示，包括：
步骤201：假设折射面水平，折射层为均匀介质，构建表征低速层厚度和三维折射初至波的旅行时之间对应关系的三维折射初至波的时距方程；
步骤202：应用二层水平介质折射模型，建立炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系；
步骤203：根据时距方程、炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系，建立三维折射初至波的旅行时与炮点静校正量、检波点静校正量和折射层慢度之间的关联关系；
步骤204：对关联关系应用最小二乘法，建立以炮检距和三维折射初至时间为已知量，以炮点静校正量、检波点静校正量和综合折射层慢度为未知量的线性方程组。
具体实施例中，假设折射面水平，折射层为均匀介质，三维折射初至波的时距方程为：
Figure BDA0002755871040000041
其中，tSR代表炮点S激发、检波点R接收的折射旅行时；hS和hR分别代表S和R处的低速层厚度；X代表S和R之间的距离；V1和V2分别代表低速层和下覆折射层的速度；α代表临界角。其中，低速层是指二层水平介质折射模型中速度相对较低的第一层；下覆折射层是指二层水平介质折射模型中速度相对较高，能产生折射波的第二层。
式(1)右端前两项之和为该直线方程的截距时，一般情况下，当V2>>V1时，cosα≈1，这时应用二层水平介质折射模型，如图3所示，式(1)右端前两项分别为炮点和检波点的延迟时间，记为：
Figure BDA0002755871040000051
其中，ST、RT分别代表炮点静校正量和检波点静校正量。炮点和检波点之间的折射层慢度P记为：
Figure BDA0002755871040000052
进行替换可以得到：tSR＝ST+RT+PX (3)，对三维排列的每一道，式(3)均成立。
其中，炮点静校正量ST、检波点静校正量RT和折射层慢度P为未知量，应用最小二乘法可建立线性方程组。
设第i炮，第j个检波点拾取的折射初至时间为tij，炮检距为Xij，炮点静校正量为STi，检波点静校正量为RTj，在给定炮检距范围内(来自同一折射层面)，该炮点附近的综合折射层慢度为Pi，应用最小二乘法使下面的目标函数达到极小：
Figure BDA0002755871040000053
由
Figure BDA0002755871040000054
可得正则方程组：
Figure BDA0002755871040000055
其中，
Figure BDA0002755871040000056
代表在给定的炮检距范围内，第n个炮点所包括的全部检波点之和，个数为Nn；
Figure BDA0002755871040000057
代表在给定的炮检距范围内，第m个检波点所包括的全部炮点之和，个数为Mm；STn代表第n个炮点的炮点静校正量；STj表示第j个炮点的炮点静校正量；RTj代表第j个检波点的检波点静校正量；RTm代表第m个检波点的检波点静校正量；Xnj代表第n个炮点与第j个检波点之间的炮检距；Xim代表第i个炮点与第m个检波点之间的炮检距；Pn代表第n个炮点附近的综合折射层慢度；Pi代表第i个炮点附近的综合折射层慢度；tnj代表第n炮、第j个检波点拾取的折射初至时间；tim代表第i炮、第m个检波点拾取的折射初至时间。
建立考虑低速层厚度的线性方程组后，将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量。具体实施过程，如图4所示，包括：
步骤401：将迭代求解方程组进行变换，形成用于求解综合折射层慢度的慢度求解迭代公式；
步骤402：对综合折射层慢度进行去高频分量处理，结合慢度求解迭代公式，确定三维观测平面上任意炮点处的综合折射层慢度值；
步骤403：根据三维观测平面上任意炮点处的综合折射层慢度值和迭代求解方程组，进行迭代求解，得到炮点静校正量和检波点静校正量。
具体实施例中，将折射层慢度Pi作为迭代求解的一个未知量，是因为利用折射波的目的，就是为了解决大剩余静校正量的问题。当静校正量很大时(尤其是检波点静校正量)，如果将检波点静校正量RTj假设为0而一次性估算折射层慢度Pi而加以应用，势必会使Pi估算误差引入到静校正量迭代计算中。但是，如果直接对线性方程组应用高斯-赛德尔迭代法求解，又保证不了该迭代法解的收敛性。为此，将线性方程组改写成下列迭代形式：
Figure BDA0002755871040000061
Figure BDA0002755871040000062
Figure BDA0002755871040000063
Figure BDA0002755871040000064
式中上角标(l)表示未知量的第L次迭代解，l＝1，2，……，L，L为给定的迭代次数且
Figure BDA0002755871040000065
利用迭代方程组形成慢度求解迭代公式，用式(5)和(6)联立求解Pn：
Figure BDA0002755871040000071
对综合折射层慢度进行去高频分量处理：
式(9)中新估算的折射层慢度Pn不是局部的，而是炮点附近一个区域内综合的折射层慢度。虽然在推导过程中假设折射层的速度(或慢度)横向是不变的，为了使本方法适应实际情况，既使折射层速度横向有一些变化，其综合折射层慢度横向也是渐变的。因此，求得全部的综合折射层慢度值后，对其去高频分量处理，即只取Pn的线性分量。
已知第i炮的坐标和求得的慢度分别为(xi，yi)和Pi，i＝1,2,…,I，I为总炮点数。设三维区域内的折射层慢度的线性函数为
P(x,y)＝P0+αx+βy (10)
其中，P0、α和β皆为常数。
下面用最小二乘拟合原理来建立三维折射层慢度函数。定义目标函数为：
Figure BDA0002755871040000072
由
Figure BDA0002755871040000073
得：
Figure BDA0002755871040000074
Figure BDA0002755871040000075
Figure BDA0002755871040000081
联立式(9)(10)(11)解方程组得：
Figure BDA0002755871040000082
Figure BDA0002755871040000083
Figure BDA0002755871040000084
计算出P0、α和β参数后，就可用式(6)的折射慢度函数，从而可求出三维观测平面上任意炮点处的综合折射层慢度值Pn (l)。
根据任意炮点处的综合折射层慢度值Pn (l)，代入式(7)和(8)求取STi和RTj新的迭代结果，作为炮点静校正量和检波点静校正量。
对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。由于拾取折射初至时间的误差、炮检距范围的限制以及坏道剔除等因素造成炮点和检波点的有效道数小于4时，为了提高静校正量计算精度，这些炮点静校正量和检波点静校正量不再有效，而是按其实际地面位置，根据相邻有效的炮点静校正量和检波点静校正量的插值所取代。
由拾取折射初至时间分解炮点静校正量和检波点静校正量，其直流分量不能唯一确定。当在炮点静校正量STi加上一个ΔT，检波点静校正量RTj上减去一个ΔT时，其折射初至时间是不变的。通过折射初至静校正方法求得的静校正量，最终直接影响反射层的构造形态。如果炮点静校正量有一个线性分量ST(x,y)＝αX+βY,检波点静校正量的线性分量RT(x,y)＝α'X+β'Y,若α'＝-α,β'＝-β，则这些分量在其中心点叠加过程中不起作用；否则会引起地下反射产状的畸变。为了使求出的静校正量对迭加剖面构造产状不产生影响，因此在最终输出的炮点、检波点静校正量中要做去低频分量处理。
已知坐标xi，yi处的炮点静校正量为STi，i＝1,2,…,n，n为该范围内的炮点数。炮点静校正量可用线性函数表示为
ST(x,y)＝ST0+αx+βy (19)
其中，ST0、α和β皆为常数。根据最小二乘原理，求炮点静校正量STi，使下列目标函数达到极小：
Figure BDA0002755871040000091
令
Figure BDA0002755871040000092
即可计算出ST0、α和β参数后，然后就可以用式(19)得到炮点静校正量线性分量。用求得的炮点静校正量STi减去平面拟合的线性分量，差就是去除平面低频分量的最终的炮点剩余静校正量。同理，把迭代求出的检波点静校正量RTj减去平面拟合的线性分量检波点静校正量作为最终的检波点剩余静校正量。
下面给出一具体实例说明本发明实施例如何进行三维初至波剩余静校正。本例应用于某三维单炮初至波剩余静校正，包括如下步骤：
1)建立线性方程组：
假设折射面水平，折射层为均匀介质，可应用最小二乘法建立炮点静校正量ST、检波点静校正量RT和折射层慢度P的线性方程组。
2)将线性方程组改写为迭代求解方程组。
3)利用迭代方程组形成慢度求解迭代公式。
4)对综合折射层慢度进行去高频分量处理。
5)求取炮点静校正量ST和检波点静校正量RT。
6)对炮点静校正量ST、检波点静校正量RT进行去低频处理，求得最终的炮点剩余静校正量、检波点剩余静校正量。
图5和图6分别为某三维单炮线性动校正初至图和某三维单炮线性动校正叠加剖面。由于静校正的影响，叠加剖面明显变差，折射初至跳动严重。图7和图8分别为应用本方法后的某三维单炮线性动校正初至图和某三维单炮线性动校正叠加剖面。线性校正的规律性明显增强，叠加剖面明显改善。由此可见，使用本发明实施例所提供的三维初至波剩余静校正方法后，有效提高了复杂地表区的静校正精度。且由于现有技术仅能在二维复杂地表资料处理中有较好的效果，对三维方法仍然处于试验探索阶段。而本发明实施例所提供的三维初至波剩余静校正方法，较好地解决了三维分解中的一些关键技术问题，取得了较好的静校正效果。
上述具体应用的实施仅为举例，其余实施方式不再一一赘述。
基于同一发明构思，本发明实施例还提供一种三维初至波剩余静校正装置，由于三维初至波剩余静校正装置所解决问题的原理与三维初至波剩余静校正方法相似，因此三维初至波剩余静校正装置的实施可以参见三维初至波剩余静校正方法的实施，重复之处不再赘述，具体结构如图9所示：
旅行时拟合模块901，用于利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；
迭代求解模块902，用于将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；
去低频模块903，用于对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。
具体实施例中，旅行时拟合模块901具体用于：
假设折射面水平，折射层为均匀介质，构建表征低速层厚度和三维折射初至波的旅行时之间对应关系的三维折射初至波的时距方程；
应用二层水平介质折射模型，建立炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系；
根据时距方程、炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系，建立三维折射初至波的旅行时与炮点静校正量、检波点静校正量和折射层慢度之间的关联关系；
对关联关系应用最小二乘法，建立以炮检距和三维折射初至时间为已知量，以炮点静校正量、检波点静校正量和综合折射层慢度为未知量的线性方程组。
具体实施例中，线性方程组为：
Figure BDA0002755871040000111
其中，
Figure BDA0002755871040000112
代表在给定的炮检距范围内，第n个炮点所包括的全部检波点之和，个数为Nn；
Figure BDA0002755871040000113
代表在给定的炮检距范围内，第m个检波点所包括的全部炮点之和，个数为Mm；STn代表第n个炮点的炮点静校正量；STj表示第j个炮点的炮点静校正量；RTj代表第j个检波点的检波点静校正量；RTm代表第m个检波点的检波点静校正量；Xnj代表第n个炮点与第j个检波点之间的炮检距；Xim代表第i个炮点与第m个检波点之间的炮检距；Pn代表第n个炮点附近的综合折射层慢度；Pi代表第i个炮点附近的综合折射层慢度；tnj代表第n炮、第j个检波点拾取的折射初至时间；tim代表第i炮、第m个检波点拾取的折射初至时间。
具体实施例中，迭代求解模块902的结构，如图10所示，包括：
慢度求解迭代单元1001，用于将迭代求解方程组进行变换，形成用于求解综合折射层慢度的慢度求解迭代公式；
慢度值确定单元1002，用于对综合折射层慢度进行去高频分量处理，结合慢度求解迭代公式，确定三维观测平面上任意炮点处的综合折射层慢度值；
迭代求解单元1003，用于根据三维观测平面上任意炮点处的综合折射层慢度值和迭代求解方程组，进行迭代求解，得到炮点静校正量和检波点静校正量。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述三维初至波剩余静校正方法。
本发明实施例还提供一种计算机可读存储介质，存储有执行上述三维初至波剩余静校正方法的计算机程序。
综上所述，本发明实施例提供的三维初至波剩余静校正方法及装置具有如下优点：
通过利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，该线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。通过建立考虑低速层厚度的线性方程组，将低速层的厚度因素考虑在内，通过利用最小二乘法进行三维折射初至波的旅行时拟合，将不同炮之间初至曲线的变化的因素考虑在内，从而提高了校正的精度；为了避免拾取折射初至时间的误差、炮检距范围的限制以及坏道剔除等因素造成炮点和检波点的有效道数小于4时带来的计算误差，通过对炮点静校正量和检波点静校正量进行去低频处理，从而提高了剩余静校正量的计算精度。
本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、装置和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种三维初至波剩余静校正方法，其特征在于，包括：
利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，所述线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；
将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；
对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。
2.如权利要求1所述的方法，其特征在于，利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组，包括：
假设折射面水平，折射层为均匀介质，构建表征低速层厚度和三维折射初至波的旅行时之间对应关系的三维折射初至波的时距方程；
应用二层水平介质折射模型，建立炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系；
根据所述时距方程、炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系，建立三维折射初至波的旅行时与炮点静校正量、检波点静校正量和折射层慢度之间的关联关系；
对所述关联关系应用最小二乘法，建立以炮检距和三维折射初至时间为已知量，以炮点静校正量、检波点静校正量和综合折射层慢度为未知量的线性方程组。
3.如权利要求1所述的方法，其特征在于，所述线性方程组为：
Figure FDA0002755871030000011
其中，
Figure FDA0002755871030000012
代表在给定的炮检距范围内，第n个炮点所包括的全部检波点之和，个数为Nn；
Figure FDA0002755871030000013
代表在给定的炮检距范围内，第m个检波点所包括的全部炮点之和，个数为Mm；STn代表第n个炮点的炮点静校正量；STj表示第j个炮点的炮点静校正量；RTj代表第j个检波点的检波点静校正量；RTm代表第m个检波点的检波点静校正量；Xnj代表第n个炮点与第j个检波点之间的炮检距；Xim代表第i个炮点与第m个检波点之间的炮检距；Pn代表第n个炮点附近的综合折射层慢度；Pi代表第i个炮点附近的综合折射层慢度；tnj代表第n炮、第j个检波点拾取的折射初至时间；tim代表第i炮、第m个检波点拾取的折射初至时间。
4.如权利要求1所述的方法，其特征在于，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量，包括：
将迭代求解方程组进行变换，形成用于求解综合折射层慢度的慢度求解迭代公式；
对综合折射层慢度进行去高频分量处理，结合慢度求解迭代公式，确定三维观测平面上任意炮点处的综合折射层慢度值；
根据三维观测平面上任意炮点处的综合折射层慢度值和迭代求解方程组，进行迭代求解，得到炮点静校正量和检波点静校正量。
5.一种三维初至波剩余静校正装置，其特征在于，包括：
旅行时拟合模块，用于利用最小二乘法进行三维折射初至波的旅行时拟合，建立考虑低速层厚度的线性方程组；其中，所述线性方程组以炮点静校正量、检波点静校正量和综合折射层慢度为未知量，炮检距和三维折射初至时间为已知量；
迭代求解模块，用于将线性方程组改写为求取综合折射层慢度的迭代求解方程组，对迭代求解方程组进行迭代求解，得到炮点静校正量和检波点静校正量；
去低频模块，用于对炮点静校正量和检波点静校正量进行去低频处理，确定炮点剩余静校正量和检波点剩余静校正量。
6.如权利要求5所述的装置，其特征在于，所述旅行时拟合模块具体用于：
假设折射面水平，折射层为均匀介质，构建表征低速层厚度和三维折射初至波的旅行时之间对应关系的三维折射初至波的时距方程；
应用二层水平介质折射模型，建立炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系；
根据所述时距方程、炮点静校正量与低速层厚度之间的对应关系和检波点静校正量与低速层厚度之间的对应关系，建立三维折射初至波的旅行时与炮点静校正量、检波点静校正量和折射层慢度之间的关联关系；
对所述关联关系应用最小二乘法，建立以炮检距和三维折射初至时间为已知量，以炮点静校正量、检波点静校正量和综合折射层慢度为未知量的线性方程组。
7.如权利要求5所述的装置，其特征在于，所述线性方程组为：
Figure FDA0002755871030000031
其中，
Figure FDA0002755871030000032
代表在给定的炮检距范围内，第n个炮点所包括的全部检波点之和，个数为Nn；
Figure FDA0002755871030000033
代表在给定的炮检距范围内，第m个检波点所包括的全部炮点之和，个数为Mm；STn代表第n个炮点的炮点静校正量；STj表示第j个炮点的炮点静校正量；RTj代表第j个检波点的检波点静校正量；RTm代表第m个检波点的检波点静校正量；Xnj代表第n个炮点与第j个检波点之间的炮检距；Xim代表第i个炮点与第m个检波点之间的炮检距；Pn代表第n个炮点附近的综合折射层慢度；Pi代表第i个炮点附近的综合折射层慢度；tnj代表第n炮、第j个检波点拾取的折射初至时间；tim代表第i炮、第m个检波点拾取的折射初至时间。
8.如权利要求5所述的装置，其特征在于，所述迭代求解模块包括：
慢度求解迭代单元，用于将迭代求解方程组进行变换，形成用于求解综合折射层慢度的慢度求解迭代公式；
慢度值确定单元，用于对综合折射层慢度进行去高频分量处理，结合慢度求解迭代公式，确定三维观测平面上任意炮点处的综合折射层慢度值；
迭代求解单元，用于根据三维观测平面上任意炮点处的综合折射层慢度值和迭代求解方程组，进行迭代求解，得到炮点静校正量和检波点静校正量。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一所述方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一所述方法的计算机程序。
CN202011202638.5A 2020-11-02 2020-11-02 三维初至波剩余静校正方法及装置 Pending CN112394413A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011202638.5A CN112394413A (zh) 2020-11-02 2020-11-02 三维初至波剩余静校正方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011202638.5A CN112394413A (zh) 2020-11-02 2020-11-02 三维初至波剩余静校正方法及装置

Publications (1)

Publication Number Publication Date
CN112394413A true CN112394413A (zh) 2021-02-23

Family

ID=74597688

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011202638.5A Pending CN112394413A (zh) 2020-11-02 2020-11-02 三维初至波剩余静校正方法及装置

Country Status (1)

Country Link
CN (1) CN112394413A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113945980A (zh) * 2021-09-30 2022-01-18 国能神东煤炭集团有限责任公司 一种初至时间校正方法及系统

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4802146A (en) * 1985-05-23 1989-01-31 Mobil Oil Corporation Method for moveout correction and stacking velocity estimation of offset VSP data
US20020075759A1 (en) * 2001-01-21 2002-06-20 Xixiang Zhou Static correction method for exploration seismic data using first arrivals of seismic waves
US6424920B1 (en) * 1999-09-17 2002-07-23 Konstantin Sergeevich Osypov Differential delay-time refraction tomography
CN101980054A (zh) * 2010-10-22 2011-02-23 中国石油化工股份有限公司 一种在高密度地震静校正处理中建立近地表速度模型的方法
CN102053275A (zh) * 2009-10-30 2011-05-11 中国石油化工股份有限公司 一种用于单点地震室内组合的相对静校正量计算方法
CN102692651A (zh) * 2011-03-24 2012-09-26 中国石油集团东方地球物理勘探有限责任公司 速度空变的初至波剩余静校正方法
US20140219054A1 (en) * 2013-02-04 2014-08-07 King Fahd University Of Petroleum And Minerals Method of first arrival picking of seismic refraction data
CN104570122A (zh) * 2015-01-14 2015-04-29 中国石油集团东方地球物理勘探有限责任公司 基于反射波的地震数据静校正方法及装置
WO2017106127A1 (en) * 2015-12-18 2017-06-22 Saudi Arabian Oil Company Automated near surface analysis by surface-consistent refraction methods
CN107219554A (zh) * 2017-05-05 2017-09-29 北京大学 陆地地震资料的剩余静校正量的自动获取方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4802146A (en) * 1985-05-23 1989-01-31 Mobil Oil Corporation Method for moveout correction and stacking velocity estimation of offset VSP data
US6424920B1 (en) * 1999-09-17 2002-07-23 Konstantin Sergeevich Osypov Differential delay-time refraction tomography
US20020075759A1 (en) * 2001-01-21 2002-06-20 Xixiang Zhou Static correction method for exploration seismic data using first arrivals of seismic waves
CN102053275A (zh) * 2009-10-30 2011-05-11 中国石油化工股份有限公司 一种用于单点地震室内组合的相对静校正量计算方法
CN101980054A (zh) * 2010-10-22 2011-02-23 中国石油化工股份有限公司 一种在高密度地震静校正处理中建立近地表速度模型的方法
CN102692651A (zh) * 2011-03-24 2012-09-26 中国石油集团东方地球物理勘探有限责任公司 速度空变的初至波剩余静校正方法
US20140219054A1 (en) * 2013-02-04 2014-08-07 King Fahd University Of Petroleum And Minerals Method of first arrival picking of seismic refraction data
CN104570122A (zh) * 2015-01-14 2015-04-29 中国石油集团东方地球物理勘探有限责任公司 基于反射波的地震数据静校正方法及装置
WO2017106127A1 (en) * 2015-12-18 2017-06-22 Saudi Arabian Oil Company Automated near surface analysis by surface-consistent refraction methods
CN107219554A (zh) * 2017-05-05 2017-09-29 北京大学 陆地地震资料的剩余静校正量的自动获取方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
崔庆辉 等: "《三维折射初至延迟时分解静校正》", 《重庆科技学院学报(自然科学版)》, vol. 11, no. 4, 31 August 2009 (2009-08-31), pages 44 - 46 *
张永峰 等: "《基于变差函数的高精度静校正融合技术及其应用》", 《岩性油气藏》, vol. 27, no. 3, 30 June 2015 (2015-06-30), pages 108 - 114 *
方勇 等: "《库车山地地震资料层析静校正方法的应用》", 《石油地球物理勘探》, vol. 52, no. 1, 30 November 2017 (2017-11-30), pages 23 - 30 *
杨华臣 等: "《地震静校正的相对时延法》", 《中国海洋大学学报(自然科学版)》, vol. 7, no. 3, 31 March 2017 (2017-03-31), pages 87 - 94 *
王克斌 等: "《折射静校正在苏里格气田三维处理中的应用》", 《石油物探》, vol. 42, no. 2, 30 June 2003 (2003-06-30), pages 248 - 251 *
穆蜀生 等: "《GeoEast静校正处理技术》", 《石油工业计算机应用》, vol. 24, no. 4, 31 December 2016 (2016-12-31), pages 25 - 35 *
苏世龙 等: "《层析静校正方法及其在东部勘探中的应用》", 《油气地球物理》, vol. 8, no. 1, 31 January 2010 (2010-01-31), pages 11 - 16 *
车建英 等: "《初至折射波静校正技术在复杂探区地震数据处理中的应用》", 《中国煤田地质》, vol. 18, no. 1, 28 February 2006 (2006-02-28), pages 54 - 56 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113945980A (zh) * 2021-09-30 2022-01-18 国能神东煤炭集团有限责任公司 一种初至时间校正方法及系统
CN113945980B (zh) * 2021-09-30 2024-02-09 国能神东煤炭集团有限责任公司 一种初至时间校正方法及系统

Similar Documents

Publication Publication Date Title
CN109917454B (zh) 基于双基准面的真地表叠前深度偏移成像方法及装置
CN109239781B (zh) 一种地震数据校正方法及装置
CN113552625B (zh) 一种用于常规陆域地震数据的多尺度全波形反演方法
CN108196305B (zh) 一种山地静校正方法
CN111638551A (zh) 地震初至波走时层析方法及装置
CN112394413A (zh) 三维初至波剩余静校正方法及装置
CN111257935B (zh) 一种加快层析反演速度收敛的速度融合方法及处理终端
CN109581494B (zh) 叠前偏移方法及装置
CN105572741A (zh) 一种计算3d高频静校正量的方法
CN109145375B (zh) 裂缝渗透率模型校正方法及装置
CN112465983A (zh) 一种高程地表模型点线状障碍物修正方法和装置
CN111435172A (zh) 层析静校正的方法和装置
CN105549075B (zh) 求取浅层速度场的方法及装置
CN112230279B (zh) 增强纵波地震数据品质的方法及装置
CN111999767B (zh) 起伏地表的偏移成像方法及装置
CN112147691B (zh) 快速编码免排序基准面校正方法及系统
CN112241021B (zh) 初至拾取方法及装置
CN110967752A (zh) 一种起伏地表偏移输入地震数据的静校正方法
CN111624655B (zh) 初至波剩余静校正量的确定方法及装置
CN112379434B (zh) 适用于沙漠区的层析反演静校正的方法及装置
CN112099089B (zh) 一种基于浮动基准面的山区地震资料静校正量计算方法
CN112711071B (zh) 地层倾角校正方法及装置
CN112444848B (zh) 一种全波形反演的方法及系统
CN112230280B (zh) 增强横波地震数据品质的方法及装置
CN110954940B (zh) 一种基于地表一致性模型估计的初至质量控制方法

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