利用ERT数据推求非均质多孔介质渗透系数初探

高密度电阻率成像法（ERT）能通过改变间距来对不同尺度上的地质体特征进行描述，近年来广泛地应用于地质、工程等领域。本文从电流场和渗流场的相似性出发来研究电阻率和渗透系数之间的关系，简要介绍了利用ERT技术所得的电阻率数据推求渗透系数的基本原理，并利用张家港大新和塘桥两处的ERT测量数据推求出渗透系数。与两地抽水试验结果和岩性取样分析的结果分别进行对比，结果较为吻合。     

Heterogeneity in the mantle—its creation, evolution and destruction

Small-scale seismic heterogeneity exists at different levels in the lower mantle, and is detected by methods that analyze scattered–not direct–energy from natural and artificial sources. Its vertical distribution, association with subduction, and its ≤ 10-km characteristic scale length strongly suggest that it is chemical/petrological in nature and originally created by melting and differentiation during mid-ocean ridge formation. What is of interest is that the scale lengths of both upper and lower mantle seismic heterogeneity are similar, which supports the view of a common origin explored here. Unlike the lower mantle however, which is broadly homogeneous in structure, the upper mantle contains things that trap and impede the dispersal and re-mixing of heterogeneity: continental crust, lithosphere and cratonic roots. These probably control the depths, the longevity and the age of heterogeneities at shallow mantle levels, and suggest that heterogeneities observed in continental mantle lithosphere are probably old, trapped by the process that grows continental roots. Alternatively, if crustal heterogeneity is controlled by the details of a magmatic process, it must either be somehow continually renewed, for which there is no recognizable surface expression, or it must be depleted over time and the present is a time when, by luck, we may still witness it.     

Effects of trench-zone scattering on receiver functions over a subduction zone: A 3-D finite-difference modelling study

As Morozov [Morozov, I. B. (2004). Crustal scattering and some artefacts in receiver function images. Bull. Seismol. Soc. Am., 94 (4), 1492–1499.] suggested, for a teleseismic array targeting subducting crust in a zone of active subduction, scattering from the strong horizontal velocity heterogeneity beneath the trench zone itself produces subhorizontally-propagating waves that should be observed as coherent dipping events in receiver functions (RF). Due to similar RF delay times and moveouts, these events could be difficult to distinguish from backscattered P- and S-wave modes. To further verify this suggestion, we performed a full-waveform, 3-D visco-elastic finite-difference modelling of teleseismic wave propagation within a simplified model of a subduction zone. The synthetics show strong scattering from the area beneath the trench, dominated by the mantle and crustal P-waves propagating at 6.2–8.1 km/s and slower. These scattered waves occupy the same time and moveout intervals as the backscattered converted modes, and also have similar amplitudes. Although their amplitude decay characters are different, the uncertainty in the knowledge of the velocity and density structure of the subduction zone could make distinguishing between these modes difficult. However, under minimal assumptions, recent observations of receiver function amplitudes decreasing away from the trench support the interpretation of (sub-) trench-zone scattering.Although still limited in its representation of crustal heterogeneity, 3-D modelling suggests that scattering from near-Moho crustal structures plays a key role in the formation of teleseismic wavefields. Recognition of scattered noise in teleseismic records could help to constrain major crustal structures, particularly those with strong horizontal velocity contrasts at near-Moho depths, such as crustal sutures, subduction fault zones, and mountain roots. Matching of the observed arrivals with wavefield synthetics could help constrain the locations and parameters of such structures and also help substantiate the interpretations.     

γ射线康普顿背散射扫描仪测量的信噪比研究

本文详细给出了ＣＢＳ测量中噪声的物理来源和降低噪声的有效措施,还对信噪比进行了定量的分析和讨论。     

Application of tomographic inversion in studying airglow in the mesopause region

It is pointed out that observations of periodic nightglow structures give excellent information on atmospheric gravity waves in the mesosphere and lower thermosphere. The periods, the horizontal wavelengths and the phase speeds of the waves can be determined from airglow images and, using several cameras, the approximate altitude of the luminous layer can also be determined by triangulation. In this paper the possibility of applying tomographic methods for reconstructing the airglow structures is investigated using numerical simulations. A ground-based chain of cameras is assumed, two-dimensional airglow models in the vertical plane above the chain are constructed, and simulated data are calculated by integrating the models along a great number of rays with different elevation angles for each camera. After addition of random noise, these data are then inverted to obtain reconstructions of the models. A tomographic analysis package originally designed for satellite radiotomography is used in the inversion. The package is based on a formulation of stochastic inversion which allows the input of a priori information to the solver in terms of regularization variances. The reconstruction is carried out in two stages. In the first inversion, constant regularization variances are used within a wide altitude range. The results are used in determining the approximate altitude range of the airglow structures. Then, in the second inversion, constant non-zero regularization variances are used inside this region and zero variances outside it. With this method reliable reconstructions of the models are obtained. The number of cameras as well as their separations are varied in order to find out the limitations of the method.     

Optical tomography of the aurora and EISCAT

