CN108845351A - A kind of VSP seismic data converted wave full waveform inversion method - Google Patents

A kind of VSP seismic data converted wave full waveform inversion method Download PDF

Info

Publication number
CN108845351A
CN108845351A CN201810666468.2A CN201810666468A CN108845351A CN 108845351 A CN108845351 A CN 108845351A CN 201810666468 A CN201810666468 A CN 201810666468A CN 108845351 A CN108845351 A CN 108845351A
Authority
CN
China
Prior art keywords
wave
longitudinal
transverse
vsp
inversion
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
CN201810666468.2A
Other languages
Chinese (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 University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201810666468.2A priority Critical patent/CN108845351A/en
Publication of CN108845351A publication Critical patent/CN108845351A/en
Pending legal-status Critical Current

Links

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/282Application of seismic models, synthetic seismograms
    • 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

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

本发明公开了一种VSP地震资料转换波全波形反演方法。通过全局优化参数反演法分离上行PS转换波;根据经验公式计算横波速度并对其进行平滑处理得到初始横波速度模型;通过求解纵横波分离的弹性波方程得到每一时间步的正传纵横波波场以及模拟VSP PS波记录;通过求解纵横波分离的伴随弹性波方程得到每一时间步的反传纵横波波场;采用正传纵波波场和反传横波波场计算横波速度梯度;采用三点抛物线拟合法求取迭代步长;通过L‑BFGS方法更新横波速度模型;运用维纳滤波从低频到高频进行多尺度反演。本发明的有益效果是通过纵横波分离策略来提高VSP地震资料转换波全波形反演的精度,提高复杂构造的反演效果。

The invention discloses a full waveform inversion method of converted wave of VSP seismic data. Separation of upgoing PS converted waves by global optimization parameter inversion method; calculation of shear wave velocity according to empirical formula and smoothing to obtain initial shear wave velocity model; solution of elastic wave equation for separation of longitudinal and transverse waves to obtain forward longitudinal and transverse waves at each time step Wave field and simulated VSP PS wave records; by solving the adjoint elastic wave equation of the separation of longitudinal and transverse waves, the reverse longitudinal and transverse wave field at each time step is obtained; the forward longitudinal wave field and the backward transverse wave field are used to calculate the shear wave velocity gradient; The three-point parabola fitting method is used to obtain the iteration step size; the shear wave velocity model is updated by the L-BFGS method; the Wiener filter is used to perform multi-scale inversion from low frequency to high frequency. The beneficial effect of the invention is to improve the accuracy of full waveform inversion of converted wave of VSP seismic data through the separation strategy of longitudinal and transverse waves, and improve the inversion effect of complex structures.

Description

一种VSP地震资料转换波全波形反演方法A Converted Wave Full Waveform Inversion Method for VSP Seismic Data

技术领域technical field

本发明属于勘探地球物理学技术领域,涉及一种全波形反演新方法,可用于VSP地震资料反演来获取高精度的转换波速度模型。The invention belongs to the technical field of exploration geophysics, and relates to a new full waveform inversion method, which can be used for VSP seismic data inversion to obtain a high-precision converted wave velocity model.

背景技术Background technique

随着油气勘探开发的深入和井中地震技术的不断提高,VSP(垂直地震剖面)勘探技术的研究与应用,也得到了快速的发展。在复杂构造或复杂介质地区,VSP技术能有效地了解井孔附近的地质结构、岩石物性等特征。通过VSP成像结合地面三维地震资料和岩性资料,能进行精细的地震属性分析,可以进一步核实关键层位、搞清岩层之间的接触关系、建立起精确的地质模型,进行地层的非均质性研究,为开发方案的调整提供技术支持。与地面地震相比,VSP资料含有更丰富的波场信息,具有分辨率高,信噪比高,且井中接收到的地震记录只穿过一次低减速带,因此相比地面资料地震波能量的衰减更小,所以对深层成像有其特有优势。其次地面地震勘探由于观测限制,很难观测到目标层盐丘体及其侧翼的和盐下信息。而VSP技术可以弥补地面资料的这些不足。因此无论在勘探领域,还是在开发领域,VSP技术都有其特有的优势和作用。With the deepening of oil and gas exploration and development and the continuous improvement of borehole seismic technology, the research and application of VSP (vertical seismic profiling) exploration technology has also been developed rapidly. In areas with complex structures or complex media, VSP technology can effectively understand the geological structure, petrophysical properties and other characteristics near the wellbore. Through VSP imaging combined with ground three-dimensional seismic data and lithological data, fine seismic attribute analysis can be carried out, and key layers can be further verified, the contact relationship between rock layers can be clarified, an accurate geological model can be established, and stratum heterogeneity Sexual research, providing technical support for the adjustment of the development plan. Compared with surface seismic data, VSP data contain richer wave field information, high resolution, high signal-to-noise ratio, and the seismic records received in the well only pass through a low deceleration zone once, so the attenuation of seismic wave energy compared with surface data Smaller, so it has its unique advantages for deep imaging. Secondly, due to observation limitations in ground seismic exploration, it is difficult to observe the salt dome body and its flanks and sub-salt information in the target layer. VSP technology can make up for these deficiencies in ground data. Therefore, no matter in the field of exploration or development, VSP technology has its unique advantages and functions.

速度是地震资料数据处理流程中的核心参数。高质量的速度建模方法对后续的偏移成像、岩性、岩相解释和储层预测起决定性作用。到目前为止,VSP纵波资料的数据处理流程已相对完善,对于纵波速度的反演精度尚且可以。但由于收诸多因素影响(如横波传播极其不对称性、纵横波耦合串扰、横波各向异性等),转换波速度反演精度仍然较低。弹性全波形反演(EFWI)具有获取精确纵波速度、横波速度和密度的潜力,因此在勘探地球物理学中变得越来越流行。但弹性波全波形反演由于多参数耦合,非线性严重、多解性强,在工业上还没能大规模推广。为消除纵波速度对横波速度反演精度的影响,有必要进一步研究转换波全波形反演方法。Velocity is the core parameter in the data processing flow of seismic data. High-quality velocity modeling methods are decisive for subsequent migration imaging, lithology, lithofacies interpretation, and reservoir prediction. So far, the data processing flow of VSP P-wave data has been relatively perfect, and the inversion accuracy of P-wave velocity is still acceptable. However, due to the influence of many factors (such as the extreme asymmetry of shear wave propagation, the coupling crosstalk of longitudinal and transverse waves, and the anisotropy of shear waves, etc.), the accuracy of converted wave velocity inversion is still low. Elastic full waveform inversion (EFWI) has become increasingly popular in exploration geophysics for its potential to obtain accurate P-wave velocities, S-wave velocities, and densities. However, due to multi-parameter coupling, elastic wave full waveform inversion has serious nonlinearity and strong multi-solution, so it has not been widely promoted in industry. In order to eliminate the influence of P-wave velocity on the accuracy of S-wave velocity inversion, it is necessary to further study the converted wave full waveform inversion method.

本发明给出一种VSP地震资料转换波全波形反演新方法。针对VSP观测方式(地表放炮、井中接收),通过全局优化参数反演法分离PS波(转换波)记录,采用分离出的转换波进行反演来缓解非线性和多解性;通过求解纵横波分离的弹性波方程和伴随方程得到不同时刻的纵横波波场,采用正向纵波波场和伴随横波波场计算目标函数关于横波速度的梯度,减弱不同波模式间的串扰效应,来提高横波速度的反演精度。The invention provides a new method for full waveform inversion of converted wave of VSP seismic data. For the VSP observation mode (surface blasting, well reception), the PS wave (converted wave) record is separated by the global optimization parameter inversion method, and the separated converted wave is used for inversion to alleviate nonlinearity and multi-solution; by solving the longitudinal and transverse waves Separate the elastic wave equation and the adjoint equation to obtain the longitudinal and transverse wave fields at different times, and use the forward longitudinal wave field and the accompanying shear wave field to calculate the gradient of the objective function with respect to the shear wave velocity, weaken the crosstalk effect between different wave modes, and increase the shear wave velocity inversion accuracy.

发明内容Contents of the invention

本发明的目的在于提供一种VSP地震资料转换波全波形反演方法。The purpose of the present invention is to provide a full waveform inversion method of converted wave of VSP seismic data.

本发明所采用的技术方案是按照以下步骤进行:The technical scheme adopted in the present invention is to carry out according to the following steps:

(1)原始VSP(垂直地震剖面)地震资料PS波(转换波)分离;(1) PS wave (converted wave) separation of original VSP (vertical seismic profile) seismic data;

通过全局优化参数反演法分离上行PS转换波;Separation of uplink PS converted waves by global optimization parameter inversion method;

(2)横波速度初始模型建立;(2) Establishment of the initial model of shear wave velocity;

输入纵波速度模型,根据经验公式获得横波速度模型,并对其进行平滑处理得到初始横波速度模型;Input the longitudinal wave velocity model, obtain the shear wave velocity model according to the empirical formula, and smooth it to obtain the initial shear wave velocity model;

(3)震源波场正向延拓,分离纵横波波场;(3) The source wave field is extended forward to separate the longitudinal and transverse wave fields;

通过求解纵横波分离的弹性波方程得到每一时间步的正传纵横波波场以及模拟的VSP PS波记录。The forward P-S wave field and the simulated VSP PS-wave record are obtained at each time step by solving the elastic wave equation for the separation of P- and S-waves.

(4)转换波残差反向延拓,分离纵横波波场;(4) Reverse continuation of the converted wave residual to separate the longitudinal and transverse wave fields;

通过求解纵横波分离的伴随弹性波方程得到每一时间步的反传纵横波波场;By solving the adjoint elastic wave equations of the separation of the longitudinal and transverse waves, the backpropagating longitudinal and transverse wave fields are obtained at each time step;

(5)计算横波速度梯度。(5) Calculate the shear wave velocity gradient.

采用第(3)步和第(4)步求得的正传纵波波场和反传横波波场计算横波速度梯度;Calculate the shear wave velocity gradient using the forward longitudinal wave field and the reverse propagating shear wave field obtained in steps (3) and (4);

(6)求取迭代步长;(6) Find the iteration step size;

采用三点抛物线拟合法求取迭代步长;Use the three-point parabola fitting method to find the iteration step size;

(7)更新横波速度模型;(7) Update the shear wave velocity model;

不满足收敛条件时,重复步骤(3)到步骤(6),满足收敛条件时,迭代结束;When the convergence condition is not met, repeat step (3) to step (6), when the convergence condition is met, the iteration ends;

(8)多尺度反演;(8) Multi-scale inversion;

运用维纳滤波对VSP PS转换波记录进行分频处理,从低频到高频逐步进行反演。Wiener filtering is used to divide the VSP PS converted wave records, and the inversion is carried out step by step from low frequency to high frequency.

具体的,VSP地震资料转换波全波形反演步骤(1)中原始VSP地震资料PS转换波分离方法如下:Specifically, the PS converted wave separation method of the original VSP seismic data in step (1) of the converted wave full waveform inversion of VSP seismic data is as follows:

参数反演法是一种基于地震波视速度和偏振特性的VSP资料波场分离方法,本发明中采用全局优化参数反演法分离PS转换波。The parameter inversion method is a method for separating wave fields of VSP data based on seismic wave apparent velocity and polarization characteristics. In the present invention, the global optimized parameter inversion method is used to separate PS converted waves.

对每个深度点,建立观测数据与模拟记录的极小化目标泛函:For each depth point, the minimized target functional of observed data and simulated records is established:

其中,ΨPdownPupSdown和ΨSup分别为下行纵波、上行纵波、下行横波和上行横波的入射角,vP和vS为纵横波速度,u和uobs为模拟和观测波场,ω为角频率,[ω11]为频率范围。Among them, Ψ Pdown , Ψ Pup , Ψ Sdown and Ψ Sup are the incident angles of downgoing P-wave, upgoing P-wave, downgoing S-wave and upgoing S-wave, respectively; v P and v S are the P- and S-wave velocities; u and u obs are the simulated and observed wave fields , ω is the angular frequency, [ω 11 ] is the frequency range.

通过全局优化遗传算法迭代可以得到最优的传播速度和传播角度。求解方程(2)实现下行纵波、上行纵波、下行横波和上行横波分离。The optimal propagation speed and propagation angle can be obtained through the iteration of the global optimization genetic algorithm. Solving equation (2) realizes the separation of downgoing longitudinal wave, upgoing longitudinal wave, downgoing shear wave and upgoing shear wave.

u(ω)=[u(-M1,ω)…u(0,ω)…u(M2,ω)]T u(ω)=[u(-M 1 ,ω)...u(0,ω)...u(M 2 ,ω)] T

=aPdownuPdown(ω)+aPupuPup(ω)+aSdownuSdown(ω)+aSupuSup(ω), (2)=a Pdown u Pdown (ω)+a Pup u Pup (ω)+a Sdown u Sdown (ω)+a Sup u Sup (ω), (2)

其中,[uPdown,uPup,uSdown,uSup]分别为下行纵波、上行纵波、下行横波和上行横波波场,[aPdown,aPup,aSdown,aup]为相应的偏振矢量,[-M1,…,0,…,M2]为空间窗内检波点的编号。Among them, [u Pdown , u Pup , u Sdown , u Sup ] are the wave fields of downgoing longitudinal wave, upgoing longitudinal wave, downgoing shear wave and upgoing shear wave respectively, and [a Pdown , a Pup , a Sdown , a up ] are the corresponding polarization vectors, [-M 1 ,…,0,…,M 2 ] are the numbers of the receiver points in the spatial window.

当震源为纵波源时,分离出的上行横波uSup(频率域)进行傅里叶反变换就得到了观测PS转换波记录(根据传播角度分解为水平分量和垂直分量)。When the seismic source is a longitudinal wave source, the separated upward shear wave u Sup (frequency domain) is inversely Fourier transformed to obtain the observed PS converted wave record (decomposed into horizontal and vertical components according to the propagation angle).

步骤(2)中横波速度初始模型建立方法如下:The establishment method of the initial model of shear wave velocity in step (2) is as follows:

纵横波速度经验公式:The empirical formula for compressional and shear wave velocity:

根据输入的纵波速度模型,采用经验公式获得横波速度模型,通过高斯窗对其进行平滑得到初始横波速度模型。According to the input P-wave velocity model, the empirical formula is used to obtain the shear-wave velocity model, which is smoothed by Gaussian window to obtain the initial shear-wave velocity model.

步骤(3)中,震源波场正向延拓,分离纵横波波场方法如下:In step (3), the source wave field is extended forward, and the method of separating the longitudinal and transverse wave fields is as follows:

一阶弹性波速度-应力方程:The first-order elastic wave velocity-stress equation:

其中,λ和μ为拉梅常数,ρ为密度,(vx,vz)为质点振动速度矢量,(τxxzzxz)为应力矢量。Among them, λ and μ are Lame constants, ρ is density, (v x , v z ) is particle vibration velocity vector, (τ xxzzxz ) is stress vector.

本发明中采用纵横波分离的速度-应力弹性波方程:Adopt the velocity-stress elastic wave equation of longitudinal and transverse wave separation in the present invention:

其中,(vxp,vzp)为纵波质点振动速度矢量,(vxs,vzs)为横波质点振动速度矢量,(τxxpzzp)为纵波应力矢量,(τxxszzsxzs)为横波应力矢量。Among them, (v xp , v zp ) is the particle vibration velocity vector of longitudinal wave, (v xs , v zs ) is the particle vibration velocity vector of shear wave, (τ xxpzzp ) is the stress vector of longitudinal wave, (τ xxszzsxzs ) is the shear wave stress vector.

采用高阶交错网格有限差分数值求解方程(6)得到每一时间步的正传纵横波波场以及模拟VSP PS波记录。The high-order staggered grid finite difference method is used to numerically solve Equation (6) to obtain the forward P- and S-wave field and simulated VSP PS wave records at each time step.

步骤(4)转换波残差反向延拓,分离纵横波波场方法如下:Step (4) Reverse continuation of the converted wave residual, and the method of separating the longitudinal and transverse wave fields is as follows:

在伴随状态法中,伴随算子满足:In the adjoint state method, the adjoint operator satisfies:

<Lu,u*>=<u,L*u*>, (7)<Lu,u * >=<u,L * u * >, (7)

其中:L为正演算子,u为质点速度和应力矢量,“*”为其伴随状态。将方程(4)带入方程(7)得到常规伴随方程:Among them: L is the forward operator, u is the particle velocity and stress vector, "*" is its adjoint state. Substituting Equation (4) into Equation (7) yields the conventional adjoint equation:

同理,基于方程(6)可以得到纵横波分离的伴随方程:Similarly, based on equation (6), the adjoint equation for the separation of P- and S-waves can be obtained:

其中,“*”表示相应波场的伴随波场,为纵横波耦合的观测和模拟VSP地震记录,为分离的观测和模拟VSP转换波地震记录。Among them, "*" indicates the accompanying wave field of the corresponding wave field, For the observation and simulation of VSP seismic records coupled with compressional and shear waves, Separated observations and simulated VSP converted wave seismic records.

以模拟和观测VSP PS波记录残差为震源,采用高阶交错有限差分数值求解方程(10)和(11)得到每一时间步的反传纵横波波场。Taking the simulated and observed VSP PS wave recording residuals as the seismic source, equations (10) and (11) are numerically solved using high-order staggered finite differences to obtain the backpropagating P-S wave field at each time step.

步骤(5)计算横波速度梯度方法如下:Step (5) The method of calculating the shear wave velocity gradient is as follows:

传统全波形反演中横波速度梯度为:The shear wave velocity gradient in traditional full waveform inversion is:

其中,E为目标泛函。Among them, E is the target functional.

本发明采用正传纵波波场和反传横波波场计算横波速度梯度:The present invention calculates the shear wave velocity gradient by using the forward longitudinal wave field and the backpropagation shear wave field:

其中,in,

通过采用新的梯度公式(方程13)可以缓解纵波速度的干扰,进而提高横波速度的反演精度。By adopting the new gradient formula (Equation 13), the interference of P-wave velocity can be alleviated, and the inversion accuracy of S-wave velocity can be improved.

步骤(6)求取迭代步长方法如下:Step (6) The method of calculating the iteration step size is as follows:

本发明中采用三点抛物线插值方法计算步长。选取三个试探点,用这三个点构建一条抛物线,这条抛物线的极值点既为最优步长。In the present invention, a three-point parabolic interpolation method is used to calculate the step size. Select three trial points, and use these three points to construct a parabola. The extreme point of this parabola is the optimal step size.

三个试探步长a0,a1,a2对应的目标泛函为Es0,Es1,Es2,则最优步长为:The target functionals corresponding to the three trial steps a 0 , a 1 , and a 2 are E s0 , E s1 , E s2 , and the optimal step size is:

步骤(7)更新横波速度模型方法如下:Step (7) The method of updating the shear wave velocity model is as follows:

通过L-BFGS方法更新横波速度:Update the shear wave velocity by the L-BFGS method:

其中,mk和mk+1分别为当前迭代和下一次迭代的横波速度模型,Hk为海森矩阵逆的近似矩阵。Among them, m k and m k+1 are the shear wave velocity models of the current iteration and the next iteration, respectively, and H k is the approximate matrix of the inverse of the Hessian matrix.

步骤8中多尺度反演方法如下:The multi-scale inversion method in step 8 is as follows:

运用维纳滤波对VSP PS转换波记录进行分频处理,从低频到高频逐步进行反演。维纳滤波器为:Wiener filtering is used to divide the VSP PS converted wave records, and the inversion is carried out step by step from low frequency to high frequency. The Wiener filter is:

其中,fwiener为wiener滤波器,Woriginal是原始主频子波,Wtarget目标主频子波,ω为角频率,ε为稳定系数,*表示复共轭。Among them, f wiener is the wiener filter, W original is the original main frequency wavelet, W target target main frequency wavelet, ω is the angular frequency, ε is the stability coefficient, and * indicates the complex conjugate.

大尺度的反演结果作为相邻小尺度反演的初始模型,重复步骤(3)到步骤(7),直到满足收敛条件。The large-scale inversion result is used as the initial model of the adjacent small-scale inversion, and steps (3) to (7) are repeated until the convergence condition is satisfied.

本发明的有益效果是通过纵横波分离策略来提高VSP地震资料转换波全波形反演的精度。发明的目的是弥补地面观测系统照明限制的缺点,提高横波速度的建模精度,特别是盐丘侧翼和陡倾角目的层的反演效果。The beneficial effect of the invention is to improve the precision of full waveform inversion of converted wave of VSP seismic data through the separation strategy of longitudinal and transverse waves. The purpose of the invention is to make up for the shortcomings of the lighting limitation of the ground observation system, improve the modeling accuracy of the shear wave velocity, especially the inversion effect of the flank of the salt dome and the target layer with a steep dip angle.

附图说明Description of drawings

图1VSP地震资料转换波全波形反演方法流程图;Fig. 1 Flowchart of converted wave full waveform inversion method of VSP seismic data;

图2复杂构造模型;Figure 2 Complex structure model;

图3第一尺度第一次迭代的横波速度梯度;Fig. 3 The shear wave velocity gradient of the first iteration at the first scale;

图4横波速度反演结果;Figure 4 S-wave velocity inversion results;

图5井旁纵侧线横波速度曲线对比图;Fig. 5 Comparison diagram of shear wave velocity curves of the longitudinal sideline beside the well;

图6转换横波目标泛函收敛曲线。Fig. 6 Convergence curve of converted shear wave target functional.

具体实施方式Detailed ways

下面结合具体实施方式对本发明进行详细说明。The present invention will be described in detail below in combination with specific embodiments.

如图1所示,为本发明实施VSP地震资料转换波全波形反演的流程图,具体包括:As shown in Figure 1, it is a flow chart of implementing VSP seismic data converted wave full waveform inversion for the present invention, specifically comprising:

(1)原始VSP地震资料PS转换波分离;(1) PS converted wave separation of original VSP seismic data;

通过全局优化参数反演法分离上行PS转换波;Separation of uplink PS converted waves by global optimization parameter inversion method;

(2)横波速度初始模型建立;(2) Establishment of the initial model of shear wave velocity;

输入纵波速度模型,根据经验公式获得横波速度模型,并对其进行平滑处理得到初始横波速度模型;Input the longitudinal wave velocity model, obtain the shear wave velocity model according to the empirical formula, and smooth it to obtain the initial shear wave velocity model;

(3)震源波场正向延拓,分离纵横波波场;(3) The source wave field is extended forward to separate the longitudinal and transverse wave fields;

通过求解纵横波分离的弹性波方程得到每一时间步的正传纵横波波场以及模拟的VSP PS波记录。The forward P-S wave field and the simulated VSP PS-wave record are obtained at each time step by solving the elastic wave equation for the separation of P- and S-waves.

(4)转换波残差反向延拓,分离纵横波波场;(4) Reverse continuation of the converted wave residual to separate the longitudinal and transverse wave fields;

通过求解纵横波分离的伴随弹性波方程得到每一时间步的反传纵横波波场;By solving the adjoint elastic wave equations of the separation of the longitudinal and transverse waves, the backpropagating longitudinal and transverse wave fields are obtained at each time step;

(5)计算横波速度梯度。(5) Calculate the shear wave velocity gradient.

采用第(3)步和第(4)步求得的正传纵波波场和反传横波波场计算横波速度梯度;Calculate the shear wave velocity gradient using the forward longitudinal wave field and the reverse propagating shear wave field obtained in steps (3) and (4);

(6)求取迭代步长;(6) Find the iteration step size;

采用三点抛物线拟合法求取迭代步长;Use the three-point parabola fitting method to find the iteration step size;

(7)更新横波速度模型;(7) Update the shear wave velocity model;

不满足收敛条件时,重复步骤(3)到步骤(6),满足收敛条件时,迭代结束;When the convergence condition is not met, repeat step (3) to step (6), when the convergence condition is met, the iteration ends;

(8)多尺度反演;(8) Multi-scale inversion;

运用维纳滤波对VSP PS转换波记录进行分频处理,从低频到高频逐步进行反演。Wiener filtering is used to divide the VSP PS converted wave records, and the inversion is carried out step by step from low frequency to high frequency.

步骤(1)中原始VSP地震资料PS转换波分离方法如下:In step (1), the PS converted wave separation method of the original VSP seismic data is as follows:

本发明中采用全局优化参数反演法分离PS转换波。In the present invention, a global optimization parameter inversion method is used to separate PS converted waves.

对每个深度点,建立观测数据与模拟记录的极小化目标泛函:For each depth point, the minimized target functional of observed data and simulated records is established:

其中,ΨPdownPupSdown和ΨSup分别为下行纵波、上行纵波、下行横波和上行横波的入射角,vP和vS为纵横波速度,u和uobs为模拟和观测波场,ω为角频率,[ω11]为频率范围。Among them, Ψ Pdown , Ψ Pup , Ψ Sdown and Ψ Sup are the incident angles of downgoing P-wave, upgoing P-wave, downgoing S-wave and upgoing S-wave, respectively; v P and v S are the P- and S-wave velocities; u and u obs are the simulated and observed wave fields , ω is the angular frequency, [ω 11 ] is the frequency range.

通过全局优化遗传算法迭代可以得到最优的传播速度和传播角度。求解方程(2)实现下行纵波、上行纵波、下行横波和上行横波分离。The optimal propagation speed and propagation angle can be obtained through the iteration of the global optimization genetic algorithm. Solving equation (2) realizes the separation of downgoing longitudinal wave, upgoing longitudinal wave, downgoing shear wave and upgoing shear wave.

u(ω)=[u(-M1,ω)…u(0,ω)…u(M2,ω)]T u(ω)=[u(-M 1 ,ω)...u(0,ω)...u(M 2 ,ω)] T

=aPdownuPdown(ω)+aPupuPup(ω)+aSdownuSdown(ω)+aSupuSup(ω), (2)=a Pdown u Pdown (ω)+a Pup u Pup (ω)+a Sdown u Sdown (ω)+a Sup u Sup (ω), (2)

其中,[uPdown,uPup,uSdown,uSup]分别为下行纵波、上行纵波、下行横波和上行横波波场,[aPdown,aPup,aSdown,aup]为相应的偏振矢量,[-M1,…,0,…,M2]为空间窗内检波点的编号。Among them, [u Pdown , u Pup , u Sdown , u Sup ] are respectively the wave fields of downgoing longitudinal wave, upgoing longitudinal wave, downgoing shear wave and upgoing shear wave, and [a Pdown , a Pup , a Sdown , a up ] are corresponding polarization vectors, [-M 1 ,…,0,…,M 2 ] are the numbers of the receiver points in the spatial window.

当震源为纵波源时,分离出的上行横波uSup(频率域)进行傅里叶反变换就得到了观测PS转换波记录(根据传播角度分解为水平分量和垂直分量)。When the seismic source is a longitudinal wave source, the separated upward shear wave u Sup (frequency domain) is inversely Fourier transformed to obtain the observed PS converted wave record (decomposed into horizontal and vertical components according to the propagation angle).

步骤(2)中横波速度初始模型建立方法如下:The establishment method of the initial model of shear wave velocity in step (2) is as follows:

纵横波速度经验公式:The empirical formula for compressional and shear wave velocity:

根据输入的纵波速度模型,采用经验公式获得横波速度模型,通过高斯窗对其进行平滑得到初始横波速度模型。According to the input P-wave velocity model, the empirical formula is used to obtain the shear-wave velocity model, which is smoothed by Gaussian window to obtain the initial shear-wave velocity model.

步骤(3)中,震源波场正向延拓,分离纵横波波场方法如下:In step (3), the source wave field is extended forward, and the method of separating the longitudinal and transverse wave fields is as follows:

一阶弹性波速度-应力方程:The first-order elastic wave velocity-stress equation:

其中,λ和μ为拉梅常数,ρ为密度,(vx,vz)为质点振动速度矢量,(τxxzzxz)为应力矢量。Among them, λ and μ are Lame constants, ρ is density, (v x , v z ) is particle vibration velocity vector, (τ xxzzxz ) is stress vector.

本发明中采用纵横波分离的速度-应力弹性波方程:Adopt the velocity-stress elastic wave equation of longitudinal and transverse wave separation in the present invention:

其中,(vxp,vzp)为纵波质点振动速度矢量,(vxs,vzs)为横波质点振动速度矢量,(τxxpzzp)为纵波应力矢量,(τxxszzsxzs)为横波应力矢量。Among them, (v xp , v zp ) is the particle vibration velocity vector of longitudinal wave, (v xs , v zs ) is the particle vibration velocity vector of shear wave, (τ xxpzzp ) is the stress vector of longitudinal wave, (τ xxszzsxzs ) is the shear wave stress vector.

采用高阶交错网格有限差分数值求解方程(6)得到每一时间步的正传纵横波波场以及模拟VSP PS波记录。The high-order staggered grid finite difference method is used to numerically solve Equation (6) to obtain the forward P- and S-wave field and simulated VSP PS wave records at each time step.

步骤(4)转换波残差反向延拓,分离纵横波波场方法如下:Step (4) Reverse continuation of the converted wave residual, and the method of separating the longitudinal and transverse wave fields is as follows:

在伴随状态法中,伴随算子满足:In the adjoint state method, the adjoint operator satisfies:

<Lu,u*>=<u,L*u*>, (7)<Lu,u * >=<u,L * u * >, (7)

其中:L为正演算子,u为质点速度和应力矢量,“*”为其伴随状态。将方程(4)带入方程(7)得到常规伴随方程:Among them: L is the forward operator, u is the particle velocity and stress vector, "*" is its adjoint state. Substituting Equation (4) into Equation (7) yields the conventional adjoint equation:

同理,基于方程(6)可以得到纵横波分离的伴随方程:Similarly, based on equation (6), the adjoint equation for the separation of P- and S-waves can be obtained:

其中,“*”表示相应波场的伴随波场,为纵横波耦合的观测和模拟VSP地震记录,为分离的观测和模拟VSP转换波地震记录。Among them, "*" indicates the accompanying wave field of the corresponding wave field, For the observation and simulation of VSP seismic records coupled with compressional and shear waves, Separated observations and simulated VSP converted wave seismic records.

以模拟和观测VSP PS波记录残差为震源,采用高阶交错有限差分数值求解方程(10)和(11)得到每一时间步的反传纵横波波场。Taking the simulated and observed VSP PS wave recording residuals as the seismic source, equations (10) and (11) are numerically solved using high-order staggered finite differences to obtain the backpropagating P-S wave field at each time step.

步骤(5)计算横波速度梯度方法如下:Step (5) The method of calculating the shear wave velocity gradient is as follows:

传统全波形反演中横波速度梯度为:The shear wave velocity gradient in traditional full waveform inversion is:

其中,E为目标泛函。Among them, E is the target functional.

本发明采用正传纵波波场和反传横波波场计算横波速度梯度:The present invention calculates the shear wave velocity gradient by using the forward longitudinal wave field and the backpropagation shear wave field:

其中,in,

通过采用新的梯度公式(方程13)可以缓解纵波速度的干扰,进而提高横波速度的反演精度。By adopting the new gradient formula (Equation 13), the interference of P-wave velocity can be alleviated, and the inversion accuracy of S-wave velocity can be improved.

步骤(6)求取迭代步长方法如下:Step (6) The method of calculating the iteration step size is as follows:

本发明中采用三点抛物线插值方法计算步长。选取三个试探点,用这三个点构建一条抛物线,这条抛物线的极值点既为最优步长。In the present invention, a three-point parabolic interpolation method is used to calculate the step size. Select three trial points, and use these three points to construct a parabola. The extreme point of this parabola is the optimal step size.

三个试探步长a0,a1,a2对应的目标泛函为Es0,Es1,Es2,则最优步长为:The target functionals corresponding to the three trial steps a 0 , a 1 , and a 2 are E s0 , E s1 , E s2 , and the optimal step size is:

步骤(7)更新横波速度模型方法如下:Step (7) The method of updating the shear wave velocity model is as follows:

通过L-BFGS方法更新横波速度:Update the shear wave velocity by the L-BFGS method:

其中,mk和mk+1分别为当前迭代和下一次迭代的横波速度模型,Hk为海森矩阵逆的近似矩阵。Among them, m k and m k+1 are the shear wave velocity models of the current iteration and the next iteration, respectively, and H k is the approximate matrix of the inverse of the Hessian matrix.

步骤(8)中多尺度反演方法如下:The multi-scale inversion method in step (8) is as follows:

运用维纳滤波对VSP PS转换波记录进行分频处理,从低频到高频逐步进行反演。维纳滤波器为:Wiener filtering is used to divide the VSP PS converted wave records, and the inversion is carried out step by step from low frequency to high frequency. The Wiener filter is:

其中,fwiener为wiener滤波器,Woriginal是原始主频子波,Wtarget目标主频子波,ω为角频率,ε为稳定系数,*表示复共轭。Among them, f wiener is the wiener filter, W original is the original main frequency wavelet, W target target main frequency wavelet, ω is the angular frequency, ε is the stability coefficient, and * indicates the complex conjugate.

大尺度的反演结果作为相邻小尺度反演的初始模型,重复步骤(3)到步骤(7),直到满足收敛条件。The large-scale inversion result is used as the initial model of the adjacent small-scale inversion, and steps (3) to (7) are repeated until the convergence condition is met.

本发明由于采取以上技术方案,其具有以下优点:1.可以提高横波速度模型的建模精度。2.可以改善复杂构造(如盐丘侧翼和陡倾角目的层)的反演精度。Due to the adoption of the above technical solutions, the present invention has the following advantages: 1. The modeling accuracy of the shear wave velocity model can be improved. 2. It can improve the inversion accuracy of complex structures (such as salt dome flanks and steep dip angle target layers).

下面通过一个例子来分析本发明中提出的VSP地震资料转换波全波形反演方法的有效性。Next, an example is used to analyze the effectiveness of the converted wave full waveform inversion method of VSP seismic data proposed in the present invention.

图2中(a)-(c)分别为某复杂构造模型的真实纵波速度,真实横波速度和初始横波速度模型。计算区域为200×200个网格点,空间间隔为10m,时间步长为1ms,总记录时间为1.7s。震源为15Hz的雷克子波,加在正应力上(纵波源)。在模型右下方有一高速体,且纵波速度和横波速度在该高速体处形态不同。井位于x=1980m处,井中布置180个检波器,从z=100m开始每隔10米一个。炮点位于地表,从x=150m处每隔70m一炮,共26炮。采用时间二阶、空间十阶精度有限差分法进行数值模拟,PML吸收边界消除边界反射。采用二个尺度进行多尺度反演,第一尺度为15Hz,第二尺度为30Hz。图3为第一尺度第一次迭代的横波速度梯度。可以看到,井旁浅层的构造清晰,深层构造,和盐丘侧翼也得到准确的刻画。图4中(a)、(b)分别为第一尺度和第二尺度横波速度的反演结果。由图可见,本发明提出的VSP地震资料转换波全波形反演方法的反演精度较高,井旁及盐丘侧翼构造清晰可见,高速异常体也得到较好恢复。图5为x=1650m处真实横波速度、初始速度和反演结果对比曲线图。反演结果和真实结果较好吻合。图6给出目标函数Es的收敛曲线,进一步表明本发明VSP地震资料转换波全波形反演方法的有效性。(a)-(c) in Fig. 2 are the true P-wave velocity, true S-wave velocity and initial S-wave velocity model of a complex structural model, respectively. The calculation area is 200×200 grid points, the spatial interval is 10m, the time step is 1ms, and the total recording time is 1.7s. The source is the 15Hz Reich wavelet, which is added to the normal stress (longitudinal wave source). There is a high-speed body at the lower right of the model, and the shape of the longitudinal wave velocity and the shear wave velocity are different at the high-speed body. The well is located at x=1980m, and 180 geophones are arranged in the well, one every 10 meters from z=100m. The shot point is located on the surface, and a shot is fired every 70m from x=150m, a total of 26 shots. The second-order time and tenth-order precision finite difference method is used for numerical simulation, and the PML absorbs the boundary to eliminate the boundary reflection. Two scales are used for multi-scale inversion, the first scale is 15Hz and the second scale is 30Hz. Figure 3 shows the shear wave velocity gradient for the first iteration of the first scale. It can be seen that the shallow structure beside the well is clear, and the deep structure and the flank of the salt dome are also accurately depicted. (a) and (b) in Fig. 4 are the inversion results of the first-scale and second-scale shear wave velocities, respectively. It can be seen from the figure that the VSP seismic data converted wave full waveform inversion method proposed by the present invention has high inversion accuracy, the wellside and salt dome flank structures are clearly visible, and the high-speed anomaly body is also well recovered. Fig. 5 is a comparative graph of true shear wave velocity, initial velocity and inversion results at x=1650m. The inversion results are in good agreement with the real results. Fig. 6 shows the convergence curve of the objective function Es, which further shows the validity of the converted wave full waveform inversion method for VSP seismic data of the present invention.

本发明是一种新的VSP地震资料转换波全波形反演方法,能够获取高精度的横波速度模型;能大大提高复杂地质构造(如盐丘侧翼和陡倾界面)的反演精度。The invention is a new converted wave full waveform inversion method of VSP seismic data, which can obtain a high-precision shear wave velocity model; and can greatly improve the inversion accuracy of complex geological structures (such as salt dome flanks and steep interfaces).

以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。The above description is only a preferred embodiment of the present invention, and does not limit the present invention in any form. Any simple modifications made to the above embodiments according to the technical essence of the present invention, equivalent changes and modifications, all belong to this invention. within the scope of the technical solution of the invention.

Claims (9)

1. A full waveform inversion method for VSP seismic data converted waves is characterized by comprising the following steps:
(1) separating PS converted waves of original VSP seismic data;
separating the uplink PS converted wave by a global optimization parameter inversion method;
(2) establishing a transverse wave velocity initial model;
inputting a longitudinal wave velocity model to obtain a transverse wave velocity model, and smoothing the transverse wave velocity model to obtain an initial transverse wave velocity model;
(3) carrying out forward continuation on a seismic source wave field, and separating a longitudinal and transverse wave field;
obtaining a positive transmission longitudinal and transverse wave field of each time step and a simulated VSPPS wave record by solving an elastic wave equation of longitudinal and transverse wave separation;
(4) converting wave residual error reverse continuation, and separating longitudinal and transverse wave fields;
obtaining a reverse transmission longitudinal and transverse wave field of each time step by solving an adjoint elastic wave equation of longitudinal and transverse wave separation;
(5) calculating the velocity gradient of the transverse wave;
calculating the shear wave velocity gradient by adopting the forward-transmission longitudinal wave field obtained in the step (3) and the backward-transmission shear wave field obtained in the step (4);
(6) calculating an iteration step length;
solving an iteration step length by adopting a three-point parabola fitting method;
(7) updating the transverse wave velocity model;
when the convergence condition is not met, repeating the steps (3) to (6), and when the convergence condition is met, finishing the iteration;
(8) performing multi-scale inversion;
and performing frequency division processing on the VSP PS converted wave record by using wiener filtering, and gradually performing inversion from low frequency to high frequency.
2. The VSP seismic data converted waveform full waveform inversion method of claim 1, wherein: the PS conversion wave separation method of the original VSP seismic data in the step (1) is as follows:
separating PS converted waves by adopting a global optimization parameter inversion method, and establishing a minimized target functional of observation data and simulation record for each depth point:
therein, ΨPdownPupSdownAnd ΨSupIncident angles v of down longitudinal wave, up longitudinal wave, down transverse wave and up transverse wave respectivelyPAnd vSIs the longitudinal and transverse wave velocity, u and uobsFor simulating and observing the wavefield, ω is the angular frequency, [ ω ]11]For the frequency range, the optimal propagation speed and propagation angle can be obtained through the iteration of a global optimization genetic algorithm; solving an equation (2) to realize separation of downlink longitudinal waves, uplink longitudinal waves, downlink transverse waves and uplink transverse waves;
u(ω)=[u(-M1,ω)…u(0,ω)…u(M2,ω)]T
=aPdownuPdown(ω)+aPupuPup(ω)+aSdownuSdown(ω)+aSupuSup(ω), (2)
wherein [ u ]Pdown,uPup,uSdown,uSup]Respectively, a down-going longitudinal wave, an up-going longitudinal wave, a down-going transverse wave and an up-going transverse wave, [ a ]Pdown,aPup,aSdown,aup]For the corresponding polarization vector, [ -M1,…,0,…,M2]The number of the detection point in the space window is the number, when the seismic source is a longitudinal wave source, the separated up-going transverse wave uSupAnd performing Fourier inversion to obtain an observation PS converted wave record.
3. The VSP seismic data converted waveform full waveform inversion method of claim 1, wherein: the method for establishing the shear wave velocity initial model in the step (2) is as follows:
longitudinal and transverse wave velocity empirical formula:
and obtaining a transverse wave velocity model by adopting an empirical formula according to the input longitudinal wave velocity model, and smoothing the transverse wave velocity model through a Gaussian window to obtain an initial transverse wave velocity model.
4. The VSP seismic data converted waveform full waveform inversion method of claim 1, wherein: the method for separating the longitudinal wave field from the transverse wave field in the step (3) comprises the following steps:
first order elastic wave velocity-stress equation:
wherein λ and μ are Lame constants, ρ is density, (v)x,vz) Is the particle vibration velocity vector (tau)xxzzxz) Is a stress vector;
adopting a speed-stress elastic wave equation of longitudinal and transverse wave separation:
wherein (v)xp,vzp) Is the vibration velocity vector of longitudinal wave particle, (v)xs,vzs) Is the transverse wave particle vibration velocity vector (tau)xxpzzp) Is a longitudinal wave stress vector (tau)xxszzsxzs) Is the transverse wave stress vector;
and (3) solving the equation (6) by adopting a high-order staggered grid finite difference numerical value to obtain a positive transmission longitudinal and transverse wave field of each time step and simulating VSP PS wave record.
5. The VSP seismic data converted waveform full waveform inversion method of claim 1, wherein: the converted wave residual error reverse continuation in the step (4) and the method for separating the longitudinal wave field and the transverse wave field are as follows:
in the adjoint state method, the adjoint operator satisfies:
<Lu,u*>=<u,L*u*>, (7)
wherein: l is the positive operator, u is the particle velocity and stress vector, "# is its companion state, substituting equation (4) into equation (7) yields the conventional companion equation:
obtaining a concomitant equation of longitudinal and transverse wave separation based on equation (6):
wherein, represents the adjoint wave field of the corresponding wave field,for observation of the coupling of longitudinal and transverse waves and for simulation of VSP seismic recordings,for separated observation and simulation VSP converted wave seismic records, simulating and observing VSP PS wave record residual errors are used as seismic sources, and equations (10) and (11) are solved by adopting high-order interleaved finite difference numerical values to obtain a reverse propagation longitudinal and transverse wave field of each time step.
6. The VSP seismic data converted waveform full waveform inversion method of claim 1, wherein: the method for calculating the shear wave velocity gradient in the step (5) comprises the following steps:
the shear wave velocity gradient in the conventional full waveform inversion is as follows:
wherein E is a target functional, and the shear wave velocity gradient is calculated by adopting a forward-transmission longitudinal wave field and a backward-transmission shear wave field:
wherein,
the new gradient formula is adopted to relieve the interference of longitudinal wave velocity, and further the inversion accuracy of transverse wave velocity is improved.
7. The VSP seismic data converted waveform full waveform inversion method of claim 1, wherein: the iterative step length calculating method in the step (6) is as follows:
selecting three trial points, and constructing a parabola by using the three points, wherein the extreme point of the parabola is the optimal step length, and the three trial step lengths a0,a1,a2Corresponding target functional is Es0,Es1,Es2Then the optimal step length is:
8. the VSP seismic data converted waveform full waveform inversion method of claim 1, wherein: the method for updating the shear wave velocity model in the step (7) is as follows:
updating the transverse wave velocity by an L-BFGS method:
wherein m iskAnd mk+1Shear wave velocity models, H, for the current iteration and the next iteration, respectivelykAn approximation matrix that is the inverse of the hessian matrix.
9. The VSP seismic data converted waveform full waveform inversion method of claim 1, wherein: the multi-scale inversion method in the step (8) is as follows:
performing frequency division processing on the VSP PS converted wave record by using wiener filtering, and gradually performing inversion from low frequency to high frequency, wherein the wiener filter is as follows:
wherein f iswienerIs a wiener filter, WoriginalIs a primary dominant frequency wavelet, WtargetAnd (3) repeating the steps (3) to (7) until a convergence condition is met, wherein omega is angular frequency, epsilon is a stable coefficient, represents complex conjugate, and the inversion result with a large scale is used as an initial model of adjacent small-scale inversion.
CN201810666468.2A 2018-06-26 2018-06-26 A kind of VSP seismic data converted wave full waveform inversion method Pending CN108845351A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810666468.2A CN108845351A (en) 2018-06-26 2018-06-26 A kind of VSP seismic data converted wave full waveform inversion method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810666468.2A CN108845351A (en) 2018-06-26 2018-06-26 A kind of VSP seismic data converted wave full waveform inversion method

