MT剖面透视新方法与托来6．0级地震

在近10余年观测研究的基础上，采用高精度MMSO2E系统，在青藏高原北缘祁连山地震危险区，建立了国内外第一个初具规模的MT透视剖面，研究区不仅电磁干扰背景较弱，而且地震活动性较强，自1986年以来，我们已经获得了11次5－7级地震前后的宝贵资料，在研究这些资料的基础上，提出了一些有意义的认识和结论，对观测资料进行了全定量化的误差分析，评估了MT响应参数和复合二元相关度和可信限，论述了MT剖面透视资料的可靠性和可信性，对介质导电性和电性横向非均匀性动态变化的原因进行了讨论，研究表明，MT剖面透社为大陆坡坏地性地震的中短期预测及震源物理发展过程提供了定量化的多参数指标体系。     

跨祁连山-海原断裂带综合观测剖面资料分析

通过对跨祁连山-海原断裂带的民乐、华藏寺和沙沟河等3个剖面2005～2009年GPS、水准及重力观测资料的分析,得到各观测量在每个时间段的变化特征,对剖面所跨主要断裂的活动特征有所了解。发现跨不同剖面的昌马-俄博断裂、毛毛山断裂、庄浪河断裂和海原断裂附近点位的位移和重力变化都较大,表明这些断层活动性较强;通过分析各剖面垂直位移与重力变化的关系,可知重力变化主要由于物质的迁移引起。GPS资料分析表明,2006～2007年水平运动基本沿断层走向,而2007～2008年水平位移偏离断层走向,以东向位移增大为主;通过进一步对剖面周围的中国地壳运动网络连续观测站坐标时间序列的去线性分析,得知2007～2008年的观测结果与汶川地震的影响有关。用2008～2009年的重力变化同重力布格异常背景相比,可知民乐剖面附近处于调整恢复阶段,而华藏寺和沙沟河剖面附近的异常仍没有完全解除。     

聊古一井数字化氡观测资料的分析与应用

对聊古一井SD－3A数字化氡观测资料的变化特征进行了分析，与该井原有的模拟观测资料进行了对比研究，并对数字化氡观测资料的应用与分析方法进行了探讨。认为尽管SD－3A自动测氡仪所测为逸出氡与模拟观测溶解氡的物理量有所不同，但两具有基本一致的变化趋势，且数字化氡具有较强的映震能力。     

赤峰井水氡异常特征与地震预报方法研究

本文使用赤峰井１９７５－１９７６年，１９８７－１９８８年的水氡观测资料，应用水文地球化学分析预报方法指南，对观测资料进行几种不同方式的处理，结果表明该井有一定的映震能力，同时研究了该井判断异常的指标。     

Argo浮标温盐剖面观测资料的质量控制技术

Argo浮标可用来监测全球大洋从海表到2000 m深层的变化,鉴于Argo浮标的剖面观测数据存在位置错误、可疑剖面、异常数据以及盐度漂移等诸多问题,必须对Argo浮标资料进行有效的质量控制.本文基于Argo观测剖面资料与法国海洋开发研究院(IFREMER)提供的可靠历史观测数据集,提出了一种Argo资料质量控制的新途径.该方法通过寻找Argo浮标不同剖面位置与其"最佳匹配"历史剖面资料对比判别的途径,可以有效地识别Argo观测误差,特别是能够将由于Argo位置环境变化引起和由Argo浮标自身漂移引起的两类Argo浮标盐度偏移现象进行有效甄别,减少了对Argo浮标盐度剖面偏移的误判,有效节约了Argo浮标质量控制时间.本文还提出基于"三倍标准差"的异常数据检测方法,并将其与传统异常数据检测法相结合进行剖面异常数据剔除,有效实现了对异常数据的剔除.基于本文提出的Argo资料质量控制方法,对中国Argo实时资料中心网站提供的全球Argo浮标剖面进行了质量控制再分析,进一步剔除和订正了其中的一些数据误差,生成了经新的质量再控制后的全球Argo浮标剖面资料集.通过将质量再控制处理前后的数据与Ishii资料进行比较发现,处理后的数据比处理前的数据误差减小,表明本文提出的方法合理有效.     

翁牛特东部地区基底折射资料的二维反演

