密度和压缩系数的散射层析成像法

本文在速度成像的基础上研究了同时对密度和压缩系数成像的散射波层析成像法.对不同散射角度的计算可以得到一系列反演图像,拟合这些图像,从而可以有效地达到对密度和压缩系数（或速度）成像的目的.与单纯的速度成像相比,增加了反演的难度.首先是对资料的方位性要求增加;其次是对资料的利用率下降.即便如此,从对较少量的炮点和检波点资料的数值计算来看,仍取得了满意的成像结果.我们对组成字母“A”的散射体结构进行了成像计算,结果能够同时再现密度和压缩系数,成像清晰,表明了方法的可行性,并能应用于复杂结构的成像问题.     

青藏高原及邻近区域的S波三维速度结构

本文收集了WWSSN台网和我国台网中13个地震台站的长周期地震记录,用140条10-90s瑞利波频散曲线和作者提出的Tarantola-Backus面波频散层析成象方法,作了青藏高原及邻区的速度反演,得出该地区岩石层速度结构的三维图象.结果表明,1.在10-110km深度范围内,速度结构出现与大地构造特征相关的分区性,显示出四个构造单元:青藏块体、柴达木-巴颜喀拉-三江块体、塔里木块体和印度块体.2.高原内部,深度为10-70km内速度较低,莫霍界面呈不对称盆形分布,藏北那曲附近地壳厚度超过70km,高原边缘壳厚为45-50km,90-110km为高速异常,表明高原内部存在上地幔盖层.3.高原北部的班公湖断裂和东部的三江断裂系是该区重要的分界线,是岩石层结构存在明显差异的重要接触部位,可能是冈瓦纳古陆与欧亚古陆的缝合带.4.柴达木-巴颜喀拉-三江块体内部速度分布不均匀,地壳厚度由北向南从45km加深到60km;在深度90-110km存在一低速层.5.塔里木地块内速度随深度均匀增加,从地壳到上地幔110km内没有发现低速层.地壳厚度约50km.     

Acoustic plane-wave scattering from a rough surface over a random fluid medium

Acoustic plane-wave scattering from a rough surface overlying a fluid half-space with a sound-speed distribution subject to a small random variation is considered. Under the assumption that the surface roughness and medium randomness are statistically independent, the scattered field may be derived by first solving for the mean field in the medium, and then incorporating with boundary-perturbation method to obtain the total mean field and the power spectral density of the scattered field. The employed algorithm is compatible to the analysis available in the existing literature so that the formulations are conveniently integrated. The results for the power spectral density have shown that the effects of medium inhomogeneities on the rough surface scattering are limited in a spectral regime where the scattered components have shallow grazing angles. The distribution of the power spectral density over the space is primarily governed by the power spectrum and correlation lengths of the rough surface.     

GPS observations of the ionospheric F2-layer behavior during the 20th November 2003 geomagnetic storm over South Korea

The ionospheric F2-layer peak density (NmF2) and its height (hmF2) are of great influence on the shape of the ionospheric electron density profile Ne (h) and may be indicative of other physical processes within the ionosphere, especially those due to geomagnetic storms. Such parameters are often estimated using models such as the semiempirical international reference ionosphere (IRI) models or are measured using moderately priced to expensive instrumentation, such as ionosondes or incoherent scatter radars. Global positioning system (GPS) observations have become a powerful tool for mapping high-resolution ionospheric structures, which can be used to study the ionospheric response to geomagnetic storms. In this paper, we describe how 3-D ionospheric electron density profiles were produced from data of the dense permanent Korean GPS network using the tomography reconstruction technique. These profiles are verified by independent ionosonde data. The responses of GPS-derived parameters at the ionospheric F2-layer to the 20th November 2003 geomagnetic storm over South Korea are investigated. A fairly large increase in the electron density at the F2-layer peak (the NmF2) (positive storm) has been observed during this storm, which is accompanied by a significant uplift in the height of the F2 layer peak (the hmF2). This is confirmed by independent ionosonde observations. We suggest that the F2-layer peak height uplift and NmF2 increase are mainly associated with a strong eastward electric field, and are not associated with the increase of the O/N₂ ratio obtained from the GUVI instruments aboard the TIMED satellite. It is also inferred that the increase in NmF2 is not caused by the changes in neutral composition, but is related to other nonchemical effects, such as dynamical changes of vertical ion motions induced by winds and E × B drifts, tides and waves in the mesosphere/lower thermosphere region, which can be dynamically coupled upward to generate ionospheric perturbations and oscillations.     

