2. 西南石油大学, 成都 610500
2. Southwest Petroleum University, Chengdu 610500, China
震源激发的能量只有一小部分投射向目的层并最终返回地表,带回了地下反射信息,其他大部分由于地层临界角度、接收排列长度限制只能在近地表附近传播,这部分能量传播过程中与近地表的非均匀介质作用形成了大量的噪声,降低了地震资料信噪比(徐峰等,2011).提升资料信噪比有两种思路:一种是绝对地提升采集数据中的反射能量,这就要增加震源初始能量,根据Sharpe(1942)的理论,爆炸压力产生的弹性波振幅正比于爆炸空腔的面积和爆炸压力,扩大爆腔和爆压就需要更大的炸药药量,而大量的生产实践表明,激发能量不会随着炸药用量线性增加,当达到一定用量后,能量基本不再变化——存在饱和激发现象(凌云,2001;杨贵祥,2005);另一种思路是相对提升能量,将激发后的能量约束聚焦在向目的层投射的角度方向上,同步压制噪声能量源,进而提升信噪比.本文探讨后一种方法.
波具有方向性,对波的方向性利用较早的是军事上的相控阵雷达(Dolph,1946),出于对空搜索目标的需要,成百上千个发射/接收单元被组合在一个天线上,通过聚焦波束增强探测距离和精度.地震波的方向性验证是20世纪70年代美国石油公司在Oklahoma的Tulsa地区做过的地震波束形成实验,该实验证明了地震波的定向方向可以通过延时时间间隔调节(Arnold,1977).2002年吉林大学林君教授课题组自主研发的相控震源系统也是基于地震波可以定向的理论(Chen and Lin.,2003).随后,姜弢等(姜弢等, 2004; 姜弢,2006)对定向地震波信号的信噪比和分辨率进行了定量分析,证明了基于相控震源的定向地震波方法在理论和实践上的可行性和优越性.
目前,无论是产生电磁波的相控阵雷达还是产生地震波的相控震源,都是研究定向波的传播,对于事先未定向的单点激发的波场数据,是否存在后期进行室内定向组合的可能,鲜有讨论.排除掉耦合性差异,波场符合线性叠加原理——即,A、B两个炮点单独激发后再叠加的波场应该与它们同时激发的波场保持一致,两个检波点单独接收的波场叠加结果应该与它们组合接收的波场保持一致.这为室内定向组合处理提供了理论依据.
进行室内定向组合首先需要知道待组合的方向,这个方向可以从波场自身的方向分解中求取.Yoon(2004, 2011)在进行逆时偏移成像研究的时候,引入坡印廷矢量来计算传播中波场的方向和倾角信息,为解决计算稳定性问题随后又进行了子波多周期叠加、时窗平滑等技术改进,他的一系列研究为室内定向组合技术提供了求解思路.秦龙等(2019)用Yoon的方法,计算了波场传播到地表的方向,推进了室内组合的实用化.
进行室内定向组合还需要确定具体的组合参数.有关组合参数设计的研究较多,Johnson(1939)从天线理论出发,推导了组合响应函数的表达式,成为组合参数设计的依据,Holzman(1963)、Rietsch E(1979)提出的切比雪夫多项式逼近方法通过调整加权系数来逼近需要的组合响应;Vermeer(2002)研究了组合形态的影响,得出四点菱形震源组合和28点圆形检波器组合在360°的范围内可均匀的获得40 dB噪声衰减的结论;Johnson(1968)研究通过时差控制来改变组合响应的方法,由检波器位置确定组合延时量;King等(1973)采用参考道与各独立地震道做互相关的方法确定延时量;秦龙等(2019)在其炮点组合的研究中也使用了时间延迟的方法来实现方向聚焦;徐峰等(2020)在检波器组合的研究中给出了组合个数、间距、时差、加权多元参数联合运用的设计方法.
本文总结前人研究成果,利用了地震波的方向性,将方向性求取技术和多元组合参数设计技术结合起来,对已采集的地震数据进行室内组合处理,达到相对提升信噪比,改善目的层成像的目的.
1 方法原理 1.1 能量定向传播机理在各向同性介质中传播的地震波不会发生方向选择——不同方向间的能量强弱关系一致.但当地震波传播到地层界面,由于透射/反射,能量会在界面两侧重新分配,这种分配与入射角度、界面两侧的弹性参数有关.广为接受的Zeoppritz方程描述了界面处能量的分配关系,式(1)给出了它较为常用的Shuey(1985)近似式:
(1) |
其中RP为纵波入射下的界面反射系数,θ为纵波入射角,公式分三部分,R0代表法向入射时纵波反射系数,中间的Gsin2θ代表中等入射角项,最后一部分F代表大角度入射项.图 1给出了VP1=2500 m·s-1, VP2=2700 m·s-1, VS1=1800 m·s-1, VS2=1900 m·s-1, ρ1=2.0 g·cm-3, ρ2=2.5 g·cm-3参数下的界面反射系数随入射角度变化曲线,可见随入射角增大,反射系数先减小后增大.
界面的存在使传播中的地震波能量在不同方向上发生了改变,这种改变会随地层的弹性参数不同而不同,当地层自身倾角较大时,方向性变化会更加明显.地震波在地层界面上发生能量的重分配,是波场具有方向性的根本原因.
下面以Marmousi模型为例来展示地震波传播的方向性.图 2为Marmousi模型的叠前时间偏移剖面,在地下R点(9250 m, 2370 m)设置震源,激发地震波,得到图 3所示的波场照明图,可以看到因复杂的地层接触关系,上传地震波场在不同方向上能量差异较大,较强的能量分布在一些弯曲的“亮线”中以特定角度出射地表.保持模型结构不变,激发点R位置不变,则地下能量传播的路径结构不会改变,出射到地表的地震波能量大小及角度也会被唯一确定.强能量意味着能量投递的高效率,在地表强能量位置,强迫震源逆向出射角度下传更多能量,强迫检波器更多地接收出射角度上传的能量,即可实现增强弱成像区照明的目的.
弹性能量在介质中的传播可以看做是弹性位能和动能的相互转化过程,系统中的弹性能量可以表示为(孙成禹,2007):
(2) |
式中,E为弹性能量,K为表征应力应变关系的介质弹性系数,ρ为介质密度,当我们研究声学介质时,p为标量声压场,V 为速质点振动速度矢量场.
定义能流密度矢量 I 表示单位时间内通过与能量传播方向垂直的单位面积内的弹性能.根据能量守恒原理,
(3) |
根据高斯定理,将面积分转化为体积分,同时考虑体积Ω的任意性,得到
(4) |
对(2)式取时间微分,并结合运动平衡微分方程,得
(5) |
因此,
(6) |
即,在声学介质中,能流密度矢量等于标量声压场与矢量质点振动速度场的乘积,方向由质点振动速度方向决定.
在直角坐标系中,能流密度矢量的倾角θ和方位角φ表达为
(7) |
其中V =vxi+vyj+vzk, θ是与Z轴正方向的夹角,φ是矢量 I 在XOY平面的投影与X轴正方向的夹角,由θ和φ就能确定空间任一点传播地震波的能流密度方向.当地震波传播到地表,θ和φ就是地震波在地表的出射角度.
图 4(a, b)为图 2所示R点激发,用一阶声波方程正演,在地表接收到的vx、vz分量模拟记录.Marmousi模型为二维模型,模拟结果中不存在vy分量,此时φ=0,θ不具有唯一性,不能用来准确刻画出射角度.重新定义一个角度α为能流密度矢量 I 与正X方向的夹角.为了计算的稳定性,统计出射角度的时候可以截取初至波的首周期计算,计算公式如(8)式所示,式中T1~T2为截取的时窗范围.
(8) |
图 5给出了依据公式(8)计算得到的出射角度α分布,可以看到在R点正上方附近地震波以90°出射地表,然后向两侧迅速改变;R点左侧地震波出射方向与X轴正方向相反,角度在100°~120°之间;R点右侧地震波出射角度由90°快速过渡到60°附近.对公式(6)取模,同样截取T1~T2时窗计算出射能量大小,得到图 6.对照图 3与图 6,后者中的高幅值与前者中的“亮线”在地表位置能够一一对应.值得注意,模型的复杂性使地表相邻位置的出射角度和能量变化较大,曲线并不光滑,尤其是在R点上方区域,高陡的断层截断了能量的通路,使透射上来的能量断续展布.
图 5的角度曲线给出了在地表做定向组合的方向,而图 6的能量曲线则可以解答在何处布设定向组合才最有效.
1.3 定向组合参数计算在得到地下目标位置上传地表的波场角度和能量之后,下一步的工作就是要通过定向组合手段将震源激发的地震波聚焦为一个窄波束,将检波器的接收口径也聚焦为一个窄波束,使这两个波束的方向指向波场角度方向.下面以炸药震源组合为例介绍定向组合参数的设计方法.
组合过程为各震源子波的线性叠加过程.如图 7,空间不同位置的几个震源Si发出的子波wi以v的速度传播了di的距离后在空间观测点P进行线性叠加.假定各震源的激发耦合性一致,激发子波的频率成份一致,子波振幅由炸药量决定(实际的情况更为复杂,为方便讨论暂且这样假设),叠加子波表达为:
(9) |
其中ai为能量加权因子(炸药量决定的系数),N为组合个数,wi(ti)为各子波的波形函数,ti是各子波参与叠加的起始时间,它由震源的激发时刻ti激发和激发后传播至观测点P的旅行时间ti旅行构成,而旅行时间又由激发点与观测点的相对距离和地震波传播速度决定.
(10) |
(11) |
(12) |
公式(9)给出了空间P点的叠加子波形式,随P点位置变化,叠加子波能量将发生改变,以震源的组合中心为圆心,P点位置为半径,绘制一个圆,利用公式(12)计算得到不同方向上的能量变化,这通常被称为能量方向因子(图 8).若组合参数的搭配使得某方向上ti近乎相等,则在该方向上子波等相位叠加,能量达到最大,实现了聚焦;在该方向之外,相位相抵,能量不同程度被削弱,实现了能量压制,这是组合定向的基本原理.
图 8a为某套组合参数下方向因子的三维空间展现,颜色值代表了归一化的能量大小,图中能量在斜向下的方向上实现了聚焦.图 8b为图 8a三维方向因子沿XOZ平面的二维切片(图中的参数1),径向为归一化能量大小,环周为空间角度,图中同时给出了能够将能量聚焦到同一角度方向的组合参数2对应的方向因子以及能量理想聚焦情况下的方向因子.
理想的定向组合是将能量只聚焦在需要的方向上,其他方向不泄露,如图 8b紫色阴影所示.而实际组合时因组合数目、组合间距、子波频率成分等因素的制约,往往在需要的压制区能量会有泄露,且越追求能量的聚焦性,压制区就会泄露越多的能量.如图 8b中,相较于红色曲线,蓝色曲线在290°附近更加聚焦,但在其他需要压制的角度上能量幅值更大,定向性反而更差.因此做定向组合时不能一味的追求聚焦性.
考虑到用公式(8)计算的角度也会存在误差,实际设计定向组合参数时推荐使用一个角度区间来代替这个单一角度,这个区间由出射角度向两侧膨胀一个角度β来实现(比如待聚焦角度α=300°,采用β=10°,则聚焦区间设定为290°~ 310°).定义一个目标函数,使聚焦区间内的能量相较于区间外最大,作为设计组合参数的依据.
(13) |
公式(13)的分子项代表聚焦区间内的能量,分母项代表聚焦区间之外其它部分的能量.式(13)是进行震源平面面积组合的参数求解目标函数,若只进行沿测线的线性组合,可以忽略φ的积分.检波器定向组合与震源组合原理一致,区别在于检波器接收的都是来自于地下半空间的信号,应用公式(13)的时候,分母中被减项θ的积分区间变为180°~360°即可.
由于叠加子波是以组合个数N、组内单元的空间位置(xi, yi, zi)、组合加权系数ai、组合激发时间ti激发为自变量的,因此定向组合问题实质就是以公式(13)为目标函数,以N、(xi, yi, zi)、ai、ti激发为自变量的多元寻优问题,一般采用数值方法求解.
2 数值模拟本文提出的定向组合技术方法,是基于某一特定地质模型的,模型结构决定了地震波的传播方向性,为验证该方法在改善弱成像区照明上的效果,采用Marmousi模型进行地表炮点的激发模拟.炮点位置的选择参照图 6,选择能量较大且曲线较平稳的区域设置两个激发点S1(7350 m, 20 m)、S2(10750 m, 20 m),从图 5中读取到S1点的出射角度113.6°,S2点65.3°,地震波在震源端是入射,根据炮检互易原理,将出射角度旋转180°可得到需要的入射角度.参照图 8b的角度坐标定义方法,S1点入射角度应为293.6°,S2点245.3°,选择膨胀角度β=20°设定两个位置的聚焦区间.Marmousi模型为二维模型,只能使用沿测线的线性组合方式,各自以S1、S2为组合中心,固定震源组合数目N=5,各子震源等深度激发zi=20 m,采用公式(13)优化5个子震源相距组合中心的位置xi、能量加权系数ai、激发时间ti激发,得到表 1.组合震源的空间展布形态及能量方向因子见图 9,能量向斜下方聚焦于指定的角度区间.
首先对S1、S2两个位置不组合进行单炮点激发照明,得到图 10(a, b),然后按照表 1的参数进行炮点定向组合激发照明,图 11a为S1炮点定向组合激发效果,图 11b为S2炮点定向组合激发效果.对比二者,可以发现定向组合后波场能量在地表聚拢后以较窄的范围倾斜下传,在激发点与地下R点的通路上能量强度明显增强,至R点所在区域,照明强度较单炮点激发提升较大.
参照图 6,地表不同位置对R点区域照明贡献程度不一,选择能量贡献较大的位置进行定向组合更有实际意义,本文采用平均值方法确定待组合位置.对图 6曲线取平均值,大于平均值的位置进行定向组合,其余部分单点激发不组合,同时保证组合激发时的总能量与单点激发的总能量一致,定向方向遵循图 5计算结果.对组合后的数据处理得到叠前时间偏移剖面(图 12),对比图 2,高陡断层下方区域的成像效果明显改善,验证了该方法的有效性.
塔里木盆地塔西南工区表层巨厚黄土覆盖,地表高程起伏较大,黄土下方覆盖一层第四系砾石,砾石层下伏老地层由于构造运动倾角普遍较大,砾石与老地层间呈角度不整合关系,高速层与低降速层间的地震波阻抗差较大,工区的地质特征非常不利于采集到高质量的地震资料,工区典型地表地貌如图 13所示.2019年,塔里木油田在该工区进行了一条线束地震采集作业,目的在于落实工区**段南北向断裂地质构造.为提高采集资料品质,现场采用了较密的空间采样,采集炮点距30 m,炮线距30 m,检波点距15 m,接收线距30 m,采用26线4炮960道奇偶观测系统,覆盖次数达到了960次.结合前期的采集试验分析,震源采用3口Inline方向线性组合,组内距5 m,激发药量3×8 kg,检波器采用3 m×7.5 m组内距的20支矩形组合模式. 观测系统及野外震源、检波器组合形态见图 14(a, b).
虽然在野外已经采用了组合激发、组合接收的方式,但受限制于极差的地表条件,采集单炮质量不高.极低信噪比的资料,复杂的构造形态,较高的空间采样密度,为应用定向组合技术提供了有利条件.
本文采用的地质模型来源于测线第一轮处理解释成果,如图 15所示.组合处理目标为加强模型中部隆起断块的成像质量.
按照前文所述的技术思路,首先对地下目标区域进行炮点向上波动方程照明,如图 16所示,统计地表接收的角度和能量,叠合显示在照明图上方.图中可以看到出射能量主体分布在隆起断块的地表位置两侧,右翼地层倾角较左翼平缓,接收能量更强;出射角度主体受控于构造两翼地层倾向和倾角大小,同时被表层黄土的非均质性复杂化,左翼均值110°左右,右翼波动较大,在黄土较厚的山体部位波场以近90°出射地表.
取图 16能量曲线均值(灰色虚线0.3的位置)以上部分对应的地表位置进行定向组合,定向方向采用对应地表位置的出射角度统计值(震源组合要进行180°的反向).以在图中的A点进行组合为例,首先根据野外观测系统确定A点附近存在哪些可以用来组合的震源和检波器,以及这些震源和检波器的空间位置关系.如图 17所示,让A点对应图中的组合中心0点,对震源采用环形7点组合,检波器采用15点矩形组合模式,此时组合个数N、组合各单元的相对空间坐标(xi, yi, zi)已经确定,读取A点在图 16中的角度110°,按此设计组合激发时间ti激发与组合能量加权系数ai,优化结果见表 2, 3,得到震源定向组合能量方向因子图 18a,检波器定向组合能量方向因子图 18b.图中X轴正方向对应生产的测线方向,Z轴负方向指向地下空间,可见下半空间能量指向性较好.
对实际数据的组合处理结果见图 19,显示的CMP点位于目标区域的地表投影位置.图 19a为原始采集数据,可见信噪比很低,记录中隐约能看到反射轴的“影子”,但被强噪声覆盖;图 19b为单独采用震源端的炮点定向组合结果,在500 ms、1000 ms、1500 ms、2000 ms,甚至更深层的3000 ms时间附近,地层反射同相轴连续性明显加强;图 19c为单独应用检波端的道间定向组合结果,反射轴改善效果更加明显,且记录整体噪声压制明显;图 19d为综合运用炮-检联合定向组合方案的效果,与图 19a原始数据相比较,信噪比改善明显.工区应用实例说明本文提出的方法对极低信噪比地区叠前数据的品质改善有效.
图 20展示了应用炮-检联合定向组合前后的水平叠加剖面,为突出对比组合效果,除静校正、速度分析、水平叠加外,未做任何去噪、提分辨率处理.取图 20a中的窗1左半侧与图 20b中的窗1右半侧数据拼接成图 21a,以清晰呈现组合前后的对比效果,图 21b做法类似.观察组合前后的两张剖面,很明显经定向组合处理后目标层位整体成像效果更好,地层同相轴的连续性得到了增强.在1000 ms以上的浅层改善效果尤为明显,这与图 19中的叠前记录效果改善能够相互印证;在中深层,定向组合压制了随机噪声,强化了被复杂构造遮挡的下方地层成像.叠加剖面的实例进一步说明本文提出的定向组合方法对改善低信噪比地区叠加数据的成像质量有效.
图 22(a, b)分别给出了应用炮-检联合定向组合前后的叠前记录(图 19)和叠加剖面(图 20)的频谱分析结果,可见组合之后频带宽度有一定程度的收窄,这是由组合本身的频率滤波特性决定的.图中20 Hz之前几乎没有变化,30 Hz~100 Hz之间存在5 dB的能量衰减,衰减趋势比较稳定,这部分频率的能量损失可以考虑在后续处理中采用Q补偿、反褶积等手段予以弥补.
本文分析了地震波场定向传播的机理,推导了地震波在地表位置的出射角度和能量大小,给出了利用组合技术实现地震波定向聚焦的计算方法,经Marmousi模型理论验证和工区采集数据实例验证,证明基于模型的定向激发/接收技术能有效改善低信噪比地区的地震资料品质.研究过程得出如下结论:
(1) 地震波在地下传播是有方向的,方向性与地质结构相关,由地层倾角、弹性参数等因素决定,地质结构越复杂,方向差异性越大;
(2) 空间任意位置处的地震波场方向和能量可以通过能流密度矢量定义,在声波介质中,该矢量方向同质点震动速度矢量方向一致,矢量大小由质点声压场值和质点速度场模量乘积决定.在地下某位置激发,计算波场传播至地表的能流密度矢量,可以用来确定进行定向组合的角度和组合位置;
(3) 定向组合可通过控制组合个数N,组内单元的空间位置(xi, yi, zi),组合加权系数ai,组合激发时间ti激发等参数实现向空间特定方向聚焦能量的目的.设计组合参数使用的目标函数应该满足使能量足够聚焦而又不产生非聚焦角度上明显能量泄露的要求;
(4) 通过基于模型的激发-接收定向组合,能够明显提升叠前数据信噪比,改善叠加剖面成像质量.
需要指出,本文展示的图片只进行了一个轮次的组合试验,由于地表的复杂性、地下目标模型的准确性等因素,无法保证计算能流密度矢量的完全准确,这将导致组合定向出现一定误差,降低信噪比的改善程度,对资料的频谱也会有一定的损伤.更为合理可靠的方式是多轮迭代逼近,即在经本文定向组合,成像质量改善的成果剖面基础上重新解释,对地质模型迭代更新,用新模型再次完成定向组合处理,视剖面质量改善程度决定退出迭代还是继续再下一轮次的更新.
致谢 本文受中国石油集团公司重大科技专项课题“复杂山地、黄土塬及大沙漠区地震关键技术研究与应用”(2018E-1807)的支持与资助,同时,感谢国家自然基金项目“炸药震源在不同岩性中激发生成弹性能量的数值模拟”(41804137)对本项目研究的资助.
Arnold M E. 1977. Beam forming with vibrator arrays. Geophysics, 42(7): 1321-1338. DOI:10.1190/1.1440796 |
Chen Z B, Lin J. 2000. Design for vibrator field experiment based-on vibrator-earth system. Journal of Geoscientific Research in Northeast Asia, 3(1): 107-113. |
Dolph C L. 1946. A Current distribution for broadside arrays which optimizes the relationship between beam width and side-lobe level. Proceedings of the Ire, 34(6): 335-348. DOI:10.1109/JRPROC.1946.225956 |
Holzman M. 1963. Chebyshev optimized geophone arrays. Geophysics, 28(2): 145-155. DOI:10.1190/1.1439160 |
Jiang T, Lin J, Chen Z B, et al. 2004. Seismic beam-forming in phased array ofvibroseis. Journal of Jilin University (Information Science Edition) (in Chinese), 22(3): 181-184. |
Jiang T. 2006. Study on seismic beam-forming method based on phased-array vibrator system[Ph. D. thesis] (in Chinese). Changchun: Jilin University.
|
Johnson C H. 1939. Steady state polar sensitivity curves. Geophysics, 4(1): 33-52. DOI:10.1190/1.1440480 |
Johnson K R. 1968. Correlation function of teleseismic noise. Bulletin of the Seismological Society of America, 58(2): 521-538. DOI:10.1785/BSSA0580020521 |
King D W, Mereu R F, Muirhead K J. 1973. The measurement of apparent velocity and azimuth using adaptive processing techniques on data from the warramunga seismic array. Geophysical Journal of the Royal Astronomical Society, 35(1-3): 137-167. |
Ling Y. 2001. Analysis of charge size and charge type. Oil Geophysical Prospecting (in Chinese), 36(5): 584-590. |
Qin L, Yin C, Cui Y F, et al. 2019. A calculation method of delay explosion time of a source array on an undulating surface. Chinese Journal of Geophysics (in Chinese), 62(4): 1492-1501. DOI:10.6038/cjg2019L0306 |
Rietsch E. 1979. Geophone sensitivities for Chebyshev optimized arrays. Geophysics, 44(6): 1142-1143. DOI:10.1190/1.1441000 |
Sharpe J A. 1942. The production of elastic waves by explosion pressures; I, Theory and empirical field observations. Geophysics, 7(2): 144-154. DOI:10.1190/1.1445002 |
Shuey R T. 1985. A simplification of the Zoeppritz equations. Geophysics, 50(4): 609-614. DOI:10.1190/1.1441936 |
Sun C Y. 2007. Theory and Methods of Seismic Waves (inChinese). Dongying: China University of Petroleum Press.
|
Vermeer G J O. 2002. 3-D Seismic Survey Design. 2nd ed. Tulsa: Society of Exploration Geophysicists.
|
Xu F, Liu F L, Liang X H. 2011. A method for source array design based on the phased-array theory. Oil Geophysical Prospecting (in Chinese), 46(2): 170-175. |
Xu F, Liu F L, Wang Y, et al. 2020. Geophone array parameter design under multi-fre-quency condition. Oil Geophysical Prospecting (in Chinese), 55(2): 235-241. |
Yang G X. 2005. The method of seismic exploration in the bare carbonate rock area. Progress in Geophysics (in Chinese), 20(4): 1108-1128. |
Yoon K, Marfurt K J, Starr W. 2004. Challenges in reverse-time migration. //74th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts, 1057-1060.
|
Yoon K, Guo M H, Cai J, et al. 2011. 3D RTM angle gathers from source wave propagation direction and dip of reflector. //81th Ann. Internat Mtg., Soc. Expi. Geophys. . Expanded Abstracts, 3136-3140.
|
姜弢, 林君, 陈祖斌, 等. 2004. 相控震源地震波定向技术. 吉林大学学报(信息科学版), 22(3): 181-184. DOI:10.3969/j.issn.1671-5896.2004.03.001 |
姜弢. 2006. 基于相控震源的地震波定向方法研究[博士论文]. 长春: 吉林大学.
|
凌云. 2001. 激发药量与药型分析. 石油地球物理勘探, 36(5): 584-590. DOI:10.3321/j.issn:1000-7210.2001.05.009 |
秦龙, 尹成, 崔永福, 等. 2019. 面向目标的起伏地表组合震源延时参数计算方法. 地球物理学报, 62(4): 1492-1501. DOI:10.6038/cjg2019L0306 |
孙成禹. 2007. 地震波理论与方法. 东营: 中国石油大学出版社.
|
徐峰, 刘福烈, 梁向豪. 2011. 基于相控理论的炮点组合设计技术. 石油地球物理勘探, 46(2): 170-175. |
徐峰, 刘福烈, 王渝, 等. 2020. 复频条件下组合接收参数设计. 石油地球物理勘探, 55(2): 235-241. |
杨贵祥. 2005. 碳酸盐岩裸露区地震勘探采集方法. 地球物理学进展, 20(4): 1108-1128. DOI:10.3969/j.issn.1004-2903.2005.04.035 |