本文介绍了采用相遇观测系统中全部走时分支曲线双向同时二维反演拟合法。对翁牛特东部地区的三条剖面基底折射资料二维反演，文中展示了各条剖面的理论走时对比图，射线追踪图，二维速度结构分布图，一维速度柱状图以及理论地震图。     

山西临汾震区地壳上地幔构造的研究

利用郑州－临汾－靖边深地震测深剖面临汾，阳城炮点所获得的太行山至靖边段的观测资料，在以往解释的基础上重新进行了对比解释。研究结果表明：临汾与其东西两侧壳幔结构与构造的差异是极其明显的，其主要特征如下：（１）对仅在临汾以西出现的强震相Ｐ２进行了解释，并在ＰＭ波之前识别出一组来自下地壳的反射波Ｐ５；（２）根据临汾以西Ｐｇ与Ｐ２波的特征，我们确认在临汾盆地下方及其西侧，中地壳的上部８－１２ｋｍ深度内存在     

华北东北部裂陷盆地与燕山隆起地壳结构

天津—北京—赤城深地震测深剖面采用密集炮点和接收点距设计构成了较为完整的观测系统,利用该剖面地震测深资料研究华北东部裂陷盆地与燕山隆起地壳结构及构造耦合.结合该区域已有研究成果,进一步分析了张渤地震带东段的地震构造环境.结果显示,华北克拉通东部裂陷盆地结晶基底构造形态和界面结构性质与新生代以来的地壳构造运动密切相关,北侧燕山隆起区高速稳定和南侧裂陷盆地低速松散截然不同的地壳结构,使宝坻－桐柏、宁河－昌黎断裂构造分区线和燕山隆起之间的中间过渡带为张渤地震带东部平原区地震的孕发提供了良好的构造环境.     

电磁辐射异常信息的分析与思考

对大丰市地震台地震电磁辐射观测仪１９９３－１１￣１１９５－１０的观测记录资料进行了整理，对出现的异常信息和地震的对应关系进行了统计分析，提出了几点思考和认识，对从事地震电磁辐射的观测人员有一定的参考意义。     

FC－1型浮充电源研制成功并投入使用

ＦＣ－１型浮充电源研制成功并投入使用ＦＣ－１型浮充电源是专为墨水记录地震仪（ＤＤ－ｌ、ＤＫ－ｌ及遥测地震记录仪）而专门设计的。该仪器的设计适应新疆地区各地震台的供电情况，保证了地震仪的正常运转，既节省了观测经费，同时又可有效地保证地震资料的连续可靠性...     

Vertical profiles of fluid velocity and suspended sediment concentration in nearshore

Based on the power function of velocity and the friction velocity,a velocity profile is obtained.By solving the Schmidt’s diffusion equation,an equilibrium suspended sediment concentration profile is further deduced.The profiles of velocity and suspended sediment concentration agree well with the field data,and the profile of suspended sediment concentration avoids the unreasonableness of the classical Rouse profile such as a zero value at the water surface.According to these profiles,an expression which is easy to use for calculating the suspended sediment transport rate is derived.     

阴山造山带及鄂尔多斯盆地北部地区上地壳的地震波属性结构及岩性推断

通过对跨越阴山造山带及鄂尔多斯块体北部的长剖面进行人工源宽角反射/折射地震探测,采集了地震记录中sg震相的走时,并进行了反演,求得了沿剖面辖区的上地壳S波速度结构模型.同时利用P波和S波速度结构求取了剖面沿线上地壳中泊松比的分布.根据所得到的地震波场属性与结构,结合该区已有的地质研究成果和钻孔信息,推断了剖面沿线不同构...     

Electrical Struture of the Crust in Manas Earthquake Area,Xinjiang,China

Magnetotelluric sounding data obtalned recently in Manas earthquake area were processed. Inthe result, curves of apparent resistivity, impedance Phase, skewness and optimum rotationangle versus period and the real magnetic induction vectors were obtained. Then the data ofall measuring points were interpreted by 2D automatic inversion. The result indicates thatalong the sounding profile the shallow crust can be divided into 5 segments and the deep crustcan be divided into 3 segments, with faults or deep-seated fault zones as the contactboundaries between them. The sedimentary cover along the profile extents down to depthabout 12 km in maximum and a low-resistivity body exists in the crust in southern section ofthe profile. The interpretation results are well consistent with geological and othergeophysical data. The Manas M7. 7 earthquake occurred near a contact zone where theelectrical structure of the crust sharply changes.     