Unified approach to amplitude attenuation and coda excitation in the randomly inhomogeneous lithosphere

A unified model is proposed for explaining the frequency dependent amplitude attenuation and the coda wave excitation on the basis of the single scattering process in the randomly inhomogeneous lithosphere. Adopting Birch's law and a direct proportion between density and wave velocity, we statistically describe the inhomogeneous medium by one random function characterized by the von Karman autocorrelation function. We calculate the amplitude attenuation from the solid angle integral of scattered wave energy on the basis of the Born approxiimation after subtracting the travel-time fluctuation effect caused by slowly varying velocity inhomogeneities. This subtraction is equivalent to neglect energy loss by scattering within a cone around the forward direction. The random inhomogeneity of the von Karman autocorrelation function of order 0.35 with the mean square fractional fluctuation of 7.2×10^–3 1.3×10^–2 and the correlation distance of 2.15.1 km well explains observed backward scattering coefficientg and the ratioQ _P ^–1 /Q _S ^–1 , and observed and partially conjecturedQ _S ^–1 for frequencies between 0.5 Hz and 30 Hz.     

电阻率图像中的地下水和断层

目的：将电阻率层析成像应用于探测潜伏断层的研究中，本文发现了断层和地下水的一些基本电阻率分布特征，这对于工程物探意义重大，一般情况下，断层两侧具有不同的电阻率特征，但是，根据电阻率层析图像中的电阻率分布，通常难以区分断层和地下水层，这是因为两者不仅都具有低电阻率值，而且还具有非常相似的电阻率异常特征。资料和方法：运用电阻率层析图像的数据，电阻率层析图像中的断层会呈现如下特征：1）由于孔隙度的加大和地下水的存在，使得断层表现出高角度的低阻线性结构。它们既可以出现在浅部盖层中，也可以存在于深部基岩中，特别是在深部区域，它们尤为明显；2）它们还呈现出高角度的线性梯度带，在该梯度带两边的电阻率结构出现整体性的差异，通常情况下，正断层的上盘表现出低阻或／和班驳状的高阻和低阻扰动区，而下盘则为完整的高阻区，这与逆冲断层正好相反；3）与断层有关的电阻率异常区常常具有良好的大尺度水平连续性，并且可以追瞎异常区附近的精细电性结构。而地下水的电阻率特征为：1）如果没有裂隙，地表水所引起的低阻区非常浅，即使存在丰富的水源以及高孔隙度的砾岩和中粗砂。一般情况下，其底端深度不超过强风化区；2）地下水的电阻率值非常低，特别在高矿化度的地区。地下水，包括岩溶水和砂岩水，的电阻率总显示出局部水平延伸或／和面团状特征；3）地下水层的深度朝某个固定方向逐渐增加，并且其电阻率图像会随季节而变；4）一般情况下，在水下渗的地区，会出现降水漏斗，其上部为高阻，而下部为低阻，从而便形成了“Y”或“V”字型的典型结构。结果：利用上述的基本特征一般可以区分断层和地下水。结论：仅依靠电阻率层析图像，可能极难准确地区分断层和裂隙水，这是因为裂隙水不但可能具有高角度的低阻线性结构，而且在一定尺度上具有很好的水平连续性，还有，由于电阻率层析成像较差的垂直分辨率，难以精确确定断层的上端点位置，所以结合其它的物探手段如钻探和浅层地震勘探是非常必要的。     

Probing mid-mantle heterogeneity using PKP coda waves

We investigate the utility of PKP coda waves for studying weak scattering from small-scale heterogeneity in the mid-mantle. Coda waves are potentially a useful probe of heterogeneity in the mid-mantle because they are not preferentially scattered near the CMB, as PKP precursors are, but are sensitive to scattering at all depths. PKP coda waves have not been used for this purpose historically because of interference with other late-arriving energy due to near-surface resonance and scattering. Any study of deep mantle scattering using coda waves requires the removal of near-surface effects from the data. We have analyzed 3624 recordings of PKP precursors and coda made by stations in the Incorporated Research Institutions for Seismology (IRIS) Global Seismographic Network (GSN). To study the range and time dependence of the scattered waves, we binned and stacked envelopes of the recordings. We have considered precursors that arrive within a 20 s window before PKP and coda waves in a 60 s window after PKP. The PKP scattered waves increase in amplitude rapidly with range as predicted by scattering theory. At ranges below 125°, we predict and observe essentially no scattered energy preceding PKP. Coda amplitudes at these ranges are independent of range and provide an estimate of energy due to near-surface effects that we can expect at all ranges. We use the average coda amplitude at ranges from 120 to 125° to correct coda amplitudes at other ranges. PKP coda waves show a strong dependence on time and range and are clearly influenced by scattering in the lower mantle. PKP coda waves, however, do not provide a tighter constraint on the vertical distribution of mantle heterogeneity than is provided by precursors. This is due, in part, to relatively large scatter in coda amplitudes as revealed by a resampling analysis. Modeling using Rayleigh–Born scattering theory and an exponential autocorrelation function shows that PKP coda amplitudes are not highly sensitive to the vertical distribution of heterogeneity in the mantle. To illustrate this we consider single-scattering in two extreme models of mantle heterogeneity. One allows heterogeneity just at the CMB; the other includes heterogeneity throughout the mantle. The amplitudes of precursors are tightly constrained by our stack and support our earlier conclusion that small-scale heterogeneity is uniformly distributed throughout the lower mantle. The best-fit model includes 8 km scale length heterogeneity with an rms velocity contrast throughout the mantle of 1%.     

