An experimental comparison of three direct methods of marine source signature estimation

Three methods for estimation of the pressure wavefield generated by a marine airgun array are tested experimentally and compared. In the trial a variety of radiation angles and array configurations were used and some large synchronization errors were deliberately introduced. The source was equipped with near-field hydrophones and a subsource ministreamer. A tethered far-field hydrophone was used so that the three estimated far-field signatures could be compared with an independent measurement. The knowledge of the source signature is important for on-board source array QC, deconvolution, multiple attenuation, stratigraphic trap prediction, modelling and inversion, AVO analysis and reservoir monitoring. The methods perform very well and give estimates whose frequency-domain spectra match the measured spectra to within a few dB and within a few tens of degrees of phase over the tested bandwidth of 3.5–110 Hz. The time-domain error-energy is typically only a few per cent of the signal energy for radiation angles within about 30° of the vertical. The third method proved to be sensitive to an experimental shortcoming leading to overloading of the ministreamer and meaningful comparison was not possible for some test configurations.     

ON THE CALIBRATION OF THE WATER GUN PRESSURE SIGNATURE*

A simple field method was proposed by the author in 1976 for measuring the absolute amplitude of the pressure pulse radiated by marine seismic sources which radiate a bubble pulse. The proposed method involves the recording of the near-field pressure signature radiated by the water gun using a wide-band hydrophone. The key feature of the proposed method is that a knowledge of the hydrophone sensitivity and its distance from the water gun are not required. It is shown that the absolute amplitudes of the pressure pulses radiated by the S80 and P400 water guns obtained using the proposed method are in agreement with those obtained using a Ref-Tek hydrophone.     

MODELLING OF GI GUN SIGNATURES1

In 1989 a new type of marine seismic source was introduced. This new air-gun, which consists of two air chambers instead of one, is called the GI gun. The main feature of this gun is that the bubble created by the gun is stabilized by an injection of extra air from the second chamber at a later time. This injection mechanism reduces the amplitude of the bubble oscillations, which also means that the acoustic signal from a GI gun shot is characterized by a very clean primary pulse followed by very small bubble oscillations. A method for calculating the acoustic signal generated by a GI gun is presented. Based on the solution of a damped Kirkwood–Bethe equation, the far-field pressure of single GI guns and of arrays of GI guns is calculated. It is shown that the optimal values for injection start time and injection period vary with injector volume and gun depth. It is also shown that the precision in the firing time for the injector should be of the order of 4 ms, while the precision of the injection period should be of the order of 8 ms. Modelled and measured far-field signatures have been compared, and the relative error energy is found to be less than 3.5% for all examples.     

DEVELOPMENT OF MORE EFFICIENT AIRGUN ARRAYS: THEORY AND EXPERIMENT*

Source strength of an airgun array may be increased by:

• — utilizing higher pressure,
• — increasing total array volume,
• — employing more guns,
• — improving gun efficiency.

One measure of gun efficiency is “specific source strength”, P_a*, defined as source strength per unit quantity of air used. Typical units are MPa m/l. Most developments are directed toward increasing gun pressure and/or gun volume to increase source strength of the array. These efforts require that more air compressors be installed onboard the ship. Consequently, a larger ship may be needed for the additional compressors, guns, and auxiliary equipment. A development program was initiated in 1976 to increase source strength of the array without using a larger ship. New guns were designed and built—one for 41.4 MPa and 7.37 liter (6000 p.s.i./450 in³) operation and another with 13.8 MPa and 4.92 liter (2000 p.s.i./300 in³) capability. Experiments were conducted with these new guns (and existing guns) over a range of pressures from 13.8 to 41.4 MPa (2000 to 6000 p.s.i.). Design of the new guns was aided by a mathematical model. The model relates physical dimensions of the airgun to acoustic pressure in the water. It consists of four nonlinear differential equations relating

• — shuttle motion,
• — bubble pressure,
• — chamber pressure,
• — bubble radius.

The last equation is the “free-bubble-oscillation equation” and represents the ideal case of a pressurized bubble released instantaneously in water. The three other equations modify this ideal case; the four equations together model an airgun of the type manufactured by Bolt Associates, Inc.     

