VERTICAL HETEROGENEITY AND MOVEOUT IN THE ONE-DIMENSIONAL MEDIUM*

Since the important contributions of Dürbaum and Dix, 30 years ago, velocity profile estimation procedures on horizontally layered and vertically heterogeneous media from seismic probing data have been based largely on hyperbolic moveout models and RMS and stacking velocity concepts. Re-examination of the fundamentals reveals that quantitative velocity heterogeneity and canonical valocity profiles have been implicit factors for moveout modelling and for profile inversion in the use of the Dix procedure. Heterogeneity h is the ratio (and v_RMS the geometric or harmonic mean) of the path-average and time-average velocities for a raypath or, in a more restricted sense, for the normal ray belonging to a velocity profile. The canonical profile for a given velocity profile or profile segment is a moveout-equivalent monotonically increasing ramp-like profile. The ramp or constant gradient in depth is the simplest velocity profile approximator which can explicitly accommodate velocity heterogeneity. A ramp model structure is detailed which facilitates moveout simulation and model parameter estimation, and the parametric effects are explored. The horizontal offset range is quantified for which this model can give good moveout approximations.     

用最大似然法进行波场分解和震相识别

本文以介质中弹性波场的传播特征和多元统计分析为基础,从数学上进行严格的推导,得到了模型与资料拟合的最大似然判据,以检验地震波列中是否存在某种特殊波型的能量。利用概率滤波可得分解后的地震图,从图中可以直接得到震相到时。用此方法处理P波资料,效果较好。对于剪切波资料,有时处理结果不十分理想。这可能是由于径向分量和横向分量之间的线性相关性造成的。对于因介质各向异性而形成的分裂剪切波,用此方法得到的S波到时为快S波到时。此方法对于单台三分向记录资料的处理是十分有效的。     

VELOCITY ANALYSIS IN THE p—x-PLANE FROM A SLANT STACK WAVEFIELD*

Classical methods of interpretation of reflection seismic data are such that interpretation and processing usually occur in the “collected” frame of reference. However, in recent times other data planes have gained increasing acceptance in seismology as a viable alternative. Through linear transformations applied to a record section, both the t—p- and p—x-planes can be produced. The r—p-domain may be obtained from the t—x-plane by a transformation known as slant stacking. Normal practice has been to do most of the data processing in the t—x-plane and then transforming to the r—p-plane. However, many of the procedures used in the t—x-domain can be modified for use in the t—p-plane to increase the coherence. Velocity inversion may be carried out either in the r—p-domain or further transformed to the p—x-plane where the modified Herglotz-Wiechert inversion may be applied. To perform the inversion, the t—p-wavefield is converted to a p—x-representation by the use of a new linear transformation technique, the cross-stack. By a simple sampling process along a particular p—x-trajectory, the Herglotz-Wiechert method can be used to reconstruct an acceptable velocity model of the subsurface. A comparison of derived velocity structures is made between that produced by the Herglotz-Wiechert technique and that of the Dix method.     

两台地磁场垂直分量日变形态的位相关系分析

本文用最大相关时间描写两台地磁场垂直分量日变形态的位相关系。通过恰当的数据处理方式,从两台（北京台和武汉台）地磁垂直分量的整点值计算出最大相关时间。它的平均值等于两台的地方时差。此外最大相关时间还有明显的以年为周期的季节性变化,变化幅度约为0.7小时。对最大相关时间作了傅氏分析,结果表明它含有周期为29.8±0.6天和13.7±0.2天,振幅约为0.06小时的周期变化。     

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.     

两台地磁场垂直分量日变形态的位相关系分析

本文用最大相关时间描写两台地磁场垂直分量日变形态的位相关系。通过恰当的数据处理方式,从两台（北京台和武汉台）地磁垂直分量的整点值计算出最大相关时间。它的平均值等于两台的地方时差。此外最大相关时间还有明显的以年为周期的季节性变化,变化幅度约为0.7小时。对最大相关时间作了傅氏分析,结果表明它含有周期为29.8±0.6天和13.7±0.2天,振幅约为0.06小时的周期变化。     

水库气枪震源产生的S波及其分裂

