CN103499835A - 初至波波形反演近地表速度模型方法 - Google Patents

初至波波形反演近地表速度模型方法 Download PDF

Info

Publication number
CN103499835A
CN103499835A CN201310476248.0A CN201310476248A CN103499835A CN 103499835 A CN103499835 A CN 103499835A CN 201310476248 A CN201310476248 A CN 201310476248A CN 103499835 A CN103499835 A CN 103499835A
Authority
CN
China
Prior art keywords
wave field
wave
utilize
equation
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
CN201310476248.0A
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.)
Institute Of Northwest Geology Of China Petroleum Group
Original Assignee
Institute Of Northwest Geology Of China Petroleum Group
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 Institute Of Northwest Geology Of China Petroleum Group filed Critical Institute Of Northwest Geology Of China Petroleum Group
Priority to CN201310476248.0A priority Critical patent/CN103499835A/zh
Publication of CN103499835A publication Critical patent/CN103499835A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了初至波波形反演近地表速度模型方法,该方法包括声波方程波场正演和最速下降法波形反演技术,其方法步骤为:①提取时间域初至波波形记录和初始模型;②利用声波方程交错网格有限差分正演模拟计算模拟波场及波场残差;③将波场残差反向传播得到回传波场;④利用回传波场和正向传播波场计算目标函数的梯度,并计算更新步长;⑤对速度模型进行更新;⑥检验是否满足迭代终止条件,满足则输出速度模型,否则返回②步骤,继续迭代更新;本发明借鉴了基于波动方程理论的全波形反演技术,利用了能量较强、波形比较稳定的初至波进行反演,降低了全波形反演的多解性,提高了反演的稳定性和计算效率;提高了静校正和浅层深度成像的精度。

Description

