CN116992616A - 一种构造尖脉冲子波提高气枪震源低频能量的方法及应用 - Google Patents
一种构造尖脉冲子波提高气枪震源低频能量的方法及应用 Download PDFInfo
- Publication number
- CN116992616A CN116992616A CN202211664560.8A CN202211664560A CN116992616A CN 116992616 A CN116992616 A CN 116992616A CN 202211664560 A CN202211664560 A CN 202211664560A CN 116992616 A CN116992616 A CN 116992616A
- Authority
- CN
- China
- Prior art keywords
- bubble
- air gun
- wavelet
- air
- pressure
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 230000008569 process Effects 0.000 claims abstract description 26
- 230000001934 delay Effects 0.000 claims abstract description 4
- 230000002706 hydrostatic effect Effects 0.000 claims description 24
- 230000008859 change Effects 0.000 claims description 23
- 230000005284 excitation Effects 0.000 claims description 22
- 238000004088 simulation Methods 0.000 claims description 20
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 10
- 230000003993 interaction Effects 0.000 claims description 10
- 238000010183 spectrum analysis Methods 0.000 claims description 10
- 238000012546 transfer Methods 0.000 claims description 9
- 230000000630 rising effect Effects 0.000 claims description 8
- 230000002708 enhancing effect Effects 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012804 iterative process Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 239000007787 solid Substances 0.000 claims 2
- 230000000977 initiatory effect Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 description 9
- 238000003384 imaging method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 230000003111 delayed effect Effects 0.000 description 3
- 239000013535 sea water Substances 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
- G01V1/133—Generating seismic energy using fluidic driving means, e.g. highly pressurised fluids; using implosion
- G01V1/137—Generating seismic energy using fluidic driving means, e.g. highly pressurised fluids; using implosion which fluid escapes from the generator in a pulsating manner, e.g. for generating bursts, airguns
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Acoustics & Sound (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于海洋地震勘探过程中的气枪震源数据识别技术领域,公开了一种构造尖脉冲子波提高气枪震源低频能量的方法及应用。根据气枪不同沉放深度与子波信号接收时间的关系,通过调节不同容量气枪的沉放深度,使阵列中部分气枪的主脉冲峰值时刻进行延迟,构造主脉冲波形较窄的尖脉冲震源子波,增大气枪震源低频能量。本发明以有规律的改变不同容量气枪的沉放深度,使子波接收信号延迟,构造尖脉冲子波,产生高分辨率的下传能力强的地震子波以适用浅海中深层地质目标勘探的相关质量要求。
Description
技术领域
本发明属于海洋地震勘探过程中的气枪震源数据识别技术领域,尤其涉及一种构造尖脉冲子波提高气枪震源低频能量的方法及应用。
背景技术
气枪震源是海洋油气资源勘探和地质调查中的关键设备,尤其是针对海洋中深层目标勘探,震源中的低频能量是至关重要的。随着海洋中深层油气资源的大量开采,复杂油气藏条件目标的勘探对气枪震源低频能量的需求越来越高。目前国内外主要通过子阵平面、立体组合的方法,利用组合子波能量,以达到拓宽震源子波频带的目的。
然而通过子阵的常规组合对于提高震源高频效果明显,但是震源频带向低频延展的效果并不理想,如何保证震源子波高频的同时延展低频,是地球物理学领域亟待解决的问题。如何在保证高频的同时,压制海面虚反射,延展低频提高地震子波的穿透能力,是本发明基于先存技术所要提高的焦点。
通过上述分析,现有技术存在的问题及缺陷为:
(1)现有技术没有以有规律的改变不同容量气枪的沉放深度,使子波接收信号延迟,构造尖脉冲子波,不能高分辨率的下传能力强的地震子波以适用于浅海中深层地质目标勘探。
(2)现有技术气枪震源针对中深层地质目标探测数据不精确,效果差。
发明内容
为克服相关技术中存在的问题,本发明公开实施例提供了一种构造尖脉冲子波提高气枪震源低频能量的方法及应用,本发明用于海洋野外宽频立体观测系统,勘探采集近海浅水中深层地震反射信号并用于海洋地质调查和油气勘探。
所述技术方案如下:通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,该方法根据气枪不同沉放深度与子波信号接收时间的关系,通过调节不同容量气枪的沉放深度,使阵列中部分气枪的主脉冲峰值时刻进行延迟,构造主脉冲波形较窄的尖脉冲震源子波,增大气枪震源低频能量;具体包括以下步骤:
S1、根据范德瓦尔斯非理想气体气枪子波模型模拟不同容量的气枪子波,设定模型初始条件;
S2、根据步骤S1设定的初始条件执行模拟过程;
S3、对模拟出的气枪子波进行分析,统计不同容量气枪子波从激发到达到主脉冲峰值的时间ti;
S4、根据模拟结果,容量越小ti越小,假设最小ti为t0,计算不同容量气枪从激发到达到主脉冲峰值的时间与t0的差值Δti;
S5、依据步骤S4得到的Δti,计算不同容量气枪的沉放深度,沉放深度使得不同容量气枪子波信号的接收延迟时间为Δti;随后进行立体阵列尖脉冲子波模拟;
S6、对步骤S5得到的气枪阵列子波做频谱分析。
在步骤S1中,范德瓦尔斯非理想气体气枪子波模型的方程表达为:
式中,a、b均为范德瓦尔斯常数,a=0.1404m6·Pa·mol-2,b=3.764×10-5m3·mol-1,Tg为有效热力学温度,Rg为普适气体常数,mg为气体质量,Vg为体积;
有效热力学温度Tg取决于气室内的高压气体:
Tg=Tw(1+Pg/Pc)(2)
式中,Pc=139MPa;
气枪激发过程中,根据能量守恒定律,气泡热传播损失和气泡物质量的传递获得的能量必须与气泡内能的变化平衡,有:
式中,T是气泡温度,P是气泡压力,mb是气泡内气体的物质量,U=CmmbT表示气泡的内能,Cm是定压比热容,dQ/dt是通过气泡壁的传热速率,传热系数k是通过对模型与试验数据拟合确定的,气泡热损失率表示为:
式中,ΔT=Tb-Tw是气泡温度Tb和周围水温Tw之间的温度差,R是气泡半径;
利用范德瓦尔斯非理想气体方程,非理想气体的内能为气体温度和体积的函数:
全微分方程表示为:
热力学第一定律转化为:
式中,Rg=Cm-Cm;
为获取气体物质的量的变化率引入节流常数τ,影响气体通过气枪端口的速率;式中Vg是气室容量,mg是气室内气体的物质量,Pg是气枪压力,Pb是气泡压力;
对于实际应用的气枪,高压气体释放到水中的速度和总量受到端口大小、端口打开时间参数的控制,影响气枪子波性能;根据幂次定律,表示为:
式中,τ0是与容量无关的端口节流常数,ζ是通过与实测数据的对比确定的节流幂律指数;根据测量和计算结果,在任意给定的时间内,通过气枪端口的气体流量取决于气枪内外的压力差,气体释放的速率表示为:
式中,mb是释放到气泡中的气体物质的量,mg|t=0是气室内气体总量,η是气泡中气体量与总量的比;
气泡壁的运动方程表示为:
式中,R是气泡半径,u和分别是气泡壁的速度和加速度,c是声波在流体介质中的速度,/>是气泡壁的焓差,ρ∞是无穷远处的静水密度,Pb是气泡压力,P∞是无穷远处静水压;气泡的静水压力在气泡由于浮力上升的过程中发生改变,气泡上升过程的气泡垂直上升速度的表达式为:
式中,z是气泡深度,g是重力加速度常数,R是气泡半径,静水压力P∞的表达式为:
式中,Patm是标准大气压,zg是气枪深度;在距离气枪1m处,气枪子波信号表示为:
在低频时,气泡间的相互作用不忽视;气泡间的这种相互作用是对流体静水压的调节;气泡间的相互作用使气泡周围的压力变化;对于地震波长,气泡是一个点,任何任意气泡周围的压力场均是静水压力加上气泡产生的时变压力场的叠加;第i各气泡处的有效静水压力为;
式中,P∞是静水压,∑k≠iΔPik是气枪阵列中所有其他气枪的压力贡献之和,ΔPik是第k个气泡引起的对第i个气泡的静水压力扰动,以及第k个气泡引起的对第i个气泡的延时和按距离标度的压力特征:
式中,rik表示第i个气泡与第k个气泡之间的气泡间距。
在步骤S1中,初始条件具体为:
步骤1.1、条件1:将气枪压力的初始值Pg|t=0设置为工作压力;
步骤1.2、条件2:气泡内初始温度设置为Tg=Tw(1+Pg/Pc);
步骤1.3、条件3:气泡初始体积Vb|t=0=Vg,初始半径为
步骤1.4、条件4:气泡壁的初始速度为u=0;
步骤1.5、条件5:气泡初始压力Pb|t=0=P∞,初始温度为水温Tw=18°,气泡内的初始物质量为
步骤1.6、条件6:设置各条气枪的放置位置(x,y,z)。
在步骤S2中,根据步骤S1设定的初始条件执行模拟过程,具体为:
步骤2.1、输入范德瓦尔斯非理想气体气枪子波模型的初始条件;
步骤2.2、开始时间循环并计算t=k时刻的气泡体积
步骤2.3、利用方程(1)计算t=k时刻的气泡压力,Pb;
步骤2.4、通过方程(4)计算气泡热损失率,
步骤2.5、通过方程(9)计算气体的释放速率,
步骤2.6、计算t=k时刻的气泡体积变化率,
步骤2.7、通过方程(7)计算t=k时刻气泡内的温度变化率,
步骤2.8、计算气泡壁的焓差,
步骤2.9、通过对方程(1)求关于时间t的微分得到气泡压力的变化率,
步骤2.10、对焓差求关于时间t的微分得到,
步骤2.11、通过方程(10)计算t=k时刻的气泡壁的速度变化率,既气泡壁的加速度;
步骤2.12、对求关于时间t的微分得到,/>
步骤2.13、由于气枪子波模拟是一个迭代过程,通过二阶泰勒级数展开,可以得到气泡壁半径、气泡壁速度、气体温度以及气泡内气体的物质量:and
步骤2.14、将气泡压力表示为焓、气泡壁速度和以及气泡半径的函数:R0为气泡中心到远场点的距离;
步骤2.15、重复步骤(2.1)至(2.14),直到t>tmax;
步骤2.16、计算气枪远场子波声压,包括海面虚反射:
Rs表示海面反射系数,D1是气枪与水听器之间的距离是,D2是气枪的海面镜像与水听器之间的距离,/>是气枪信号经过D1和D2的时间延迟。
在步骤S3中,统计不同容量气枪子波从激发到达到主脉冲峰值的时间ti包括:统计45cu.in\70cu.in\100cu.in\150cu.in\250cu.in等容量气枪激发后,达到主脉冲峰值的时间ti。
在步骤S4中,计算不同容量气枪从激发到达到主脉冲峰值的时间与t0的差值Δti为计算70cu.in\100cu.in\150cu.in\250cu.in容量的气枪激发后达到主脉冲峰值的时间与45cu.in容量的气枪达到主脉冲峰值的时间差值Δti。
在步骤S6中,对步骤S5得到的气枪阵列子波做频谱分析具体包括:计算气枪阵列子波的主脉冲峰值、虚反射值、气泡脉冲峰值,通过傅里叶变换对子波做频谱分析,以最大振幅-6dB为判断有效带宽的标准,求得子波主脉冲的有效带宽,并得到子波主脉冲主频。
本发明的另一目的在于提供一种气枪阵列,实施所述的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,所述气枪阵列包括工作枪33条,空枪6条,总容量4040cu.in;单枪容量与条数分别包括:
6条45cu.in;
4条70cu.in;
10条100cu.in,其中包含2条空枪;
11条150cu.in,其中包含2条空枪;
8条250cu.in,其中包含2条空枪。
本发明的另一目的在于提供一种所述气枪阵列在海洋油气资源勘探气枪震源设备上的应用。
本发明的另一目的在于提供一种所述气枪阵列在地质调查中气枪震源设备上的应用。
结合上述的所有技术方案,本发明所具备的优点及积极效果为:
第一、针对上述现有技术存在的技术问题以及解决该问题的难度,紧密结合本发明的所要保护的技术方案以及研发过程中结果和数据等,详细、深刻地分析本发明技术方案如何解决的技术问题,解决问题之后带来的一些具备创造性的技术效果,具体描述如下:本发明利用气枪沉放深度与子波信号接收时间的关系,通过有规律的调整不同容量气枪的沉放深度,使子波接收信号延迟,构造尖脉冲子波,实现气枪阵列的频带在保证高频的基础上向低频拓展的目的,最终获得高低频丰富,能量下传能力强的针对中深层地质目标的气枪震源。
第二、把技术方案看作一个整体或者从产品的角度,本发明所要保护的技术方案具备的技术效果和优点,具体描述如下:本发明以有规律的改变不同容量气枪的沉放深度,使子波接收信号延迟,构造尖脉冲子波,产生高分辨率的下传能力强的地震子波以适用浅海中深层地质目标勘探的相关质量要求。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本公开的实施例,并与说明书一起用于解释本公开的原理;
图1是本发明实施例提供的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法流程图;
图2是本发明实施例提供的模拟出来的不同容量的气枪单枪子波图;图中:bar m为气枪激发产生的气泡的腔室内的压力,cu.in为气枪容量,立方英寸;
图3是本发明实施例提供的所用的气枪阵列中各容量气枪单枪模拟所得从激发到达到主脉冲峰值时间统计图;
图4是本发明实施例提供的所用气枪阵列的平面图;
图5(a)是本发明实施例提供的模拟出来的尖脉冲气枪子波图;
图5(b)是本发明实施例提供的模拟出来的频谱图;
图6(a)是本发明实施例提供的常规阵列模拟出来的子波图;
图6(b)是本发明实施例提供的利用本发明模拟出来的尖脉冲气枪子波图;
图7为marmousi模型的速度场模型;
图8(a)是本发明实施例提供的利用现有技术雷克子波做的波场成像图;
图8(b)是本发明实施例提供的利用本发明构造的尖脉冲子波做的波场成像和图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其他方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似改进,因此本发明不受下面公开的具体实施的限制。
一、解释说明实施例:
本发明实施例1提供一种通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法包括:
根据气枪不同沉放深度与子波信号接收时间的关系,通过有规律的调节不同容量气枪的沉放深度,使得阵列中部分气枪的主脉冲峰值时刻适当的延迟,构造主脉冲波形较窄的尖脉冲震源子波,实现提高气枪震源低频能量的目的。
实施例1
如图1所示,本发明实施例提供地通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法包括以下步骤:
S101、根据范德瓦尔斯非理想气体气枪子波模型模拟不同容量的气枪子波,设定模型初始条件;
S102、根据步骤S101设定的初始条件执行模拟过程;
S103、对模拟出的气枪子波进行分析,统计不同容量气枪子波从激发到达到主脉冲峰值的时间ti;
S104、根据模拟结果,容量越小ti越小,假设最小ti为t0,计算不同容量气枪从激发到达到主脉冲峰值的时间与t0的差值Δti;
S105、依据步骤S104得到的Δti,计算不同容量气枪的沉放深度,该深度恰好使得不同容量气枪子波信号的接收延迟时间为Δti。随后进行立体阵列尖脉冲子波模拟。
S106、对步骤S105得到的气枪阵列子波做频谱分析。
实施例2
基于本发明实施例1提供地通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,所述步骤S101中的范德瓦尔斯非理想气体气枪子波模型的方程表达为:
式中,a=0.1404m6·Pa·mol-2,b=3.764×10-5m3·mol-1是范德瓦尔斯常数,Tg为有效热力学温度,Rg为普适气体常数,mg为气体质量,Vg为体积。
Laws等认为有效温度Tg取决于气室内的高压气体:
Tg=Tw(1+Pg/Pc)(2)
式中,Pc=139MPa。
气枪激发过程中,高压气体从腔室内喷出,形成气泡,这个过程中热量通过气泡壁向外进行传递,符合开放的热力学系统特征。根据能量守恒定律,气泡热传播损失和气泡物质量的传递获得的能量必须与气泡内能的变化平衡,因此有:
式中,T是气泡温度,P是气泡压力,mb是气泡内气体的物质量,U=CmmbT表示气泡的内能,Cm和Cm分别是定压比热容和定体比热容,dQ/dt是通过气泡壁的传热速率,传热系数k是通过对模型与试验数据拟合确定的,气泡热损失率可以表示为:
式中,ΔT=Tb-Tw是气泡温度Tb和周围水温Tw之间的温度差,R是气泡半径。气体分子的动能受到温度的影响,气体分子的势能受变化的体积的影响。因此,利用范德瓦尔斯非理想气体方程,非理想气体的内能是气体温度和体积的函数:
全微分方程可表示为:
此外,热力学第一定律可以转化为:
式中,Rg=Cm-Cm。
为了推导气体物质的量的变化率引入节流常数τ,其决定了气体通过气枪端口的速率;式中Vg是气室容量,mg是气室内气体的物质量,Pg是气枪压力,Pb是气泡压力。
对于实际应用的气枪,高压气体释放到水中的速度和总量受到端口大小、端口打开时间等参数的控制,进而影响了气枪子波性能。由于端口面积固定,节流常数τ提高了模型与实测数据的一致性,假设不同容量气枪的节流常数只与气室大小有关。根据幂次定律,可以表示为:
式中,τ0是与容量无关的端口节流常数,ζ是通过与实测数据的对比确定的节流幂律指数。根据测量和计算结果,气体逸入水中的气泡可以持续几毫秒。在任意给定的时间内,通过气枪端口的气体流量取决于气枪内外的压力差,从而气体释放的速率可以表示为:
式中,mb是释放到气泡中的气体物质的量,mg|t=0是气室内气体总量,η是气泡中气体量与总量的比。
气泡壁的运动方程可以表示为:
式中,R是气泡半径,u和分别是气泡壁的速度和加速度,c是声波在流体介质中的速度,/>是气泡壁的焓差,ρ∞是无穷远处的静水密度,Pb是气泡压力,P∞是无穷远处静水压。气泡的静水压力在气泡由于浮力上升的过程中发生改变,因此必须考虑气泡的上升。考虑气泡上升过程的气泡垂直上升速度的表达式为:
式中,z是气泡深度,g是重力加速度常数,R是气泡半径,因此,静水压力P∞的表达式为:
式中,Patm是标准大气压,zg是气枪深度。在距离气枪1m处,气枪子波信号可表示为:
在低频时,气泡间的相互作用不可忽视。气泡间的这种相互作用可以看作是对流体静水压的调节。气泡间的相互作用会使气泡周围的压力变化。相对于地震波长,气泡可以看作是一个点,所以,任何任意气泡周围的压力场均是静水压力加上气泡产生的时变压力场的叠加。因此,第i各气泡处的有效静水压力为;
式中,P∞是静水压,∑k≠iΔPik是气枪阵列中所有其他气枪的压力贡献之和,ΔPik是第k个气泡引起的对第i个气泡的静水压力扰动,以及第k个气泡引起的对第i个气泡的延时和按距离标度的压力特征:
式中,rik表示第i个气泡与第k个气泡之间的气泡间距。
实施例2
基于本发明实施例1提供地通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,所述步骤S101中初始条件具体为:
步骤1.1、条件1:将气枪压力的初始值Pg|t=0设置为工作压力;
步骤1.2、条件2:气泡内初始温度设置为Tg=Tw(1+Pg/Pc);
步骤1.3、条件3:气泡初始体积Vb|t=0=Vg,初始半径为
步骤1.4、条件4:气泡壁的初始速度为u=0;
步骤1.5、条件5:气泡初始压力Pb|t=0=P∞,初始温度为水温Tw=18°,气泡内的初始物质量为
步骤1.6、条件6:设置各条气枪的放置位置(x,y,z)。
实施例3
基于本发明实施例1提供地通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,所述步骤S102中根据步骤S101设定的初始条件进行的模拟过程具体为:
步骤2.1、输入范德瓦尔斯非理想气体气枪子波模型的初始条件;
步骤2.2、开始时间循环并计算t=k时刻的气泡体积
步骤2.3、利用方程(1)计算t=k时刻的气泡压力,Pb;
步骤2.4、通过方程(4)计算气泡热损失率,
步骤2.5、通过方程(9)计算气体的释放速率,
步骤2.6、计算t=k时刻的气泡体积变化率,
步骤2.7、通过方程(7)计算t=k时刻气泡内的温度变化率,
步骤2.8、计算气泡壁的焓差,
步骤2.9、通过对方程(1)求关于时间t的微分得到气泡压力的变化率,
步骤2.10、对焓差求关于时间t的微分得到,
步骤2.11、通过方程(10)计算t=k时刻的气泡壁的速度变化率,既气泡壁的加速度;
步骤2.12、对求关于时间t的微分得到,/>
步骤2.13、由于气枪子波模拟是一个迭代过程,通过二阶泰勒级数展开,可以得到气泡壁半径、气泡壁速度、气体温度以及气泡内气体的物质量:and
步骤2.14、将气泡压力表示为焓、气泡壁速度和以及气泡半径的函数:R0为气泡中心到远场点的距离;
步骤2.15、重复步骤(2.1)至(2.14),直到t>tmax;
步骤2.16、计算气枪远场子波声压,包括海面虚反射:
Rs表示海面反射系数,D1是气枪与水听器之间的距离是,D2是气枪的海面镜像与水听器之间的距离,/>是气枪信号经过D1和D2的时间延迟。
实施例4
基于本发明实施例1提供地通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,所述步骤S103中分析步骤S102模拟的气枪阵列子波,统计不同容量气枪从激发到达到主脉冲峰值的时间ti,即统计45cu.in\70cu.in\100cu.in\150cu.in\250cu.in等容量气枪激发后,达到主脉冲峰值的时间ti。
实施例5
基于本发明实施例1提供地通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,所述步骤S104中根据步骤S103统计所得ti计算Δti,即计算70cu.in\100cu.in\150cu.in\250cu.in容量的气枪激发后达到主脉冲峰值的时间与45cu.in容量的气枪达到主脉冲峰值的时间差值Δti。
实施例6
基于本发明实施例1提供地通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,所述步骤S105中气枪阵列具体为:气枪阵列中共有39条枪,其中工作枪33条。空枪6条,总容量4040cu.in,单枪容量与条数分别为45cu.in(6条)、70cu.in(4条)、100cu.in(10条,包含2条空枪)、150cu.in(11条,包含2条空枪)、250cu.in(8条,包含2条空枪)。根据Δti,计算不同容量气枪沉放深度,该深度使得不同容量气枪的子波信号接收延迟时间恰好为Δti。各容量气枪沉放深度具体为,45cu.in/70cu.in容量气枪沉放5m,100/150cu.in容量气枪沉放6m,250cu.in容量气枪沉放7m。缆的沉放深度8m。
实施例7
基于本发明实施例1提供地通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,所述步骤S106中对步骤S105所述气枪阵列的子波进行频谱分析。具体包括:计算气枪阵列子波的主脉冲峰值、虚反射值、气泡脉冲峰值,通过傅里叶变换对子波做频谱分析,以最大振幅-6dB为判断有效带宽的标准,所以本发明以-6dB来求得子波主脉冲的有效带宽,并得到子波主脉冲主频。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。
上述装置/单元之间的信息交互、执行过程等内容,由于与本发明方法实施例基于同一构思,其具体功能及带来的技术效果,具体可参见方法实施例部分,此处不再赘述。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。上述系统中单元、模块的具体工作过程,可以参考前述方法实施例中的对应过程。
二、应用实施例:
应用实施例1
本发明实施例提供的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法包括以下步骤:
步骤1、根据范德瓦尔斯非理想气体条件气枪子波模型模拟气枪子波,首先设定模型初始条件。具体包括气枪沉放深度7m,缆的沉放深度为8m,气枪容量和条数分别为45cu.in、70cu.in、100cu.in、150cu.i、250cu.in。采样间隔为0.0005s,海水密度为1.03g/cm3,海水速度为1500m/s,海水温度为293.15开尔文,海面反射系数为-0.9等。
步骤2、根据步骤1所设的初始条件来进行模拟过程。具体执行过程为:
a)输入气枪子波模型所有的初始条件。
b)开始时间循环并计算t=k时刻的气泡体积
c)利用方程(1)计算t=k时刻的气泡压力,Pb;
d)通过方程(4)计算气泡热损失率,
e)通过方程(9)计算气体的释放速率,
f)计算t=k时刻的气泡体积变化率,
g)通过方程(7)计算t=k时刻气泡内的温度变化率,
h)计算气泡壁的焓差,
i)通过对方程(1)求关于时间t的微分得到气泡压力的变化率,
j)对焓差求关于时间t的微分得到,
k)通过方程(10)计算t=k时刻的气泡壁的速度变化率,既气泡壁的加速度;
l)对求关于时间t的微分得到,/>
m)由于气枪子波模拟是一个迭代过程,通过二阶泰勒级数展开,可以得到气泡壁半径、气泡壁速度、气体温度以及气泡内气体的物质量: and/>
n)将气泡压力表示为焓、气泡壁速度和以及气泡半径的函数: R0为气泡中心到远场点的距离;
o)重复步骤(a)至(n),直到t>tmax;
p)计算气枪远场子波声压,包括海面虚反射:
q)Rs表示海面反射系数,D1是气枪与水听器之间的距离是,D2是气枪的海面镜像与水听器之间的距离,/>是气枪信号经过D1和D2的时间延迟。经过步骤(a)到(q)得到不同容量的气枪子波,并统计得到如图2所示的不同容量气枪子波。
步骤3、统计不同容量气枪子波从激发到达到主脉冲峰值所用时间ti,如图3所示。
步骤4、根据步骤3统计所得ti计算Δti;
步骤5、根据设定的初始条件,进行尖脉冲子波气枪阵列模拟,如图4本发明实施例提供的所用气枪阵列的平面图所示;
具体过程为:气枪阵列中共有39条枪,其中工作枪33条。空枪6条,总容量4040cu.in,单枪容量与条数分别为45cu.in(6条)、70cu.in(4条)、100cu.in(10条,包含2条空枪)、150cu.in(11条,包含2条空枪)、250cu.in(8条,包含2条空枪)。根据Δti,计算不同容量气枪沉放深度,该深度使得不同容量气枪的子波信号接收延迟时间恰好为Δti。各容量气枪沉放深度具体为,45cu.in/70cu.in容量气枪沉放5m,100/150cu.in容量气枪沉放6m,250cu.in容量气枪沉放7m。缆的沉放深度8m。
步骤6、根据傅里叶变换对尖脉冲子波进行频谱分析,得到图5(a)是本发明实施例提供的模拟出来的尖脉冲气枪子波图;图5(b)是本发明实施例提供的模拟出来的频谱图;
为了更好的展示本发明的技术优越性,与渤海某区块实际气枪阵列子波进行对比,得到图6(a)本发明实施例提供的常规阵列模拟出来的子波图;图6(b)本发明实施例提供的利用本发明模拟出来的尖脉冲气枪子波图;
通过图6(a)-图6(b)中子波对比可以看出,采用本发明方法之后,所得到的气枪子波的频带宽度得到拓宽,频谱中低频能量得到增强。
应用实施例2
本发明实施例还提供了一种利用获得的尖脉冲子波进行高精度的全波形波场成像。该全波形波场成像的具体方法包括marmousi模型的速度场数据、本发明构造的尖脉冲子波以及进行波场成像的计算方法。
所述marmousi模型的速度场数据如图7所示,该模型基本包含了包括裂缝、断层、凹陷、向斜、背斜、潜山等各种地质构造,通过对该模型的成像可以较全面的反应地震子波的质量优劣。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。
三、实施例相关效果的证据:
图8所示为通过本发明构造的尖脉冲子波与常规子波进行波场成像的效果对比图。其中常规阵列为渤海某一区块所用实际气枪阵列,阵列编号为718_4040-7-66-air,本发明所用阵列是在该阵列的基础上进行优化所得,阵列编号为718_4040_7_66-air-depth_narrow-1。
通过图8(a)-图8(b)可以看出,利用本发明构造的尖脉冲子波完成的波场成像效果更好,具体为能量更强,下传深度更深,宽度更广,中深层潜山和内幕地区构造刻画更清晰。
以上所述,仅为本发明较优的具体的实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,其特征在于,该方法根据气枪不同沉放深度与子波信号接收时间的关系,通过调节不同容量气枪的沉放深度,使阵列中部分气枪的主脉冲峰值时刻进行延迟,构造主脉冲波形较窄的尖脉冲震源子波,增大气枪震源低频能量;具体包括以下步骤:
S1、根据范德瓦尔斯非理想气体气枪子波模型模拟不同容量的气枪子波,设定模型初始条件;
S2、根据步骤S1设定的初始条件执行模拟过程;
S3、对模拟出的气枪子波进行分析,统计不同容量气枪子波从激发到达到主脉冲峰值的时间ti;
S4、根据模拟结果,容量越小ti越小,假设最小ti为t0,计算不同容量气枪从激发到达到主脉冲峰值的时间与t0的差值Δti;
S5、依据步骤S4得到的Δti,计算不同容量气枪的沉放深度,沉放深度使得不同容量气枪子波信号的接收延迟时间为Δti;随后进行立体阵列尖脉冲子波模拟;
S6、对步骤S5得到的气枪阵列子波做频谱分析。
2.根据权利要求1所述的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,其特征在于,在步骤S1中,范德瓦尔斯非理想气体气枪子波模型的方程表达为:
式中,a、b均为范德瓦尔斯常数,a=0.1404m6·Pa·mol-2,b=3.764×10-5m3·mol-1,Tg为有效热力学温度,Rg为普适气体常数,mg为气体质量,Vg为体积;
有效热力学温度Tg取决于气室内的高压气体:
Tg=Tw(1+Pg/Pc) (2)
式中,Pc=139MPa;
气枪激发过程中,根据能量守恒定律,气泡热传播损失和气泡物质量的传递获得的能量必须与气泡内能的变化平衡,有:
式中,T是气泡温度,P是气泡压力,mb是气泡内气体的物质量,U=CmmbT表示气泡的内能,Cm是定压比热容,dQ/dt是通过气泡壁的传热速率,传热系数k是通过对模型与试验数据拟合确定的,气泡热损失率表示为:
式中,ΔT=Tb-Tw是气泡温度Tb和周围水温Tw之间的温度差,R是气泡半径;
利用范德瓦尔斯非理想气体方程,非理想气体的内能为气体温度和体积的函数:
全微分方程表示为:
热力学第一定律转化为:
式中,Rg=Cm-Cm;
为获取气体物质的量的变化率引入节流常数τ,影响气体通过气枪端口的速率;式中Vg是气室容量,mg是气室内气体的物质量,Pg是气枪压力,Pb是气泡压力;
对于实际应用的气枪,高压气体释放到水中的速度和总量受到端口大小、端口打开时间参数的控制,影响气枪子波性能;根据幂次定律,表示为:
式中,τ0是与容量无关的端口节流常数,ζ是通过与实测数据的对比确定的节流幂律指数;根据测量和计算结果,在任意给定的时间内,通过气枪端口的气体流量取决于气枪内外的压力差,气体释放的速率表示为:
式中,mb是释放到气泡中的气体物质的量,mg|t=0是气室内气体总量,η是气泡中气体量与总量的比;
气泡壁的运动方程表示为:
式中,R是气泡半径,u和分别是气泡壁的速度和加速度,c是声波在流体介质中的速度,/>是气泡壁的焓差,ρ∞是无穷远处的静水密度,Pb是气泡压力,P∞是无穷远处静水压;气泡的静水压力在气泡由于浮力上升的过程中发生改变,气泡上升过程的气泡垂直上升速度的表达式为:
式中,z是气泡深度,g是重力加速度常数,R是气泡半径,静水压力P∞的表达式为:
式中,Patm是标准大气压,zg是气枪深度;在距离气枪1m处,气枪子波信号表示为:
在低频时,气泡间的相互作用不忽视;气泡间的这种相互作用是对流体静水压的调节;气泡间的相互作用使气泡周围的压力变化;对于地震波长,气泡是一个点,任何任意气泡周围的压力场均是静水压力加上气泡产生的时变压力场的叠加;第i各气泡处的有效静水压力为;
式中,P∞是静水压,∑k≠iΔPik是气枪阵列中所有其他气枪的压力贡献之和,ΔPik是第k个气泡引起的对第i个气泡的静水压力扰动,以及第k个气泡引起的对第i个气泡的延时和按距离标度的压力特征:
式中,rik表示第i个气泡与第k个气泡之间的气泡间距。
3.根据权利要求2所述的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,其特征在于,在步骤S1中,初始条件具体为:
步骤1.1、条件1:将气枪压力的初始值Pg|t=0设置为工作压力;
步骤1.2、条件2:气泡内初始温度设置为Tg=Tw(1+Pg/Pc);
步骤1.3、条件3:气泡初始体积Vb|t=0=Vg,初始半径为
步骤1.4、条件4:气泡壁的初始速度为u=0;
步骤1.5、条件5:气泡初始压力Pb|t=0=P∞,初始温度为水温Tw=18°,气泡内的初始物质量为
步骤1.6、条件6:设置各条气枪的放置位置(x,y,z)。
4.根据权利要求2所述的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,其特征在于,在步骤S2中,根据步骤S1设定的初始条件执行模拟过程,具体为:
步骤2.1、输入范德瓦尔斯非理想气体气枪子波模型的初始条件;
步骤2.2、开始时间循环并计算t=k时刻的气泡体积
步骤2.3、利用方程(1)计算t=k时刻的气泡压力,Pb;
步骤2.4、通过方程(4)计算气泡热损失率,
步骤2.5、通过方程(9)计算气体的释放速率,
步骤2.6、计算t=k时刻的气泡体积变化率,
步骤2.7、通过方程(7)计算t=k时刻气泡内的温度变化率,
步骤2.8、计算气泡壁的焓差,
步骤2.9、通过对方程(1)求关于时间t的微分得到气泡压力的变化率,
步骤2.10、对焓差求关于时间t的微分得到,
步骤2.11、通过方程(10)计算t=k时刻的气泡壁的速度变化率,既气泡壁的加速度;
步骤2.12、对求关于时间t的微分得到,/>
步骤2.13、由于气枪子波模拟是一个迭代过程,通过二阶泰勒级数展开,可以得到气泡壁半径、气泡壁速度、气体温度以及气泡内气体的物质量:and
步骤2.14、将气泡压力表示为焓、气泡壁速度和以及气泡半径的函数:R0为气泡中心到远场点的距离;
步骤2.15、重复步骤(2.1)至(2.14),直到t>tmax;
步骤2.16、计算气枪远场子波声压,包括海面虚反射:
Rs表示海面反射系数,D1是气枪与水听器之间的距离是,D2是气枪的海面镜像与水听器之间的距离,/>是气枪信号经过D1和D2的时间延迟。
5.根据权利要求2所述的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,其特征在于,在步骤S3中,统计不同容量气枪子波从激发到达到主脉冲峰值的时间ti包括:统计45cu.in\70cu.in\100cu.in\150cu.in\250cu.in等容量气枪激发后,达到主脉冲峰值的时间ti。
6.根据权利要求2所述的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,其特征在于,在步骤S4中,计算不同容量气枪从激发到达到主脉冲峰值的时间与t0的差值Δti为计算70cu.in\100cu.in\150cu.in\250cu.in容量的气枪激发后达到主脉冲峰值的时间与45cu.in容量的气枪达到主脉冲峰值的时间差值Δti。
7.根据权利要求2所述的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,其特征在于,在步骤S6中,对步骤S5得到的气枪阵列子波做频谱分析具体包括:计算气枪阵列子波的主脉冲峰值、虚反射值、气泡脉冲峰值,通过傅里叶变换对子波做频谱分析,以最大振幅-6dB为判断有效带宽的标准,求得子波主脉冲的有效带宽,并得到子波主脉冲主频。
8.一种气枪阵列,实施权利要求1~7任意一项所述的通过立体阵列构造尖脉冲子波提高气枪震源低频能量的方法,其特征在于,所述气枪阵列包括工作枪33条,空枪6条,总容量4040cu.in;单枪容量与条数分别包括:
6条45cu.in;
4条70cu.in;
10条100cu.in,其中包含2条空枪;
11条150cu.in,其中包含2条空枪;
8条250cu.in,其中包含2条空枪。
9.一种如权利要求8所述气枪阵列在海洋油气资源勘探气枪震源设备上的应用。
10.一种如权利要求8所述气枪阵列在地质调查中气枪震源设备上的应用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211664560.8A CN116992616A (zh) | 2022-12-23 | 2022-12-23 | 一种构造尖脉冲子波提高气枪震源低频能量的方法及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211664560.8A CN116992616A (zh) | 2022-12-23 | 2022-12-23 | 一种构造尖脉冲子波提高气枪震源低频能量的方法及应用 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116992616A true CN116992616A (zh) | 2023-11-03 |
Family
ID=88520169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211664560.8A Pending CN116992616A (zh) | 2022-12-23 | 2022-12-23 | 一种构造尖脉冲子波提高气枪震源低频能量的方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116992616A (zh) |
-
2022
- 2022-12-23 CN CN202211664560.8A patent/CN116992616A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101241192B (zh) | 一种消除气枪近场子波虚反射的方法 | |
CN108549100B (zh) | 基于非线性高次拓频的时间域多尺度全波形反演方法 | |
CN108108331B (zh) | 一种基于拟空间域弹性波方程的有限差分计算方法 | |
CN110687617B (zh) | 一种地震勘探气枪阵列远场子波模拟方法、脸谱评价方法及装置 | |
Zhang et al. | Numerical study on attenuation of bubble pulse through tuning the air-gun array with the particle swarm optimization method | |
CN110007340A (zh) | 基于角度域直接包络反演的盐丘速度密度估计方法 | |
CN105911584A (zh) | 一种隐式交错网格有限差分弹性波数值模拟方法及装置 | |
CN111665556B (zh) | 地层声波传播速度模型构建方法 | |
CN115980830A (zh) | 通过尖脉冲子波提高气枪震源低频能量的方法及应用 | |
CN111257930B (zh) | 一种黏弹各向异性双相介质区域变网格求解算子 | |
CN116992616A (zh) | 一种构造尖脉冲子波提高气枪震源低频能量的方法及应用 | |
CN115980829A (zh) | 延迟激发构造宽脉冲子波拓宽气枪震源频带的方法及应用 | |
CN115877440A (zh) | 立体排列和延迟激发拓宽气枪震源子波频带的方法及应用 | |
CN111665546B (zh) | 用于可燃冰探测的声学参数获取方法 | |
Wang et al. | Simulating the signature produced by a single airgun under real gas conditions | |
Chunchuzov et al. | Atmospheric boundary layer as a laboratory for modeling infrasound propagation and scattering in the atmosphere | |
CN110231102A (zh) | 一种基于声波速率测试松散煤体温度的方法 | |
CN111665550A (zh) | 地下介质密度信息反演方法 | |
Wang et al. | Pressure waves from air gun bubbles: A numerical analysis based on the finite volume method | |
CN111665551B (zh) | 用于桥梁基底探测的声学参数获取方法 | |
Zhao et al. | Analysis of factors affecting the performance of air gun source under shallow water conditions: Insights from the non-ideal gas air gun wavelet model | |
Huaibang et al. | Seismic acceleration signal theory and application onshore | |
Fei et al. | Design of broadband sweep signal using the statistical near field wavelet | |
Ortiz-Osornio et al. | Physical modeling of the reflectivity and transmissivity dependence on tilt and azimuth of a material with orthorhombic symmetry | |
Silva et al. | Azimuthal variations of P-wave seismic attributes over a HTI medium: A physical modeling study |
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 |