人工气枪震源在陆地水库可以有效激发S波,S波能量较强,与ML1.6天然地震相当。气枪可用于S波分裂研究,对布置在燕山隆起带的流动地震台的气枪信号进行了S波分裂参数分析,结果表明,快剪切波偏振优势方向为NWW和NNE向,偏振方向和断裂的性质密切相关。气枪是高度可重复性人工震源,利用气枪定点激发和定点接收有可能精确获取S波分裂参数随时间的变化规律,为地震预测探索实践提供可靠的物理途径     

双相介质中的地震波运动学特征

以 Biot理论为基础 ,推导了双相介质中弹性波的本构关系、运动方程及解。结果表明 :在均匀双相介质中可获得第一 P波、第二 P波和 S波。第一 P波近似于弹性介质中的 P波 ,第二 P波称为慢波 ,其速度一般小于 S波 ,但在个别介质模型中大于 S波。选择了合理的参数 ,对均匀双相介质中的波速进行了数值模拟     

PHASE DISTORTION DUE TO ABSORPTION IN SEISMOGRAMS AND VSP*

Absorption of seismic energy in the earth reduces amplitudes and changes phases of the propagating seismic waves. Amplitudes are usually recovered according to an estimated exponential decay curve, while phase distortions are generally disregarded. Therefore, accurate processing of seismic data requires a careful investigation of the relationship between absorption and phases. In this paper a procedure is suggested to achieve this goal, and some related topics are worked out. A method is outlined for computing synthetic seismograms and vertical seismic profiles with phase distortion due to absorption. The algorithm works in the frequency domain, and it provides for absorption according to the usual model of exponential decay of amplitude with distance. The absorption coefficient is a linear function of frequency and is related to the quality factor Q of the rocks. Complex seismic velocities are introduced and minimum-phase delay due to absorption is assumed for all cases considered. Methods for estimating Q profiles from seismic well surveys and seismic data are described. Comparison between field and synthetic data shows the effectiveness and benefits of the procedure. Some applications of the method to phase distortion recovery and wavelet processing are presented.     

ORIGIN OF THE ANOMALIES OF GRAVITY AND ITS VERTICAL GRADIENT OVER CAVITIES IN BRITTLE ROCK*

Investigations show that the depth range of the gravity method for detecting cavities is considerably greater than expected on the basis of theoretical calculations that consider only the depth and shape of the cavity. The cavity generates its own gravity field that is the sum of the activity of the cavity itself and the density changes caused in the surrounding rocks by the process of their destruction as a result of the cavity's expansion. The latter factor is decisive for the cavity's detection by the gravity method in several cases. The extent of the zones of change density in the surrounding rocks and their actual shape depend on many parameters. Thus, the determination of the depth of the cavity cannot be made by comparing measured and theoretical curves of gravity anomalies. The location of the center of gravity of the area disturbed by the cavity or upper boundary of its development towards the ground surface can be determined by denoting the singular point of gravity as a result of the downward analytical continuation of Ag. The Strakhov-Grigorieva-Lapina (1977) method fulfils the necessary conditions in the field.     

THE AMPLITUDE AND PHASE RESPONSE OF A SEISMIC VIBRATOR*

The amplitude and phase response of a simple model is compared with the performance of a real vibrator working in the field. The field results show a characteristic phase response which confirms that the real drive force applied to the baseplate and its load impedance is faithfully represented by the acceleration of the reaction mass. It follows that all the parameters necessary to calculate the load impedance and the true power dissipated in the earth can be measured at the output of the vibrator. It also follows that the current method of baseplate phase compensation should be reconsidered.     

VERTICAL SEISMIC PROFILING IN COAL*

