CN112305590B - 粘声介质叠前时间偏移计算方法及装置 - Google Patents

粘声介质叠前时间偏移计算方法及装置 Download PDF

Info

Publication number
CN112305590B
CN112305590B CN202011006297.4A CN202011006297A CN112305590B CN 112305590 B CN112305590 B CN 112305590B CN 202011006297 A CN202011006297 A CN 202011006297A CN 112305590 B CN112305590 B CN 112305590B
Authority
CN
China
Prior art keywords
constant
variable
seismic trace
trace data
sequence
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
CN202011006297.4A
Other languages
English (en)
Other versions
CN112305590A (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 CN202011006297.4A priority Critical patent/CN112305590B/zh
Publication of CN112305590A publication Critical patent/CN112305590A/zh
Application granted granted Critical
Publication of CN112305590B publication Critical patent/CN112305590B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/282—Application of seismic models, synthetic seismograms
    • 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/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G01V1/305—Travel times
    • 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/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

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)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种粘声介质叠前时间偏移计算方法及装置，该方法包括：根据地震波传播的最大旅行时、最小等效Q因子以及预设的变量R的计算公式计算出变量R的最大值；对所述变量R的最大值进行离散化处理生成常数序列；得到时间域地震道序列，其中，所述常数序列中的各常数与所述时间域地震道序列中的各时间域地震道数据一一对应；确定成像样点位置对应的变量R在所述常数序列中的位置，确定所述常数序列中与该位置相邻的两个常数，并确定该相邻的两个常数各自对应的时间域地震道数据；进行横向加权Kirchhoff积分法求和，以完成所有地震道数据的偏移计算。本发明实现了快速计算粘声介质叠前时间偏移的有益效果。

Description

