利用强震记录分析场地的地震反应

利用唐山近场强震记录对如何分析场地反应进行研究，分析比较了多种方法:线性反演法、传统的谱比法以及接收函数法．通过比较研究，发现各种方法都能给出土层场地的卓越周期，但每种方法给出的场地反应不完全相同．其中，传统的谱比法与线性反演法的结果较为接近，而接收函数法的结果则同其它两种方法的结果相差较大．在线性反演法中，给出了Ｓ波的品质因子QS，QS 在0.5～32Ｈz范围内与频率有关:QS（f）＝67 f 1.1．      

电磁波层析成像在绵阳污水处理工程场地勘探中的应用

电磁波层析成像(CT)是一种新的地球物理勘探方法，不仅勘探精度高，而且还具有快速、省钱的优点，本文叙述了电磁波层析成像技术的基本原理、工作方法及污水处理工程场地的勘探实例。最终以二维图像的形式直观、清晰地重现地下介质的构造形态，为工程场地地基安全稳定性评价提供重要依据。     

Preliminary study of crust-upper mantle structure of the Tibetan Plateau by using broadband teleseismic body waveforms

As part of a joint Sino-U.S. research project to study the deep structure of the Tibetan Plateau, 11 broadband digital seismic recorders were deployed on the Plateau for one year of passive seismic recording. In this report we use teleseimic P waveforms to study the seismic velocity structure of crust and upper mantle under three stations by receiver function inversion. The receiver function is obtained by first rotating two horizontal components of seismic records into radial and tangential components and then deconvolving the vertical component from them. The receiver function depends only on the structure near the station because the source and path effects have been removed by the deconvolution. To suppress noise, receiver functions calculated from events clustered in a small range of back-azimuths and epicentral distances are stacked. Using a matrix formalism describing the propagation of elastic waves in laterally homogeneous stratified medium, a synthetic receiver function and differential receiver functions for the parameters in each layer can be calculated to establish a linearized inversion for one-dimensional velocity structure. Preliminary results of three stations, Wen-quan, Golmud and Xigatze (Coded as WNDO, TUNL and XIGA), located in central, northern and southern Plateau are given in this paper. The receiver functions of all three stations show clear P-S converted phases. The time delays of these converted phases relative to direct P arrivals are: WNDO 7.9s (for NE direction) and 8.3s (for SE direction), TUNL 8.2s, XIGA 9.0s. Such long time delays indicate the great thickness of crust under the Plateau. The differences between receiver function of these three station shows the tectonic difference between southern and north-central Plateau. The waveforms of the receiver functions for WNDO and TUNL are very simple, while the receiver function of XIGA has an additional midcrustal converted phase. The S wave velocity structures at these three stations are estimated from inversions of the receiver function. The crustal shear wave velocities at WNDO and TUNL are vertically homogeneous, with value between 3.5–3.6 km/s down to Moho. This value in the lower crust is lower than the normal value for the lower crust of continents, which is consistent with the observed strong Sn attenuation in this region. The velocity structure at XIGA shows a velocity discontinuity at depth of 20 km and high velocity value of 4.0 km/s in the midcrust between 20–30 km depth. Similar results are obtained from a DSS profile in southern Tibet. The velocity under XIGA decreases below a depth of 30 km, reaching the lowest value of 3.2 km/s between 50–55 km. depth. This may imply that the Indian crust underthrusts the low part of Tibetan crust in the southern Plateau, forming a “double crust”. The crustal thickness at each of these sites is: WNDO, 68 km; TUNL, 70 km; XI-GA, 80 km. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, Supp., 581–592, 1992.     

西昆仑接收函数反演与构造解析

中法合作布设的青藏高原西昆仑叶城-噶尔(狮泉河)天然地震探测剖面穿越了西昆仑构造的主要构造单元, 通过对36个台站128个远震P波事件的接收函数分析,取得了各台站的接收函数随方位角和震中距变化的趋势, 分析了变化趋势与地壳构造的关系。通过台站平均接收函数的反演拟合,取得了120km深度范围以上的横波速度分布,综合地质构造解释,发现康西瓦断裂地表倾斜渐变为直立并略向南拖曳的现象,塔里木前陆盆地沉积与逆冲构造有比较明显的显示,同时推测了盆地下的滑脱构造。班公湖-怒江断裂也有向北俯冲的显示,得出了印藏碰撞远距离效应在西昆仑受阻,应力通过地表次级断裂释放,同时形成盆地沉积的前陆构造模式。     