Publications (1)

Publication Number Publication Date
CN108845351A true CN108845351A (en) 2018-11-20

Family

ID=64201961

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810666468.2A Pending CN108845351A (en) 2018-06-26 2018-06-26 A kind of VSP seismic data converted wave full waveform inversion method

Country Status (1)

Country Link
CN (1) CN108845351A (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109459787A (en) * 2018-10-09 2019-03-12 中国地质大学(武汉) Underground coal mine structure imaging method and system based on SEISMIC CHANNEL WAVE METHOD full waveform inversion
CN109459789A (en) * 2018-12-18 2019-03-12 吉林大学 Time-domain full waveform inversion method based on amplitude decaying and linear interpolation
CN111751876A (en) * 2020-07-07 2020-10-09 中油奥博(成都)科技有限公司 A method and device for variable offset VSP converted shear wave one-way wave prestack depth migration
EP3761067A1 (en) * 2019-07-01 2021-01-06 CGG Services SAS Full waveform inversion approach to building an s-wave velocity model using ps data
CN112773396A (en) * 2021-01-13 2021-05-11 佟小龙 Medical imaging method based on full waveform inversion, computer equipment and storage medium
CN113156498A (en) * 2021-02-26 2021-07-23 中海石油(中国)有限公司 Pre-stack AVO three-parameter inversion method and system based on homotopy continuation
CN113376695A (en) * 2021-06-11 2021-09-10 中国矿业大学 Full waveform inversion method suitable for complex collapse column of coal seam floor
CN113589366A (en) * 2020-04-30 2021-11-02 中国石油化工股份有限公司 Broadband fusion modeling method based on full waveform inversion
CN113987888A (en) * 2021-11-16 2022-01-28 华中科技大学 A Numerical Separation Method of Bedrock Incident-Down Waves Based on Array Observation
CN114660659A (en) * 2022-03-30 2022-06-24 中国矿业大学 Longitudinal and transverse wave decoupling high-precision microseism positioning method
CN115166823A (en) * 2022-08-10 2022-10-11 中国石油大学(华东) Seismic longitudinal and transverse wave data separation method, model construction method and system
CN115436995A (en) * 2022-08-09 2022-12-06 中煤科工集团西安研究院有限公司 Multi-source mixed seismic wave field separation method based on first-arrival wave field constraint
CN116540297A (en) * 2023-05-06 2023-08-04 中国科学院地质与地球物理研究所 Method, system and equipment for full waveform inversion of elastic wave seismic data

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103499835A (en) * 2013-10-13 2014-01-08 中国石油集团西北地质研究所 Method for inverting near-surface velocity model by utilizing preliminary waveforms
CN106842295A (en) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 The waveform inversion method of logging information constrained
WO2018035370A1 (en) * 2016-08-19 2018-02-22 Halliburton Energy Services, Inc. Full waveform inversion of vertical seismic profile data for anisotropic velocities using pseudo-acoustic wave equations

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103499835A (en) * 2013-10-13 2014-01-08 中国石油集团西北地质研究所 Method for inverting near-surface velocity model by utilizing preliminary waveforms
CN106842295A (en) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 The waveform inversion method of logging information constrained
WO2018035370A1 (en) * 2016-08-19 2018-02-22 Halliburton Energy Services, Inc. Full waveform inversion of vertical seismic profile data for anisotropic velocities using pseudo-acoustic wave equations

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
LEONARDO GO ´MEZ ET AL.: "Wavefield separation methods for gradient filtering in time-domain full-waveform inversion", 《FIFTEENTH INTERNATIONAL CONGRESS OF THE BRAZILIAN GEOPHYSICAL SOCIETY》 *
MIN LOU ET AL.: "An improved parametric inversion methodology to separate P and Sv wavefields from VSP data", 《SEG HOUSTON 2013 ANNUAL MEETING》 *
ZHIMING REN ET AL.: "A hierarchical elastic full-waveform inversion scheme based on wavefield separation and the multistep-length approach", 《GEOPHYSICS》 *
任志明: "声波和弹性波波动方程有限差分正反演方法研究", 《中国博士学位论文全文数据库 基础科学辑》 *
孙思宇: "时间域弹性波全波形反演", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
郭振波: "弹性介质波形反演方法研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109459787B (en) * 2018-10-09 2019-12-06 中国地质大学(武汉) Coal Mine Underground Structural Imaging Method and System Based on Seismic Channel Wave Full Waveform Inversion
CN109459787A (en) * 2018-10-09 2019-03-12 中国地质大学(武汉) Underground coal mine structure imaging method and system based on SEISMIC CHANNEL WAVE METHOD full waveform inversion
CN109459789A (en) * 2018-12-18 2019-03-12 吉林大学 Time-domain full waveform inversion method based on amplitude decaying and linear interpolation
EP3761067A1 (en) * 2019-07-01 2021-01-06 CGG Services SAS Full waveform inversion approach to building an s-wave velocity model using ps data
US11215720B2 (en) 2019-07-01 2022-01-04 Cgg Services Sas Full waveform inversion approach to building an S-wave velocity model using PS data
CN113589366A (en) * 2020-04-30 2021-11-02 中国石油化工股份有限公司 Broadband fusion modeling method based on full waveform inversion
CN113589366B (en) * 2020-04-30 2023-10-20 中国石油化工股份有限公司 Broadband fusion modeling method based on full waveform inversion
CN111751876A (en) * 2020-07-07 2020-10-09 中油奥博(成都)科技有限公司 A method and device for variable offset VSP converted shear wave one-way wave prestack depth migration
CN111751876B (en) * 2020-07-07 2022-05-20 中油奥博(成都)科技有限公司 A method and device for variable offset VSP converted shear wave one-way wave prestack depth migration
CN112773396A (en) * 2021-01-13 2021-05-11 佟小龙 Medical imaging method based on full waveform inversion, computer equipment and storage medium
CN112773396B (en) * 2021-01-13 2023-06-16 佟小龙 Medical imaging method based on full waveform inversion, computer equipment and storage medium
CN113156498A (en) * 2021-02-26 2021-07-23 中海石油(中国)有限公司 Pre-stack AVO three-parameter inversion method and system based on homotopy continuation
CN113156498B (en) * 2021-02-26 2024-01-26 中海石油(中国)有限公司 Pre-stack AVO three-parameter inversion method and system based on homotopy continuation
CN113376695A (en) * 2021-06-11 2021-09-10 中国矿业大学 Full waveform inversion method suitable for complex collapse column of coal seam floor
CN113376695B (en) * 2021-06-11 2022-07-05 中国矿业大学 A full-waveform inversion method suitable for complex collapse columns of coal seam floor
CN113987888A (en) * 2021-11-16 2022-01-28 华中科技大学 A Numerical Separation Method of Bedrock Incident-Down Waves Based on Array Observation
CN113987888B (en) * 2021-11-16 2024-07-09 华中科技大学 Bedrock incidence-downlink wave numerical value separation method based on array observation
CN114660659A (en) * 2022-03-30 2022-06-24 中国矿业大学 Longitudinal and transverse wave decoupling high-precision microseism positioning method
CN115436995A (en) * 2022-08-09 2022-12-06 中煤科工集团西安研究院有限公司 Multi-source mixed seismic wave field separation method based on first-arrival wave field constraint
CN115166823A (en) * 2022-08-10 2022-10-11 中国石油大学(华东) Seismic longitudinal and transverse wave data separation method, model construction method and system
CN116540297A (en) * 2023-05-06 2023-08-04 中国科学院地质与地球物理研究所 Method, system and equipment for full waveform inversion of elastic wave seismic data
CN116540297B (en) * 2023-05-06 2024-03-08 中国科学院地质与地球物理研究所 Full waveform inversion method, system and equipment for elastic wave seismic data

Similar Documents

Publication Publication Date Title
CN108845351A (en) A kind of VSP seismic data converted wave full waveform inversion method
CN108873066B (en) Elastic medium wave equation reflected wave travel time inversion method
CN108037526B (en) Reverse-time migration method based on all-wave wave field VSP/RVSP seismic data
CN102508293B (en) A Pre-stack Inversion Method for Identifying Oil and Gas in Thin Layers
CN110187382B (en) Traveling time inversion method for wave equation of reverse wave and reflected wave
WO2011146161A1 (en) Method and system for checkpointing during simulations
CN107462924B (en) A kind of absolute wave impedance inversion method independent of well-log information
CN106842295A (en) The waveform inversion method of logging information constrained
CN101201409B (en) Method for revising earthquake data phase
CN113740901A (en) Land seismic data full-waveform inversion method and apparatus based on complex undulating surface
CN111025387A (en) Pre-stack earthquake multi-parameter inversion method for shale reservoir
CN110579795A (en) Joint Velocity Inversion Method Based on Passive Source Seismic Waveform and Its Reverse Time Imaging
CN107229066B (en) VSP data full waveform inversion modeling method based on surface seismic structure constraint
CN110737018B (en) Method for modeling anisotropy of VSP seismic data
Qu et al. Joint acoustic and decoupled-elastic least-squares reverse time migration for simultaneously using water-land dual-detector data
CN111025388A (en) Multi-wave combined prestack waveform inversion method
Li et al. Waveform inversion of seismic first arrivals acquired on irregular surface
CN115236730A (en) A seismic wavefield migration imaging method based on interlayer multiple Fourier finite difference
CN113536638B (en) High-precision seismic wave field simulation method based on discontinuous finite element and staggered grid
Jin et al. 2D multiscale non‐linear velocity inversion
CN108680957B (en) Local cross-correlation time-frequency domain Phase-retrieval method based on weighting
CN102466818A (en) Method for imaging anisotropic medium by using borehole seismic data
Mittet et al. Imaging of offset VSP data with an elastic iterative migration scheme
CN111175822B (en) Improved Direct Envelope Inversion and Perturbation Decomposition Inversion Method for Strong Scattering Media
CN109901221B (en) Seismic data anisotropy modeling method based on dynamic correction velocity parameter

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20181120