Crustal structure of the continental margin of Korea in the East Sea (Japan Sea) from deep seismic sounding data: evidence for rifting affected by the hotter than normal mantle

Despite the various opening models of the southwestern part of the East Sea (Japan Sea) between the Korean Peninsula and the Japan Arc, the continental margin of the Korean Peninsula remains unknown in crustal structure. As a result, continental rifting and subsequent seafloor spreading processes to explain the opening of the East Sea have not been adequately addressed. We investigated crustal and sedimentary velocity structures across the Korean margin into the adjacent Ulleung Basin from multichannel seismic (MCS) reflection and ocean bottom seismometer (OBS) data. The Ulleung Basin shows crustal velocity structure typical of oceanic although its crustal thickness of about 10 km is greater than normal. The continental margin documents rapid transition from continental to oceanic crust, exhibiting a remarkable decrease in crustal thickness accompanied by shallowing of Moho over a distance of about 50 km. The crustal model of the margin is characterized by a high-velocity (up to 7.4 km/s) lower crustal (HVLC) layer that is thicker than 10 km under the slope base and pinches out seawards. The HVLC layer is interpreted as magmatic underplating emplaced during continental rifting in response to high upper mantle temperature. The acoustic basement of the slope base shows an igneous stratigraphy developed by massive volcanic eruption. These features suggest that the evolution of the Korean margin can be explained by the processes occurring at volcanic rifted margins. Global earthquake tomography supports our interpretation by defining the abnormally hot upper mantle across the Korean margin and in the Ulleung Basin.     

用TRMM资料研究江淮、华南降水的微波特性

热带测雨卫星TRMM (Tropical Rainfall Measuring Mission) 于1997年11月发射成功, 其首次携带了空载雷达, 有关资料已在网上对公众发布。利用热带测雨卫星上的微波成像仪TMI (TRMM Microw ave Imager) 资料以及其和测雨雷达TRMM/PR (Precipitation Radar) 资料联合反演的地面瞬时降水产品, 采用散射指数 (Is) 法从理论上探讨了我国江淮、华南降水尤其是暴雨的微波特性, 其中I_s表达式通过江淮、华南晴空TMI资料统计回归得到。以联合反演的地面瞬时降水产品为真值, 用面积相当法对14个降水个例求I_s降水阈值, 研究了阈值和降水面积以及85.5 GHz垂直通道最低亮温的关系, 并寻求了I_s和降水的相关特征。研究表明:I_s降水阈值随降水面积的增大或85.5 GHz垂直通道最低亮温的降低有增高的趋势; I_s与强对流性降水瞬时雨强对应很好, I_s≥60 K是一个好的暴雨指标。最后进行了初步的雨强反演试验研究, 由于TMI资料分辨率的提高以及时空配合较好的真值, 反演的地面瞬时降水与真值相关效果大大提高。     

大别造山带上部地壳结构的有限差分层析成像

将有限差分层析反演方法用于大别造山带人工地震测深剖面上部地壳初至波的走时层析成像，得到这一剖面上部地壳横向不均匀结构图像表明，在南大别地区的超高压变质带发现榴辉岩相岩石，其下方3km深度以内的基底仍具有大陆地壳构造中正常的结晶基底速度（6．00km／s左右）；北大别地区在整个上地壳保持正常的速度，而在超高压带下方3km深度以下为6．20-630km／s的相对高速异常区，这一现象可能与超高压变质岩的含量增大有关．同时还表明，南大别和北大别之间至少在上地壳已有明显差异，它们之间的水吼一五河断裂可能是大别造山带内部的主要构造分界线．