初至波波形反演近地表速度模型方法
技术领域
本发明属于石油地震勘探地震资料处理技术领域,旨在提高地震资料处理中的近地表速度模型精度,进而提高静校正精度和浅层深度成像精度,具体地说是一种基于初至波波形(有别于传统的初至走时)反演,建立高精度近地表速度模型的新技术、新方法。
背景技术
目前反演近地表速度模型,绝大部分是基于射线理论,利用初至波旅行时间反演近地表速度模型,如折射静校正、层析反演等等,由于射线理论是对地震波场的高频近似,只能反演出近地表速度模型的低波数成分(模型尺度远大于地震波长),反演的速度模型分辨率受到限制。基于波动理论的全波形反演方法,利用全波场信息,反演地下速度模型,虽理论比较完善,但受限于数据的信噪比、反演方法本身的多解性和计算效率等问题,实际应用受到很大限制。
为了提高近地表速度模型的反演精度,借鉴全波形反演方法的技术路线,研究了利用初至波波形(旅行时间、振幅)信息反演近地表速度模型方法。
发明内容
鉴于已有技术的不足,本发明的目的是提供一种基于初至波波形反演,建立高精度近地表速度模型的新技术、新方法;该方法基于波动方程理论,和传统的基于射线理论的走时层析反演相比,能够反演出近地表速度模型的高波数成分(模型尺度接近或小于地震波长),提高了近地表速度模型的分辨率;和全波形反演方法相比,该方法只利用能量较强、波形比较稳定的初至波,降低了全波形反演方法的多解性,提高了计算精度和计算效率;该方法具有提高静校正精度和浅层深度成像精度的特点。
实现解决技术问题的技术方案是:初至波波形反演近地表速度模型方法,该方法包括声波方程交错网格有限差分正演模拟技术和最速下降法波形反演技术,其实现过程的方法步骤如下:
(1)、提取时间域初至波波形记录Pobs和初始速度模型;
(2)、利用声波方程交错网格有限差分正演模拟计算模拟波场Pcal,并计算波场残差δp=Pobs-Pcal
(3)、将波场残差反向传播得到回传波场P′;
(4)、利用回传波场和正向传播波场计算目标函数的梯度,即速度的更新方向,并计算更新步长;
(5)、对速度模型进行更新;
(6)、检验是否满足迭代终止条件,满足则输出速度模型的反演结果,否则返回步骤(2),继续迭代更新。
所述声波方程交错网格有限差分正演模拟技术,利用一阶声波方程有限差分方法,数值模拟炮点激发的正向传播波场;非均匀介质一阶声波方程:应力-质点速度方程;时间2阶、空间12阶的高阶有限差分,PML边界条件,交错网格。
所述最速下降法波形反演技术,⑴地面观测记录加时窗函数—初至波窗函数,作为一阶声波方程的初值条件,类似于炮点激发正向传播波场正演技术,数值模拟初至波反向传播的地震波场和数值模拟的正向传播波场求残差,利用L2模建立反演目标函数;⑵利用正向传播波场和波场残差函数,计算目标函数的最速下降方向即梯度,确定出速度扰动的最大变化方向;⑶利用梯度和Frechet导数计算出速度扰动的迭代步长;⑷利用速度更新方向和更新步长更新速度。
本发明相比现有技术的显著效果是:
1)本发明基于波动方程理论而没有高频近似假设,能够反演出近地表速度模型的高波数成分(模型尺度接近或小于地震波长),提高了速度模型的分辨率,进而提高静校正的精度,并为近地表和地下深度偏移奠定了良好的基础;
2)本发明借鉴了基于波动方程理论的全波形反演技术,利用了能量较强、波形比较稳定的初至波进行反演,降低了全波形反演的多解性,提高了反演的稳定性和计算效率;
3)本发明采用时间2阶、空间12阶的高精度交错网格有限差分正演模拟方法,并利用PML(完全匹配层)边界条件消除边界效应,提高了正演模拟精度,有效地减小了人工边界反射对有效波场的影响;
4)本发明采用声波方程正演模拟技术,和传统的标量方程规则网格有限差分相比,该技术可以降低频散,减少边界效应,从而提高数值模拟的精度和稳定性;
5)本发明采用最速下降法速度反演技术,并通过正演和反演的多次迭代,直至目标函数完全收敛,就完成了初至波波形反演近地表速度模型;为了提高反演的效率和精度,初始模型可采用常规射线理论层析反演的近地表模型;
6)本发明采用波形反演的同时利用振幅和走时信息,比传统基于射线理论的旅行时层析反演方法具有更高的反演精度,能够获取较为精确的境地表速度模型,为静校正和浅层深度域成像精度创造了条件。
附图说明
图1是本发明近地表速度异常体模型的示意图
图2是本发明声波方程有限差分正演模拟结果示意图
图2a为规则网格有限差分正演模拟得到的单炮波形记录图
图2c为规则网格有限差分正演模拟得到的200ms波场快照图
图2b为本方法正演模拟得到的单炮波形记录示意图
图2d为本方法正演模拟得到的200ms波场快照示意图
图3是本发明实际速度模型示意图
图4是本发明传统旅行时层析反演结果示意图
图5是本发明初至波波形反演结果示意图
图5a为第1次迭代反演结果示意图
图5b为第100次迭代反演结果示意图
图5c为第500次迭代反演结果示意图
图6是本发明初至波波形反演目标函数梯度示意图
上述图1-图6均为计算机绘画软件suffer绘制的附图
具体实施方式
下面结合附图和实例对该技术方案作进一步描述
初至波波形反演近地表速度模型方法,该方法包括声波方程交错网格有限差分正演模拟和最速下降法波形反演两大技术,其实现过程的方法步骤如下:
⑴提取时间域初至波波形记录Pobs和初始速度模型;
⑵利用声波方程交错网格有限差分正演模拟计算模拟波场Pcal,并计算波场残差δp=Pobs-Pcal
⑶将波场残差反向传播得到回传波场P′;
⑷利用回传波场和正向传播波场计算目标函数的梯度,即速度的更新方向,并计算更新步长;
⑸对速度模型进行更新;
⑹检验是否满足迭代终止条件,满足则输出速度模型的反演结果,否则返回步骤2,继续迭代更新。
所述1、声波方程交错网格有限差分正演模拟技术
利用一阶声波方程有限差分方法,数值模拟炮点激发的正向传播波场;
技术要点:非均匀介质一阶声波方程(应力-质点速度方程,公式(1));高阶有限差分(时间2阶、空间12阶),PML边界条件,交错网格;与传统的标量方程规则网格有限差分相比,该技术可以降低频散,减少边界效应,从而提高数值模拟的精度和稳定性。
∂ u ( x , y , z , t ) ∂ t = - k ( ∂ v x ∂ x + ∂ v y ∂ y + ∂ v z ∂ z ) ∂ v x ∂ x = - 1 ρ ∂ u ( x , y , z , t ) ∂ x ∂ v y ∂ y = - 1 ρ ∂ u ( x , y , z , t ) ∂ y ∂ v z ∂ z = - 1 ρ ∂ u ( x , y , z , t ) ∂ z - - - ( 1 )
公式中,k为体变模量,ρ为密度,vx、vy、vz为质点振动速度,u为波场函数。
所述2、最速下降法波形反演技术
1)地面观测记录加时窗函数(初至波窗函数),作为一阶声波方程的初值条件,类似于炮点激发正向传播波场正演技术,数值模拟初至波反向传播的地震波场和数值模拟的正向传播波场求残差,利用L2模建立反演目标函数(公式(2));
E = 1 2 Σ s Σ g ∫ dtδ p gs 2 ( t ) δ p gs ( r → g , t ; r → s ) = p ( r → g , t ; r → s ) obs - p ( r → g , t ; r → s ) cal - - - ( 2 )
公式中,
Figure BDA0000395017190000052
为检波点位置,
Figure BDA0000395017190000053
为炮点位置,t为时间变量,
Figure BDA0000395017190000054
为实际记录波场,
Figure BDA0000395017190000055
为初始模型正演模拟计算波场,
Figure BDA0000395017190000056
为波场残差,E为波场残差对所有炮点、检波点和记录时间的积分,成为总的波场残差,为全波形反演的目标函数。
2)利用正向传播波场和波场残差函数,计算目标函数的最速下降方向(梯度),确定出速度扰动的最大变化方向(公式(3))。
γ ( r → ) = 1 c ( r → ) 3 Σ s ∫ dt p · ( r → , t ; r → s ) p · ′ ( r → , t ; r → s ) p ′ ( r → , t ; r → s ) = Σ g g ( r → , - t ; r → g , 0 ) * δp ( r → g , t ; r → s ) - - - ( 3 )
公式中,
Figure BDA0000395017190000058
为速度模型,波场函数p上的小点代表波场函数p对时间的导数。
3)利用梯度和Frechet导数计算出速度扰动的迭代步长(公式(4))。
α k = [ γ k ( r → ) ] t [ γ k ( r → ) ] [ F γ k ( r → ) ] t [ F γ k ( r → ) ] [ F γ k ( r → ) ] = p [ c k ( r → ) + ϵ γ k ( r → ) ] - p [ c k ( r → ) ] ϵ - - - ( 4 )
公式中,
Figure BDA00003950171900000510
为第k次迭代的初始速度模型,
Figure BDA00003950171900000511
为速度更新方向,ε为速度扰动步长,
Figure BDA00003950171900000512
为Frechet导数,上标t表示向量的转置,αk为速度更新步长。
4)利用速度更新方向和更新步长更新速度(公式(5))。
c k + 1 ( r → ) = c k ( r → ) + α k · γ k ( r → ) - - - ( 5 )
公式中,
Figure BDA00003950171900000514
为第k次迭代的初始速度模型,为速度更新方向,αk为速度更新步长,
Figure BDA00003950171900000516
为更新后的速度模型,将作为第k+1次迭代的初始速度模型。
至此,完成了速度模型的一次反演;通过正演和反演的多次迭代,直至目标函数完全收敛,就完成了初至波波波形反演近地表速度模型。为了提高反演的效率和精度,初始模型可以采用常规射线理论层析反演的近地表模型。
初至波波形反演近地表速度模型方法的流程框图如下:
Figure BDA0000395017190000061
理论模型测试
为了进一步说明本方法的实现思路及实现过程并证明方法的有效性效果,用一个近地表具有6个不同尺度异常体地质模型进行试验,并和常规层析反演的结果进行比较。
1)建立一个近地表速度异常体地质模型(如图1所示),模型分为上下两层,上层速度为2000m/s,下层速度为2800m/s,上层中还包含六个尺寸不同的异常体,异常体速度均为3000m/s;
模型尺寸:深度2Km,宽度7Km;网格尺寸20*20m;模型速度:400m以下,2800m/s,400m以上,背景速度2000m/s,具有6个不同尺度的地质模型,分别为50*50m,50*100m,100*100m,100*200m,200*200m,200*400m,异常体顶界面深度100m。
2)利用声波方程有限差分正演模拟技术,得到炮点激发的正向传播波场;观测系统:接收点距20m,炮点距40m,滚动放炮,接收道数350道,激发炮数160炮。震源子波:Ricker子波,主频30Hz;与传统的正演模拟技术相比,由于采用了时间2阶,空间12阶的高精度有限差分,交错网格和PML(完美匹配层)边界条件等技术,正演模拟的地震记录和传播波场,频散效应降低,边界效应减少。如图2a、c所示,用传统交错网有限差分正演模拟得到的单炮波形记录和200ms的波场快照,可以看出边界反射明显;本发明采用时间2阶,空间12阶交错网格有限差分,并采用PML边界条件消除边界反射。如图2b、d所示,用本方法模拟的单炮记录和200ms的波场快照,边界效应明显降低,基本看不到边界反射。
3)由真实的实际速度模型模拟地震波场(如图3所示),地面观测波场作为单炮记录,进行初至波旅行时间拾取,用传统射线理论层析反演技术,反演近地表速度模型。如图4a所示,是旅行时层析反演近地表速度模型,与图3相比,只反演出大尺度异常体的大概位置和轮廓,没有反演出小尺度异常体,精度较低;但可以利用这个速度模型作为波形反演的初始速度。
4)利用层析反演得到的初始速度模型,模拟正向传播波场;对真实速度模型模拟的单炮记录进行时窗函数处理,只保留初至波和异常体绕射波信息,并作为初值条件,进行反向传播波场地震模拟(速度用初始速度);根据最速下降法求取速度扰动量,修改速度模型,完成一次初至波波形反演迭代;如图5a所示是第一次反演的近地表速度模型,6个不同尺度的近地表异常体都有一定响应,但精度较低。
5)把反演的结果作为初始模型,重复2)和4)步骤,逐级提高近地表速度模型的精度,直至结果满意为止。如图5b、图5c所示,分别是100次、500次迭代反演的结果,可以看到,随着迭代次数的增加,模型精度逐级提高;500次反演的结果,6个异常体边界和400m处的水平层清晰,速度值十分接近真实的速度值。如图6所示,是第1次、第100次、第300次、第500次迭代时的梯度剖面(归一化显示),可以看到随着速度精度的提高,梯度剖面幅值逐渐减少,从另一方面也说明该方法的稳定性、收敛性和有效性。