Quantification of Hydrophone Records of the 2004 Sumatra Tsunami

The 2004 Sumatra-Andaman tsunami was recorded by hydrophones of the International Monitoring System at Site H08 near Diego Garcia, notably in frequency bands extending outside the range of the Shallow Water Approximation. Despite the severe high-pass filtering involved in this instrumentation, we show that the spectral amplitudes recovered around T = 87 s can be successfully explained by generation from the seismic source, in the framework of the normal mode theory of tsunami excitation. At the lower frequencies characteristic of more conventional tsunami waves (800 to 3200 s), the signal is probably present in the hydrophone records, but reliable deconvolution of its spectral amplitude is precluded by the fact that the instrumental filters lowered the tsunami signal to the level of resolution of the instrument digitizer. In the context of distant tsunami warning, hydrophone records could provide useful insight into high-frequency tsunami components, and even at lower, more conventional, frequencies, provided that an unfiltered channel could be recorded routinely.     

CALIBRATION OF MARINE SEISMIC SOURCES USING A HYDROPHONE OF UNKNOWN SENSITIVITY*

The absolute amplitude of the pressure pulse radiated by a marine seismic source is one of the prime criteria used in evaluating its performance. A technique for this measurement is proposed which is applicable to all sources which radiate a bubble pulse. The technique is described with reference to an air-gun array in which the pulse radiated from a single gun is compared with that radiated by the full array. The advantage of the method is that absolute values of pressure are obtained without any need for a calibrated hydrophone. In theory this may seem a trivial advantage but in practice sensitivity factors for the hydrophone channel cannot always be relied upon. The proposed technique is illustrated by an example.     

PERFORMANCE OF 2000 AND 6000 PSI AIR GUNS: THEORY AND EXPERIMENT*

Air guns have been used in various applications for a number of years. They were first used in coal-mining operations and were operated at up to 16000 psi charge pressures. Later, single air guns, operated at 2000 psi, found application as an oceanographic survey tool. Air gun arrays were first used in offshore seismic exploration in the mid-1960's. These early arrays were several hundred cubic inches in total volume and were operated at 2000 psi; they were either tuned arrays or several large guns of the same size with wave-shape kits. Today's arrays have total volumes greater than 5000 cu in. and are typically operated at 2000 psi. Recently, higher-pressure, lower-volume arrays operated at 4000–5000 psi have been introduced; guns used in these arrays are descendants of the coal-mining gun. On first thought one would equate increased gun pressure linearly with the amplitude of the initial pulse. This is approximately true for the signature radiated by a “free-bubble” (no confining vessel) and recorded broadband. The exact relation depends on the depth at which the gun is operated; from solution of the free-bubble oscillation equation, the relation is If P_c,1= 6014.7 psia, P_c,2= 2014.7 psia and P_{O, 1}=P_{O, 2}= 25.8 psia (corresponding to absolute pressure at 25 ft water depth), then Experiments were conducted offshore California in deep water to determine the performance of several models of air guns at pressures ranging from 2000 to 6000 psi and gun volumes ranging from 5 to 300 cu in. At a given gun pressure, the initial acoustic pulse P_a correlated with gun volume V_c according to the classical relation For 1 ms sampled data the ratio varied between 4.5 and 5.5 dB depending on gun model. Pulse width of the 2000 psi signatures indicated they are compatible with 2 ms sample-rate recording while pulse width of the 6000 psi signatures was greater, indicating they are less compatible with 2 ms sample-rate recording. Conclusions reached were that 2000 psi air guns are more efficient than higher pressure guns and are more compatible with 2 ms sample-rate requirements.     

气枪震源激发模式及应用

气枪震源是一种重要的人工地震震源。气枪阵列理论的提出,使得气枪阵列设计技术日趋成熟,并能够在石油勘探和地球物理探测中得到更加广泛的运用。气枪震源在不同领域中应用时,需要不同的组合和激发模式,以适应不同的探测要求。加强主脉冲和加强气泡脉冲,是目前两种主要的气枪激发模式。通过比较研究两种气枪激发模式,讨论各种激发模式在激发时间、气枪间距、频率、分辨率等方面的差异,为气枪震源的广泛运用提供依据。     