Single GPS receiver time-relative positioning with loop misclosure corrections

Time-relative positioning makes use of observations taken at two different epochs and stations with a single global positioning system (GPS) receiver to determine the position of the unknown station with respect to the known station. The limitation of this method is the degradation over time of the positioning accuracy due to the temporal variation of GPS errors (ionospheric delay, satellite clock corrections, satellite ephemerides, and tropospheric delay). The impact of these errors is significantly reduced by adding to the one-way move from the known to the unknown station, a back move to the known station. A loop misclosure is computed from the coordinates obtained at the known station at the beginning and at the end of the loop, and is used to correct the coordinates of the unknown station. The field tests, presented in this paper, show that using the loop misclosure corrections, time-relative positioning accuracy can be improved by about 60% when using single frequency data, and by about 40% with dual frequency data. For a 4-min processing interval (an 8-min loop) and a 95% probability level, errors remain under 20 cm for the horizontal components and 36 cm for the vertical component with single frequency data; and under 11 cm for the horizontal components and 29 cm for the vertical component with dual frequency data.     

Imaging Offsets in the Moho: Synthetic Tests using Gaussian Beams with Teleseismic Waves

We carry out a sequence of numerical tests to understand conditions under which rapid changes in crustal thickness can be reliably imaged by teleseismic body waves. Using the finite-difference method over a 2-D grid, we compute synthetic seismograms resulting from a planar P-wavefield incident below the grid. We then image the Moho using a migration scheme based on the Gaussian beam representation of the wavefield. The use of Gaussian beams for the downward propagation of the wavefield is particularly advantageous in certain geologically critical cases such as overthrusting of continental lithosphere, resulting in the juxtaposition of high-velocity mantle material over crustal rocks. In contrast to ray-based methods, Gaussian beam migration requires no special treatment to handle such heterogeneities. Our results suggest that with adequate station spacing and signal-to-noise ratios, offsets of the Moho, on the order of 10 km in height, can be reliably imaged beneath thickened crust at depths of about 50 km. Furthermore, even sharp corners and edges are faithfully imaged when precise values of seismic wave speeds are available. Our tests also demonstrate that flexibility in choices of different types of seismic phases is important, because any single phase has trade-offs in issues such as spatial resolution, array aperture, and amplitude of signals.     

Thickness estimates of the upper-mantle transition zone from bootstrapped velocity spectrum stacks of receiver functions

We modify the receiver-functions stacking technique known as velocity spectrum stacking (VSS) so as to estimate combinations of velocity model ( V_P and V_S ) and depth that stack the Ps conversion from upper-mantle discontinuities most coherently. We find that by estimating the differences in the depths to the 660 and 410 km discontinuities using velocities that maximize the stacked amplitudes of P410s and P660s phases we can estimate the thickness of the transition zone more accurately than the depths to either of these discontinuities. We present two examples indicating that the transition zone beneath Obninsk, Russia, is 252±6 km thick and that beneath Pasadena, California, is only 220±6 km thick.     

略论海洋测高卫星GPS厘米级定轨的实施基础

基于对多种飞机的机载GPS测量实践,认为我国HY-2海洋二号卫星实现厘米级精度的星载GPS定轨测量的基础是:选择适合天线,捕获多颗在视GPS卫星;注重天线安装位置,减弱多路径效应影响;选择适合GPS信号接收机,确保星载GPS测量数据优质。     

基于波束控制的卫星导航接收机抗干扰算法研究

卫星导航接收机抗干扰技术是接收机设计的一项关键技术。在研究基于功率反演的空域自适应抗干扰算法基础上，提出了采用基于波束控制的抗干扰算法，并进行了仿真计算，结果表明：该算法明显提高了接收机抗干扰性能。