Claims (3)

1.一种初至波波形反演近地表速度模型方法,该方法包括声波方程交错网格有限差分正演模拟和最速下降法波形反演技术,其特征在于:实现过程的方法步骤如下:
(1)、提取时间域初至波波形记录Pobs和初始速度模型;
(2)、利用声波方程交错网格有限差分正演模拟计算模拟波场Pcal,并计算波场残差δp=Pobs-Pcal
(3)、将波场残差反向传播得到回传波场P′;
(4)、利用回传波场和正向传播波场计算目标函数的梯度,即速度的更新方向,并计算更新步长;
(5)、对速度模型进行更新;
(6)、检验是否满足迭代终止条件,满足则输出速度模型的反演结果,否则返回步骤(2),继续迭代更新。
2.根据权利要求1所述的初至波波形反演近地表速度模型方法,其特征在于,所述声波方程交错网格有限差分正演模拟技术,利用一阶声波方程有限差分方法,数值模拟炮点激发的正向传播波场;非均匀介质一阶声波方程:应力-质点速度方程;时间2阶、空间12阶的高阶有限差分,PML边界条件,交错网格。
3.根据权利要求1所述的初至波波形反演近地表速度模型方法,其特征在于,所述最速下降法波形反演技术,⑴地面观测记录加时窗函数—初至波窗函数,作为一阶声波方程的初值条件,类似于炮点激发正向传播波场正演技术,数值模拟初至波反向传播的地震波场和数值模拟的正向传播波场求残差,利用L2模建立反演目标函数;⑵利用正向传播波场和波场残差函数,计算目标函数的最速下降方向即梯度,确定出速度扰动的最大变化方向;⑶利用梯度和Frechet导数计算出速度扰动的迭代步长;⑷利用速度更新方向和更新步长更新速度。
CN201310476248.0A 2013-10-13 2013-10-13 初至波波形反演近地表速度模型方法 Pending CN103499835A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310476248.0A CN103499835A (zh) 2013-10-13 2013-10-13 初至波波形反演近地表速度模型方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310476248.0A CN103499835A (zh) 2013-10-13 2013-10-13 初至波波形反演近地表速度模型方法

