Modelling and modal analysis of seismic vibrator baseplate

The vibroseis method must be extended to its limits as the search for oil and gas continues on land. To successfully improve vibroseis data quality, it is crucial to evaluate each element in the vibroseis data acquisition system and ensure that the contribution from each element is successful. Vibroseis systems depend greatly upon the ability of vibrators to generate synchronous, repeatable ground-force sweeps over a broad frequency range. This requires that the reaction mass and the baseplate of the vibrator move as rigid bodies. However, rigid-body motion is not completely true for high- frequency vibrations, especially for the vibrator baseplate. In order to accurately understand the motion of the vibrator baseplate, a finite element analysis model of the vibrator baseplate and the coupled ground has been developed. This model is useful for simulating the vibrator baseplate dynamics, evaluating the impact of the baseplate on the coupled ground and vibrator baseplate design. Model data demonstrate that the vibrator baseplate and its stilt structure are subject to six significant resonant frequencies in the range of 10–80 Hz. Due to the low rigidity of the baseplate, the baseplate stilt structure experiences severe rocking motions at lower frequencies and the baseplate pad experiences severe flexing motions at higher frequencies. Flexing motions cause partial decoupling, which gives rise to increased levels of harmonic distortion and less useable signal energy. In general, the baseplate pad suffers more bending and flexing motions at high frequencies than low frequencies, leading less efficiency in transmitting the useable energy into the ground.     

Break through the limits of vibroseis data quality

The vibroseis method has become the principal data acquisition method in land seismic exploration. It seems that this method has been extended to its limits as the search for energy resources continues. Many practical issues arising from field operations have remained theoretically unexplained, for example, variations in wavelet arrival time, inaccurate wavelet estimation and harmonics in the wavelet itself. The focus of this paper is the proposal of a new model, which is referred to as the vibrator‐coupled ground model, to simulate the filtering effects of a complex coupling system consisting of the coupling between the baseplate and the ground as well as the coupling between the captured ground mass near the vibrator baseplate and the surrounding earth. With this vibrator‐coupled ground model many of the practical issues mentioned above were reasonably addressed. Furthermore, it was demonstrated from experimental tests that both the pilot sweep and the weighted‐sum groundforce, when filtered by the vibrator‐coupled ground model, are proportional to the far‐field particle velocity whereas the unfiltered signals are not. The harmonics on the filtered weighted‐sum groundforce successfully maintain a proportional relationship with the harmonics seen in the far‐field signal.     

Imaging shallow structure with active-source surface wave signal recorded by distributed acoustic sensing arrays

Distributed acoustic sensing (DAS) is one recently developed seismic acquisition technique that is based on fiber-optic sensing. DAS provides dense spatial spacing that is useful to image shallow structure with surface waves. To test the feasibility of DAS in shallow structure imaging, the PoroTomo team conducted a DAS experiment with the vibroseis truck T-Rex in Brady’s Hot Springs, Nevada, USA. The Rayleigh waves excited by the vertical mode of the vibroseis truck were analyzed with the Multichannel Analysis of Surface Waves (MASW) method. Phase velocities between 5 and 20 Hz were successfully extracted for one segment of cable and were employed to build a shear-wave velocity model for the top 50 meters. The dispersion curves obtained with DAS agree well with the ones extracted from co-located geophones data and from the passive source Noise Correlation Functions (NCF). Comparing to the co-located geophone array, the higher sensor density that DAS arrays provides help reducing aliasing in dispersion analysis, and separating different surface wave modes. This study demonstrates the feasibility and advantage of DAS in imaging shallow structure with surface waves.     

The effects of inter‐vibrator interference on vibrator performance

The output from the hydraulic vibrators typically used for land seismic surveys is controlled by monitoring the acceleration measured by accelerometers mounted on the reaction mass and baseplate. The considerable energy output by such vibrators, which are coupled with the sensitivity of the accelerometers used, results in crosstalk if more than one vibrator is being used. In this paper, we present the results of a field experiment in which we measured the crosstalk between two adjacent vibrators. We found that the level of crosstalk was approximately ‐20 dB when the vibrators were adjacent but decreased with increasing frequency and separation. This result has implications for measurements of vibrator performance, source‐signature deconvolution, and in particular, estimates of the total energy output by a fleet of vibrators.     

相控震源对地震信号信噪比的改善研究