The intellection of seismic wave propagation in coal measures demands direct observation of the wavefield progression. Two vertical seismic profiles with high spatial and temporal sampling, were recently recorded in the Sydney Basin coalfields as part of an experimental coal seismic program. Static corrections and interval velocities were obtained by an automated system to determine first kicks and pulse rise times. Upgoing and downgoing waves were separated in the f—k-plane using a novel technique of contour slice filtering. The isolated upgoing waves clearly display reflections from the major coal seams within the stratigraphic sequence. The downgoing wave spectra were subjected to attenuation analysis. The deduced specific quality factor Q for Permian coal measure rocks lies in the range 20–70. Similar estimates were obtained in the time domain from measurements of pulse broadening. Synthetic VSP seismograms, computed using an exact recursive formulation, are an indispensable aid to interpretation. They illustrate the filtering effects of coal seams and sequences, and the effects of the contribution of internal and free-surface multiple reflections in the recorded wavetrains.     

SOME IMPROVEMENTS IN THE INTERPRETATION OF VERTICAL ELECTRIC SOUNDING CURVES*

The process of VES interpretation is discussed, including the following points. (a) Preliminary interpretation by means of master curves. It is shown that the positions of the auxiliary points K and Q depend on the resistivity of the substratum. The interpretation is improved if the auxiliary points are determined separately for each master curve. (b) The individual parts of the measured curves are shifted in overlapping MN electrode positions so that the total sum of squares of the shifts is minimal. (c) The ambiguity may be reduced by means of supplementary information or assumptions on the resistivities. Fixing the resistivities is not always possible because discrepancies may arise between the ground measurement and the well-logging data. The simultaneous interpretation of several VES curves is recommended assuming constant resistivities. This assumption may be subsequently verified by means of the F-test. (d) A nonlinear algorithm is proposed for the determination of confidence intervals. As the multi-dimensional confidence intervals are often very complicated, it is recommended to construct only one-dimensional confidence intervals for the estimable parametric functions. (e) A ‘double-least-squares’ optimization technique is presented. The optimization is performed on the estimable parametric functions, and the individual parameters are determined so that the solution remains near the initial guess. This technique is faster than the singular value decomposition.     

THE INTERPRETATION OF THE VERTICAL FAULT PROBLEM IN GEOELECTRICS USING HANKEL INTEGRAL*

It is shown how to interpret, without curve-matching, Schlumberger resistivity soundings carried out with the array parallel to the surface trace of a vertical boundary plane separating two media with different resistivities (the vertical fault problem). To this end, it is demonstrated that the apparent resistivity function can be expressed as an Hankel-transformable integral that allows to invert the apparent resistivity curve into an associated resistivity transform curve. The study of the asymptotical properties of this function and of some mathematical properties of a related reduced transform function allows to realize a simple procedure for deriving the parameters of the model. In practice, this procedure consists in fitting a straight-line to a semi-logarithmic plot of the reduced transform function and in evaluating the intercept along the vertical axis and the slope.     

人造双行星与地球组成的太阳和行星际观测网

A Solar and interplanetary observational system composed of two artificial planets and the earth is suggested for the solar-terrestrial study and prediction.The two planets should move along the earth's orbit around the sun. The angle distance between the two planets and between them with the earth is 120°. This system can be used to improve greatly the short-term,midiu-mterm solar activity predictions and monitor an earth-toward coronal mass ejection,giving an accurate warning of a solar-terrestrial disturbance. Observational data obtained by this system would be very useful in a model work of solar wind and an evolution study of structures in solar atmosphere such as solar active regions and magnetic fields with various scales.     

利用小波分析研究河西地区地壳垂直形变特征

利用小波分析对河西地区1971—1979、1979—1983年的地壳垂直形变场进行计算,结果表明:三阶以内的小波分解结果与地质构造走向基本一致,说明河西地区地壳垂直运动的短波主要反应的是断层活动;四阶小波变换结果为地壳垂直形变的低频部分,反映的是大区域地质构造运动所引起的地壳垂直形变。1971—1979年,门源南部的地壳垂直运动相对剧烈,而到了1979—1983年,门源地区地壳垂直运动主要活动区域由门源南部地区转移到门源附近区域,这可能与门源6.4级地震孕育的区域构造活动有关。     

THE DISPLACEMENT AND STRESS FIELD GENERATED BY THE COLLAPSE IN PINGYI COUNTY,SHANGDONG PROVINCE,ON DECEMBER 25, 2015