粘声介质叠前时间偏移计算方法及装置
技术领域
本发明涉及反射波地震数据处理过程中的粘声介质叠前时间域成像技术，具体而言，涉及一种粘声介质叠前时间偏移计算方法及装置。
背景技术
传统的地震波传播理论一般假设地下介质为均匀各向同性和完全弹性，而真实的地球内部介质的黏滞性广泛存在。当地震波在黏滞性介质中传播时，会引起地震波的振幅衰减及相位畸变，且高频成分吸收衰减尤为严重。地下介质黏滞性的存在会导致地震数据频带宽度的减小以及分辨率的降低。对这些数据进行处理时，地震成像往往会遇到深部反射层成像精度低的问题，不利于深层数据的地震解释。因此，发展基于衰减补偿的偏移方法可有效提高地震资料的成像精度。
和传统叠前时间偏移方法相同，粘声介质叠前时间偏移仍然沿用Kirchhoff积分法，其频率域的表达式为：
其中：
I(x,y,T)---成像结果，(x,y)为成像点的空间位置，T为成像时间；
F(ω)---输入地震道在频率域的表达式；
τg和τs---分别是上、下行波的旅行时；
Qeff---等效Q因子；
ω0---地震数据的主频；
令t＝τs+τr，并对上述公式进行整理可以得到：
上述公式(1)中和Qeff相关的两项分别代表振幅衰减项和频散项，每一项是频率、旅行时间、等效Q三个变量的函数，对以上两项的补偿必须对上述三个变量分别离散、按照三重循环的方式完成，所以其计算量及其庞大，较常规叠前时间偏移有三个数量级的增加。因此现有技术缺少一种能够快速计算粘声介质叠前时间偏移的方法。
发明内容
本发明为了解决上述背景技术中的至少一个技术问题，提出了一种粘声介质叠前时间偏移计算方法及装置。
为了实现上述目的，根据本发明的一个方面，提供了一种粘声介质叠前时间偏移计算方法，该方法包括：
根据地震波传播的最大旅行时、最小等效Q因子以及预设的变量R的计算公式计算出变量R的最大值；
对所述变量R的最大值进行离散化处理生成常数序列；
通过快速傅里叶变换把地震道数据转化到频率域得到地震道数据频率域表达式；
分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式；
采用快速反傅里叶变换把每个所述补偿函数表达式转化到时间域，得到时间域地震道序列，其中，所述常数序列中的各常数与所述时间域地震道序列中的各时间域地震道数据一一对应；
从用户提供的Q场数据中获取成像样点位置的等效Q因子和成像样点位置的旅行时，并计算出成像样点位置对应的变量R；
确定成像样点位置对应的变量R在所述常数序列中的位置，确定所述常数序列中与该位置相邻的两个常数，并确定该相邻的两个常数各自对应的时间域地震道数据；
根据该相邻的两个常数以及该相邻的两个常数各自对应的时间域地震道数据进行横向加权Kirchhoff积分法求和，以完成所有地震道数据的偏移计算。
可选的，所述变量R为地震波传播的旅行时与等效Q因子的商。
可选的，所述分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式，具体为：
将各常数分别代入以下公式对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式：
其中，G′1(ω)为补偿函数表达式，F(ω)为地震道数据频率域表达式，R表示所述常数序列中的各常数。
可选的，所述对所述变量R的最大值进行离散化处理生成常数序列，包括：
对所述变量R的最大值按照步长8进行离散化处理，生成常数序列。
为了实现上述目的，根据本发明的另一方面，提供了一种粘声介质叠前时间偏移计算装置，该装置包括：
变量R最大值计算单元，用于根据地震波传播的最大旅行时、最小等效Q因子以及预设的变量R的计算公式计算出变量R的最大值；
常数序列生成单元，用于对所述变量R的最大值进行离散化处理生成常数序列；
快速傅里叶变换单元，用于通过快速傅里叶变换把地震道数据转化到频率域得到地震道数据频率域表达式；
补偿单元，用于分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式；
快速反傅里叶变换单元，用于采用快速反傅里叶变换把每个所述补偿函数表达式转化到时间域，得到时间域地震道序列，其中，所述常数序列中的各常数与所述时间域地震道序列中的各时间域地震道数据一一对应；
成像样点位置数据获取单元，用于从用户提供的Q场数据中获取成像样点位置的等效Q因子和成像样点位置的旅行时，并计算出成像样点位置对应的变量R；
确定单元，用于确定成像样点位置对应的变量R在所述常数序列中的位置，确定所述常数序列中与该位置相邻的两个常数，并确定该相邻的两个常数各自对应的时间域地震道数据；
偏移计算单元，用于根据该相邻的两个常数以及该相邻的两个常数各自对应的时间域地震道数据进行横向加权Kirchhoff积分法求和，以完成所有地震道数据的偏移计算。
为了实现上述目的，根据本发明的另一方面，还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述粘声介质叠前时间偏移计算方法中的步骤。
为了实现上述目的，根据本发明的另一方面，还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序在计算机处理器中执行时实现上述粘声介质叠前时间偏移计算方法中的步骤。
本发明的有益效果为：本发明实施例通过引入变量R，达到了减少粘声介质叠前时间偏移计算公式变量的目的，并用常数扫描代替变量计算，使地震道在频率域的振幅和相位补偿的计算量减少了3个数量级，能够有效的减少粘声介质叠前时间偏移的计算量，实现了对粘声介质叠前时间偏移快速计算。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1是本发明实施例粘声介质叠前时间偏移计算方法的流程图；
图2是叠前时间偏移旅行时计算示意图；
图3是常数离散粘弹性衰减补偿示意图；
图4是二维粘声介质模型示意图；
图5是图4模型数据叠前时间偏移剖面对比示意图；
图6是图4模型Q异常体正下方第二层底界面的频谱对比示意图；
图7是实际数据叠前时间偏移剖面示意图；
图8是图7剖面应用时窗均衡后结果示意图；
图9是实际数据两种叠前时间偏移方法频谱对比示意图；
图10是本发明实施例计算机设备示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
图1是本发明实施例粘声介质叠前时间偏移计算方法的流程图，如图1所示，本发明实施例的粘声介质叠前时间偏移计算方法包括步骤S1至步骤S8。
步骤S1，根据地震波传播的最大旅行时、最小等效Q因子以及预设的变量R的计算公式计算出变量R的最大值。
在本发明实施例中，本发明首先定义变量R、确定其取值范围、并进行离散化处理。在本发明实施例中，所述变量R为地震波传播的旅行时与等效Q因子的商，其表达式为：
R＝t/Qeff
其中，t＝ts+tr为地震波传播的旅行时间，tr和ts分别是上、下行波旅行时，Qeff是等效Q因子。
本发明通过引入变量R，可以达到减少粘声介质叠前时间偏移计算公式变量的目的。
在本发明实施例中，根据地震记录的长度确定地震波传播的最大旅行时tmax。
在本发明实施例中，粘声介质Kirchhoff积分法叠前时间偏移公式(1)中t为旅行时，由两段组成：炮点到成像点ts和成像点到检波点tr，如图2所示。公式(2)是经典的双平方根旅行时计算公式。
公式(2)中变量可以参看图2，含义如下：
T：公式(1)中某一样点的成像时间；
hs：炮点到成像点的距离OS；
hr：检波点到成像点的距离OR；
vrms：成像点位置处的均方根速度；
公式(2)旅行时的最大值应为输入地震道的记录长度tmax。
在本发明实施例中，本发明可以遍历用户提供的等效Q场数据体，统计最小等效Q因子：
本发明需要的等效Q场由用户以外部文件的形式提供，其本质是一个三维数据体，三个维度分别是Inline线号、Xline线号、时间，每个样点对应一正实数值，该数据体中必然存在最大、最小值。本发明需要得到等效Q场的最小值。
在本发明实施例中，本发明用地震波传播的最大旅行时tmax除以最小等效Q因子得到变量R的最大值Rmax，具体公式可以如下：
步骤S2，对所述变量R的最大值进行离散化处理生成常数序列。
在本发明实施例中，本步骤将步骤S1得到的Rmax按照步长8进行离散化处理，形成常数序列R1、R2…Rn。
步骤S3，通过快速傅里叶变换把地震道数据转化到频率域得到地震道数据频率域表达式。
在本发明实施例中，本步骤先读取一个地震道数据，可以用f(t)表示，然后用快速傅里叶变换(FFT)把地震道数据f(t)转化到频率域，得到地震道数据频率域表达式F(ω)。
步骤S4，分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式。
在本发明一个实施例中，本步骤具体可以为，将所述常数序列中的各常数分别代入以下公式对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式：
其中，G′1(ω)为补偿函数表达式，F(ω)为地震道数据频率域表达式，R表示所述常数序列中的各常数，ω为频率为变量，ω0为初始频率为常亮。
在函数G′1(ω)表达式中，除了ω是变量外，离散值Ri、ω0都是常量，所以粘滞介质叠前时间偏移公式(1)可以表示为：
公式(3)和声波弹性介质Kirchhoff叠前时间偏移表达式完全相同。利用该公式将可使粘滞介质叠前时间偏计算量减少到常规叠前时间偏移的级别。
步骤S5，采用快速反傅里叶变换把每个所述补偿函数表达式转化到时间域，得到时间域地震道序列，其中，所述常数序列中的各常数与所述时间域地震道序列中的各时间域地震道数据一一对应。
在本发明实施例中，本步骤用快速反傅里叶变换(IFFT)把每个所述常数各自对应的补偿函数表达式G′1(ω)转化到时间域，得到时间域地震道序列f1(t)、f2(t)…fn(t)。其中，常数序列R1、R2…Rn中的各常数与时间域地震道序列f1(t)、f2(t)…fn(t)中的各时间域地震道数据一一对应。
图3是本步骤的物理意义，图示为一个二维地震道集，横坐标对应发明中变量R，纵坐标表示地震样点时间，从左到右每一道对应一个R的离散值。
本发明通过以上步骤，用常数扫描代替变量计算，使地震道在频率域的振幅和相位补偿的计算量减少了3个数量级。
步骤S6，从用户提供的Q场数据中获取成像样点位置的等效Q因子和成像样点位置的旅行时，并计算出成像样点位置对应的变量R。
步骤S7，确定成像样点位置对应的变量R在所述常数序列中的位置，确定所述常数序列中与该位置相邻的两个常数，并确定该相邻的两个常数各自对应的时间域地震道数据。
步骤S8，根据该相邻的两个常数以及该相邻的两个常数各自对应的时间域地震道数据进行横向加权Kirchhoff积分法求和，以完成所有地震道数据的偏移计算。
步骤S6至步骤S8对应粘声介质Kirchhoff积分法求和，具体的步骤包括：
1、根据地震道数据f(t)所在中心点位置的Inline线号和Xline线号、用户提供的偏移孔径参数值计算出本地震道涉及的积分法求和的Inline线号范围inlinemin、inlinemax和Xnline线号范围xlinemin、xlinemax；
2、根据用户提供的参数确定偏移结果最小时间TOmin、最大时间TOmax和增量TOinc
3、Kirchhoff积分法求和过程用计算机程序语言表示为：
本发明通过采用公式(4)的加权线性求和，保障计算结果的精度。本发明采用横向两个地震道求和，两个地震道的权值由两个离散R和实际R的差值确定。
在积分法求和三重循环中，黑体字为本专利发明新增加内容，其意义在图3中用双箭头标记，即在积分法求和过程中，要计算每个样点处变量R的瞬时值Rk，根据瞬时值Rk找到与之相邻的两个常数Ri和Ri+1，再用与这两个常数对应的地震道上的样点经线性插值计算偏移结果的值。
下面结合实例来介绍本发明粘声介质叠前时间偏移计算方法的实际应用效果。
本发明技术已经在我国大型地震数据处理解释一体化系统(GeoEast)成功开发为工业化程序，并分别用二维模型和三维实际数据进行了应用效果验证。
图4是二维粘声介质模型，从左到右分别是层速度、Q和等效Q场，模型共三层，深度依次为1300米、2300米和3000米，层速度分别是2300米/秒、3000米/秒、4000米/秒。在模型横向中间位置第一层底部设计了一个Q异常体，其值为10(强粘声介质)，其它部位Q值是10000(表示声波弹性介质)。等效Q是根据实际Q利用理论公式计算得到的。图5是常规叠前时间偏移偏移剖面(图5上)和Q叠前时间偏移偏移剖面(图5下)，在常规叠前时间偏移剖面中，Q异常正下方，两个层位能量相比两端同相轴有明显的变弱，在Q叠前时间偏移剖面中弱能量得到明显的加强。图6是Q异常体正下方第二层底界面的频谱：图6左边为常规叠前时间偏移、图6右边是Q叠前时间偏移，左边频谱主频不到10HZ，右边频谱主频接近14HZ，频率提高约4HZ。
图7是三维实际数据不同处理方法的叠前时间偏移结果对比，左边是常规叠前时间偏移剖面，从浅到深地震振幅有明显差异，右边为Q叠前时间偏移成果，地震振幅相对一致。由于两种叠前时间偏移结果地震振幅量级存在差异，很难进行详细对比，为此对两个偏移结果分别应用相同参数的时窗均衡，图8是图7数据时窗均衡后剖面，同样左边是常规叠前时间偏移剖面，右边为Q叠前时间偏移剖面。在红色椭圆标记处，同相轴连续性得到明显改善，在红色箭头标记处，Q叠前时间偏移结果波组被打开，分辨率得到明显提升。图9是实际三维数据叠前时间偏移的频谱：左上对应常规叠前时间偏移浅层、右上对应Q叠前时间偏移浅层、左下对应常规叠前时间偏移深层、右下对应Q叠前时间偏移深层。从频谱图上可以看出，两种方法浅层频谱基本一致，但是深层频谱有明显差别，常规叠前时间偏移深层频谱主频约为23HZ，Q叠前时间偏移深层频谱主频提升到40HZ，频率提高约17HZ。
需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
基于同一发明构思，本发明实施例还提供了一种粘声介质叠前时间偏移计算装置，可以用于实现上述实施例所描述的粘声介质叠前时间偏移计算方法，如下面的实施例所述。由于粘声介质叠前时间偏移计算装置解决问题的原理与粘声介质叠前时间偏移计算方法相似，因此粘声介质叠前时间偏移计算装置的实施例可以参见粘声介质叠前时间偏移计算方法的实施例，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
在本发明实施例中，本发明的粘声介质叠前时间偏移计算装置包括：
变量R最大值计算单元，用于根据地震波传播的最大旅行时、最小等效Q因子以及预设的变量R的计算公式计算出变量R的最大值；
常数序列生成单元，用于对所述变量R的最大值进行离散化处理生成常数序列；
快速傅里叶变换单元，用于通过快速傅里叶变换把地震道数据转化到频率域得到地震道数据频率域表达式；
补偿单元，用于分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式；
快速反傅里叶变换单元，用于采用快速反傅里叶变换把每个所述补偿函数表达式转化到时间域，得到时间域地震道序列，其中，所述常数序列中的各常数与所述时间域地震道序列中的各时间域地震道数据一一对应；
成像样点位置数据获取单元，用于从所述地震道数据获取成像样点位置的等效Q因子和成像样点位置的旅行时，并计算出成像样点位置对应的变量R；
确定单元，用于确定成像样点位置对应的变量R在所述常数序列中的位置，确定所述常数序列中与该位置相邻的两个常数，并确定该相邻的两个常数各自对应的时间域地震道数据；
偏移计算单元，用于根据该相邻的两个常数以及该相邻的两个常数各自对应的时间域地震道数据进行Kirchhoff积分法求和，以完成所有地震道数据的偏移计算。
在本发明一个实施例中，所述变量R为地震波传播的旅行时与等效Q因子的商。
在本发明一个实施例中，所述补偿单元，具体用于将各常数分别代入以下公式对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式：
其中，G′1(ω)为补偿函数表达式，F(ω)为地震道数据频率域表达式，R表示所述常数序列中的各常数。
在本发明一个实施例中，所述常数序列生成单元，具体用于对所述变量R的最大值按照步长8进行离散化处理，生成常数序列。
为了实现上述目的，根据本申请的另一方面，还提供了一种计算机设备。如图10所示，该计算机设备包括存储器、处理器、通信接口以及通信总线，在存储器上存储有可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述实施例方法中的步骤。
处理器可以为中央处理器(Central Processing Unit，CPU)。处理器还可以为其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片，或者上述各类芯片的组合。
存储器作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及单元，如本发明上述方法实施例中对应的程序单元。处理器通过运行存储在存储器中的非暂态软件程序、指令以及模块，从而执行处理器的各种功能应用以及作品数据处理，即实现上述方法实施例中的方法。
存储器可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储处理器所创建的数据等。此外，存储器可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中，存储器可选包括相对于处理器远程设置的存储器，这些远程存储器可以通过网络连接至处理器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
所述一个或者多个单元存储在所述存储器中，当被所述处理器执行时，执行上述实施例中的方法。
上述计算机设备具体细节可以对应参阅上述实施例中对应的相关描述和效果进行理解，此处不再赘述。
为了实现上述目的，根据本申请的另一方面，还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序在计算机处理器中执行时实现上述粘声介质叠前时间偏移计算方法中的步骤。本领域技术人员可以理解，实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)、随机存储记忆体(RandomAccessMemory，RAM)、快闪存储器(Flash Memory)、硬盘(HardDisk Drive，缩写：HDD)或固态硬盘(Solid-State Drive，SSD)等；所述存储介质还可以包括上述种类的存储器的组合。
显然，本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种粘声介质叠前时间偏移计算方法，其特征在于，包括：
根据地震波传播的最大旅行时、最小等效Q因子以及预设的变量R的计算公式计算出变量R的最大值，变量R为地震波传播的旅行时与等效Q因子的商，其计算公式为：
R＝t/Qeff
其中，t＝ts+tr为地震波传播的旅行时，tr和ts分别是上、下行波旅行时，Qeff是等效Q因子；
对所述变量R的最大值进行离散化处理生成常数序列；
通过快速傅里叶变换把地震道数据转化到频率域得到地震道数据频率域表达式；
分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式；
采用快速反傅里叶变换把每个所述补偿函数表达式转化到时间域，得到时间域地震道序列，其中，所述常数序列中的各常数与所述时间域地震道序列中的各时间域地震道数据一一对应；
从用户提供的Q场数据体中获取成像样点位置的等效Q因子和成像样点位置的旅行时，并计算出成像样点位置对应的变量R；
确定成像样点位置对应的变量R在所述常数序列中的位置，确定所述常数序列中与该位置相邻的两个常数，并确定该相邻的两个常数各自对应的时间域地震道数据；
根据该相邻的两个常数以及该相邻的两个常数各自对应的时间域地震道数据进行横向加权Kirchhoff积分法求和，以完成所有地震道数据的偏移计算；
所述分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式，具体为：
将各常数分别代入以下公式对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式：
其中，G'1(w)为补偿函数表达式，F(w)为地震道数据频率域表达式，R表示所述常数序列中的各常数，w为频率，w0为初始频率，j表示复数的虚部。
2.根据权利要求1所述的粘声介质叠前时间偏移计算方法，其特征在于，所述变量R为地震波传播的旅行时与等效Q因子的商。
3.根据权利要求1所述的粘声介质叠前时间偏移计算方法，其特征在于，所述对所述变量R的最大值进行离散化处理生成常数序列，包括：
对所述变量R的最大值按照步长8进行离散化处理，生成常数序列。
4.一种粘声介质叠前时间偏移计算装置，其特征在于，包括：
变量R最大值计算单元，用于根据地震波传播的最大旅行时、最小等效Q因子以及预设的变量R的计算公式计算出变量R的最大值，变量R为地震波传播的旅行时与等效Q因子的商，其计算公式为：
R＝t/Qeff
其中，t＝ts+tr为地震波传播的旅行时，tr和ts分别是上、下行波旅行时，Qeff是等效Q因子；
常数序列生成单元，用于对所述变量R的最大值进行离散化处理生成常数序列；
快速傅里叶变换单元，用于通过快速傅里叶变换把地震道数据转化到频率域得到地震道数据频率域表达式；
补偿单元，用于分别采用所述常数序列中的各常数对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式；
快速反傅里叶变换单元，用于采用快速反傅里叶变换把每个所述补偿函数表达式转化到时间域，得到时间域地震道序列，其中，所述常数序列中的各常数与所述时间域地震道序列中的各时间域地震道数据一一对应；
成像样点位置数据获取单元，用于从用户提供的Q场数据体中获取成像样点位置的等效Q因子和成像样点位置的旅行时，并计算出成像样点位置对应的变量R；
确定单元，用于确定成像样点位置对应的变量R在所述常数序列中的位置，确定所述常数序列中与该位置相邻的两个常数，并确定该相邻的两个常数各自对应的时间域地震道数据；
偏移计算单元，用于根据该相邻的两个常数以及该相邻的两个常数各自对应的时间域地震道数据进行横向加权Kirchhoff积分法求和，以完成所有地震道数据的偏移计算；
所述补偿单元，具体用于将各常数分别代入以下公式对所述地震道数据频率域表达式进行补偿，得到每个所述常数各自对应的补偿函数表达式：
其中，G1'(w)为补偿函数表达式，F(w)为地震道数据频率域表达式，R表示所述常数序列中的各常数，w为频率，w0为初始频率，j表示复数的虚部。
5.根据权利要求4所述的粘声介质叠前时间偏移计算装置，其特征在于，所述变量R为地震波传播的旅行时与等效Q因子的商。
6.根据权利要求4所述的粘声介质叠前时间偏移计算装置，其特征在于，所述常数序列生成单元，具体用于对所述变量R的最大值按照步长8进行离散化处理，生成常数序列。
7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至3任一项所述的方法。
8.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序在计算机处理器中执行时实现如权利要求1至3任意一项所述的方法。
CN202011006297.4A 2020-09-23 2020-09-23 粘声介质叠前时间偏移计算方法及装置 Active CN112305590B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011006297.4A CN112305590B (zh) 2020-09-23 2020-09-23 粘声介质叠前时间偏移计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011006297.4A CN112305590B (zh) 2020-09-23 2020-09-23 粘声介质叠前时间偏移计算方法及装置