LOW-FREQUENCY STANDARD HYDROPHONE * FOR CALIBRATING LINE HYDROPHONE ARRAYS (SEISMIC STREAMERS) AND MARINE SEISMIC SOURCES **

The design of a standard hydrophone with a maximally flat (Butterworth) response in the frequency range 8.0 Hz-1.0 kHz is described. The standard hydrophone has been developed primarily for calibrating line hydrophone arrays (seismic streamers) and marine seismic sources. The standard hydrophone has been used successfully during the past eight years for monitoring the output of a single air gun. It can be used for the calibration of a marine seismic streamer.     

SYSTEM APPROACH TO AIR-GUN ARRAY DESIGN *

To specify intelligently a nondynamite source in a marine seismic data-collection system, it is important to use all known parameters of the system—source, receiver, and recording-system characteristics. A technique has been developed to design the far-field pressure pulse of an air-gun array by taking these parameters into account. Important source variables to consider are interaction among guns in the array and the depth of the array. Near-field pressure signatures of individual guns, which are relatively unaffected by boundaries, have been used to‘construct’the far-field pressure pulse of the array by considering these variables. Comparison between constructed pulses and measured far-field pulses shows substantial agreement. Streamer depth and recording-system bandpass should also be considered when designing an air-gun array. Comparison of far-field pressure pulses for several bandpasses clearly shows the importance of considering this variable; e.g., the initial pulse is severely attenuated when a high-cut filter is used. Likewise, an additional filtering effect due to the streamer's proximity to the surface should be taken into account. Design of an air-gun array using the principles just outlined are illustrated by an example.     

Prediction of signatures of airgun arrays using dual near-field hydrophones

Successful estimation of airgun-array signatures from near-field measurements depends on the accuracy of poorly controlled model parameters such as the effective sea surface reflection coefficient and source depth. We propose a method for prediction of robust source signatures, which are insensitive to fluctuations of the latter parameters. The method uses vertical pairs of near-field hydrophones to measure near-field pressure and its vertical gradient, combination of which eliminates sea surface reflections from the near-field data. This excludes the uncertainty related to the fluctuating sea state and source depth from the process of inversion of the near-field data for source signature. The method explicitly separates the recorded near-field pressure into its up and down going components, which allows one to measure the effective frequency- and angle-dependent sea surface reflection coefficient right at the source, as well as to estimate the actual source depth. Tests on synthetic and field data demonstrate robust performance of the method.     

Marine seismic wavefield measurement to remove sea-surface multiples

We propose a new method for removing sea-surface multiples from marine seismic reflection data in which, in essence, the reflection response of the earth, referred to a plane just above the sea-floor, is computed as the ratio of the plane-wave components of the upgoing wave and the downgoing wave. Using source measurements of the wavefield made during data acquisition, three problems associated with earlier work are solved: (i) the method accommodates source arrays, rather than point sources; (ii) the incident field is removed without simultaneously removing part of the scattered field; and (iii) the minimum-energy criterion to find a wavelet is eliminated. Pressure measurements are made in a horizontal plane in the water. The source can be a conventional array of airguns, but must have both in-line and cross-line symmetry, and its wavefield must be measured and be repeatable from shot to shot. The problem is formulated for multiple shots in a two-dimensional configuration for each receiver, and for multiple receivers in a two-dimensional configuration for each shot. The scattered field is obtained from the measurements by subtracting the incident field, known from measurements at the source. The scattered field response to a single incident plane wave at a single receiver is obtained by transforming the common-receiver gather to the frequency–wavenumber domain, and a single component of this response is obtained by Fourier transforming over all receiver coordinates. Each scattered field component is separated into an upgoing wave and a downgoing wave using the zero-pressure condition at the water-surface. The upgoing wave may then be expressed as a reflection coefficient multiplied by the incident downgoing wave plus a sum of scattered downgoing plane waves, each multiplied by the corresponding reflection coefficient. Keeping the upgoing scattered wave fixed, and using all possible incident plane waves for a given frequency, yields a set of linear simultaneous equations for the reflection coefficients which are solved for each plane wave and for each frequency. To create the shot records that would have been measured if the sea-surface had been absent, each reflection coefficient is multiplied by complex amplitude and phase factors, for source and receiver terms, before the five-dimensional Fourier transformation back to the space–time domain.     