A collapse happened in Pingyi County, Shandong Province, on December 25, 2015. The displacement field, stress field and Coulomb failure stress change on the Mengshan frontal fault generated by the collapse are calculated by using point collapse model in isotropic medium. The result shows that: (1) The maximum horizontal displacement is located at the center of the collapse with value of~18mm. The horizontal displacements are greater than 1mm within~5km of the collapse with its direction pointing to the collapse center. The maximum subsidence is located at the center of the collapse with the value of 4mm. The subsidence is greater than 1mm within ~3km of the collapse. The displacement field decays so rapidly that can be ignored at far away from the collapse for the shallow source, which caused local displacement field. (2) Influenced by the free surface, the contraction area stress within ~5km of the collapse with the order of 1000Pa and expansion area stress in farther away areas at depth of 2km are estimated. the expansion area stress of 1000Pa is estimated at the~5km from the collapse center. Then the expansion area stress decays to 100Pa at the distance of ~10km from the collapse. The maximum compressive and extensional principal stresses are estimated as 10000Pa at the depth of 2km. The compressive stress axes present radical direction pointing to the collapse within ~5km of the center. In farther away from the collapse, The extensional principal stress axes present radical direction pointing to the center of the collapse. With farther distance to the collapse, the compressive and extensional stress decay rapidly to the order of 100Pa. (3) The Coulomb failure stress on the northwestern part of the Mengshan frontal fault, which is known as active segment of the Mengshan frontal fault, is decreased by the collapse with maximum value of 2500Pa. Whereas, the Coulomb failure stress on the southeastern part of the Mengshan frontal fault, which is known as left-lateral normal slip fault segment in Quaternary period, is increased by the collapse with maximum of 2400Pa, to which attention would be paid in seismic hazard analysis.     

VERTICAL FREQUENCY EFFECT SOUNDING IN INDUCED POLARIZATION AND GALVANIC RESISTIVITY METHODS *

It has been shown that for a polarizable layer as a transmission medium there is an indirect proportionality between the frequency effect (fe) in induced polarization (IP) and the wave number of the electrical induced field. Making use of this relationship for a two layered earth a polarization transform function has been obtained. Since the mathematical expression for the polarization transform function is the same as that of the resistivity transform function, it is possible to make direct interpretation for IP frequency effect field curves. Thus, AA or QQ type resistivity sequences can be interpreted from induced polarization response of a horizontally stratified earth without resistivity extrema. A depth factor has been defined in order to obtain the true depth using the apparent depth. In this way, some electromagnetic effects between horizontal layers with different polarizabilities can partly be eliminated.     

THE USE OF VERTICAL LINE SOURCES IN ELECTRICAL PROSPECTING FOR HYDROCARBON*

Conventional electrical prospecting can be extended to the search for deep-seated hydrocarbon deposits, by using the steel casings of drill-holes as vertical line sources. These sources produce at depth a density of current higher than the density created by point sources located at the ground surface. Several tests have shown that the contrast of conductivity between resistive hydrocarbon deposits and the surrounding salt water produces relevant anomalies on a resistivity map obtained with vertical line sources, especially where there exists a superficial masking effect caused by a highly resistive layer. In a survey carried out in the USSR, combined measurements were performed, both with line source and with surface point sources. The detected residual resistivity anomaly roughly delineates the contours of the known hydrocarbon deposit.     

双相各向异性介质中弹性波方程的有限元解法及波场模拟

基于双相各向异性介质模型,首先推导了双相各向异性介质中弹性波传播的动力学方程及其Galerkin变分方程和有限元运动方程,然后给出了孔隙弹性波方程的有限元数值解法以及二维双相PTL介质中波场模拟的人为吸收边界条件. 最后,利用本文给出的有限元方法对双相PTL介质和双相各向同性介质中的弹性波传播进行了数值模拟. 结果表明:有限元方法和吸收边界条件有效、可行,在理想相界条件下,不论是从固体位移,还是从流体位移的波场快照都能看到明显的慢速拟P波;在黏滞相界情况下,能否观察到慢速拟P波,与含流体地层介质的耗散性质有关.对实际含流体介质,从流体位移分量的波场快照比从固体位移波场快照更容易观察到慢速拟P波.