Testing of observation operators designed to estimate profile soil moisture from surface measurements

The accurate estimation of profile soil moisture is usually difficult due to the associated costs, strong spatiotemporal variability, and nonlinear relationship between surface and profile moisture. Here, we used data sets from the Soil and Climate Analysis Network to test how reliably observation operators developed based on the cumulative distribution function matching method can predict profile soil moisture from surface measurements. We first analysed how temporal resolution (hourly, daily, and weekly) and data length (half year, 1 year, 2 years, and 4 years) affected the performance of observation operators. The results showed that temporal resolution had a negligible influence on the performance of observation operators. However, a leave‐one‐year‐out cross‐validation showed that data length affected the performance of observation operators; a 2‐year interval was identified as the most suitable duration due to its low uncertainty in prediction accuracy. The robustness of the observation operators was then tested in three primary climates (humid continental, humid subtropical, and semiarid) of the continental United States, with the exponential filter employed as an independent method. The results indicated that observation operators reliably predicted profile soil moisture for most of the tested stations and performed almost equally well as the exponential filter method. The presented results verified the feasibility of using the cumulative distribution function matching method to predict profile soil moisture from surface measurements.     

Lithology and fluid prediction from prestack seismic data using a Bayesian model with Markov process prior

We invert prestack seismic amplitude data to find rock properties of a vertical profile of the earth. In particular we focus on lithology, porosity and fluid. Our model includes vertical dependencies of the rock properties. This allows us to compute quantities valid for the full profile such as the probability that the vertical profile contains hydrocarbons and volume distributions of hydrocarbons. In a standard point wise approach, these quantities can not be assessed. We formulate the problem in a Bayesian framework, and model the vertical dependency using spatial statistics. The relation between rock properties and elastic parameters is established through a stochastic rock model, and a convolutional model links the reflectivity to the seismic. A Markov chain Monte Carlo (MCMC) algorithm is used to generate multiple realizations that honours both the seismic data and the prior beliefs and respects the additional constraints imposed by the vertical dependencies. Convergence plots are used to provide quality check of the algorithm and to compare it with a similar method. The implementation has been tested on three different data sets offshore Norway, among these one profile has well control. For all test cases the MCMC algorithm provides reliable estimates with uncertainty quantification within three hours. The inversion result is consistent with the observed well data. In the case example we show that the seismic amplitudes make a significant impact on the inversion result even if the data have a moderate well tie, and that this is due to the vertical dependency imposed on the lithology fluid classes in our model. The vertical correlation in elastic parameters mainly influences the upside potential of the volume distribution. The approach is best suited to evaluate a few selected vertical profiles since the MCMC algorithm is computer demanding.     

The geological structure background and the crustal structure in the northeastern margin of the Qinghai-Tibetan plateau