Moment magnitude (Mw) from hydrophone records of low energy volcanic quakes

Earthquake magnitude calibration using hydrophone records has been carried out at Campi Flegrei caldera, an active area close to the highly populated area of Naples city, partly undersea. Definite integrals of the hydrophone records amplitude spectra, between the limits of 1 and 20 Hz, were calculated on a set of small volcano-tectonic earthquakes with moment magnitudes ranging from 1 to 3.3. The coefficients of a linear relationship between the logarithm of these integrals and the magnitude were obtained by linear optimization, thus defining a useful equation to calculate the moment magnitude from the hydrophone record spectra. This method could be easily exported to other volcanic areas, where submerged volcanoes are monitored by networks of hydrophones and seismic sensors on land. The proposed approach allows indeed magnitude measurements of small magnitude earthquakes occurring at sea, thus adding useful information to the seismicity of these volcanoes.

     

大容量气枪震源子波激发特性分析

大容量气枪水库激发作为陆地震源的可行性与有效性已经得到成功验证.为进一步提高气枪震源激发效果,本文通过水库气枪激发试验对单枪容量为2000 in³的气枪震源激发子波特征及规律进行了研究.依据近场水听器和远场短周期地震仪记录数据,分析气枪震源沉放深度、工作压力等不同激发条件对压力脉冲和气泡脉冲的影响.有助于人们根据不同尺度地下结构探测对震源激发信号的要求,调整气枪激发参数和激发环境,获得最佳激发效果.试验结果表明：（1）沉放深度对压力脉冲波形影响较小,其优势频率不随沉放深度而改变;（2）随着沉放深度从5 m增加到11 m,气泡脉冲的优势频率由5 Hz增加至7 Hz,其最大振幅亦近线性递增;（3）工作压力越大,激发压力脉冲能量越强,而对气泡脉冲的影响主要体现在主频降低.适合远距离深穿透地下结构探测的低频信号主要来自大容量气枪所激发气泡的反复振荡,由于气枪振荡过程非常复杂,本文通过较为简洁的数学和物理模型进行了解释.     

Sea surface shape derivation above the seismic streamer

The rough sea surface causes perturbations in the seismic data that can be significant for time‐lapse studies. The perturbations arise because the reflection response of the non‐flat sea perturbs the seismic wavelet. In order to remove these perturbations from the received seismic data, special deconvolution methods can be used, but these methods require, as input, the time varying wave elevation above each hydrophone in the streamer. In addition, the vertical displacement of the streamer itself must also be known at the position of each hydrophone and at all times. This information is not available in conventional seismic acquisition. However, it can be obtained from the hydrophone measurements provided that the hydrophones are recorded individually (not grouped), that the recording bandwidth is extended down to 0.05 Hz and that data are recorded without gaps between the shot records. The sea surface elevation, and also the wave‐induced vertical displacement of the streamer, can be determined from the time‐varying pressure that the sea waves cause in the hydrophone measurements. When this was done experimentally, using a single sensor seismic streamer without a conventional low cut filter, the wave induced pressure variations were easily detected. The inversion of these experimental data gives results for the sea surface elevation that are consistent with the weather and sea state at the time of acquisition. A high tension approximation allows a simplified solution of the equations that does not demand a knowledge of the streamer tension. However, best results at the tail end of the streamer are obtained using the general equation.     

西南印度洋中脊三维地震探测中炮点与海底地震仪的位置校正