Publications (2)

Publication Number Publication Date
CN112305590A CN112305590A (zh) 2021-02-02
CN112305590B true CN112305590B (zh) 2023-12-22

Family

ID=74489092

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011006297.4A Active CN112305590B (zh) 2020-09-23 2020-09-23 粘声介质叠前时间偏移计算方法及装置

Country Status (1)

Country Link
CN (1) CN112305590B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102590862A (zh) * 2012-01-19 2012-07-18 中国科学院地质与地球物理研究所 补偿吸收衰减的叠前时间偏移方法
CN103424777A (zh) * 2013-07-01 2013-12-04 中国科学院地质与地球物理研究所 一种提高地震成像分辨率的方法
CN106842304A (zh) * 2017-01-03 2017-06-13 中国石油天然气集团公司 一种叠前深度偏移方法及装置
CN108919354A (zh) * 2018-09-27 2018-11-30 中国科学院地质与地球物理研究所 近地表q偏移方法及装置
CN109085648A (zh) * 2018-07-16 2018-12-25 中国科学院地质与地球物理研究所 叠前深度偏移方法和装置
CN110031898A (zh) * 2019-05-14 2019-07-19 中国科学院地质与地球物理研究所 数据优化方法及积分法叠前深度偏移方法
CN110873900A (zh) * 2018-09-04 2020-03-10 中国石油化工股份有限公司 频率域叠前地震道q补偿方法及系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9829592B2 (en) * 2014-12-16 2017-11-28 Pgs Geophysical As Seismic imaging with visco-acoustic reverse-time migration using pseudo-analytical method
EP3076205B1 (en) * 2015-03-31 2023-06-14 CGG Services SAS Method for survey data processing compensating for visco-acoustic effects in tilted transverse isotropy reverse time migration
CA2998522A1 (en) * 2015-10-02 2017-04-06 Exxonmobil Upstream Research Company Q-compensated full wavefield inversion
US11294087B2 (en) * 2018-03-27 2022-04-05 Exxonmobil Upstream Research Company Directional Q compensation with sparsity constraints and preconditioning

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102590862A (zh) * 2012-01-19 2012-07-18 中国科学院地质与地球物理研究所 补偿吸收衰减的叠前时间偏移方法
CN103424777A (zh) * 2013-07-01 2013-12-04 中国科学院地质与地球物理研究所 一种提高地震成像分辨率的方法
CN106842304A (zh) * 2017-01-03 2017-06-13 中国石油天然气集团公司 一种叠前深度偏移方法及装置
CN109085648A (zh) * 2018-07-16 2018-12-25 中国科学院地质与地球物理研究所 叠前深度偏移方法和装置
CN110873900A (zh) * 2018-09-04 2020-03-10 中国石油化工股份有限公司 频率域叠前地震道q补偿方法及系统
CN108919354A (zh) * 2018-09-27 2018-11-30 中国科学院地质与地球物理研究所 近地表q偏移方法及装置
CN110031898A (zh) * 2019-05-14 2019-07-19 中国科学院地质与地球物理研究所 数据优化方法及积分法叠前深度偏移方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Improving Imaging and Inversion Through Least-Squares Q-Kirchhoff APSDM: A Case Study from Offshore China;Chenghai Jiao等;IPTC;全文 *
积分法Ｑ成像技术的应用研究;赵玲芝等;石油地球物理勘探;全文 *
粘性稳相叠前时间偏移方法的研究及应用;袁安龙等;中国石油学会2019年物探技术研讨会;全文 *
近地表Q值模型建立方法及其地震叠前补偿应用;陈志德;王成;刘国友;丁吉丰;;石油学报(第02期);全文 *
黏滞声学介质地震波吸收补偿叠前时间偏移方法;陈志德;赵忠华;王成;;石油地球物理勘探(第02期);全文 *