IntroductionThestudyingareaislocatedwithin32(~40(N,100(~108(E,consistingofseveralgeotectonicsubdivisions.Thisarea,locatedonthenortheasternmarginoftheuprisingregionoftheQinghai-Tibetanplateau,isaveryactiveregionofnewtectonicmovement.Moreover,afewoflargedeepfaultsdistributethere,mostofwhicharenotonlyboundaryfaultsofimportantgeotectonicdivisions,butalsoactivefaultscontrollingcontemporarystrongseismicactivities.AsshowninFigure1,manystrongearthquakeshaveoccurredinthisresearchareasuchassouthernTi…     

井约束非稳态相位校正方法

在地震勘探资料处理中,子波的零相位化有助于提高地震资料的分辨率、改善叠加剖面的质量.常规的相位校正方法是利用测井合成记录对井旁地震记录进行相位估计,然后对整条剖面进行常相位校正,该方法没有考虑地震子波相位的非稳态性（相位随时间和空间变化）.虽然通过局部相似度方法利用最大方差模准则或包络最大相似度准则可以估计出随时间和空间变化的相位属性,但是由于零相位判别准则本身具有一定的局限性,因此精度有限.针对这一问题,本文在局部地震属性和局部平面波模型下,提出了一种井约束的非稳态相位校正方法,该方法不仅考虑了子波相位的非稳态性,而且充分利用了测井合成地震记录进行相位校正精度高的优点.理论模型和实际资料处理表明,本文方法可以有效实现信号的零相位化,有利于改善叠加效果,提高资料的分辨率.     

Imaging of offset VSP data acquired in complex areas with modified reverse-time migration

A modified reverse-time migration algorithm for offset vertical seismic profiling data is proposed. This algorithm performs depth imaging of target areas in the borehole vicinity without taking into account the overburden. Originally recorded seismograms are used; reliable results can be obtained using only the velocity profile obtained along the well. The downgoing wavefield emitted from a surface source is approximated in the target area using the transmitted P-wave, recorded by the receivers deployed in the well. This is achieved through a reverse-time extrapolation of the direct transmitted P-wave into the target area after its separation in offset vertical seismic profiling seismograms generated using a finite-difference scheme for the solution of the scalar wave equation.
The proposed approach produces 'kinematically' reliable images from reflected PP- and PS-waves and, furthermore, can be applied as a salt proximity tool for salt body flank imaging based on the transmitted PS-waves. Our experiments on synthetic data demonstrate that the modified reverse-time migration provides reliable depth images based on offset vertical seismic profiling data even if only the velocity profile obtained along the borehole is used.     

Mapping the thickness of sediments in the Ljubljana Moor basin (Slovenia) using microtremors

The Ljubljana Moor basin is characterized by moderate bedrock topography and thicknesses of Quaternary lacustrine and fluvial sediments ranging from 0 to 200 m. More than 65 boreholes which reached the bedrock were drilled in the area, but their distribution in the basin is very uneven and some data from the boreholes uncertain. There are also no data on S-velocity distribution within the basin, but seismic refraction measurements pointed out a rather uniform increase of P-velocity with depth, great impedance contrast with the bedrock and relatively small lateral velocity variations. The microtremor horizontal-to-vertical spectral ratio (HVSR) method was therefore applied as a complementary tool to seismic refraction survey to map the thickness of sediments. First, microtremors were measured at the locations of boreholes which reached the bedrock and the resonance frequencies determined. The inverse power relationship between the resonance frequency and the thickness of sediments was then determined from 53 data pairs. The quality of the correlation is moderate due to possible heterogeneities in sediments and possible 3D effects in some minor areas, but the obtained parameters correspond well to the values obtained in six other European basins. Secondly, a 16 km-long discontinuous seismic refraction profile was measured across the whole basin, leaving uncovered some larger segments where active seismic measurements were not possible. Microtremors were then measured at 64 locations along the same profile, using 250 m point spacing, without leaving any gaps. The frequency–thickness relationship was used to invert resonance frequencies to depths. These were first validated using the results of the seismic refraction survey, which showed good agreement, and finally used for interpolation in the segments of missing refraction data to obtain a continuous depth profile of the bedrock. The study has shown that the microtremor method can be used as a complementary tool for mapping the thickness of unconsolidated sediments also in areas characterized by moderate bedrock topography. As the input data are always to some extent uncertain, it is important to have a sufficiently large number of borehole data to establish a frequency–thickness relationship, as well as some additional independent geophysical information for its validation.     

LONG-OFFSET TRANSIENT EM SOUNDING NORTH OF THE RHINE-RUHR COAL DISTRICT,GERMANY1

During the past two decades electromagnetic methods have been increasingly applied to problems in hydrocarbon exploration. It has been shown that EM methods achieve satisfactory results which can be supplementary to results of the seismic method. We extend the application of the Long-Offset Transient ElectroMagnetic Sounding (lotem ) to exploration problems in the coal mining industry. A survey was conducted north of the Rhine-Ruhr coal district using a grounded dipole source and receivers laid out along several profile lines giving areal coverage. The strong noise from cultural as well as natural electromagnetic sources was reduced using prestack data processing techniques resulting in reliable data with a good signal-to-noise ratio. The interpretation is self-consistent and matches the well logs and the geological information of the area quite well. The high station density permitted the mapping of the top of the Carboniferous formation to a depth of between 800 and 1000 m below surface.