在海上实施三维地震探测过程中,人工震源枪阵中心与船上GPS的距离及地震探测作业中的船行方向造成炮点实际位置与预设位置有一定偏差;自由落体投放的OBS由于海流的影响会偏离原定设计位置(投放点),因此,炮点与海底地震仪(OBS)的位置校正是三维地震结构研究中的基本环节.本文利用艏向信息校正了炮点位置;采用蒙特卡洛和最小二乘法方法对海底地震仪的位置进行了校正,并探讨了直达水波曲线特征.结果表明 OBS位置一般偏离设计点1 km左右,其误差范围在20 m以内,校正后的OBS记录剖面展示了真实的记录情况.该研究结果为下一步西南印度洋的三维层析成像研究提供了坚实数据基础,同时为今后南海的三维深部地壳结构探测提供经验与借鉴.     

震源和水听器沉放深度对单道地震勘探的影响

为研究震源和水听器电缆沉放深度对单道地震勘探的影响,建立鬼波影响模型,采用雷克子波模拟震源子波,分析因震源及水听器电缆沉放深度产生的鬼波对一次有效反射地震波的影响,并采用不同震源和水听器电缆沉放深度的数据对本次研究结论进行验证。研究表明：地震剖面分辨率和震源及水听器电缆沉放的总深度呈反比;随着震源或水听器电缆沉放深度的增加,地震剖面的分辨率递减,接收地震波能量从低值递增到最大值,再从最大值缓慢降低到稳定值;当水听器电缆沉放深度小于168.8/f(m)时,接收地震波能量随震源沉放深度变化的波动相对平缓,水听器电缆沉放深度浅有利于地震波的一致性,但沉放深度过浅时接收地震波能量将被极大削弱,会明显降低地震波的有效穿透深度;当震源沉放深度约为292.3/f(m)时,一次有效反射波的峰值和震源鬼波的峰值重叠,接收地震波能量最大;当震源沉放深度大于584.6/f(m)时,不考虑地震波能量在水中衰减的情况下,震源沉放深度继续增加对接收地震波能量的增加无效。     

甘肃祁连山主动源不同枪组激发效果对比

大容量气枪震源激发探测因具绿色环保、传播距离远和高度重复性等特点成为地下介质探测的新方法,利用其进行地下介质物性变化高精度探测,则对不同激发条件下气枪传播信号的一致性提出要求。本文选取了祁连山主动源不同枪组（4支和3支气枪）在激发环境一致情况下分别激发时各远台波形记录,通过信号叠加、互相关及频谱分析等方法,对比分析了各台站记录的不同枪组激发信号,及其与震中距的关系,并讨论了不同枪组激发对地下介质走时变化观测的影响。结果表明：（1）不同枪组的激发信号波形相关度随震中距增大而降低,信号震相随激发能量有微弱的变化,震中距超过50 km后其变化几乎分辨不出;（2）远台接收信号的能量大小随主动源激发能量变化,各台站能量变化主要体现在优势频率范围内;（3）接收台站台基越坚固不同激发能量均方根振幅比越小;（4）叠加一定数量的波形可以有效提高远台记录信号的信噪比,增加气枪激发能量则相对更容易提高信号信噪比。     

ATTENUATION OF COHERENT NOISE IN MARINE SEISMIC EXPLORATION USING VERY LONG ARRAYS *

The use of arrays to separate primary reflections from unwanted coherent seismic events is common practice in land seismic surveys. Very long source and receiver arrays have been used recently to reduce the effects of waterbottom multiples on marine seismic data. The source array consists of five uniformly spaced identical subarrays, each with five different airguns, where the distance between the subarrays may vary from 20 m56 m. The volume of each subarray is 10.3 1 (630 cu.in.) which gives a total volume of the array of 51.5 1 (3150 cu.in.) operated at a pressure of 14 MPa (2000 psi). In order to have a flexible receiver system it was decided to implement the extended receiver array in data processing by computing a weighted sum of two to five traces. The hydrophone cable consists of fifty-four channels with a group length of 50 m. Data shot with the superlong airgun array are processed by a combination of standard techniques and special procedures. In particular, the quality of the stack section is improved by using a weighted stack. The stack weights are computed by a program which takes into account the primary-to-multiple ratio. Comparisons with conventional data show significant improvements in data quality obtained by using the superlong airgun array. Examples show that the waterbottom multiples have been strongly attenuated and the deep seismic events have been enhanced. The combined array response function for dipping events is given in an appendix.