当测区噪声很强,采用组合震源工作仍然不能满足信噪比要求时,为进一步提高地震勘探信噪比,引入能形成定向地震波的相控震源.针对水平层状介质模型,对组合震源、相控震源地震进行了数值模拟,并定量计算了采用4,8,15,20单元的组合与相控震源地震数据的信噪比,同组合地震相比,相控地震使来自不同反射层的反射波信号的信噪比平均提高了154～990 dB.结果表明,相控震源地震得到的反射波信号信噪比高于组合地震,并且随着激震器数目增加,相控震源合成地震记录中的反射波地震数据信噪比明显提高.     

ON THE DETERMINATION OF THE DOWNGOING P-WAVES RADIATED BY THE VERTICAL SEISMIC VIBRATOR*

It is well recognized that in order to realize the full potential of the Vibroseis technique, one needs to ensure accurate phase locking and a meaningful cross-correlation. To achieve these two important objectives we require an accurate estimate of the compressional stress wave radiated by the vibrator into the ground. In this paper a simple method (subject of a patent application) is developed for predicting the compressional stress waves radiated by a vertical vibrator. The main feature of the proposed method is that it involves the field measurement of the acceleration of the reaction mass and the baseplate, respectively. The method is illustrated by computing the compressional stress waves generated by a typical vertical vibrator radiating into ice, chalk, sand, and mud. It is shown that for a seismic vibrator radiating into hard ground the pressure of the downgoing P-wave is 180° out of phase with the baseplate velocity. It is also shown that when the driving force of the seismic vibrator has a flat amplitude spectrum, the amplitude spectrum of the downgoing P-wave falls off by 6 dB/octave towards low frequencies.     

SIMULTANEOUS VIBROSEIS RECORDING1

An experiment was undertaken at BP's Fulbeck Geophysical test site to compare the viability of various simultaneous vibroseis recording techniques, which are often recommended as a means of improving data acquisition production rates for 3D seismic surveys. Of particular interest were: (a) the ability to separate the signals from each source during processing, (b) the generation and suppression of harmonics and (c) the effects of any source interaction. Two vibrators were deployed with a baseplate separation of 10 m, about a borehole containing a vertical array of geophones. Our analysis concentrated on the groundforce signals measured at each vibrator and the far-field signatures measured using a vertical geo-phone at a depth of 204 m. By comparing single vibrator records with similar but separated records from a simultaneous recording sequence, signal separability, harmonic suppression and vibrator interaction could be fully studied. Separated far-field signatures from simultaneous vibroseis methods using combinations of up and downsweeps exhibited unsuppressed harmonics and substantial energy from the undesired source which leaked through the correlation process. The ‘up/down’ method was capable of separating the signal from each source by only 12.7 dB, and is therefore unsuitable as a field technique. The variphase simultaneous vibroseis methods studied afforded some harmonic suppression and gave signal separations of about 30.0 dB. Use of variphase simultaneous vibroseis methods will compromise the quality of the data recorded, when compared with single-source acquisition methods. None of the simultaneous vibroseis methods tested provided adequate signal separation and, therefore, cannot be recommended as data acquisition techniques. The ‘alternate sweeping’ method coupled with multispread recording will give the desired improvement in data acquisition rates, while preserving the necessary quality of our seismic data.     

地震作用下均质土坡动力特性的振动台试验研究

设计并完成比例尺1∶100的边坡振动台模型试验,讨论模型的相似关系、传感器的布置及模型的建造,并编制相应的动荷载加载方案。通过输入不同类型、不同幅值、频率的动荷载,分析模型边坡在地震作用下的动力响应规律以及地震动参数对边坡动力特性和动力响应的影响。试验结果表明,在坡体的表面和坡内的竖直方向上,加速度放大系数均随着高程增加而明显增大。当输入不同压缩比的地震波时,压缩比越大坡体的动力响应越明显,即随着输入动荷载的频率增加,越接近土体的自振频率加速度放大效应越明显;在坡体的同一高程处,坡面的加速度放大系数略大于坡体内的加速度放大系数,表现出一定的趋表效应,同时随着输入地震波振幅的增加,加速度放大系数整体出现递减的趋势。实验结果分析有助于揭示土质边坡在地震作用下的失稳破坏机制,为今后边坡工程的抗震设计提供积极的参考。     

Evaluation of hard rock spectral models for the Taiwan region on the basis of the 1999 Chi–Chi earthquake data

We analyze the ability of different spectral models to describe the frequency content of ground motion during the 1999 Chi–Chi earthquake (M_W=7.6, Taiwan) and two large (M_L=6.8) aftershocks. The spectral models evaluated include the one-corner model of Brune applied with various key parameters (seismic moment and stress drop), and the two-corner-frequency models proposed for eastern North America [Bull. Seismol. Soc. Am. 83 (1993) 1778] and California [Bull. Seismol. Soc. Am. 90 (2000) 255]. The ground-motion spectra predicted by these spectral models for hypothetical very hard rock site were compared with the Chi–Chi earthquake data obtained on rock (class B) and soft rock or very dense soil (class C) sites. The approach also allows us evaluating the generalized empirical amplification function for class B and C sites in the region.

It has been found that, the amplitude spectra of recorded ground acceleration (the mainshock and aftershocks) for frequencies larger than 0.3–0.4 Hz agree with the modelled two-corner-frequency spectra calculated using the model proposed for California. The single-corner-frequency model also provides a good agreement with the observations when using so-called ‘short-period seismic moment’ [Phys. Earth Planet. Interiors 37 (1985) 108] instead of the reported values obtained from long-period waves. The key parameters used in the single-corner model coincide with parameters of subsources evaluated for the complicated mainshock source. Therefore, it is possible to confirm the suggestion that the short-period seismic waves, at least for the thrust earthquakes, are generated mainly from the fracture of small-scale heterogeneities. The use of two-corner-frequency source model for earthquake spectrum that is based on long-period seismic moment value is equivalent, for frequencies larger than 0.3–0.4 Hz, to the use of single-corner-frequency model that is based on the parameters of major subsource.     

RESEARCH ON PASSIVE SERVO SEISMIC ROTATIONAL ACCELERATION SENSOR WITH LARGE DAMPING

Many strong motion records show that under the strong seismic vibration of, the torsional disfigurement of building structures is a common and serious damage. At present, there are no special sensors for measuring seismic rotation in the world. Most of the experts obtain rotational components through observing deformation, theoretical analysis and calculation. The theory of elastic wave and source dynamics also prove the conclusion that the surface of the earth will rotate when an earthquake occurs. Based on a large number of investigations and experiments, a rotational acceleration sensor was developed for the observation of the rotational component of strong ground motions. This acceleration sensor is a double-pendulum passive servo large-damped seismic rotational acceleration sensor with the moving coil transducer. When an earthquake occurs, the seismic rotational acceleration acts on the bottom plate at the same time. The magnetic circuit system and the middle shaft fixedly connected to the bottom plate follow the bottom plate synchronous vibration, and the moving part composed of the mass ring, the swing frame and the moving ring produces relative corners to the central axis. The two working coils mounted on the two pendulums produce the same relative motion with respect to the magnetic gaps of the two magnetic circuits. Both working coils at this time generate an induced electromotive force by cutting magnetic lines of force in the respective magnetic gaps. The generated electromotive forces are respectively input to respective passive servo large damper dynamic ring transducer circuits and angular acceleration adjusting circuits, and the signals are simultaneously input to the synthesizing circuit after conditioning. Finally, the composite circuit outputs a voltage signal proportional to the seismic rotational acceleration to form a seismic rotational acceleration sensor. The paper presents the basic principles of the rotational acceleration sensor, including its mechanical structure diagram, circuit schematic diagram and mathematical models. The differential equation of motion and its circuit equation are derived to obtain the expressions of the main technical specifications, such as the damping ratio and sensitivity. The calculation shows that when the damping ratio is much larger than 1, the output voltage of the passive servo large damping dynamic coil transducer circuit is proportional to the ground rotation acceleration, and the frequency characteristic of bandpass is wider when the damping ratio is larger. Based on the calibration test, the dynamic range is greater than or equal to 100dB and the linearity error is less than 0.05%. The amplitude-frequency characteristics, the phase-frequency characteristics and their corresponding curves of the passive servo rotational acceleration sensor are acquired through the calculations. Based on the accurate measurement of the micro-vibration of the precision rotating vibration equipment, the desired result is obtained. The measured data are presented in the paper, which verify the correctness of the calculation result. The passive servo large damping rotational acceleration sensor has simple circuit design, convenient operation and high resolution, and can be widely applied to seismic acceleration measurement of earthquake or structure.     

SEISMIC VIBRATOR MODELLING1

The wavefield in, and at the surface of, a homogeneous, isotropic, perfectly elastic half-space, excited by a traction distribution at the surface of the medium is investigated. The emitted wavefield is a spatial convolution of the surface tractions and the spatial impulse response. The properties of the wavefield in the far-field of the medium are derived and it is shown that the far-field particle velocity is essentially equal to a weighted sum of the time derivative of the integrated surface tractions, that is, of the components of the ‘ground force’. The theory is valid for an arbitrary geometry and orientation of the surface tractions, and is independent of the boundary conditions at the surface of the medium. The surface tractions are related to a source that consists of a mass distribution with an arbitrary force distribution imposed upon it. A boundary condition is introduced that accounts for the mass load and the forces applied to it but neglects vibrations within the mass. The boundary condition follows from the equation of motion of the surface mass load. The theory is applied to the Vibroseis configuration, using a P-wave vibrator model with a uniformly distributed force imposed on top of the baseplate, and assuming that horizontal surface traction components are absent. The distribution of displacement and stress directly underneath the baseplate of a single vibrator and an array of vibrators is investigated. Three different boundary conditions are used: (1) assuming uniform pressure, (2) assuming uniform displacement, (3) using the equation of motion of the baseplate as a boundary condition. The calculations of the distribution of stress and displacement over the plate for different elastic media and several frequencies of operation show that only the results obtained with the mixed boundary condition agree with measurements made in the field. The accuracy of three different phase-feedback signals is compared using synthetic data. Baseplate velocity phase-feedback leads to huge deviations in the determination of the far-field wavelet; reaction mass acceleration phase-feedback looks stable but neglects the differentiating earth filter; and phase-feedback to a weighted sum of baseplate and reaction mass accelerations becomes unstable with increasing frequency. The instability can be overcome using measurements over the whole baseplate. The model can be extended to a lossy layered earth.     

Investigation into dynamic response of regional sites to seismic waves using shaking table testing

This study addresses the changes in acceleration,pore water pressure and Fourier spectrums of different types of seismic waves with various amplitudes via large-scale shaking table tests from two sites:a sand-containing regional site and an all-clay site.Comparative analyses of the test results show that the pore water pressures in sand-soil layers of the regional site initially increase and then decrease as the amplitudes of the seismic accelerations increase.The actions of the vertical and vibrational seismic waves contribute to greater pore water pressures.The amplification coefficient of the sand-layer regional site becomes smaller as the seismic waves grow stronger,so that both sites are capable of filtering high frequencies and amplifying low frequencies of seismic waves.This is more apparent with the increase in the peak value of the acceleration,and the natural vibration frequencies of both sites decrease with the transmission of the seismic waves from the basement to the ground surface.The decreasing frequency value of the sand-containing regional site is smaller than that of the all-clay site.     

地震SV波斜入射下沉管隧道的地震响应分析

在ABAQUS黏弹性人工边界时域波动方法的基础上,首先运用等效应力输入方法实现地震SV波倾斜入射,半空间算例验证该方法具有较好的计算精度,进而基于所建立的斜入射方法研究地震波斜入射对海河沉管隧道地震响应的影响。计算结果表明:SV波斜入射情况下,沉管隧道的地震响应规律与垂直入射时具有明显差异;随入射角增加,沉管隧道结构应力增大,应力较大点出现在沉管隧道的四个角点及隔墙与底板、顶板的连接处,其中中隔墙为最薄弱点;随入射角增加,侧墙和隔墙的相对最大水平位移增大,其中中隔墙位移最大;随入射角增加,沉管隧道结构竖向加速度峰值明显增大。因此在沉管隧道结构抗震设计中应考虑地震波斜入射的影响。     

Retrieval of local surface wave velocities from traffic noise – an example from the La Barge basin (Wyoming)

In regions where active source seismic exploration is constrained by limitations of energy penetration and recovery, cost and logistical concerns, or regulatory restrictions, analysis of natural source seismic data may provide an alternative. In this study, we investigate the feasibility of using locally‐generated seismic noise in the 2–6 Hz band to obtain a subsurface model via interferometric analysis. We apply this technique to three‐component data recorded during the La Barge Passive Seismic Experiment, a local deployment in south‐western Wyoming that recorded continuous seismic data between November 2008 and June 2009. We find traffic noise from a nearby state road to be the dominant source of surface waves recorded on the array and observe surface wave arrivals associated with this source up to distances of 5 kms. The orientation of the road with respect to the deployment ensures a large number of stationary points, leading to clear observations on both in‐line and cross‐line virtual source‐receiver pairs. This results in a large number of usable interferograms, which in turn enables the application of standard active source processing methods like signal processing, common offset stacking and traveltime inversion. We investigate the dependency of the interferograms on the amount of data, on a range of processing parameters and on the choice of the interferometry algorithm. The obtained interferograms exhibit a high signal‐to‐noise ratio on all three components. Rotation of the horizontal components to the radial/transverse direction facilitates the separation of Rayleigh and Love waves. Though the narrow frequency spectrum of the surface waves prevents the inversion for depth‐dependent shear‐wave velocities, we are able to map the arrival times of the surface waves to laterally varying group and phase velocities for both Rayleigh and Love waves. Our results correlate well with the known geological structure. We outline a scheme for obtaining localized surface wave velocities from local noise sources and show how the processing of passive data benefits from a combination with well‐established exploration seismology methods. We highlight the differences with interferometry applied to crustal scale data and conclude with recommendations for similar deployments.     

天津滨海场地土动力学参数研究

天津滨海地区是典型的软土区,属晚第四系以来的浅海相沉积软土。本文收集了该地区地震安全性评价报告中的土动力学参数的实验资料,并补充了15个钻孔的测试实验,统计分析了淤泥质粘土、淤泥质粉质粘土、粉土、粉质粘土、粘土、粉砂、细砂等7种土类的实测土动力学参数,给出了它们在不同深度下的动剪切模量比和阻尼比的统计值。然后,选取了2个典型钻孔并建立了土层地震反应分析模型,分别运用本文统计值、94规范值和袁晓铭等(2000)的推荐值进行土层地震反应计算,从反应谱形状、地表峰值加速度和反应谱特征周期等方面说明了本文统计值的适用性和针对性。本文的研究成果对天津滨海地区各类工程建设的场地地震安全性评价工作具有一定的借鉴和参考价值,且有利于该地区工程场地地震安全性评价工作的开展。     

地下综合管廊土-结构接触面参数对地震动力响应的影响数值分析

为研究土-结构接触面参数对地下综合管廊地震动力响应特征的影响,建立动力有限元数值模型,模型边界采用激励侧固定边界、远离激励侧黏性边界、其余侧自由场边界的优化组合动力边界,土体本构采用HSS模型,接触面采用改进Goodman单元,动力荷载考虑三种情况(Rayleigh波的作用、底部激励了美国加利福尼亚Upland地震波以及前两者的共同作用),分别研究不同地震动输入、接触面折减系数的改变对综合管廊内力及加速度的影响。研究结果表明:在相同的折减系数条件下,与静力作用相比,动力作用下的结构内力明显增大,综合管廊设计时应考虑地震荷载作用下内力增大的情况;随着界面折减系数的增加,正弯矩极值减小,负弯矩极值增大,加速度峰值增大;在相同接触面折减系数条件下,底部地震波输入产生的结构内力极值显著高于仅有Rayleigh波输入的情况;考虑Rayleigh波和地震波共同作用条件下,引起的管廊结构内力极值与仅考虑底部地震波输入时的结构内力极值差异不大。研究成果可供地下综合管廊结构地震响应精细化数值模拟及抗震设计参考。     

Geophone‐ground coupling with flat bases

Seismic acquisition can be costly and inefficient when using spiked geophones. In most cases, such as the desert, the most practical solution is the use of flat bases, where geophone‐ground coupling is based on an optimal choice of the mass and area of contact between the receiver and the ground. This optimization is necessary since areas covered by sand are loose sediments and poor coupling occurs. Other cases include ground coupling in stiff pavements, for instance urban areas and ocean‐bottom nodes. We consider three different approaches to analyse coupling and model the geophone with a flat base (plate) resting on an elastic half‐space. Two existing models, based on the full‐wave theory, which we refer to as the Wolf and Hoover‐O'Brien models, predict a different behaviour with respect to the novel method introduced in this work. This method is based on the transmission coefficient of upgoing waves impinging in the geophone‐ground contact, where the ground is described as an anelastic half‐space. The boundary conditions at the contact have already been used to model fractures and are shown here to provide the equation of the damped oscillator. This fracture‐contact model depends on the stiffness characteristic of the contact between the geophone base plate and the ground. The transmission coefficient from the ground to the plate increases for increasing weight and decreasing base plate area. The new model predicts that the resonant frequency is independent of the geophone weight and plate radius, while the recorded energy increases with increasing weight and decreasing base plate area (as shown from our own experiments and measurements by Krohn) which is contrary to the theories developed by Wolf and Hoover‐O'Brien. The transient response is obtained by an inverse Fourier transform. Optimal geophone‐ground coupling and energy transmission are required, the first concept meaning that the geophone is following the motion of the ground and the second one that the signal is detectable. As a final example, we simulate seismic acquisition based on the novel theory, showing the differences between optimal and poor ground‐to‐geophone energy transmission.     

地震动非一致性对隧道动力响应的影响分析

针对隧道的抗震安全性,借助有限元软件平台对其地震时程响应进行非线性有限元分析。利用弹性模型考虑混凝土隧道在循环荷载作用下的拉、压应力-应变关系,同时对土体采用MohrCoulomb模型,选用人工边界作为边界条件,以人工合成的多点地震加速度时程作为地震动输入,分别考虑一致输入、行波效应、相干效应及行波加相干效应对隧道的影响。分析结果表明,与一致激励的计算结果相比,非一致地震激励会显著增加隧道结构的内力和位移响应,从而对隧道抗震产生不利影响。研究结论可为长隧道结构的抗震设计和分析提供科学依据。     

Effective elimination of subharmonic ghost events from vibroseis data

Harmonic or subharmonic noise is often present in vibroseis data as reverberation‐like, laterally coherent bands occurring parallel to and before or after, the main events. Such periodic noise is typically generated during the standard correlation process when the actual source signal travelling through the subsurface is, for whatever reason, different from the desired source signal, i.e., the pilot‐sweep controlling the baseplate and used for correlation. A typical cause can be that harmonic or subharmonic frequency partials are generated in addition to the vibroseis sweep's desired fundamental frequencies. These harmonics produce strong ‘ghost events’ during correlation of the geophone trace with the pilot‐sweep, originating from additional correlations between the fundamental and harmonic frequencies. Especially subharmonic ‘ghosts’ will overlap with ‘good’ fundamental signals, since for typically used up‐sweeps they are folded to later traveltimes, where the signal/noise‐ratio is already lower, thus aggravating or preventing a reliable interpretation of possible later reflections. Here, a method is introduced to remove these unwanted noise trains (with only negligible impact on the fundamental signal) by transforming the seismogram traces into a so‐called ‘(sub)harmonic domain’. In this domain, the respective harmonic noise portions are focused and separated from the fundamental signals, enabling easier detection and appropriate suppression. After back‐transformation to the x‐T domain, the records are free from the corresponding harmonic contamination and can then be processed as usual. The method operates in a data‐driven fashion, i.e., the traces are not uniformly processed but are processed depending upon their actual (sub)harmonic content. The decontamination procedure can be applied universally, i.e., to uncorrelated/correlated and/or vertically unstacked/stacked data either in a manual, semiautomated or fully automated manner. The method works perfectly for synthetic vibroseis traces with or without harmonic/subharmonic portions. The application to real, crustal‐scale vibroseis records that were acquired in 2006 in the Dead Sea region, Israel and that were severely contaminated by subharmonic ground‐roll ghosts covering reflectivity from the basement to the Moho, shows the robustness and success of the presented method.     

Seismic airgun exploration of continental crust structures

The active sources generate seismic waves transmitting appropriate through the deep is underground key and can be used to image Abstract high-resolution subsurface structures.Therefore,an seismic source the factor to active source exploration.In order to study the structure of continental crust and its temporal variations,we selected an artificial seismic source generated from releasing air bubbles in water(airgun source hereinafter)out of a variety and of artificial sources like the is explosion,new electronic sparkers,source hammering,eccentric proven vibration,be heavy-duty train vibration,vibroseis etc.Airgun Three source Fixed a type of artificial that have been to environmentally friendly,safe,and highly efficient.Airgun western Signal China Transmission and Stations(FASTS)have been for built a few years ago in Yunnan,Xinjiang,and Gansu provinces in have been continuously them running several years.Seismic waves generated away by the the airgun sources are highly seismic reproducible waves and stacking in of can produce can good seismograms on 1300 stations km far from source,for instance,an produced Xinjiang FASTS be well about recorded 60 nearly away after 5000 stacking,China covering area of 6 million km2 and penetrating down to of a depth of km.Establishing about 10 FASTSs in would enable long-term illuminate continuous subsurface underground structures,monitoring can all 9.6 million km2 of land area.Treating from airgun sky active sources as lanterns to we achieve the situation with"Beidou surveys the and lantern illuminates underground".