CN104062678B - 一种优化***阵列的方法和装置 - Google Patents
一种优化***阵列的方法和装置 Download PDFInfo
- Publication number
- CN104062678B CN104062678B CN201410300715.9A CN201410300715A CN104062678B CN 104062678 B CN104062678 B CN 104062678B CN 201410300715 A CN201410300715 A CN 201410300715A CN 104062678 B CN104062678 B CN 104062678B
- Authority
- CN
- China
- Prior art keywords
- bubble
- peak
- particle
- particles
- air
- 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
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种优化***阵列的方法和装置,包括:初始化***阵列的所有粒子的参数;根据初始化的参数计算所有粒子的远场压力子波;根据计算得到的所有粒子的远场压力子波计算气泡比和/或峰峰值;根据计算得到的所有粒子的气泡比和/或峰峰值更新所有粒子的状态向量,继续计算更新后的气泡比和/或峰峰值,直到计算得到的气泡比和/或峰峰值与上一次迭代计算得到的气泡比和/或峰峰值的差值小于等于预设阈值;从迭代结果中选择气泡比和/或峰峰值最大的粒子作为最优***阵列。本发明能够自动优化***阵列,工作量小,效率高。
Description
技术领域
本发明涉及海洋石油地球物理勘探领域,尤指一种优化***阵列的方法和装置。
背景技术
***阵列因为环保,稳定等特点在海洋资源勘探中得到了广泛应用。***阵列的远场压力子波是考察***阵列优劣的重要指标,同时***阵列的远场压力子波也是后续地震资料处理的重要参数之一。例如,海上地震勘探大多使用***或***阵列作为震源激发地震波,而***阵列的远场子波的气泡比越大,其激发的地震波的低频信号越强。在***设计过程中,由于***阵列由多枪组合而成,根据勘探目的的需要往往需要将不同容积的***参与到组合中。
现有的优化***阵列的方法是:首先,将不同容积的***在***阵列中布局好后,根据油气地球物理服务公司(PGS,Petroleum Geo-Services)开发的Nucleus软件对布局好的***阵列进行模拟仿真得到该***阵列对应的远场压力子波,根据得到的远场压力子波计算气泡比和/或峰峰值;其次,改变***阵列的布局重复进行模拟仿真,直到得到所有可能的***阵列的布局对应的远场压力子波的气泡比和/或峰峰值;最后,从得到的所有远场压力子波中选择气泡比和/或峰峰值最大的远场压力子波对应的***阵列即为最优的***阵列。
现有的优化***阵列的方法要得到满足要求的最优的***阵列必须遍历由所有***种类构成的所有组合,工作量较大,效率较低。
发明内容
为了解决上述问题,本发明提出了一种优化***阵列的方法和装置,能够提高优化***阵列的效率。
为了达到上述目的,本发明提出了一种优化***阵列的方法,包括:
初始化***阵列的所有粒子的参数;
根据初始化的参数计算所有粒子的远场压力子波;
根据计算得到的所有粒子的远场压力子波计算气泡比和/或峰峰值;
根据计算得到的所有粒子的气泡比和/或峰峰值更新所有粒子的状态向量,继续计算更新后的气泡比和/或峰峰值,直到计算得到的气泡比和/或峰峰值与上一次迭代计算得到的气泡比和/或峰峰值的差值小于等于预设阈值;
从迭代结果中选择气泡比和/或峰峰值最大的粒子作为最优***阵列。
优选地,所述粒子的参数至少包括:***数目、粒子数、每个粒子的状态向量,每个粒子的初始速度;
所述每个粒子的状态向量是由所述***阵列中所有***的三维空间坐标位置、压力和气泡体积等效容积构成的多维空间向量。
优选地,根据所述初始化的参数采用Ziolkowski模型计算所述所有粒子的远场压力子波。
优选地,根据公式计算所述远场压力子波;
其中,t为时间,Pi(t)为第i个气泡在预设距离处产生的压力,Pi∞为第i个气泡在无穷远处产生的压力,ρi为第i个气泡处海水的密度,r为预设距离,f′i(t)为与第i个气泡表面的参数有关的参数,Ri(t)为第i个气泡的半径,Hi(t)为第i个气泡的气体的焓,P(t)为预设距离处的压力(即远场压力子波),M为***数。
优选地,根据公式Pi∞=P0+ρghi计算所述第i个气泡处海水的密度;其中,P0为标准大气压,g为重力加速度,hi为所述第i个气泡在水中的深度。
优选地,根据公式计算所述第i个气泡的半径;其中,ΔPib为除了所述第i个气泡外其他所有气泡对所述第i个气泡形成的辐射动压力,Pib为所述第i个气泡内的压力,Pjb为第j个气泡在预设距离处产生的压力,Pj∞为所述第j个气泡在无穷远处产生的压力,η为相干枪影响因子。
优选地,根据公式计算第i个气泡的气体的焓,其中,t为时间,Pi(t)为第i个气泡在预设距离处产生的压力,Pi∞为第i个气泡在无穷远处产生的压力,ρi为第i个气泡处海水的密度,r为预设距离,f′i(t)为与第i个气泡表面的参数有关的参数,Ri(t)为第i个气泡的半径,Hi(t)为第i个气泡的气体的焓,P(t)为预设距离处的压力(即远场压力子波)。
优选地,根据公式计算所述第i个气泡内的压力;其中,mi为第i个***内气体的物质的量,Cv为理想气体的等体比热容,Ti气泡为形成气泡后第i个气泡内气体的温度,Rg为气体常数值,Ti***为第i个***内气体的初始温度,Vi为第i个气泡的体积,为第i个气泡吸收的热量功率,τi为第i个***的发射时间,ψ为传热系数,Si为第i个气泡的表面积,ΔTi为第i个气泡和海水的温度之差。
优选地,根据公式更新所述所有粒子的状态向量;
其中,x′id为更新后的粒子的状态向量,xid为更新前的粒子的状态向量,v′id为更新后的粒子的速度,vid为更新前的粒子的速度,w为权重因子,c1和c2为学习因子,r1和r2为[0,1]范围内的均匀随机数,pid为迭代过程中各个粒子的最优值对应的粒子的状态向量,pgd为整个粒子群迄今为止搜索到全局极值对应的粒子的状态向量。
优选地,根据公式更新所述所有粒子的状态向量;
其中,x′id为更新后的粒子的状态向量,xid为更新前的粒子的状态向量,v′id为更新后的粒子的速度,vid为更新前的粒子的速度,w为权重因子,c1和c2为学习因子,r1和r2为[0,1]范围内的均匀随机数,pid为上一次迭代过程中气泡比和/或峰峰值最大的粒子的状态向量,pgd为整个粒子群迄今为止气泡比和/或峰峰值最大的粒子的状态向量。
优选地,初始化后的所有粒子的同一参数的所有取值满足正态分布。
优选地,根据公式A=A0+a0randn1初始化所述粒子的参数中的优化参数;
其中,A为所有***的三维空间坐标位置,压力,气泡体积等效容积中的任何一个的初始值,A0为A的取值范围的最大值和最小值的平均值,a0为A变化的梯度值,randn1为符合正态分布,在[-1,1]范围内的随机数。
优选地,根据公式v=b0randn2初始化所述粒子的速度,其中,v为粒子的初始速度,b0为粒子的速度变化的梯度值,randn2为符合正态分布,在[-1,1]范围内的随机数。
本发明还提出了一种优化***阵列的装置,至少包括:
初始化模块,用于初始化***阵列的所有粒子的参数;
计算模块,用于根据初始化的参数计算所有粒子的远场压力子波;根据计算得到的所有粒子的远场压力子波计算气泡比和/或峰峰值;根据计算得到的所有粒子的气泡比和/或峰峰值更新所有粒子的状态向量,继续计算更新后的气泡比和/或峰峰值,直到计算得到的气泡比和/或峰峰值与上一次迭代计算得到的气泡比和/或峰峰值的差值小于等于预设阈值;
选择模块,用于从迭代结果中选择气泡比和/或峰峰值最大的粒子作为最优***阵列。
与现有技术相比,本发明包括:初始化***阵列的所有粒子的参数;根据初始化的参数计算所有粒子的远场压力子波;根据计算得到的所有粒子的远场压力子波计算气泡比和/或峰峰值;根据计算得到的所有粒子的气泡比和/或峰峰值更新所有粒子的状态向量,继续计算更新后的气泡比和/或峰峰值,直到计算得到的气泡比和/或峰峰值与上一次迭代计算得到的气泡比和/或峰峰值的差值小于等于预设阈值;从迭代结果中选择气泡比和/或峰峰值最大的粒子作为最优***阵列。通过本发明的方案,能够根据计算得到的粒子群的气泡比和/或峰峰值自动优化***阵列,工作量小,效率高。
附图说明
下面对本发明实施例中的附图进行说明,实施例中的附图是用于对本发明的进一步理解,与说明书一起用于解释本发明,并不构成对本发明保护范围的限制。
图1为本发明的优化***阵列的方法流程图;
图2为本发明的优化***阵列的装置结构组成示意图。
具体实施方式
为了便于本领域技术人员的理解,下面结合附图对本发明作进一步的描述,并不能用来限制本发明的保护范围。
参见图1,本发明提出了一种优化***阵列的方法,包括:
步骤100、初始化***阵列的所有粒子的参数。
本步骤中,所有粒子的参数至少包括:***数目M、粒子数N、每个粒子的状态向量,每个粒子的初始速度。
其中,每个粒子的状态向量是由***阵列中所有***的三维空间坐标位置(X,Y,Z)、压力P***和气泡体积等效容积V构成的多维空间向量。
其中,***的压力是指***发射的气泡的初始压力。
其中,可以将***的三维空间坐标位置、***的压力、***的气泡体积等效容积中的一个或多个参数设为优化参数。
其中,粒子数N为大于等于2的正整数,基于计算速度的考虑,粒子数N可以取值为适当的值,例如,可以取值为各个优化参数具体取值的个数的乘积,或者优化参数的个数和每个优化参数具体取值的个数的乘积。
其中,坐标位置(X,Y,Z)所在的坐标系的Z轴垂直于海平面,XoY平面平行于海平面。
在初始化X,Y,Z,P***,V时,如果其为非优化参数,则取为固定值,该固定值根据实际需求进行取值,例如,当优化等效容积V,而不优化X,Y,Z,P***时,可以将P***取为常用的2000英镑每平方英寸(psi,Plasma Surface Interaction),X取为2米,Y取为2米,Z取为2米。
如果其为优化参数,则根据公式(1)进行初始化。
A=A0+a0randn1 (1)
其中,A可以是所有***的X,Y,Z,P***,V中的任何一个的初始值,A0为A的取值范围的最大值和最小值的平均值,a0为A变化的梯度值,a0可以是A的取值范围的最大值和最小值的差值的公约数,randn1为符合正态分布,在[-1,1]范围内的随机数。
根据公式(2)初始化粒子的初始速度。
v=b0randn2 (2)
其中,v为粒子的初始速度,b0为粒子的速度变化的梯度值,b0可以是v的取值范围的最大值和最小值的差值的公约数,randn2为符合正态分布,在[-1,1]范围内的随机数。
步骤101、根据初始化的参数计算所有粒子的远场压力子波。
本步骤中,根据初始化的参数采用Ziolkowski模型计算所有粒子的远场压力子波。
本步骤中,根据公式(3)计算所有气泡在距离***阵列设计中心位置垂直于海平面方向预设距离r(如9千米(km))处产生的远场压力子波(即预设距离处的压力随时间的变化)。
其中,t为时间,Pi(t)为第i个气泡在预设距离处产生的压力,Pi∞为第i个气泡在无穷远处产生的压力,ρi为第i个气泡处海水的密度,r为预设距离,f′i(t)为与第i个气泡表面的参数有关的参数,Ri(t)为第i个气泡的半径,Hi(t)为第i个气泡的气体的焓,P(t)为预设距离处的压力(即远场压力子波)。
其中,***阵列设计中心位置的三维空间坐标位置(X中心,Y中心,Z中心)分别为所有***的三维空间坐标位置的均值,即
其中,(Xi,Yi,Zi)为第i把***的三维空间坐标位置。
根据公式(5)计算Pi∞。
Pi∞=P0+ρighi (5)
其中,P0为标准大气压,g为重力加速度,hi为第i个气泡在水中的深度。
根据公式(6)计算Ri(t)。
其中,c为声速。
根据公式(7)计算Hi(t)。
其中,ΔPib为除了第i个气泡外其他所有气泡对第i个气泡形成的辐射动压力,Pib为第i个气泡内的压力(即绝对压力),Pjb为第j个气泡在预设距离处产生的压力,Pj∞为第j个气泡在无穷远处产生的压力,η为相干枪影响因子(可取值为0.6)。
根据公式(8)计算Pib。
其中,mi为第i个***内气体的物质的量,Cv为理想其他的等体比热容(空气的等比热容为2到3倍Rg),Ti气泡为形成气泡后第i个气泡内气体的温度,Rg为气体常数值(8.314焦耳每摩尔每开(J/(mol·K)),Ti***为第i个***内气体的初始温度,Vi为第i个气泡的体积,为第i个气泡吸收的热量功率,τi为第i个***的发射时间(经过对试验数据的研究,τ的范围在10~15毫秒左右),ψ为传热系数(其值一般在6000J/(K·m2·s)左右),Si为第i个气泡的表面积,ΔTi为第i个气泡和海水的温度之差。
步骤102、根据计算得到的所有粒子的远场压力子波计算气泡比和/或峰峰值。
如何根据计算得到的远场压力子波计算气泡比和/或峰峰值属于现有技术,并不用于限定本发明的保护范围。
例如,对远场压力子波上主峰值的最大能量和气泡最大能量进行比值得到气泡比,对主峰值的最大能量和最小能量求差得到峰峰值。
步骤103、根据计算得到的所有粒子的气泡比和/或峰峰值更新所有粒子的状态向量,继续计算更新后的气泡比和/或峰峰值,直到计算得到的气泡比和/或峰峰值与上一次迭代计算得到的气泡比和/或峰峰值的差值小于等于预设阈值。
本步骤中,根据公式(9)或公式(10)更新每个粒子的参数。
其中,x′id为更新后的粒子的状态向量,xid为更新前的粒子的状态向量,v′id为更新后的粒子的速度,vid为更新前的粒子的速度,v′id,vid∈[-vmax,vmax],vmax为常数,由用户设定来限制粒子的速度,w为权重因子,c1和c2为学习因子,也称加速常数(accelerationconstant),r1和r2为[0,1]范围内的均匀随机数,pid为上一次迭代过程中气泡比和/或峰峰值最大的粒子的状态向量,pgd为整个粒子群迄今为止气泡比和/或峰峰值最大的粒子的状态向量。
公式v′id=wvid+c1r1(pid-xid)+c2r2(pgd-xid)右边由三部分组成,第一部分为“惯性(inertia)”或“动量(momentum)”部分,反映了粒子的运动“习惯(habit)”,代表粒子有维持自己先前速度的趋势;第二部分为“认知(cognition)”部分,反映了粒子对自身历史经验的记忆(memory)或回忆(remembrance),代表粒子有向自身历史最佳位置逼近的趋势;第三部分为“社会(social)”部分,反映了粒子间协同合作与知识共享的群体历史经验,代表粒子有向群体或邻域历史最佳位置逼近的趋势,根据经验,通常c1=c2=2。
根据公式(9)更新粒子的速度时,如果c1的值较大,则会使得粒子过多的在局部范围内徘徊,而如果c2的值较大,则会促使粒子过早的收敛到局部最小值,从而影响迭代的结果。而采用公式(10)更新粒子的速度时能够有效控制粒子的速度使得算法达到全局探测与局部开采两者间的有效平衡,提高了求取最优解的收敛速度。
为了保证算法的顺利求解,c1+c2必须大于4。
步骤104、从迭代结果中选择气泡比和/或峰峰值最大的粒子作为最优***阵列。
参见图2,本发明还提出了一种优化***阵列的装置,至少包括:
初始化模块,用于初始化***阵列的所有粒子的参数;
计算模块,用于根据初始化的参数计算所有粒子的远场压力子波;根据计算得到的所有粒子的远场压力子波计算气泡比和/或峰峰值;根据计算得到的所有粒子的气泡比和/或峰峰值更新所有粒子的状态向量,继续计算更新后的气泡比和/或峰峰值,直到计算得到的气泡比和/或峰峰值与上一次迭代计算得到的气泡比和/或峰峰值的差值小于等于预设阈值;
选择模块,用于从迭代结果中选择气泡比和/或峰峰值最大的粒子作为最优***阵列。
需要说明的是,以上所述的实施例仅是为了便于本领域的技术人员理解而已,并不用于限制本发明的保护范围,在不脱离本发明的发明构思的前提下,本领域技术人员对本发明所做出的任何显而易见的替换和改进等均在本发明的保护范围之内。
Claims (13)
1.一种优化***阵列的方法,其特征在于,包括:
初始化***阵列的所有粒子的参数;
根据初始化的参数计算所有粒子的远场压力子波;
根据计算得到的所有粒子的远场压力子波计算气泡比和/或峰峰值;
根据计算得到的所有粒子的气泡比和/或峰峰值更新所有粒子的状态向量,继续计算更新后的气泡比和/或峰峰值,直到计算得到的气泡比和/或峰峰值与上一次迭代计算得到的气泡比和/或峰峰值的差值小于等于预设阈值;
从迭代结果中选择气泡比和/或峰峰值最大的粒子作为最优***阵列;
根据公式A=A0+a0randn1初始化所述所有粒子的参数中的优化参数;
其中,A为所有***的三维空间坐标位置、压力、气泡体积等效容积中的任何一个的初始值,A0为A的取值范围的最大值和最小值的平均值,a0为A变化的梯度值,randn1为在[-1,1]范围内符合正态分布的随机数。
2.根据权利要求1所述的方法,其特征在于,所述粒子的参数至少包括:***数目、粒子数、每个粒子的状态向量、每个粒子的初始速度;
所述每个粒子的状态向量是由所述***阵列中所有***的三维空间坐标位置、压力和气泡体积等效容积构成的多维空间向量。
3.根据权利要求1所述的方法,其特征在于,根据所述初始化的参数采用Ziolkowski模型计算所述所有粒子的远场压力子波。
4.根据权利要求3所述的方法,其特征在于,
根据公式
计算所述远场压力子波;
其中,t为时间,Pi(t)为第i个气泡在预设距离处产生的压力,Pi∞为第i个气泡在无穷远处产生的压力,ρi为第i个气泡处海水的密度,r为预设距离,fi”(t)为与第i个气泡表面的参数有关的参数,Ri(t)为第i个气泡的半径,Hi(t)为第i个气泡的气体的焓,P(t)为远场压力子波,M为***数。
5.根据权利要求4所述的方法,其特征在于,根据公式Pi∞=P0+ρghi计算所述第i个气泡处海水的密度;其中,P0为标准大气压,g为重力加速度,hi为所述第i个气泡在水中的深度。
6.根据权利要求4所述的方法,其特征在于,根据公式
计算第i个气泡的半径;
c为声速。
7.根据权利要求4所述的方法,其特征在于,根据公式计算第i个气泡的气体的焓,其中,t为时间,Pi∞为第i个气泡在无穷远处产生的压力,ρi为第i个气泡处海水的密度,Hi(t)为第i个气泡的气体的焓;
其中,ΔPib为除了所述第i个气泡外其他所有气泡对所述第i个气泡形成的辐射动压力,Pib为所述第i个气泡内的压力,Pjb为第j个气泡在预设距离处产生的压力,Pj∞为所述第j个气泡在无穷远处产生的压力,η为相干枪影响因子。
8.根据权利要求4所述的方法,其特征在于,根据公式
计算第i个气泡内的压力;其中,mi为第i个***内气体的物质的量,Cv为理想气体的等体比热容,Ti气泡为形成气泡后第i个气泡内气体的温度,Rg为气体常数值,Ti***为第i个***内气体的初始温度,Vi为第i个气泡的体积,为第i个气泡吸收的热量功率,τi为第i个***的发射时间,ψ为传热系数,Si为第i个气泡的表面积,△Ti为第i个气泡和海水的温度之差。
9.根据权利要求1所述的方法,其特征在于,根据公式
更新所述所有粒子的状态向量;
其中,xi′d为更新后的粒子的状态向量,xid为更新前的粒子的状态向量,v′id为更新后的粒子的速度,vid为更新前的粒子的速度,w为权重因子,c1和c2为学习因子,r1和r2为[0,1]范围内的均匀随机数,pid为上一次迭代过程中气泡比和/或峰峰值最大的粒子的状态向量,pgd为整个粒子群迄今为止气泡比和/或峰峰值最大的粒子的状态向量。
10.根据权利要求1所述的方法,其特征在于,根据公式
更新所述所有粒子的状态向量;
其中,x′id为更新后的粒子的状态向量,xid为更新前的粒子的状态向量,v′id为更新后的粒子的速度,vid为更新前的粒子的速度,c1和c2为学习因子,r1和r2为[0,1]范围内的均匀随机数,pid为上一次迭代过程中气泡比和/或峰峰值最大的粒子的状态向量,pgd为整个粒子群迄今为止气泡比和/或峰峰值最大的粒子的状态向量。
11.根据权利要求1~10任意一项所述的方法,其特征在于,初始化后的所有粒子的同一参数的所有取值满足正态分布。
12.根据权利要求11所述的方法,其特征在于,根据公式v=b0randn2初始化所述所有粒子的速度,其中,v为粒子的初始速度,b0为粒子的速度变化的梯度值,randn2为在[-1,1]范围内符合正态分布的随机数。
13.一种优化***阵列的装置,其特征在于,至少包括:
初始化模块,用于初始化***阵列的所有粒子的参数;
计算模块,用于根据初始化的参数计算所有粒子的远场压力子波;根据计算得到的所有粒子的远场压力子波计算气泡比和/或峰峰值;根据计算得到的所有粒子的气泡比和/或峰峰值更新所有粒子的状态向量,继续计算更新后的气泡比和/或峰峰值,直到计算得到的气泡比和/或峰峰值与上一次迭代计算得到的气泡比和/或峰峰值的差值小于等于预设阈值;
选择模块,用于从迭代结果中选择气泡比和/或峰峰值最大的粒子作为最优***阵列;
根据公式A=A0+a0randn1初始化所述所有粒子的参数中的优化参数;
其中,A为所有***的三维空间坐标位置、压力、气泡体积等效容积中的任何一个的初始值,A0为A的取值范围的最大值和最小值的平均值,a0为A变化的梯度值,randn1为在[-1,1]范围内符合正态分布的随机数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410300715.9A CN104062678B (zh) | 2014-06-27 | 2014-06-27 | 一种优化***阵列的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410300715.9A CN104062678B (zh) | 2014-06-27 | 2014-06-27 | 一种优化***阵列的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104062678A CN104062678A (zh) | 2014-09-24 |
CN104062678B true CN104062678B (zh) | 2017-05-24 |
Family
ID=51550476
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410300715.9A Active CN104062678B (zh) | 2014-06-27 | 2014-06-27 | 一种优化***阵列的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104062678B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104849747B (zh) * | 2015-05-19 | 2017-06-30 | 中国海洋石油总公司 | 一种优化***阵列的方法和装置 |
CN110162907B (zh) * | 2019-05-29 | 2023-04-07 | 包头钢铁(集团)有限责任公司 | 一种利用数值模拟研究得到表征板材成形性各项参数窗口值的方法 |
CN110197029A (zh) * | 2019-05-29 | 2019-09-03 | 包头钢铁(集团)有限责任公司 | 一种仿真模拟用材料参数的分析方法 |
CN113655519B (zh) * | 2021-08-23 | 2023-10-13 | 中海石油(中国)有限公司 | ***节流作用系数和气体释放效率参数获取方法及*** |
CN114089411B (zh) * | 2021-11-05 | 2022-10-04 | 青岛海洋地质研究所 | 一种改进的***阵列声波传输建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6788618B2 (en) * | 2002-03-14 | 2004-09-07 | Input/Output, Inc. | Method and apparatus for marine source diagnostics |
GB2414314A (en) * | 2004-05-20 | 2005-11-23 | Pgs Exploration | Method of seismic source monitoring using modelled source signatures with calibration functions |
CN103163567A (zh) * | 2013-02-27 | 2013-06-19 | 中国石油集团东方地球物理勘探有限责任公司 | 基于能量的***子波气泡比分析方法及装置 |
CN103852782A (zh) * | 2012-11-30 | 2014-06-11 | 中国石油天然气集团公司 | 一种确定最佳***阵列的方法 |
-
2014
- 2014-06-27 CN CN201410300715.9A patent/CN104062678B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6788618B2 (en) * | 2002-03-14 | 2004-09-07 | Input/Output, Inc. | Method and apparatus for marine source diagnostics |
CN1748380B (zh) * | 2002-03-14 | 2010-05-26 | 离子地球物理公司 | 测试声源的方法和装置 |
GB2414314A (en) * | 2004-05-20 | 2005-11-23 | Pgs Exploration | Method of seismic source monitoring using modelled source signatures with calibration functions |
CN103852782A (zh) * | 2012-11-30 | 2014-06-11 | 中国石油天然气集团公司 | 一种确定最佳***阵列的方法 |
CN103163567A (zh) * | 2013-02-27 | 2013-06-19 | 中国石油集团东方地球物理勘探有限责任公司 | 基于能量的***子波气泡比分析方法及装置 |
Non-Patent Citations (2)
Title |
---|
Measurement of air-gun bubble oscillations;Anton Ziolkowski;《GEOPHYSICS》;19981231;第63卷(第6期);第2009-2024页 * |
一种PID参数整定的粒子群优化算法;张兴华 等;《计算机工程与应用》;20071130;第43卷(第33期);第227-229页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104062678A (zh) | 2014-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104062678B (zh) | 一种优化***阵列的方法和装置 | |
CN106646645B (zh) | 一种重力正演加速方法 | |
CN102053258B (zh) | 基于复杂地质构造的自适应三维射线追踪方法 | |
CN106650046B (zh) | 一种舰船空气流场的非定常特性获取方法 | |
Blaise et al. | A dynamic hp-adaptive discontinuous Galerkin method for shallow-water flows on the sphere with application to a global tsunami simulation | |
NO343091B1 (no) | Elektromagnetisk avbildning av undergrunnen ved fire-dimensjonal parallell databehandling | |
Zhang et al. | Numerical study on attenuation of bubble pulse through tuning the air-gun array with the particle swarm optimization method | |
CN111460597B (zh) | 一种基于改进的多目标粒子群优化算法的雷达布站方法 | |
Shen et al. | Three‐dimensional MHD simulation of the evolution of the April 2000 CME event and its induced shocks using a magnetized plasma blob model | |
CN114970213B (zh) | 一种孤子内波起伏效应下的海杂波动态模拟方法 | |
CN110287535A (zh) | 高海况下水下拖曳***运动响应的仿真计算方法 | |
CN106503365B (zh) | 一种用于sph算法的分区搜索方法 | |
Mo et al. | Tracing analytic ray curves for light and sound propagation in non-linear media | |
CN106569270B (zh) | 规则网格速度模型的自适应非结构三角网格化方法 | |
Salapaka et al. | Constraints on locational optimization problems | |
CN115841076A (zh) | 基于bp神经网络模型的浅海分层海底地声参数反演方法 | |
CN105093279A (zh) | 针对山前带的三维地震初至波菲涅尔体层析反演方法 | |
CN105184010A (zh) | 基于快速多极间接边界元法的高频地震波散射模拟方法 | |
Ding et al. | Regional seismic waveform inversion using swarm intelligence algorithms | |
WO2014040338A1 (zh) | 一种优化系数获取方法、装置及相关波场模拟方法、装置 | |
CN110764154B (zh) | 基于改进粒子群算法的时间域航空电磁一维反演方法 | |
Springel | Hydrodynamic simulations on a moving Voronoi mesh | |
CN104849747A (zh) | 一种优化***阵列的方法和装置 | |
CN111256946A (zh) | 波浪模拟装置及方法 | |
CN114647990A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Co-patentee after: China Oilfield Services Limited Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Co-patentee before: China Oilfield Services Limited Patentee before: China National Offshore Oil Corporation |