Also Published As

Publication number Publication date
CN112305590A (zh) 2021-02-02

Similar Documents

Publication Publication Date Title
KR101092668B1 (ko) 파형 역산을 이용한 지하 구조의 영상화 장치와 방법
Groos et al. The role of attenuation in 2D full-waveform inversion of shallow-seismic body and Rayleigh waves
Liu et al. Seismic data interpolation beyond aliasing using regularized nonstationary autoregression
CN104133241B (zh) 波场分离方法和装置
US8990053B2 (en) Method of wavelet estimation and multiple prediction in full wavefield inversion
CN107153216B (zh) 确定地震波场的坡印廷矢量的方法、装置以及计算机存储介质
Wu et al. Efficient scattering-angle enrichment for a nonlinear inversion of the background and perturbations components of a velocity model
US10185046B2 (en) Method for temporal dispersion correction for seismic simulation, RTM and FWI
Liu et al. Adaptive prediction filtering in txy domain for random noise attenuation using regularized nonstationary autoregression
Chong et al. Joint inversion of body‐wave receiver function and Rayleigh‐wave ellipticity
Mittet Second-order time integration of the wave equation with dispersion correction procedures
Wang et al. Separation and imaging of seismic diffractions using a localized rank-reduction method with adaptively selected ranks
WO2018102813A2 (en) Seismic acquisition geometry full-waveform inversion
CN105425298B (zh) 一种消除有限差分正演过程中数值频散的方法和装置
Hawkins A spectral element method for surface wave dispersion and adjoints
CN111060961B (zh) 基于多信息约束反演的品质因子确定方法、装置及系统
Tatlock et al. An assessment of the roller approach for wave breaking in a hybrid finite-volume finite-difference Boussinesq-type model for the surf-zone
CN108445539B (zh) 一种消除地震子波旁瓣干扰的方法、设备及系统
de Lucena et al. Rayleigh wave modeling: A study of dispersion curve sensitivity and methodology for calculating an initial model to be included in an inversion algorithm
CN112305590B (zh) 粘声介质叠前时间偏移计算方法及装置
CN106547019A (zh) 一种确定地层品质因子的方法
CN113866823A (zh) 一种粘声各向异性介质中的正演成像方法
You et al. Q-compensated wavefield depth extrapolation-based migration using a viscoacoustic wave equation
Zhou et al. Value of nonstationary wavelet spectral balancing in mapping a faulted fluvial system, Bohai Gulf, China
CN111665550A (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
GR01 Patent grant
GR01 Patent grant