Publications (1)

Publication Number Publication Date
CN103499835A true CN103499835A (zh) 2014-01-08

Family

ID=49865058

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310476248.0A Pending CN103499835A (zh) 2013-10-13 2013-10-13 初至波波形反演近地表速度模型方法

Country Status (1)

Country Link
CN (1) CN103499835A (zh)

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977607A (zh) * 2014-04-09 2015-10-14 中国石油集团东方地球物理勘探有限责任公司 利用变步长网格声波波场模拟的时间域全波形反演方法
CN104977608A (zh) * 2014-04-09 2015-10-14 中国石油集团东方地球物理勘探有限责任公司 利用固定网格声波波场模拟的时间域全波形反演方法
CN105445798A (zh) * 2014-08-21 2016-03-30 中国石油化工股份有限公司 一种基于梯度处理的全波形反演方法和系统
CN106054244A (zh) * 2016-06-16 2016-10-26 吉林大学 截断时窗的低通滤波多尺度全波形反演方法
CN106772596A (zh) * 2016-12-08 2017-05-31 中国石油天然气集团公司 一种确定叠前时间偏移速度场的方法及装置
CN107179547A (zh) * 2017-06-06 2017-09-19 中海石油(中国)有限公司 一种地震波阻抗反演低频模型建立方法
CN107203002A (zh) * 2017-06-12 2017-09-26 中国科学院地质与地球物理研究所 反演速度模型及其建立方法和地下结构的像的获得方法
CN107256290A (zh) * 2017-05-19 2017-10-17 四川腾盾科技有限公司 一种基于残差摄动法的边界条件变化对流场影响差量的计算方法
CN107450102A (zh) * 2017-07-28 2017-12-08 西安交通大学 基于分辨率可控包络生成算子的多尺度全波形反演方法
CN108369289A (zh) * 2015-12-18 2018-08-03 埃克森美孚上游研究公司 使用全波场反演点扩展函数分析设计地球物理勘测的方法
CN108646288A (zh) * 2018-02-27 2018-10-12 中国石油天然气集团有限公司 近地表模型的建立方法和装置
CN108845351A (zh) * 2018-06-26 2018-11-20 中国石油大学(华东) 一种vsp地震资料转换波全波形反演方法
CN109725345A (zh) * 2018-11-15 2019-05-07 中国石油天然气集团有限公司 一种初至波正演模拟方法及装置
WO2019105173A1 (zh) * 2017-11-30 2019-06-06 中国石油天然气集团有限公司 检波点定位准确度评价方法和装置
CN110426739A (zh) * 2019-08-02 2019-11-08 中铁第四勘察设计院集团有限公司 一种地质勘探检测方法、装置和存储介质
WO2020057286A1 (zh) * 2018-09-20 2020-03-26 中国科学院地质与地球物理研究所 波场正演模拟方法及装置
CN110942090A (zh) * 2019-11-11 2020-03-31 北京迈格威科技有限公司 模型训练、图像处理方法、装置、电子设备及存储介质
CN111007565A (zh) * 2019-12-24 2020-04-14 清华大学 三维频率域全声波成像方法及装置
CN111190224A (zh) * 2020-01-09 2020-05-22 中国石油大学(华东) 一种基于三维地震波反向照明的动态采样全波形反演系统及方法
CN114966826A (zh) * 2021-02-18 2022-08-30 中国石油化工股份有限公司 基于野外近地表调查资料的高低频融合近地表建模方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
ALBERT TARANTOLA: "Inversion of seismic reflection data in the acoustic approximation", 《GEOPHYSICS》 *
D.V.VIGH,ET AL.: "3D Prestack Time Domain Full Waveform Inversion", 《71ST EAGE CONFERENCE & EXHIBITION》 *
HASSAN KHANIANI,ET AL.: "Time domain full waveform inversion algorithm using common scatter point gathers", 《CREWES RESEARCH REPORT》 *
JOSEPH MCNEELY,ET AL.: "3D acoustic waveform inversion of land data: a case study from Saudi Arabia", 《SEG LAS VEGAS 2012 ANNUAL MEETING》 *
徐小云: "过井地震数据的声波方程走时和波形反演", 《图书与石油科技信息》 *
潘艳梅,等: "初至波波形反演方法及其数值模拟试验", 《石油物探》 *
潘艳梅,等: "起伏地表速度建模初步研究", 《新疆石油地质》 *
陈章: "叠前全波形反演方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977608A (zh) * 2014-04-09 2015-10-14 中国石油集团东方地球物理勘探有限责任公司 利用固定网格声波波场模拟的时间域全波形反演方法
CN104977607A (zh) * 2014-04-09 2015-10-14 中国石油集团东方地球物理勘探有限责任公司 利用变步长网格声波波场模拟的时间域全波形反演方法
CN104977608B (zh) * 2014-04-09 2017-07-07 中国石油集团东方地球物理勘探有限责任公司 利用固定网格声波波场模拟的时间域全波形反演方法
CN104977607B (zh) * 2014-04-09 2017-07-07 中国石油集团东方地球物理勘探有限责任公司 利用变步长网格声波波场模拟的时间域全波形反演方法
CN105445798A (zh) * 2014-08-21 2016-03-30 中国石油化工股份有限公司 一种基于梯度处理的全波形反演方法和系统
CN105445798B (zh) * 2014-08-21 2018-08-07 中国石油化工股份有限公司 一种基于梯度处理的全波形反演方法和系统
CN108369289A (zh) * 2015-12-18 2018-08-03 埃克森美孚上游研究公司 使用全波场反演点扩展函数分析设计地球物理勘测的方法
CN106054244B (zh) * 2016-06-16 2017-11-28 吉林大学 截断时窗的低通滤波多尺度全波形反演方法
CN106054244A (zh) * 2016-06-16 2016-10-26 吉林大学 截断时窗的低通滤波多尺度全波形反演方法
CN106772596A (zh) * 2016-12-08 2017-05-31 中国石油天然气集团公司 一种确定叠前时间偏移速度场的方法及装置
CN107256290A (zh) * 2017-05-19 2017-10-17 四川腾盾科技有限公司 一种基于残差摄动法的边界条件变化对流场影响差量的计算方法
CN107179547A (zh) * 2017-06-06 2017-09-19 中海石油(中国)有限公司 一种地震波阻抗反演低频模型建立方法
CN107203002A (zh) * 2017-06-12 2017-09-26 中国科学院地质与地球物理研究所 反演速度模型及其建立方法和地下结构的像的获得方法
CN107203002B (zh) * 2017-06-12 2019-05-24 中国科学院地质与地球物理研究所 反演速度模型的建立方法和地下结构的像的获得方法
CN107450102A (zh) * 2017-07-28 2017-12-08 西安交通大学 基于分辨率可控包络生成算子的多尺度全波形反演方法
US11493595B2 (en) 2017-11-30 2022-11-08 China National Petroleum Corporation Method and apparatus for evaluating accuracy in positioning a receiver point
WO2019105173A1 (zh) * 2017-11-30 2019-06-06 中国石油天然气集团有限公司 检波点定位准确度评价方法和装置
CN108646288A (zh) * 2018-02-27 2018-10-12 中国石油天然气集团有限公司 近地表模型的建立方法和装置
CN108845351A (zh) * 2018-06-26 2018-11-20 中国石油大学(华东) 一种vsp地震资料转换波全波形反演方法
WO2020057286A1 (zh) * 2018-09-20 2020-03-26 中国科学院地质与地球物理研究所 波场正演模拟方法及装置
CN109725345A (zh) * 2018-11-15 2019-05-07 中国石油天然气集团有限公司 一种初至波正演模拟方法及装置
CN109725345B (zh) * 2018-11-15 2020-08-11 中国石油天然气集团有限公司 一种初至波正演模拟方法及装置
CN110426739A (zh) * 2019-08-02 2019-11-08 中铁第四勘察设计院集团有限公司 一种地质勘探检测方法、装置和存储介质
CN110426739B (zh) * 2019-08-02 2021-07-16 中铁第四勘察设计院集团有限公司 一种地质勘探检测方法、装置和存储介质
CN110942090A (zh) * 2019-11-11 2020-03-31 北京迈格威科技有限公司 模型训练、图像处理方法、装置、电子设备及存储介质
CN110942090B (zh) * 2019-11-11 2024-03-29 北京迈格威科技有限公司 模型训练、图像处理方法、装置、电子设备及存储介质
CN111007565A (zh) * 2019-12-24 2020-04-14 清华大学 三维频率域全声波成像方法及装置
CN111190224A (zh) * 2020-01-09 2020-05-22 中国石油大学(华东) 一种基于三维地震波反向照明的动态采样全波形反演系统及方法
CN114966826A (zh) * 2021-02-18 2022-08-30 中国石油化工股份有限公司 基于野外近地表调查资料的高低频融合近地表建模方法