Tomographic reconstruction of the three-dimensional auroral are emission is used to obtain vertical and horizontal distributions of the optical auroral emission. Under the given experimental conditions with a very limited angular range and a small number of observers, algebraic reconstruction methods generally yield better results than transform techniques. Different algebraic reconstruction methods are tested with an auroral are model and the best results are obtained with an iterative least-square method adapted from emission-computed tomography. The observation geometry used during a campaign in Norway in 1995 is tested with the are model and root-mean-square errors, to be expected under the given geometrical conditions, are calculated. Although optimum geometry was not used, root-mean-square errors of less than 2% for the images and of the order of 30% for the distribution could be obtained. The method is applied to images from real observations. The correspondence of original pictures and projections of the reconstructed volume is discussed, and emission profiles along magnetic field lines through the three-dimensionally reconstructed arc are calibrated into electron density profiles with additional EISCAT measurements. Including a background profile and the temporal changes of the electron density due to recombination, good agreement can be obtained between measured profiles and the time-sequence of calculated profiles. These profiles are used to estimate the conductivity distribution in the vicinity of the EISCAT site. While the radar can only probe the ionosphere along the radar beam, the three-dimensional tomography enables conductivity estimates in a large area around the radar site.Former address: MPE Garching     

二维速度随机分布逐步迭代射线追踪方法

基于Ｓｎｅｌｌ定理，研究一种新的射线追踪方法──逐步迭代射线追踪法。首先，从一端出发，根据射线路径上任意连续三点均满足Ｓｎｅｌｌ定律，利用一个近似公式逐段迭代，求取中间折射点，从而实现逐步迭代计算。该方法的追踪路径结果符合射线追踪要求，计算速度快，方法精度高，可以克服射线追踪的路径非唯一性；并且在追踪过程中可以给出透射走时，为层析成像方法提供了一种有效的射线追踪过程。     

Vector-scattering of elastic waves by directional structural space gradients

A nonstochastic and noniterative theory of vector scattering in inhomogeneous media is presented. The elastodynamic vector wave-equation for 3D inhomogeneous media is solved for a weak heterogeneity at the high-frequency region. It is shown that there exists a forward scattered field which decays slowly along the source-receiver path. Its rate of attenuation depends on the azimuth of the path relative to the direction of the inhomogeneity, but is independent of frequency. The Green's tensor for the above regime is derived in closed form and leads to the quantification of fields of dipolar sources in weak inhomogeneous media. The inhomogeneity at the source creates a source-induced scattering (in addition to path-scattering) having a radiation-pattern that bears the signature of the source. The availability of the analytic Green's tensor, in conjunction with the Huygens-Kirchhoff-Helmholtz formalism, opens new ways to calculate the scattered fields due to various structural inhomogeneities applicable to exploration and earthquake seismology. The theoretical results of this study point to the conclusion that the scalar wave approximation may not always be valid for the propagation of seismic waves in the earth's lithosphere.     

Stratigraphic filtering of seismic waves: The influence of mode conversion multiples onP-wave andS-wave phase and attenuation

We demonstrate how multiples, generated at the interfaces of plane parallel beds, modify the propagation characteristics of an originally coherent seismic wave. For waves propagating at an angle to the bedding plane we find that theSV andP-waves couple so that neither is a pure mode. TheSH-wave, while modified in its propagation characteristics by multiples, remains a pure mode. The coupling ofSV-multiples into the quasi-P-mode appears weaker than the coupling ofP-wave multiples into the quasi-SV mode; at least this is so for the two simple cases of (a) density fluctuations only and (b) correlatedV _p andV _s fluctuations which conserve Poisson's ratio.We also find that the coupling is sensitive to both the angle of propagation and frequency. In addition there is a cut-off angle forP-wave multiples influencing the quasi-SV mode. Propagation angles larger than the cut-off permit theP-multiples to modify the phase of the quasi-SV mode, but not its effective attenuation. No such cut-off effect is found for SV-multiples influencing the quasi-P mode, whose angle-dependent and frequency-dependent phase distortion and effective attenuation are influenced both byP-wave multiples andSV-multiples.In view of the mathematical complexity of the expressions describing the phase, and effective attenuation of modes when allowance is made forP-andS-wave multiples, we strongly advocate numerical coding of the major mathematical formulae. By so doing a systematic study can be undertaken of the frequency and offset dependence of seismic waves as a function of seismic source input and power spectral behavior of the fluctuations in density and elastic constants of beds. It is our opinion that the full mathematical expressions are too involved to permit an analytic, systematic investigation to be given of the phase and attenuation of seismic waves with any degree of sophistication or generality.     

中国大陆地壳上地幔S波品质因子三维层析成像

根据欧亚大陆及西太平洋地区58个数字地震台站2861个长周期波形记录，从中挑选出符合条件的223对经过中国大陆的双台路径记录，采用相匹配滤波和频率域维纳滤波相结合的方法，计算出双台之间的混合路径衰减系数. 对中国大陆区域进行网格划分，通过网格化反演得到中国大陆区域每个网格单元的纯路径衰减系数. 根据速度和密度模型，对每个网格单元进行深度反演，得到10～350km的S波品质因子Qβ，进而得到整个中国大陆的S波品质因子分布. 结果显示：研究区域的Qβ横向和纵向变化非常显著，能够明显区分软流圈结构，各个大构造单元的Qβ结构呈现出显著的变化. 这为解释中国大陆区域地球内部介质性质和探讨地球动力学提供了丰富的信息.