CN103018779A - 一种海上地震勘探***震源子波模拟方法及*** - Google Patents
一种海上地震勘探***震源子波模拟方法及*** Download PDFInfo
- Publication number
- CN103018779A CN103018779A CN201210507082XA CN201210507082A CN103018779A CN 103018779 A CN103018779 A CN 103018779A CN 201210507082X A CN201210507082X A CN 201210507082XA CN 201210507082 A CN201210507082 A CN 201210507082A CN 103018779 A CN103018779 A CN 103018779A
- Authority
- CN
- China
- Prior art keywords
- bubble
- unit
- initial
- gas
- air gun
- 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.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种海上地震勘探***震源子波模拟方法及***:测量海水及***装置的基本参数;计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径;确定气泡振动过程中气泡内温度随时间的变化函数,气泡内气体的摩尔质量随时间的变化函数,气泡内的气体体积随时间的变化函数;模拟随时间变化的海水粘度;考虑海水粘度和虚反射压力场,模拟气泡产生的压力波;由气泡产生的压力波,模拟其虚反射压力波;由气泡产生的压力波及虚反射压力波叠加合成***子波。本发明增强了模拟子波与实测子波的一致性。
Description
技术领域
本发明涉及油气地球物理勘探中的地震资料采集领域,尤其涉及一种海上地震勘探***震源子波模拟方法及***。
背景技术
海上地震勘探是一种利用人工地震技术探测地下构造的勘探方法。它按照一定的方式在海面附近激发地震波,产生称之为地震子波的振动信号,地震子波由震源开始向下传播,遇到地质界面之后,在界面处发生透射和反射,透射的地震子波继续向下传播,而反射子波在界面位置向上传播。来自不同深度界面的地震子波以不同的时间到达海水表面,通过布设在海面的一种称之为检波器的接收装置,接收来自不同深度地质界面的反射波,接收到的数字信号称为地震记录。
海上地震采集技术来源于陆上地震勘探,因此,早期的海上勘探也使用陆地勘探的***震源。但很快暴露出其致命弱点:其一,自动化程度差,人工操作危险性大;其二,有悖于环保理念,对海洋及其海洋生物危害很大。有鉴于此,空***、蒸***、水枪等非***震源应运而生。其中,空***以其性能稳定、自动化程度高、成本较低等诸多优点逐渐占据主导地位。
***震源所激发的地震子波是海上地震勘探的基本地震信号,其形态特征和频率特征直接关系到海上地震勘探的最终质量,因此,***震源子波数值模拟和特征分析对海上地震勘探具有十分重要的作用。
Keller和Kolodner(1956)认为***产生的高压气体是一个自由振荡的球形气泡,并提出了“自由气泡振荡原理”。基于这一理论,Ziolkowski(1970)提出了用于模拟***子波的基本数学模型,但是与实测子波相比,Ziolkowski模型所模拟的子波初始脉冲值太大,而且气泡振荡衰减太慢。Schulze-Gattermann(1972)提出了假设枪体为刚性球体的新的气泡震荡模型;Ziolkowski(1984)对枪体节流作用和气泡壁热传导作用做出了详细分析,总结出气体在枪室中并非瞬间释放,而且气体的释放过程会影响***子波的形状;Langhammer(1993)通过实验性研究量化了海水粘度对***子波的影响。李国发等(2010)综合考虑枪体节流作用、气泡壁的热传导作用、气泡上浮以及流体粘度等因素对***子波的影响,完成了考虑多种影响因素的单枪子波模拟。
然而,以上模型存在两个基本缺陷,首先,上述模型假设流体粘度为一常数,忽略了气泡在振动过程中温度和压力的变化对流体粘度的影响。再者,上述模型虽然也考虑了虚反射对远场子波的影响,但只是简单地将虚反射叠加在气泡振动之上,忽略了虚反射产生的压力场对气泡振动的再次影响。
发明内容
本发明实施例提供一种海上地震勘探***震源子波模拟方法及***,以建立一种更加完善的***子波模拟***,使得其所模拟的子波与野外实测子波具有更高的一致性。
一方面,本发明实施例提供了一种海上地震勘探***震源子波模拟方法,所述海上地震勘探***震源子波模拟方法包括:
(1)测量并获取海水及其***装置的基本参数;
(2)利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径;
(3)按照Ziolkowski气泡振动理论,确定气泡振动过程中气泡内温度随时间的变化函数T(iΔt),i=0,1,…n,气泡内气体的摩尔质量随时间的变化函数m(iΔt),i=0,1,…n,气泡内的气体体积随时间的变化函数V(iΔt),i=0,1,…n,其中,i时间是采样序号,Δt是采样间隔,单位是ms,n是时间采样个数;
(4)由气泡内的温度T(iΔt),i=0,1,…n,模拟随时间变化的海水粘度μ(iΔt),i=0,1…n,单位为Kg/(m·s);
μ(iΔt)=100+1012×T-6.08(iΔt);
(5)考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n;
其中,R为海水的反射系数,R=-1,无量纲,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
可选的,在本发明一实施例中,所述(1)中,测量并获取海水及其***装置的基本参数,包括:
(1.1)测量海水温度Tw,单位为K,和海水密度ρ,单位为Kg/m3;
(1.2)测量***装置的下列参数:
a.枪室的工作压强pg,单位为psi;
b.***枪室的容积Vg,单位为in3;
c.枪体的释放效率η,无量纲;
d.与枪体容量无关的常数τ0和节流指数β,无量纲;
e.***沉放深度h,单位为m。
可选的,在本发明一实施例中,所述(2)中,利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径,包括:
(2.1)由下式计算枪体节流常数τ,无量纲;
τ=τ0(Vg)β;
(2.2)由海水温度Tw和***枪室的初始压强pg,计算***枪室的初始温度Tg,单位为K,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol,
其中RG为普氏常量,其值为8.2;
m′g=mgη;
(2.4)由海水密度ρ和***沉放深度h,计算初始静水压强p∞,单位为Pa,
p∞=p0+ρgh,
其中大气压强p0=101300Pa,重力加速度g=9.8m/s2;
(2.5)由初始静水压强p∞,计算气泡内气体的初始摩尔质量m(0),单位为mol;
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p∞,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
rg=r0×0.23;
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3;
可选的,在本发明一实施例中,所述(5)中,考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n,具体步骤包括:
(5.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(iΔt),i=0,1,…n,其中,R(iΔt)是随时间变化的气泡半径,初值为单位为m,是气泡半径随时间的变化率,即气泡壁的速度,初值为单位为m/s;
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
其中,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.4)由以上结果计算i+1时刻的气泡壁加速度i=0,1,…n,单位为m/s2,
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w[(i+1)Δt],由此迭代方式得到所有时刻的压力波w(iΔt),i=0,1,…n。
另一方面,本发明实施例提供了一种海上地震勘探***震源子波模拟***,所述海上地震勘探***震源子波模拟***包括:
测量单元,用于测量并获取海水及其***装置的基本参数;
计算单元,用于利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径;
函数单元,用于按照Ziolkowski气泡振动理论,确定气泡振动过程中气泡内温度随时间的变化函数T(iΔt),i=0,1,…n,气泡内气体的摩尔质量随时间的变化函数m(iΔt),i=0,1,…n,气泡内的气体体积随时间的变化函数V(iΔt),i=0,1,…n,其中,i时间是采样序号,Δt是采样间隔,单位是ms,n是时间采样个数;
模拟单元,用于由气泡内的温度T(iΔt),i=0,1,…n,模拟随时间变化的海水粘度μ(iΔt),i=0,1…n,单位为Kg/(m·s);
μ(iΔt)=100+1012×T-6.08(iΔt);
用于考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n;
其中,R为海水的反射系数,R=-1,无量纲,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms
可选的,在本发明一实施例中,所述测量单元,进一步用于测量并获取海水及其***装置的基本参数,包括:
(1.1)测量海水温度Tw,单位为K,和海水密度ρ,单位为Kg/m3;
(1.2)测量***装置的下列参数:
a.枪室的工作压强pg,单位为psi;
b.***枪室的容积Vg,单位为in3;
c.枪体的释放效率η,无量纲;
d.与枪体容量无关的常数τ0和节流指数β,无量纲;
e.***沉放深度h,单位为m。
可选的,在本发明一实施例中,所述计算单元,进一步用于利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径,包括:
(2.1)由下式计算枪体节流常数τ,无量纲;
τ=τ0(Vg)β;
(2.2)由海水温度Tw和***枪室的初始压强pg,计算***枪室的初始温度Tg,单位为K,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol,
其中RG为普氏常量,其值为8.2;
mg′=mgη;
(2.4)由海水密度ρ和***沉放深度h,计算初始静水压强p∞,单位为Pa,
p∞=p0+ρgh,
其中大气压强p0=101300Pa,重力加速度g=9.8m/s2;
(2.5)由初始静水压强p∞,计算气泡内气体的初始摩尔质量m(0),单位为mol;
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p∞,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
rg=r0×0.23;
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3;
可选的,在本发明一实施例中,所述模拟单元,进一步用于考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n,具体包括:
(5.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(iΔt),i=0,1,…n,其中,R(iΔt)是随时间变化的气泡半径,初值为单位为m,是气泡半径随时间的变化率,即气泡壁的速度,初值为单位为m/s;
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
首先由气泡半径R(iΔt)和气泡壁的速度得到中间函数f(iΔt);
其中,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w[(i+1)Δt],由此迭代方式得到所有时刻的压力波w(iΔt),i=0,1,…n
上述技术方案具有如下有益效果:因为采用综合考虑流体粘度随温度变化以及虚反射对气泡周围压力的影响的技术手段,所以达到了如下的技术效果:建立了一种更加完善的***子波模拟***,增强模拟子波与实测子波的一致性,为海上地震勘探提供基础数据,提高海上地震勘探精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一种海上地震勘探***震源子波模拟方法流程图;
图2为本发明实施例一种海上地震勘探***震源子波模拟***结构示意图;
图3为本发明应用实例本***所模拟的气泡内温度随时间的变化;
图4为本发明应用实例本***所模拟气泡中气体质量随时间的变化;
图5为本发明应用实例本***所模拟气泡中气体体积随时间的变化;
图6为本发明应用实例本***所模拟海水粘度随温度的变化;
图7为本发明应用实例本***所模拟的气泡产生的压力波随时间的变化;
图8为本发明应用实例本***所模拟的气泡半径随时间的变化;
图9为本发明应用实例本***所模拟的气泡壁速度随时间的变化;
图10为本发明应用实例本***所模拟的气泡壁加速度随时间的变化;
图11为本发明应用实例本***所模拟的气泡壁压力随时间的变化;
图12为本发明应用实例本***所模拟的气泡虚反射压力波随时间的变化;
图13为本发明应用实例本***所模拟的***子波随时间的变化;
图14为本发明应用实例某海上工区实际激发的***子波与按照李国发等所提供的方法模拟的***子波对比示意图;
图15为本发明应用实例某海上工区实际激发的***子波与本发明所模拟的***子波对比示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种海上地震勘探***震源子波模拟方法及***,以建立一种更加完善的***子波模拟***,使得其所模拟的子波与野外实测子波具有更高的一致性。
如图1所示,为本发明实施例一种海上地震勘探***震源子波模拟方法流程图,所述海上地震勘探***震源子波模拟方法包括:
(1)测量并获取海水及其***装置的基本参数;
(2)利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径;
(3)按照Ziolkowski气泡振动理论,确定气泡振动过程中气泡内温度随时间的变化函数T(iΔt),i=0,1,…n,气泡内气体的摩尔质量随时间的变化函数m(iΔt),i=0,1,…n,气泡内的气体体积随时间的变化函数V(iΔt),i=0,1,…n,其中,i时间是采样序号,Δt是采样间隔,单位是ms,n是时间采样个数;
(4)由气泡内的温度T(iΔt),i=0,1,…n,模拟随时间变化的海水粘度μ(iΔt),i=0,1…n,单位为Kg/(m·s);μ(iΔt)=100+1012×T-6.08(iΔt);
(5)考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n;
(6)由气泡产生的压力波w(iΔt),i=0,1,…n,模拟其虚反射压力波i=0,1,…n,单位为Pa,
其中,R为海水的反射系数,R=-1,无量纲,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
可选的,在本发明一实施例中,所述(1)中,测量并获取海水及其***装置的基本参数,包括:
(1.1)测量海水温度Tw,单位为K,和海水密度ρ,单位为Kg/m3;
(1.2)测量***装置的下列参数:
a.枪室的工作压强pg,单位为psi;
b.***枪室的容积Vg,单位为in3;
c.枪体的释放效率η,无量纲;
d.与枪体容量无关的常数τ0和节流指数β,无量纲;
e.***沉放深度h,单位为m。
可选的,在本发明一实施例中,所述(2)中,利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径,包括:
(2.1)由下式计算枪体节流常数τ,无量纲;
τ=τ0(Vg)β;
(2.2)由海水温度Tw和***枪室的初始压强pg,计算***枪室的初始温度Tg,单位为K,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol,
其中RG为普氏常量,其值为8.2;
mg′=mgη;
(2.4)由海水密度ρ和***沉放深度h,计算初始静水压强p∞,单位为Pa,
p∞=p0+ρgh,
其中大气压强p0=101300Pa,重力加速度g=9.8m/s2;
(2.5)由初始静水压强p∞,计算气泡内气体的初始摩尔质量m(0),单位为mol;
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p∞,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
rg=r0×0.23;
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3;
可选的,在本发明一实施例中,所述(5)中,考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n,具体步骤包括:
(5.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(iΔt),i=0,1,…n,其中,R(iΔt)是随时间变化的气泡半径,初值为单位为m,是气泡半径随时间的变化率,即气泡壁的速度,初值为单位为m/s;
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
再利用中间函数f(iΔt)及其对时间的导数f′(iΔt),模拟气泡产生的压力波w(iΔt),i=0,1,…n和虚反射压力场
其中,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w「(i+1)Δt]由此迭代方式得到所有时刻的压力波w(iΔt),i=0,1,…n。
如图2所示,为本发明实施例一种海上地震勘探***震源子波模拟***结构示意图,所述海上地震勘探***震源子波模拟***包括:
测量单元21,用于测量并获取海水及其***装置的基本参数;
计算单元22,用于利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径;
函数单元23,用于按照Ziolkowski气泡振动理论,确定气泡振动过程中气泡内温度随时间的变化函数T(iΔt),i=0,1,…n,气泡内气体的摩尔质量随时间的变化函数m(iΔt),i=0,1,…n,气泡内的气体体积随时间的变化函数V(iΔt),i=0,1,…n,其中,i时间是采样序号,Δt是采样间隔,单位是ms,n是时间采样个数;
模拟单元24,用于由气泡内的温度T(iΔt),i=0,1,…n,模拟随时间变化的海水粘度μ(iΔt),i=0,1…n,单位为Kg/(m·s);
μ(iΔt)=100+1012×T-6.08(iΔt);
用于考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n;
其中,R为海水的反射系数,R=-1,无量纲,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
可选的,在本发明一实施例中,所述测量单元21,进一步用于测量并获取海水及其***装置的基本参数,包括:
(1.1)测量海水温度Tw,单位为K,和海水密度ρ,单位为Kg/m3;
(1.2)测量***装置的下列参数:
a.枪室的工作压强pg,单位为psi;
b.***枪室的容积Vg,单位为in3;
c.枪体的释放效率η,无量纲;
d.与枪体容量无关的常数τ0和节流指数β,无量纲;
e.***沉放深度h,单位为m。
可选的,在本发明一实施例中,所述计算单元22,进一步用于利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径,包括:
(2.1)由下式计算枪体节流常数τ,无量纲;
τ=τ0(Vg)β;
(2.2)由海水温度Tw和***枪室的初始压强pg,计算***枪室的初始温度Tg,单位为K,
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol,
其中RG为普氏常量,其值为8.2;
mg′=mgη;
(2.4)由海水密度ρ和***沉放深度h,计算初始静水压强p∞,单位为Pa,
p∞=p0+ρgh,
其中大气压强p0=101300Pa,重力加速度g=9.8m/s2;
(2.5)由初始静水压强p∞,计算气泡内气体的初始摩尔质量m(0),单位为mol;
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p∞,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
rg=r0×0.23;
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3;
可选的,在本发明一实施例中,所述模拟单元24,进一步用于考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n,具体包括:
(5.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(iΔt),i=0,1,…n,其中,R(iΔt)是随时间变化的气泡半径,初值为单位为m,是气泡半径随时间的变化率,即气泡壁的速度,初值为单位为m/s;
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
其中,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w[(i+1)Δt]由此迭代方式得到所有时刻的压力波w(iΔt),i=0,1,…n。
本发明上述方法或***实施例上述技术方案具有如下有益效果:因为采用综合考虑流体粘度随温度变化以及虚反射对气泡周围压力的影响的技术手段,所以达到了如下的技术效果:建立了一种更加完善的***子波模拟***,增强模拟子波与实测子波的一致性,为海上地震勘探提供基础数据,提高海上地震勘探精度。
以下举应用实例进行详细说明:
本发明应用实例通过综合考虑流体粘度随温度的变化以及虚反射对气泡周围压力的影响,建立一种更完善的模拟***子波的模型,使其所模拟的子波与野外实测子波具有更高的一致性。
下面以某***震源海上勘探的实际数据为例。如图3所示,为本发明应用实例本***所模拟的气泡内温度随时间的变化;如图4所示,为本发明应用实例本***所模拟气泡中气体质量随时间的变化;如图5所示,为本发明应用实例本***所模拟气泡中气体体积随时间的变化;如图6所示,为本发明应用实例本***所模拟海水粘度随温度的变化;如图7所示,为本发明应用实例本***所模拟的气泡产生的压力波随时间的变化;如图8所示,为本发明应用实例本***所模拟的气泡半径随时间的变化;如图9所示,为本发明应用实例本***所模拟的气泡壁速度随时间的变化;如图10所示,为本发明应用实例本***所模拟的气泡壁加速度随时间的变化;如图11所示,为本发明应用实例本***所模拟的气泡壁压力随时间的变化;如图12所示,为本发明应用实例本***所模拟的气泡虚反射压力波随时间的变化;如图13所示,为本发明应用实例本***所模拟的***子波随时间的变化。以下说明具体实施过程:
(1)测量海水及其***装置的基本参数:
(1.1)测量海水温度Tw=300K和海水密度ρ=1000Kg/m3;
(1.2)测量***装置的下列参数:
a.枪室的工作压强pg=2000psi;
b.***枪室的容积Vg=40in3;
c.枪体的释放效率η=0.8317,无量纲;
d.与枪体容量无关的常数τ0=0.12和节流指数β=0.52,无量纲;
e.***沉放深度h=5m;
(2)利用上述测量参数计算枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径。
(2.1)由下式计算枪体节流常数τ,无量纲;
τ=τ0(Vg)β=0.0027
(2.2)由海水温度Tw和***枪室的初始压强pg,计算***枪室的初始温度Tg,单位为K
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol
其中RG为普氏常量,其值为8.2;
mg′=mgη=2.7784mol
(2.4)由海水密度ρ和***沉放深度h,计算初始静水压强p∞,单位为Pap∞=p0+ρgh=135900Pa
其中大气压强p0=101300Pa,重力加速度g=9.8m/s2;
(2.5)由初始静水压强p∞,计算气泡内气体的初始摩尔质量m(0),单位为mol;
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p∞,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
rg=r0×0.23=0.0509m
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3;
(3)按照Ziolkowski气泡振动理论,确定气泡振动过程中如图3所示气泡内温度随时间的变化函数T(iΔt),i=0,1,…n,如图4所示的气泡内气体摩尔质量随时间的变化函数m(iΔt),i=0,1,…n,如图5所示的气泡内气体体积随时间的变化函数V(iΔt),i=0,1,…n,其中,i时间是采样序号,Δt是采样间隔,单位是ms,n是时间采样个数,本例中,采样间隔Δt=0.1ms,采样个数n=3001;
(4)由气泡内的温度T(iΔt),i=0,1,…n,模拟如图6所示随时间变化的海水粘度μ(iΔt),i=0,1…n,单位为Kg/(m·s);
μ(iΔt)=100+1012×T-6.08(iΔt)
(5)考虑海水粘度和虚反射压力场对气泡振动的影响,模拟如图7所示气泡产生的压力波w(iΔt),i=0,1,…n,具体步骤为:
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
首先由气泡半径R(iΔt)和气泡壁的速度得到中间函数f(iΔt);
其中,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.6)返回步骤(5.2),计算i+1时刻气泡产生的压力波w[(i+1)Δt],由此迭代方式得到所有时刻的压力波w(iΔt),i=0,1,…n。在此过程中产生的中间成果有,如图8所示的气泡半径R(iΔt),i=0,1,…n随时间的变化情况,如图9所示的气泡半径随时间的变化率即气泡壁速度i=0,1,…n随时间的变化情况,如图10所示的气泡壁加速度i=0,1,…n随时间的变化情况,如图11所示的气泡壁压力p(iΔt),i=0,1,…n随时间的变化情况。
(6)由气泡产生的压力波w(iΔt),i=0,1,…n,模拟如图12所示虚反射压力波i=0,1,…n,单位为Pa
其中,R为海水的反射系数,R=-1,无量纲,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
例如:本实施例为某***震源海上地震勘探工区的应用实例,其中***容量Vg=80in3,陈放深度h=8m,***工作压力pg=2000psi;海水中声波速度为c=1.483m/ms,海水温度Tw=20℃。
本实施例***子波数据采样间隔为0.1ms,记录长度为300ms。如图14所示,为本发明应用实例某海上工区实际激发的***子波与按照李国发等(2010)所提供的方法(具体方法参见Li Guofa,Cao Mingqiang,et al.2010,Modeling of the Signature of Air Gun in Marine SeismicExploration Considering the effects of multiple Practical Physics.Applied Geophysics,7(2):158-165)模拟的***子波对比示意图,其中,实线是该海上工区实际激发的***子波,虚线是按照李国发等(2010)所提供的方法在没有考虑流体粘度随温度的变化及其虚反射压力场对气泡振动影响情况下模拟的***子波。如图15所示,为本发明应用实例某海上工区实际激发的***子波与本发明所模拟的***子波对比示意图,其中,实线是该海上工区实际激发的***子波,虚线是本发明的方法在考虑流体粘度随温度的变化及其虚反射压力场对气泡振动影响情况下模拟的***子波。与图14进行对比可以看出,本发明所模拟的子波与实际子波具有更高的一致性。
本发明综合考虑流体粘度随温度的变化以及虚反射对气泡周围压力的影响,对传统***子波模拟模型进行修正,建立了一种更加完善的***子波模拟***,所模拟的子波与野外实测子波具有更高的一致性,对海上地震勘探具有重要的指导意义和参考价值。
本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrativelogical block),单元,和步骤可以通过电子硬件、电脑软件,或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability),上述的各种说明性部件(illustrativecomponents),单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个***的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明实施例中所描述的各种说明性的逻辑块,或单元都可以通过通用处理器,数字信号处理器,专用集成电路(ASIC),现场可编程门阵列(FPGA)或其它可编程逻辑装置,离散门或晶体管逻辑,离散硬件部件,或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器,可选地,该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现,例如数字信号处理器和微处理器,多个微处理器,一个或多个微处理器联合一个数字信号处理器核,或任何其它类似的配置来实现。
本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地,存储媒介可以与处理器连接,以使得处理器可以从存储媒介中读取信息,并可以向存储媒介存写信息。可选地,存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中,ASIC可以设置于用户终端中。可选地,处理器和存储媒介也可以设置于用户终端中的不同的部件中。
在一个或多个示例性的设计中,本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现,这些功能可以存储与电脑可读的媒介上,或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如,这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置,或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外,任何连接都可以被适当地定义为电脑可读媒介,例如,如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘,磁盘通常以磁性复制数据,而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种海上地震勘探***震源子波模拟方法,其特征在于,所述海上地震勘探***震源子波模拟方法包括:
(1)测量并获取海水及其***装置的基本参数;
(2)利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径;
(3)按照Ziolkowski气泡振动理论,确定气泡振动过程中气泡内温度随时间的变化函数T(iΔt),i=0,1,…n,气泡内气体的摩尔质量随时间的变化函数m(iΔt),i=0,1,…n,气泡内的气体体积随时间的变化函数V(iΔt),i=0,1,…n,其中,i时间是采样序号,Δt是采样间隔,单位是ms,n是时间采样个数;
(4)由气泡内的温度T(iΔt),i=0,1,…n,模拟随时间变化的海水粘度μ(iΔt),i=0,1…n,单位为Kg/(m·s);
μ(iΔt)=100+1012×T-608(iΔt);
(5)考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n;
其中,R为海水的反射系数,R=-1,无量纲,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
2.如权利要求1所述海上地震勘探***震源子波模拟方法,其特征在于,所述(1)中,测量并获取海水及其***装置的基本参数,包括:
(1.1)测量海水温度Tw,单位为K,和海水密度ρ,单位为Kg/m3;
(1.2)测量***装置的下列参数:
a.枪室的工作压强pg,单位为psi;
b.***枪室的容积Vg,单位为in3;
c.枪体的释放效率η,无量纲;
d.与枪体容量无关的常数τ0和节流指数β,无量纲;
e.***沉放深度h,单位为m。
3.如权利要求1所述海上地震勘探***震源子波模拟方法,其特征在于,所述(2)中,利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径,包括:
(2.1)由下式计算枪体节流常数τ,无量纲;
τ=τ0(Vg)β;
(2.2)由海水温度Tw和***枪室的初始压强pg,计算***枪室的初始温度Tg,单位为K,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol,
其中RG为普氏常量,其值为8.2;
mg=mgη;(2.4)由海水密度ρ和***沉放深度h,计算初始静水压强p∞,单位为Pa,p∞=p0+ρgh,
其中大气压强p0=101300Pa,重力加速度g=9.8m/s2;
(2.5)由初始静水压强p∞,计算气泡内气体的初始摩尔质量m(0),单位为mol;
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p∞,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
rg=r0×0.23;
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3;
4.如权利要求1所述海上地震勘探***震源子波模拟方法,其特征在于,所述(5)中,考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n,具体步骤包括:
(5.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(iΔt),i=0,1,…n,其中,R(iΔt)是随时间变化的气泡半径,初值为 单位为m,是气泡半径随时间的变化率,即气泡壁的速度,初值为单位为m/s;
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡 内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
再利用中间函数f(iΔt)及其对时间的导数f′(iΔt),模拟气泡产生的压力波w(iΔt),i=0,1,…n和虚反射压力场
其中,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w[(i+1)Δt],由此迭代方式得到所有时刻的压力波w(iΔt),i=0,1,…n。
5.一种海上地震勘探***震源子波模拟***,其特征在于,所述海上地震勘探***震源子波模拟***包括:
测量单元,用于测量并获取海水及其***装置的基本参数;
计算单元,用于利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径;
函数单元,用于按照Ziolkowski气泡振动理论,确定气泡振动过程中气泡内温度随时间的变化函数T(iΔt),i=0,1,…n,气泡内气体的摩尔质量随时间的变化函数m(iΔt),i=0,1,…n,气泡内的气体体积随时间的变化函数V(iΔt),i=0,1,…n,其中,i时间是采样序号,Δt是采样间隔,单位是ms,n是时间采样个数;
模拟单元,用于由气泡内的温度T(iΔt),i=0,1,…n,模拟随时间变化的海水粘度μ(iΔt),i=0,1…n,单位为Kg/(m·s);
μ(iΔt)=100+1012×T-6.08(iΔt);
用于考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n;
其中,R为海水的反射系数,R=-1,无量纲,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
6.如权利要求5所述海上地震勘探***震源子波模拟***,其特征在于,所述测量单元,进一步用于测量并获取海水及其***装置的基本参数,包括:
(1.1)测量海水温度Tw,单位为K,和海水密度ρ,单位为Kg/m3;
(1.2)测量***装置的下列参数:
a.枪室的工作压强pg,单位为psi;
b.***枪室的容积Vg,单位为in3;
c.枪体的释放效率η,无量纲;
d.与枪体容量无关的常数τ0和节流指数β,无量纲;
e.***沉放深度h,单位为m。
7.如权利要求5所述海上地震勘探***震源子波模拟***,其特征在于,所述计算单元,进一步用于利用所述基本参数计算:枪体节流常数,***枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径,包括:
(2.1)由下式计算枪体节流常数τ,无量纲;
τ=τ0(Vg)β;
(2.2)由海水温度Tw和***枪室的初始压强pg,计算***枪室的初始温度Tg,单位为K,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol,
其中RG为普氏常量,其值为8.2;mg′=mgη;
(2.4)由海水密度ρ和***沉放深度h,计算初始静水压强p∞,单位为Pa, p∞=p0+ρgh,
其中大气压强p0=101300Pa,重力加速度g=9.8m/s2;
(2.5)由初始静水压强p∞,计算气泡内气体的初始摩尔质量m(0),单位为mol;
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p∞,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
rg=r0×0.23;
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3;
8.如权利要求5所述海上地震勘探***震源子波模拟***,其特征在于,所述模拟单元,进一步用于考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n,具体包括:
(5.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(iΔt),i=0,1,…n,其中,R(iΔt)是随时间变化的气泡半径,初值为 单位为m,是气泡半径随时间的变化率,即气泡壁的速 度,初值为单位为m/s;
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
首先由气泡半径R(iΔt)和气泡壁的速度得到中间函数f(iΔt);
其中,h为***沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w[(i+1)Δt],由此迭代方式得到所有时刻的压力波w(iΔt),i=0,1,…n。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210507082.XA CN103018779B (zh) | 2012-11-30 | 2012-11-30 | 一种海上地震勘探***震源子波模拟方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210507082.XA CN103018779B (zh) | 2012-11-30 | 2012-11-30 | 一种海上地震勘探***震源子波模拟方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103018779A true CN103018779A (zh) | 2013-04-03 |
CN103018779B CN103018779B (zh) | 2015-10-21 |
Family
ID=47967591
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210507082.XA Active CN103018779B (zh) | 2012-11-30 | 2012-11-30 | 一种海上地震勘探***震源子波模拟方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103018779B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107024714A (zh) * | 2017-03-24 | 2017-08-08 | 中国石油天然气集团公司 | 一种实现***震源质控的处理方法及装置 |
CN108061915A (zh) * | 2016-11-09 | 2018-05-22 | 中国石油化工股份有限公司 | 一种***子波模拟方法 |
CN109239769A (zh) * | 2018-11-01 | 2019-01-18 | 国家***第二海洋研究所 | 利用深度和容量组合压制残留气泡的***震源设计方法 |
CN110687617A (zh) * | 2019-08-26 | 2020-01-14 | 中国海洋大学 | 一种地震勘探***阵列远场子波模拟方法、脸谱评价方法及装置 |
CN113655519A (zh) * | 2021-08-23 | 2021-11-16 | 中海石油(中国)有限公司 | ***节流作用系数和气体释放效率参数获取方法及*** |
CN114647990A (zh) * | 2022-02-22 | 2022-06-21 | 哈尔滨工程大学 | 一种高效高精度计算***震源压力子波的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6081765A (en) * | 1995-12-20 | 2000-06-27 | Ziolkowski; Antoni Marjan | Signatures of marine seismic sources |
CN101241192A (zh) * | 2007-02-06 | 2008-08-13 | 中国石油集团东方地球物理勘探有限责任公司 | 一种消除***近场子波虚反射的方法 |
CN102298156A (zh) * | 2010-06-10 | 2011-12-28 | 地球物理维里达斯集团公司 | 用于反虚反射地震数据的方法和装置 |
-
2012
- 2012-11-30 CN CN201210507082.XA patent/CN103018779B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6081765A (en) * | 1995-12-20 | 2000-06-27 | Ziolkowski; Antoni Marjan | Signatures of marine seismic sources |
CN101241192A (zh) * | 2007-02-06 | 2008-08-13 | 中国石油集团东方地球物理勘探有限责任公司 | 一种消除***近场子波虚反射的方法 |
CN102298156A (zh) * | 2010-06-10 | 2011-12-28 | 地球物理维里达斯集团公司 | 用于反虚反射地震数据的方法和装置 |
Non-Patent Citations (5)
Title |
---|
唐杰 等: "大容量***震源的实验与模拟研究", 《中国地震》 * |
李国发 等: "枪口节流作用和热传导作用对海上勘探***子波的影响", 《特种油气藏》 * |
李国发 等: "海上地震勘探GI枪子波数值模拟", 《吉林大学学报(地球科学版)》 * |
李国发 等: "考虑多种因素影响的海上地震勘探***震源单枪子波数值模拟", 《应用地球物理》 * |
陈浩林 等: "***震源单枪子波计算机模拟", 《物探装备》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108061915A (zh) * | 2016-11-09 | 2018-05-22 | 中国石油化工股份有限公司 | 一种***子波模拟方法 |
CN107024714A (zh) * | 2017-03-24 | 2017-08-08 | 中国石油天然气集团公司 | 一种实现***震源质控的处理方法及装置 |
CN107024714B (zh) * | 2017-03-24 | 2019-02-15 | 中国石油天然气集团公司 | 一种实现***震源质控的处理方法及装置 |
CN109239769A (zh) * | 2018-11-01 | 2019-01-18 | 国家***第二海洋研究所 | 利用深度和容量组合压制残留气泡的***震源设计方法 |
CN110687617A (zh) * | 2019-08-26 | 2020-01-14 | 中国海洋大学 | 一种地震勘探***阵列远场子波模拟方法、脸谱评价方法及装置 |
CN113655519A (zh) * | 2021-08-23 | 2021-11-16 | 中海石油(中国)有限公司 | ***节流作用系数和气体释放效率参数获取方法及*** |
CN113655519B (zh) * | 2021-08-23 | 2023-10-13 | 中海石油(中国)有限公司 | ***节流作用系数和气体释放效率参数获取方法及*** |
CN114647990A (zh) * | 2022-02-22 | 2022-06-21 | 哈尔滨工程大学 | 一种高效高精度计算***震源压力子波的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103018779B (zh) | 2015-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103018779B (zh) | 一种海上地震勘探***震源子波模拟方法及*** | |
Chouet et al. | Seismic source mechanism of degassing bursts at Kilauea Volcano, Hawaii: Results from waveform inversion in the 10–50 s band | |
CN108181652B (zh) | 一种海底节点地震资料上下行波场数值模拟方法 | |
CN103713315B (zh) | 一种地震各向异性参数全波形反演方法及装置 | |
CN111562611B (zh) | 基于波动方程驱动的半监督深度学习地震数据反演方法 | |
WO2023087451A1 (zh) | 基于观测数据自编码的多尺度无监督地震波速反演方法 | |
CN105629298A (zh) | 确定地震节点位置的方法和计算机*** | |
CN1630827A (zh) | 用于倾斜断裂的层剥转换反射波形 | |
CN103513277B (zh) | 一种地震地层裂隙裂缝密度反演方法及*** | |
CN102253415A (zh) | 基于裂缝等效介质模型的地震响应模式建立方法 | |
CN107817516B (zh) | 基于初至波信息的近地表建模方法及*** | |
CN102239428A (zh) | 用于内部多次反射预测的波场外推建模 | |
NO20140261A1 (no) | Systemer og fremgangsmåter for å fjerne innsamlingsrelaterte effekter fra seismiske data | |
CN103630931A (zh) | 从近场测量和建模假想特征计算假想源特征的方法和*** | |
CN106569262B (zh) | 低频地震数据缺失下的背景速度模型重构方法 | |
MX2014008500A (es) | Fuente de onda de cizalla para vsp y exploracion sismica de superficie. | |
CN111352160A (zh) | 一种海底地震仪自动重定位装置及方法 | |
US11209559B2 (en) | Method and system for analyzing a borehole using passive acoustic logging | |
CN104199088A (zh) | 一种提取入射角道集的方法及*** | |
Bird et al. | Topographically scattered infrasound waves observed on microbarometer arrays in the lower stratosphere | |
Tylor-Jones et al. | A Practical guide to seismic reservoir characterization | |
King | Air-gun bubble-ghost interactions | |
CN104483702A (zh) | 一种适用于非均匀运动水体的地震正演模拟方法 | |
CN109581484A (zh) | 地震纵波分量多次波总照明度指标分析方法及装置 | |
US20220026593A1 (en) | Implicit property modeling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |