Phase-array decomposition of a seismic section

A seismic re fraction/wide-angle reflection profile is analysed for the presence of correlated events ('phases'). The correlation problem is formulated in terms of temporally, spatially and frequency-local complex covariances. For robustness, the method concentrates on phase rather than amplitude information. This allows a computationally efficient algorithm that can make allowance for signal correlation length and can model curved wavefronts. A statistical test based on residual phase misfit across the analysed subarray is used to assess the probability that a detected event represents a real correlated signal.
With our chosen analysis parameters and confidence level (over 99.9 per cent). 1222 events were detected in the data. Using simple techniques based on 1-D earth models, detected events are associated with a small number of particular wave types. In this way, we have succeeded in classifying almost 95 per cent of the detected events. Those that remain describe those components of the data that are inconsistent with our simple ray paths in the 1-D assumption and with our prescribed tolerance. These include reverberations, near-surface guided waves and reflected waves from strongly laterally inhomogeneous structures. According to our modelling, about 25 per cent of the detected events are consistent with simple P -wave reflected energy, and these are to a very large extent (over 85 per cent) distinct from all the other wave-type models we have used. A direct mapping of the detected events into the offset-depth domain reveals dear internal and external consistencies among the detections for the various wave types. Estimated earth structure is consistent with models from previous analyses based on much larger data sets.
We have thus succeeded in extracting correlated events from the data and decomposing these, approximately but meaningfully, into distinct classes (ray paths)     

Noise reduction and detection of weak, coherent signals through phase-weighted stacks

We present a new tool for efficient incoherent noise reduction for array data employing complex trace analysis. An amplitude-unbiased coherency measure is designed based on the instantaneous phase, which is used to weight the samples of an ordinary, linear stack. The result is called the phase-weighted stack (PWS) and is cleaned from incoherent noise. PWS thus permits detection of weak but coherent arrivals. The method presented can easily be extended to phase-weighted cross-correlations or be applied in the τ p domain. We illustrate and discuss the advantages and disadvantages of PWS in comparison with other coherency measures and present examples. We further show that our non-linear stacking technique enables us to detect a weak lower-mantle P -to- S conversion from a depth of approximately 840 km on array data. Hints of an 840 km discontinuity have been reported; however, such a discontinuity is not yet established due to the lack of further evidence.     

Site effects at a vertical accelerometer array near Paducah, Kentucky

The downhole vertical accelerometer array VSAP near Paducah, KY, consists of three-component accelerometers at the surface, the top of the McNairy Formation (−41 m), and the top of the Paleozoic bedrock (—99 m). The array is at the northern end of the Mississippi Embayment, and it was installed to verify the ground-motion modeling for the site, assuming a significant earthquake in the New Madrid Seismic Zone. Accelerograms from 4.2 and 2.0 m_b earthquakes were used to check aspects of the modeling pertaining to linear behavior of the soil column, and to review the soil column models derived by drilling and geotechnical methods and through the use of high-resolution P- and SH-wave seismic refraction in reflection techniques. Results of the study indicate that for the linear case the soil column models derived by the two techniques are equivalent, and that the most important boundary in the soil column, with respect to amplification of the ground motions, is the interface between the limestone bedrock and soil.     

Applicability of a Backprojection Algorithm to Reconstruct Images of Subsurface Horizontal Planes for Laboratory Experiments in Electrical Resistance Tomography

This paper describes the use of a backprojection algorithm to reconstruct subsurface images of the electrical resistivity in horizontal planes parallel to the surface. The algorithm can be applied to detect buried objects such as tanks or pipes and possible leakages from them. Two imaging strategies are compared: juxtaposition of vertical planes, and 3D reconstruction from the sensitivity matrix corresponding to the entire volume whose surface is explored. The electrode arrays used for voltage measurement are the dipole–dipole array and a modified Schlumberger array. A personal computer controls current injection, electrode switching, and voltage detection. The system injects 1 kHz, 20 V peak-to-peak square waveforms, thus avoiding electrode polarization effects. Experimental laboratory measurements show that the algorithm detects localized objects such as an insulating sphere and a conductive cylinder immersed in water. Furthermore, covering half of the cylinder by a rubber sleeve to simulate a nonconductive leak, yields a distinct image for the leak. The backprojection algorithm does not need any regularization parameter and it is very fast in inverting the sensitivity matrix because it approximates the inverse matrix by its transposed. The dipole–dipole array usually yields a lower overall pixel error than the modified Schlumberger array but both allow the detection of simulated underground leaks.     

Simultaneous identification of time and space variant dynamic soil properties during the 1995 Hyogoken-Nanbu earthquake

Time and space variant soil properties at a liquefied site were simultaneously identified in the time domain by using borehole array strong motion records. During soil liquefaction at a site, soils usually show a wide variety of non-linear behavior along the depth as well as non-stationary behavior. Strong ground motion records were obtained at Port Island borehole array observatory, Kobe, during the 1995 Hyogoken-Nanbu earthquake. In this study, the instrumented soil was modeled by the equivalent linear MDOF system, and an extended Kalman filter with local iteration was employed for the identification of the soils. The identification process was successfully conducted, and the stress–strain relationships of the soils at the liquefied site were obtained from different depths all at once.     

Observed Evolution of Linear and Nonlinear Effects at the Dahan Downhole Array,Taiwan: Analysis of the September 21, 1999 M 7.3 Chi-Chi Earthquake Sequence

— In this paper, the site characteristics of the Dahan downhole array are studied by analyzing the September 21, 1999 M 7.3 Chi-Chi earthquake sequence including the main shock and some aftershocks. The four-level array (0 m, 50 m, 100 m and 200 m) is located to the north of Hualien City in eastern Taiwan. Polarization analysis is used to check the orientation errors of the seismometers at different levels of depth. If the surface instrument is chosen as reference, the angle between the major polarization axes of the surface and any downhole records is the orientation error that must be corrected for the downhole accelerographs. The orientation errors at depths of 50 m, 100 m and 200 m are 32°, 120° and –84°. After the corrections, the coherency between the surface and downhole records is substantially improved. Spectral ratio analysis shows that the predominant frequency of the Chi-Chi main shock shifts to a lower frequency. We also simulate ground motions at different depths by using the Haskell method with a linear velocity structure model. The record at surface is chosen as the input motion. Compared with the observed data, ground acceleration can be well reproduced for the aftershocks (weak-motion events) of the September 21, 1999 M 7.3 Chi-Chi earthquake. However, for the Chi-Chi main shock, the synthetic waveform cannot match well with the observation neither in amplitude nor in phase. This indicates that large ground shaking probably induced the nonlinear site effect at that time, and the model used cannot support it.Acknowledgement. The authors would like to express their thanks to Dr. L.F. Bonilla and one anonymous reviewer for their valuable suggestions. This research was supported by the National Science Council under grant number NSC 89-2921-M-194-007. The Institute of Earth Sciences, Academia Sinica supplied the strong-motion data. The support of these organizations is gratefully acknowledged.     

“探宝号”船可变阵列震源的实现及应用

2001年，广州海洋地质调查局对“探宝号”船原有震源更新升级为阵列可变的BOLT震源系统。随着更新改造工程的成功，该船已成为具备承担高分辨率、常规、低频深水油气地震勘探的多功能地震调查船。本文介绍了地震勘探震源的基本工作原理以及三种阵列的分布、波形、频谱等特性。     

Interaction of cnoidal waves with an array of vertical concentric porous cylinders

An array of large concentric porous cylinder arrays is mounted in shallow water exposed to cnoidal waves. The interactions between waves and cylinders are studied theoretically using an eigenfunction expansion approach. Semi-analytical solutions of hydrodynamic loads and wave run-up on each cylinder are obtained using first approximation to cnoidal waves. The square array configuration of four-legged identical concentric porous cylinder is investigated in present study. Numerical results reveal the variation of dimensionless wave force and wave run-up on individual cylinder with angle of incidence, porosity parameter, spacing between outer and inner cylinders, spacing between concentric porous cylinders and wave parameter. Different mechanism of wave force is found under different range of scattering parameter.     

利用上海台阵对常熟2001年2002年震群的精定位

2001年6月上海台阵正式运行,上海台阵的建立提高了上海及周边地区的地震监测能力.本文着重在定位方法和震相判读上,充分利用上海台阵整体观测系统的优势和台阵独特的数据处理方法,对常熟2001年和2002年二次小地震序列进行了精确定位.对精确定位的结果进行了对比与分析,并以此为基础,对常熟地区及周边的地质构造等进行了解释与推断.