CN103018779B - 一种海上地震勘探气枪震源子波模拟方法及系统 - Google Patents

一种海上地震勘探气枪震源子波模拟方法及系统 Download PDF

Info

Publication number
CN103018779B
CN103018779B CN201210507082.XA CN201210507082A CN103018779B CN 103018779 B CN103018779 B CN 103018779B CN 201210507082 A CN201210507082 A CN 201210507082A CN 103018779 B CN103018779 B CN 103018779B
Authority
CN
China
Prior art keywords
bubble
delta
unit
initial
gas
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201210507082.XA
Other languages
English (en)
Other versions
CN103018779A (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 Beijing
Original Assignee
China University of Petroleum Beijing
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 Beijing filed Critical China University of Petroleum Beijing
Priority to CN201210507082.XA priority Critical patent/CN103018779B/zh
Publication of CN103018779A publication Critical patent/CN103018779A/zh
Application granted granted Critical
Publication of CN103018779B publication Critical patent/CN103018779B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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;
(6)由气泡产生的压力波w(iΔt),i=0,1,…n,模拟其虚反射压力波 w g ( i Δ t - 2 h c ) , i = 0 , 1 , ... n , 单位为Pa,
w g ( i &Delta; t - 2 h c ) = 0 i &Delta; t < 2 h c R &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h c
其中,R为海水的反射系数,R=-1,无量纲,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(7)由气泡产生的压力波w(iΔt),i=0,1,…n及虚反射压力波叠加合成气枪子波W(iΔt),i=0,1,…n,单位Bar·m;
所述(5)中,考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n,具体步骤包括:
(5.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(iΔt),i=0,1,…n,其中,R(iΔt)是随时间变化的气泡半径,初值为 R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 , 单位为m,是气泡壁的速度,初值为 R &CenterDot; ( 0 ) = 0 , 单位为m/s;Vg为气枪枪室的容积,单位为in3;rg为枪体半径,单位为m;
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa,
p ( i &Delta; t ) = m ( i &Delta; t ) R G T ( i &Delta; t ) V ( i &Delta; t ) + P V - 2 &sigma; R ( i &Delta; t ) - 4 &mu; ( i &Delta; t ) R ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
(5.3)考虑虚反射压力场对气泡振动的影响,由第i时刻的气泡半径R(iΔt)和气泡壁的速度模拟第i时刻气泡产生的压力波w(iΔt)和虚反射压力场单位均为Pa,具体过程为:
首先由气泡半径R(iΔt)和气泡壁的速度得到中间函数f(iΔt);
f ( i &Delta; t ) = R 2 ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
再利用中间函数f(iΔt)及其对时间导数f′(iΔt),模拟气泡产生的压力波w(iΔt),i=0,1,…n和虚反射压力场
w ( i &Delta; t ) = &rho; &lsqb; - ( - f &prime; ( i &Delta; t ) ) - 1 2 ( f &prime; ( i &Delta; t ) c + f ( i &Delta; t ) ) 2 &rsqb; ,
p g ( i &Delta; t - 2 h c ) = 0 , i &Delta; t < 2 h v - 1 1 + 2 h &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h v ,
其中,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;ρ为海水密度,单位为Kg/m3
(5.4)由以上结果计算i+1时刻的气泡壁加速度单位为m/s2
R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) = - 1 &rho; &lsqb; p &infin; + p g ( i &Delta; t - 2 h c ) - p ( i &Delta; t ) - R ( i &Delta; t ) c p &CenterDot; ( i &Delta; t ) &rsqb; - 3 2 R &CenterDot; 2 ( i &Delta; t ) &lsqb; 1 - 4 R &CenterDot; ( i &Delta; t ) 3 c &rsqb; R ( i &Delta; t ) &lsqb; 1 - 2 R &CenterDot; ( i &Delta; t ) c &rsqb; ,
其中是气泡壁压力p(iΔt)的变化率,ρ为海水密度,单位为Kg/m3,p为初始静水压强,单位为Pa;
(5.5)由i+1时刻的气泡壁的加速度计算该时刻的气泡壁速度和气泡半径
R &CenterDot; ( ( i + 1 ) &Delta; t ) = R &CenterDot; ( i &Delta; t ) + R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ,
R ( ( i + 1 ) &Delta; t ) = R ( i &Delta; t ) + R &CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ;
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w[(i+1)Δt],由此迭代方式得到所有时刻的压力波w(iΔt),i=0,1,…n。
可选的,在本发明一实施例中,所述(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)β;Vg为气枪枪室的容积,单位为in3,τ0为与枪体容量无关的常数;
(2.2)由海水温度Tw和气枪枪室的初始压强pg,计算气枪枪室的初始温度Tg,单位为K,
T g = T w ( 1 + p g p c ) ,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中气体的初始摩尔质量mg和气泡平衡时气泡内气体的摩尔质量mg′,单位均为mol,
m g = p g V g T g R G ,
其中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;
m ( 0 ) = p &infin; V g T w R G ;
(2.6)由气泡平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
r 0 = ( 3 m g &prime; R G T w 4 &pi;p &infin; ) 1 3 ,
rg=r0×0.23;
(2.7)由气枪枪室的容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 ;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3
V ( 0 ) = 4 3 &pi; R ( 0 ) 3 - 4 3 &pi;r g 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是时间采样个数;
模拟单元,用于由气泡内的温度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;
并用于由气泡产生的压力波w(iΔt),i=0,1,…n,模拟其虚反射压力波 w g ( i &Delta; t - 2 h c ) , i = 0 , 1 , ... n , 单位为Pa,
w g ( i &Delta; t - 2 h c ) = 0 i &Delta; t < 2 h c R &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h c
其中,R为海水的反射系数,R=-1,无量纲,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms
合成单元,用于由气泡产生的压力波w(iΔt),i=0,1,…n及虚反射压力波叠加合成气枪子波W(iΔt),i=0,1,…n,单位Bar·m;
所述模拟单元,进一步用于考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(iΔt),i=0,1,…n,具体包括:
(5.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(iΔt),i=0,1,…n,其中,R(iΔt)是随时间变化的气泡半径,初值为 R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 , 单位为m,是气泡壁的速度,初值为 R &CenterDot; ( 0 ) = 0 , 单位为m/s;Vg为气枪枪室的容积,单位为in3;rg为枪体半径,单位为m;
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa,
p ( i &Delta; t ) = m ( i &Delta; t ) R G T ( i &Delta; t ) V ( i &Delta; t ) + P V - 2 &sigma; R ( i &Delta; t ) - 4 &mu; ( i &Delta; t ) R ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
(5.3)考虑虚反射压力场对气泡振动的影响,由第i时刻的气泡半径R(iΔt)和气泡壁的速度模拟第i时刻气泡产生的压力波w(iΔt)和虚反射压力场单位均为Pa,具体过程为:
首先由气泡半径R(iΔt)和气泡壁的速度得到中间函数f(iΔt);
f ( i &Delta; t ) = R 2 ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
再利用中间函数f(iΔt)及其对时间的导数f′(iΔt),模拟气泡产生的压力波w(iΔt),i=0,1,…n和虚反射压力场
w ( i &Delta; t ) = &rho; &lsqb; - ( - f &prime; ( i &Delta; t ) ) - 1 2 ( f &prime; ( i &Delta; t ) c + f ( i &Delta; t ) ) 2 &rsqb; ,
p g ( i &Delta; t - 2 h c ) = { 0 , i &Delta; t < 2 h v - 1 1 + 2 h &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h v ,
其中,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;ρ为海水密度,单位为Kg/m3
(5.4)由以上结果计算i+1时刻的气泡壁加速度单位为m/s2
R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) = - 1 &rho; &lsqb; p &infin; + p g ( i &Delta; t - 2 h c ) - p ( i &Delta; t ) - R ( i &Delta; t ) c p &CenterDot; ( i &Delta; t ) &rsqb; - 3 2 R &CenterDot; 2 ( i &Delta; t ) &lsqb; 1 - 4 R &CenterDot; ( i &Delta; t ) 3 c &rsqb; R ( i &Delta; t ) &lsqb; 1 - 2 R &CenterDot; ( i &Delta; t ) c &rsqb; ,
其中是气泡壁压力p(iΔt)的变化率,ρ为海水密度,单位为Kg/m3,p为初始静水压强,单位为Pa;
(5.5)由i+1时刻的气泡壁的加速度计算该时刻的气泡壁速度和气泡半径R((i+1)Δt);
R &CenterDot; ( ( i + 1 ) &Delta; t ) = R &CenterDot; ( i &Delta; t ) + R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ,
R ( ( i + 1 ) &Delta; t ) = R ( i &Delta; t ) + R &CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ;
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w[(i+1)Δt],由此迭代方式得到所有时刻的压力波w(iΔt),i=0,1,…n。
可选的,在本发明一实施例中,所述测量单元,进一步用于测量并获取海水及其气枪装置的基本参数,包括:
(1.1)测量海水温度Tw,单位为K,和海水密度ρ,单位为Kg/m3
(1.2)测量气枪装置的下列参数:
a.枪室的工作压强pg,单位为psi;
b.气枪枪室的容积Vg,单位为in3
c.枪体的释放效率η,无量纲;
d.与枪体容量无关的常数τ0和节流指数β,无量纲;
e.气枪沉放深度h,单位为m。
可选的,在本发明一实施例中,所述计算单元,进一步用于利用所述基本参数计算:枪体节流常数,气枪枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径,包括:
(2.1)由下式计算枪体节流常数τ,无量纲;
τ=τ0(Vg)β;Vg为气枪枪室的容积,单位为in3,τ0为与枪体容量无关的常数;
(2.2)由海水温度Tw和气枪枪室的初始压强pg,计算气枪枪室的初始温度Tg,单位为K,
T g = T w ( 1 + p g p c ) ,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中气体的初始摩尔质量mg和气泡平衡时气泡内气体的摩尔质量mg′,单位均为mol,
m g = p g V g T g R G ,
其中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;
m ( 0 ) = p &infin; V g T w R G ;
(2.6)由气泡平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
r 0 = ( 3 m g &prime; R G T w 4 &pi;p &infin; ) 1 3 ,
rg=r0×0.23;
(2.7)由气枪枪室的容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 ;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3
V ( 0 ) = 4 3 &pi; R ( 0 ) 3 - 4 3 &pi;r g 3 .
上述技术方案具有如下有益效果:因为采用综合考虑流体粘度随温度变化以及虚反射对气泡周围压力的影响的技术手段,所以达到了如下的技术效果:建立了一种更加完善的气枪子波模拟系统,增强模拟子波与实测子波的一致性,为海上地震勘探提供基础数据,提高海上地震勘探精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图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,模拟其虚反射压力波 w g ( i &Delta; t - 2 h c ) , i = 0 , 1 , ... n , 单位为Pa,
w g ( i &Delta; t - 2 h c ) = 0 i &Delta; t < 2 h c R &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h c
其中,R为海水的反射系数,R=-1,无量纲,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(7)由气泡产生的压力波w(iΔt),i=0,1,…n及虚反射压力波叠加合成气枪子波W(iΔt),i=0,1,…n,单位Bar·m。
可选的,在本发明一实施例中,所述(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,
T g = T w ( 1 + p g p c ) ,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol,
m g = p g V g T g R G ,
其中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;
m ( 0 ) = p &infin; V g T w R G ;
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
r 0 = ( 3 m g &prime; R G T w 4 &pi;p &infin; ) 1 3 ,
rg=r0×0.23;
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 ;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3
V ( 0 ) = 4 3 &pi; R ( 0 ) 3 - 4 3 &pi;r g 3 .
可选的,在本发明一实施例中,所述(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,
p ( i &Delta; t ) = m ( i &Delta; t ) R G T ( i &Delta; t ) V ( i &Delta; t ) + P V - 2 &sigma; R ( i &Delta; t ) - 4 &mu; ( i &Delta; t ) R ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
(5.3)考虑虚反射压力场对气泡振动的影响,由第i时刻的气泡半径R(iΔt)和气泡壁的速度模拟第i时刻气泡产生的压力波w(iΔt)和虚反射压力场单位均为Pa,具体过程为:
首先由气泡半径R(iΔt)和气泡壁的速度得到中间函数f(iΔt);
f ( i &Delta; t ) = R 2 ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
再利用中间函数f(iΔt)及其对时间的导数f′(iΔt),模拟气泡产生的压力波w(iΔt),i=0,1,…n和虚反射压力场
w ( i &Delta; t ) = &rho; &lsqb; - ( - f &prime; ( i &Delta; t ) ) - 1 2 ( f &prime; ( i &Delta; t ) c + f ( i &Delta; t ) ) 2 &rsqb; ,
p g ( i &Delta; t - 2 h c ) = 0 , i &Delta; t < 2 h v - 1 1 + 2 h &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h v ,
其中,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.4)由以上结果计算i+1时刻的气泡壁加速度单位为m/s2
R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) = - 1 &rho; &lsqb; p &infin; + p g ( i &Delta; t - 2 h c ) - p ( i &Delta; t ) - R ( i &Delta; t ) c p &CenterDot; ( i &Delta; t ) &rsqb; - 3 2 R &CenterDot; 2 ( i &Delta; t ) &lsqb; 1 - 4 R &CenterDot; ( i &Delta; t ) 3 c &rsqb; R ( i &Delta; t ) &lsqb; 1 - 2 R &CenterDot; ( i &Delta; t ) c &rsqb; ,
其中是气泡壁压力p(iΔt)的变化率;
(5.5)由i+1时刻的气泡壁的加速度计算该时刻的气泡壁速度和气泡半径R((i+1)Δt);
R &CenterDot; ( ( i + 1 ) &Delta; t ) = R &CenterDot; ( i &Delta; t ) + R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ,
R ( ( i + 1 ) &Delta; t ) = R ( i &Delta; t ) + R &CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ;
(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;
并用于由气泡产生的压力波w(iΔt),i=0,1,…n,模拟其虚反射压力波 w g ( i &Delta; t - 2 h c ) , i = 0 , 1 , ... n , 单位为Pa,
w g ( i &Delta; t - 2 h c ) = 0 i &Delta; t < 2 h c R &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h c
其中,R为海水的反射系数,R=-1,无量纲,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
合成单元25,用于由气泡产生的压力波w(iΔt),i=0,1,…n及虚反射压力波叠加合成气枪子波W(iΔt),i=0,1,…n,单位Bar·m。
可选的,在本发明一实施例中,所述测量单元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,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol,
m g = p g V g T g R G ,
其中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;
m ( 0 ) = p &infin; V g T w R G ;
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
r 0 = ( 3 m g &prime; R G T w 4 &pi;p &infin; ) 1 3 ,
rg=r0×0.23;
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 ;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3
V ( 0 ) = 4 3 &pi; R ( 0 ) 3 - 4 3 &pi;r g 3 .
可选的,在本发明一实施例中,所述模拟单元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,
p ( i &Delta; t ) = m ( i &Delta; t ) R G T ( i &Delta; t ) V ( i &Delta; t ) + P V - 2 &sigma; R ( i &Delta; t ) - 4 &mu; ( i &Delta; t ) R ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
(5.3)考虑虚反射压力场对气泡振动的影响,由第i时刻的气泡半径R(iΔt)和气泡壁的速度模拟第i时刻气泡产生的压力波w(iΔt)和虚反射压力场单位均为Pa,具体过程为:
首先由气泡半径R(iΔt)和气泡壁的速度得到中间函数f(iΔt);
f ( i &Delta; t ) = R 2 ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
再利用中间函数f(iΔt)及其对时间的导数f′(iΔt),模拟气泡产生的压力波w(iΔt),i=0,1,…n和虚反射压力场
w ( i &Delta; t ) = &rho; &lsqb; - ( - f &prime; ( i &Delta; t ) ) - 1 2 ( f &prime; ( i &Delta; t ) c + f ( i &Delta; t ) ) 2 &rsqb; ,
p g ( i &Delta; t - 2 h c ) = { 0 , i &Delta; t < 2 h v - 1 1 + 2 h &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h v ,
其中,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.4)由以上结果计算i+1时刻的气泡壁加速度单位为m/s2
R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) = - 1 &rho; &lsqb; p &infin; + p g ( i &Delta; t - 2 h c ) - p ( i &Delta; t ) - R ( i &Delta; t ) c p &CenterDot; ( i &Delta; t ) &rsqb; - 3 2 R &CenterDot; 2 ( i &Delta; t ) &lsqb; 1 - 4 R &CenterDot; ( i &Delta; t ) 3 c &rsqb; R ( i &Delta; t ) &lsqb; 1 - 2 R &CenterDot; ( i &Delta; t ) c &rsqb; ,
其中是气泡壁压力p(iΔt)的变化率;
(5.5)由i+1时刻的气泡壁的加速度计算该时刻的气泡壁速度和气泡半径R((i+1)Δt);
R &CenterDot; ( ( i + 1 ) &Delta; t ) = R &CenterDot; ( i &Delta; t ) + R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ,
R ( ( i + 1 ) &Delta; t ) = R ( i &Delta; t ) + R &CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ;
(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
T g = T w ( 1 + p g p c ) = 330.2190 K
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中初始气体的摩尔质量mg和平衡时气泡内的摩尔质量mg′,单位均为mol
m g = p g V g T g R G = 3.3406 m o 1
其中RG为普氏常量,其值为8.2;
mg′=mgη=2.7784mol
(2.4)由海水密度ρ和气枪沉放深度h,计算初始静水压强p,单位为Pa
p=p0+ρgh=135900Pa
其中大气压强p0=101300Pa,重力加速度g=9.8m/s2
(2.5)由初始静水压强p,计算气泡内气体的初始摩尔质量m(0),单位为mol;
m ( 0 ) = p &infin; V g T w R G = 0.04 m o 1
(2.6)由平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
r 0 = ( 3 m g &prime; R G T w 4 &pi;p &infin; ) 1 3 = 0.2214 m
rg=r0×0.23=0.0509m
(2.7)由枪体容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 = 0.0661 m
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3
V ( 0 ) = 4 3 &pi; R ( 0 ) 3 - 4 3 &pi;r g 3 = 0.00065548 m 3
(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.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(iΔt),i=0,1,…n,其中,R(iΔt)是随时间变化的气泡半径,初值为 R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 = 0.0661 m , 是气泡半径随时间的变化率,即气泡壁的速度,初值为 R &CenterDot; ( 0 ) = 0 m / s ;
(5.2)由第i时刻气泡内气体的摩尔质量m(iΔt)、气泡内的气体体积V(iΔt)、气泡内的温度T(iΔt)和海水粘度μ(iΔt),模拟第i时刻的气泡壁压力p(iΔt),单位为Pa
p ( i &Delta; t ) = m ( i &Delta; t ) R G T ( i &Delta; t ) V ( i &Delta; t ) + P V 2 &sigma; R ( i &Delta; t ) - 4 &mu; ( i &Delta; t ) R ( i &Delta; t ) R &CenterDot; ( i &Delta; t )
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
(5.3)考虑虚反射压力场对气泡振动的影响,由第i时刻的气泡半径R(iΔt)和气泡壁的速度模拟第i时刻气泡产生的压力波w(iΔt)和虚反射压力场单位均为Pa,具体过程为:
首先由气泡半径R(iΔt)和气泡壁的速度得到中间函数f(iΔt);
f ( i &Delta; t ) = R 2 ( i &Delta; t ) R &CenterDot; ( i &Delta; t )
再利用中间函数f(iΔt)及其对时间的导数f′(iΔt),模拟气泡产生的压力波w(iΔt)和虚反射压力场
w ( i &Delta; t ) = &rho; &lsqb; - ( - f &prime; ( i &Delta; t ) ) - 1 2 ( f &prime; ( i &Delta; t ) c + f ( i &Delta; t ) ) 2 &rsqb;
p g ( i &Delta; t - 2 h c ) = 0 , i &Delta; t < 2 h v - 1 1 + 2 h &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h v
其中,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(5.4)由以上结果计算i+1时刻气泡壁的加速度单位为m/s2
R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) = - 1 &rho; &lsqb; p &infin; + p g ( i &Delta; t - 2 h c ) - p ( i &Delta; t ) - R ( i &Delta; t ) c p &CenterDot; ( i &Delta; t ) &rsqb; - 3 2 R &CenterDot; 2 ( i &Delta; t ) &lsqb; 1 - 4 R &CenterDot; ( i &Delta; t ) 3 c &rsqb; R ( i &Delta; t ) &lsqb; 1 - 2 R &CenterDot; ( i &Delta; t ) c &rsqb;
其中是气泡壁压力p(iΔt)的变化率;
(5.5)由第i+1时刻的气泡壁加速度计算该时刻的气泡壁速度和气泡半径R((i+1)Δt);
R &CenterDot; ( ( i + 1 ) &Delta; t ) = R &CenterDot; ( i &Delta; t ) + R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t
R ( ( i + 1 ) &Delta; t ) = R ( i &Delta; t ) + R &CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t
(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所示的气泡半径随时间的变化率即气泡壁速度随时间的变化情况,如图10所示的气泡壁加速度随时间的变化情况,如图11所示的气泡壁压力p(iΔt),i=0,1,…n随时间的变化情况。
(6)由气泡产生的压力波w(iΔt),i=0,1,…n,模拟如图12所示虚反射压力波 w g ( i &Delta; t - 2 h c ) , i = 0 , 1 , ... n , 单位为Pa
w g ( i &Delta; t - 2 h c ) = 0 i &Delta; t < 2 h c R &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h c
其中,R为海水的反射系数,R=-1,无量纲,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(7)由气泡产生的压力波w(iΔt),i=0,1,…n及虚反射压力波叠加合成如图13所示气枪子波W(iΔt),i=0,1,…n,单位Bar·m。
例如:本实施例为某气枪震源海上地震勘探工区的应用实例,其中气枪容量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 (6)

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,模拟其虚反射压力波 w g ( i &Delta; t - 2 h c ) , i = 0 , 1 , ... n , 单位为Pa,
w g ( i &Delta; t - 2 h c ) = 0 i &Delta; t < 2 h c R &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h c
其中,R为海水的反射系数,R=-1,无量纲,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
(7)由气泡产生的压力波w(i△t),i=0,1,…n及虚反射压力波叠加合成气枪子波W(i△t),i=0,1,…n,单位Bar·m;
所述(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;Vg为气枪枪室的容积,单位为in3;rg为枪体半径,单位为m;
(5.2)由第i时刻气泡内气体的摩尔质量m(i△t)、气泡内的气体体积V(i△t)、气泡内的温度T(i△t)和海水粘度μ(i△t),模拟第i时刻的气泡壁压力p(i△t),单位为Pa,
p ( i &Delta; t ) = m ( i &Delta; t ) R G T ( i &Delta; t ) V ( i &Delta; t ) + P V - 2 &sigma; R ( i &Delta; t ) - 4 &mu; ( i &Delta; t ) R ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
(5.3)考虑虚反射压力场对气泡振动的影响,由第i时刻的气泡半径R(i△t)和气泡壁的速度模拟第i时刻气泡产生的压力波w(i△t)和虚反射压力场单位均为Pa,具体过程为:
首先由气泡半径R(i△t)和气泡壁的速度得到中间函数f(i△t);
f ( i &Delta; t ) = R 2 ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
再利用中间函数f(i△t)及其对时间的导数f′(i△t),模拟气泡产生的压力波w(i△t),i=0,1,…n和虚反射压力场
w ( i &Delta; t ) = &rho; &lsqb; - ( - f &prime; ( i &Delta; t ) ) - 1 2 ( f &prime; ( i &Delta; t ) c + f ( i &Delta; t ) ) 2 &rsqb; ,
p g ( i &Delta; t - 2 h c ) = 0 , i &Delta; t < 2 h v - 1 1 + 2 h &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h v ,
其中,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;ρ为海水密度,单位为Kg/m3
(5.4)由以上结果计算i+1时刻的气泡壁加速度单位为m/s2
R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) = - 1 &rho; &lsqb; p &infin; + p g ( i &Delta; t - 2 h c ) - p ( i &Delta; t ) - R ( i &Delta; t ) c p &CenterDot; ( i &Delta; t ) &rsqb; - 3 2 R &CenterDot; 2 ( i &Delta; t ) &lsqb; 1 - 4 R &CenterDot; ( i &Delta; t ) 3 c &rsqb; R ( i &Delta; t ) &lsqb; 1 - 2 R &CenterDot; ( i &Delta; t ) c &rsqb; ,
其中是气泡壁压力p(i△t)的变化率,ρ为海水密度,单位为Kg/m3,p为初始静水压强,单位为Pa;
(5.5)由i+1时刻的气泡壁的加速度计算该时刻的气泡壁速度和气泡半径
R &CenterDot; ( ( i + 1 ) &Delta; t ) = R &CenterDot; ( i &Delta; 1 ) + R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ,
R ( ( i + 1 ) &Delta; t ) = R ( i &Delta; 1 ) + R &CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ;
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w[(i+1)△t],由此迭代方式得到所有时刻的压力波w(i△t),i=0,1,…n。
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)β;Vg为气枪枪室的容积,单位为in3,τ0为与枪体容量无关的常数;
(2.2)由海水温度Tw和气枪枪室的初始压强pg,计算气枪枪室的初始温度Tg,单位为K,
T g = T w ( 1 + p g p c ) ,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中气体的初始摩尔质量mg和气泡平衡时气泡内气体的摩尔质量mg′,单位均为mol,
m g = p g V g T g R G ,
其中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;
m ( 0 ) = p &infin; V g T w R G ;
(2.6)由气泡平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
r 0 = ( 3 m g &prime; R G T w 4 &pi;p &infin; ) 1 3 ,
rg=r0×0.23;
(2.7)由气枪枪室的容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 :
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3
V ( 0 ) = 4 3 &pi; R ( 0 ) 3 - 4 3 &pi; r g 3 .
4.一种海上地震勘探气枪震源子波模拟系统,其特征在于,所述海上地震勘探气枪震源子波模拟系统包括:
测量单元,用于测量并获取海水及其气枪装置的基本参数;
计算单元,用于利用所述基本参数计算:枪体节流常数,气枪枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径;
函数单元,用于按照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;
并用于由气泡产生的压力波w(i△t),i=0,1,…n,模拟其虚反射压力波 w g ( i &Delta; t - 2 h c ) , i = 0 , 1 , ... n , 单位为Pa,
w g ( i &Delta; t - 2 h c ) = 0 i &Delta; t < 2 h c R &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h c
其中,R为海水的反射系数,R=-1,无量纲,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;
合成单元,用于由气泡产生的压力波w(i△t),i=0,1,…n及虚反射压力波叠加合成气枪子波W(i△t),i=0,1,…n,单位Bar·m;
所述模拟单元,进一步用于考虑海水粘度和虚反射压力场对气泡振动的影响,模拟气泡产生的压力波w(i△t),i=0,1,…n,具体包括:
(5.1)由初始时刻i=0开始,以递归循环方式模拟气泡产生的压力波w(i△t),i=0,1,…n,其中,R(i△t)是随时间变化的气泡半径,初值为单位为m,是气泡壁的速度,初值为单位为m/s;Vg为气枪枪室的容积,单位为in3;rg为枪体半径,单位为m;
(5.2)由第i时刻气泡内气体的摩尔质量m(i△t)、气泡内的气体体积V(i△t)、气泡内的温度T(i△t)和海水粘度μ(i△t),模拟第i时刻的气泡壁压力p(i△t),单位为Pa,
p ( i &Delta; t ) = m ( i &Delta; t ) R G T ( i &Delta; t ) V ( i &Delta; t ) + P V - 2 &sigma; R ( i &Delta; t ) - 4 &mu; ( i &Delta; t ) R ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
其中,水蒸气压力PV=2337Pa,气泡表面张力σ=0.07061N/m;
(5.3)考虑虚反射压力场对气泡振动的影响,由第i时刻的气泡半径R(i△t)和气泡壁的速度模拟第i时刻气泡产生的压力波w(i△t)和虚反射压力场单位均为Pa,具体过程为:
首先由气泡半径R(i△t)和气泡壁的速度得到中间函数f(i△t);
f ( i &Delta; t ) = R 2 ( i &Delta; t ) R &CenterDot; ( i &Delta; t ) ,
再利用中间函数f(i△t)及其对时间的导数f′(i△t),模拟气泡产生的压力波w(i△t),i=0,1,…n和虚反射压力场
w ( i &Delta; t ) = &rho; &lsqb; - ( - f &prime; ( i &Delta; t ) ) - 1 2 ( f &prime; ( i &Delta; t ) c + f ( i &Delta; t ) ) 2 &rsqb; ,
p g ( i &Delta; t - 2 h c ) = 0 , i &Delta; t < 2 h v - 1 1 + 2 h &times; w ( i &Delta; t ) , i &Delta; t &GreaterEqual; 2 h v ,
其中,h为气枪沉放深度,单位为m,c为海水中的声波速度,单位为m/ms,c=1.483m/ms;ρ为海水密度,单位为Kg/m3
(5.4)由以上结果计算i+1时刻的气泡壁加速度单位为m/s2
R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) = - 1 &rho; &lsqb; p &infin; + p g ( i &Delta; t - 2 h c ) - p ( i &Delta; t ) - R ( i &Delta; t ) c p &CenterDot; ( i &Delta; t ) &rsqb; - 3 2 R &CenterDot; 2 ( i &Delta; t ) &lsqb; 1 - 4 R &CenterDot; ( i &Delta; t ) 3 c &rsqb; R ( i &Delta; t ) &lsqb; 1 - 2 R &CenterDot; ( i &Delta; t ) c &rsqb; ,
其中是气泡壁压力p(i△t)的变化率,ρ为海水密度,单位为Kg/m3,p为初始静水压强,单位为Pa;
(5.5)由i+1时刻的气泡壁的加速度计算该时刻的气泡壁速度和气泡半径R((i+1)△t);
R &CenterDot; ( ( i + 1 ) &Delta; t ) = R &CenterDot; ( i &Delta; t ) + R &CenterDot;&CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ,
R ( ( i + 1 ) &Delta; t ) = R ( i &Delta; t ) + R &CenterDot; ( ( i + 1 ) &Delta; t ) &Delta; t ;
(5.6)返回步骤(5.2),计算i+1时刻的气泡产生的压力波w[(i+1)△t],由此迭代方式得到所有时刻的压力波w(i△t),i=0,1,…n。
5.如权利要求4所述海上地震勘探气枪震源子波模拟系统,其特征在于,所述测量单元,进一步用于测量并获取海水及其气枪装置的基本参数,包括:
(1.1)测量海水温度Tw,单位为K,和海水密度ρ,单位为Kg/m3
(1.2)测量气枪装置的下列参数:
a.枪室的工作压强pg,单位为psi;
b.气枪枪室的容积Vg,单位为in3
c.枪体的释放效率η,无量纲;
d.与枪体容量无关的常数τ0和节流指数β,无量纲;
e.气枪沉放深度h,单位为m。
6.如权利要求4所述海上地震勘探气枪震源子波模拟系统,其特征在于,所述计算单元,进一步用于利用所述基本参数计算:枪体节流常数,气枪枪室的初始温度,枪室中气体的初始摩尔质量,气泡内气体的初始摩尔质量,气泡平衡时气泡内气体的摩尔质量,初始静水压强,气泡的初始半径,气泡中气体的初始体积,气泡平衡半径,包括:
(2.1)由下式计算枪体节流常数τ,无量纲;
τ=τ0(Vg)β;Vg为气枪枪室的容积,单位为in3,τ0为与枪体容量无关的常数;
(2.2)由海水温度Tw和气枪枪室的初始压强pg,计算气枪枪室的初始温度Tg,单位为K,
T g = T w ( 1 + p g p c ) ,
其中常量pc=139Mpa;
(2.3)由枪室初始温度Tg和枪体的释放效率η,计算枪室中气体的初始摩尔质量mg和气泡平衡时气泡内气体的摩尔质量mg′,单位均为mol,
m g = p g V g T g R G ,
其中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;
m ( 0 ) = p &infin; V g T w R G ;
(2.6)由气泡平衡时气泡内气体的摩尔质量mg′、海水温度Tw和初始静水压强p,计算气泡的平衡半径r0以及枪体半径rg,单位均为m;
r 0 = ( 3 m g &prime; R G T w 4 &pi;p &infin; ) 1 3 ,
rg=r0×0.23;
(2.7)由气枪枪室的容积Vg和枪体半径rg,计算气泡的初始半径R(0),单位为m;
R ( 0 ) = &lsqb; 3 ( V g + 4 3 &pi;r g 3 ) 4 &pi; &rsqb; 1 3 ;
(2.8)由气泡初始半径R(0)和枪体半径rg计算气泡中气体的初始体积V(0),单位为m3
V ( 0 ) = 4 3 &pi; R ( 0 ) 3 - 4 3 &pi;r g 3 .
CN201210507082.XA 2012-11-30 2012-11-30 一种海上地震勘探气枪震源子波模拟方法及系统 Active CN103018779B (zh)

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 CN103018779A (zh) 2013-04-03
CN103018779B true 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)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108061915A (zh) * 2016-11-09 2018-05-22 中国石油化工股份有限公司 一种气枪子波模拟方法
CN107024714B (zh) * 2017-03-24 2019-02-15 中国石油天然气集团公司 一种实现气枪震源质控的处理方法及装置
CN109239769A (zh) * 2018-11-01 2019-01-18 国家海洋局第二海洋研究所 利用深度和容量组合压制残留气泡的气枪震源设计方法
CN110687617B (zh) * 2019-08-26 2021-06-29 中国海洋大学 一种地震勘探气枪阵列远场子波模拟方法、脸谱评价方法及装置
CN113655519B (zh) * 2021-08-23 2023-10-13 中海石油(中国)有限公司 气枪节流作用系数和气体释放效率参数获取方法及系统
CN114647990A (zh) * 2022-02-22 2022-06-21 哈尔滨工程大学 一种高效高精度计算气枪震源压力子波的方法

Citations (3)

* Cited by examiner, † Cited by third party
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 地球物理维里达斯集团公司 用于反虚反射地震数据的方法和装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
大容量气枪震源的实验与模拟研究;唐杰 等;《中国地震》;20090331;第25卷(第1期);第1-10页 *
李国发 等.海上地震勘探GI枪子波数值模拟.《吉林大学学报(地球科学版)》.2011,第41卷(第1期),第277-282页. *
枪口节流作用和热传导作用对海上勘探气枪子波的影响;李国发 等;《特种油气藏》;20110430;第18卷(第2期);第26-28页 *
气枪震源单枪子波计算机模拟;陈浩林 等;《物探装备》;20021231;第12卷(第4期);第241-244页 *
考虑多种因素影响的海上地震勘探气枪震源单枪子波数值模拟;李国发 等;《应用地球物理》;20100630;第7卷(第2期);第158-165页 *

Also Published As

Publication number Publication date
CN103018779A (zh) 2013-04-03

Similar Documents

Publication Publication Date Title
CN103018779B (zh) 一种海上地震勘探气枪震源子波模拟方法及系统
CN103713315B (zh) 一种地震各向异性参数全波形反演方法及装置
Chouet et al. Seismic source mechanism of degassing bursts at Kilauea Volcano, Hawaii: Results from waveform inversion in the 10–50 s band
CA3122509A1 (en) Machine learning-augmented geophysical inversion
CN103959099B (zh) 使用梯度传感器的波场分离
CN111562611B (zh) 基于波动方程驱动的半监督深度学习地震数据反演方法
CN106526674A (zh) 一种三维全波形反演能量加权梯度预处理方法
CA2946711A1 (en) System and method of building complex earth models
CN106569262B (zh) 低频地震数据缺失下的背景速度模型重构方法
CN103499835A (zh) 初至波波形反演近地表速度模型方法
CN102253415A (zh) 基于裂缝等效介质模型的地震响应模式建立方法
CN107817516B (zh) 基于初至波信息的近地表建模方法及系统
CN103630931B (zh) 从近场测量和建模假想特征计算假想源特征的方法和系统
Orescanin et al. Unsteady high‐pressure flow experiments with applications to explosive volcanic eruptions
CN105093278A (zh) 基于激发主能量优化算法的全波形反演梯度算子的提取方法
MX2014008500A (es) Fuente de onda de cizalla para vsp y exploracion sismica de superficie.
CN104975851B (zh) 用于振幅随炮检距变化道集分析的油藏模型优化方法
CN104965222A (zh) 三维纵波阻抗全波形反演方法及装置
US11209559B2 (en) Method and system for analyzing a borehole using passive acoustic logging
CN104483702B (zh) 一种适用于非均匀运动水体的地震正演模拟方法
King Air-gun bubble-ghost interactions
King et al. Boundary conditions for simulations of oscillating bubbles using the non-linear acoustic approximation
CN109581484A (zh) 地震纵波分量多次波总照明度指标分析方法及装置
CN102012517B (zh) 地下介质成像方法和装置
CN109001804A (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
C14 Grant of patent or utility model
GR01 Patent grant