Similar Documents

Publication Publication Date Title
CN103499835A (zh) 初至波波形反演近地表速度模型方法
CN103713315B (zh) 一种地震各向异性参数全波形反演方法及装置
CN108345031B (zh) 弹性介质主动源和被动源混采地震数据全波形反演方法
CN104977618B (zh) 一种评价页岩气储层及寻找甜点区的方法
CN103454685B (zh) 利用测井约束波阻抗反演预测砂体厚度的方法和装置
CN103913774B (zh) 基于微地震事件的储层地质力学参数反演方法
CN104570082B (zh) 一种基于格林函数表征的全波形反演梯度算子的提取方法
CN102183790A (zh) 基于时空双变网格的弹性波正演模拟技术
CN106842295A (zh) 测井信息约束的波形反演方法
CN103293553B (zh) 一种复杂海底上下缆地震采集数据边界元延拓校正方法
CN103513277B (zh) 一种地震地层裂隙裂缝密度反演方法及系统
CN108508482A (zh) 一种地下裂缝地震散射响应特征模拟方法
CN104749617A (zh) 一种多尺度裂缝储层正演模型建立方法
CN102901985B (zh) 一种适用于起伏地表的深度域层速度修正方法
CN105549080B (zh) 一种基于辅助坐标系的起伏地表波形反演方法
CN105549079A (zh) 一种地球物理参数的全波形反演模型的建立方法和装置
CN102116869A (zh) 高精度叠前域最小二乘偏移地震成像技术
CN105093278A (zh) 基于激发主能量优化算法的全波形反演梯度算子的提取方法
CN107462924A (zh) 一种不依赖于测井资料的绝对波阻抗反演方法
CN103576200A (zh) 低信噪比区浅层波阻抗界面静校正方法
CN103913768A (zh) 基于地震波资料对地表中浅层进行建模的方法及装置
CN107894618A (zh) 一种基于模型平滑算法的全波形反演梯度预处理方法
CN104122581A (zh) 一种叠后声波阻抗反演方法
CN110531410A (zh) 一种基于直达波场的最小二乘逆时偏移梯度预条件方法
CN107179547A (zh) 一种地震波阻抗反演低频模型建立方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140108