Use of microtremor measurement for assessing site effects in Northern Belgium – interpretation of the observed intensity during the MS = 5.0 June 11 1938 earthquake

To provide quantitative information on site effects in the northern part of Belgium, forty-seven H/V microtremor measurements were performed with 5 second seismometers over an area of about 15.000 km². Most of the results show a northward regular increase of the fundamental period in agreement with the augmentation of the Mesozoic and Cainozoic soft sedimentsthickness from a few meters 40 km south of Brussels to 900 m at the Netherlands-Belgium border. The measured resonance frequency values wereconsistent with theoretical computations performed at different sites onthe basis of existing information and shallow seismic experiments. At oneparticular site (Uccle) where borehole data were available, microtremor measurements using an array of four seismological stations with differentapertures allowed to obtain the low frequency part of the Rayleigh wave dispersion curve, extending the range covered by the analysis of surfacewaves artificially generated. The Vs profile derived from the surface waveinversion corroborates the 1 Hz natural frequency of the site. Comparison of these results with the macroseismic information concerning the M_S = 5.0 1938 earthquake which occurred 50 km west of Brussels,confirmed the hypothesis that the geological structure of the Brabant massifis likely to control damage distribution during such an earthquake. Comparisonbetween the intensity map of the 1938 earthquake and the resonance period ofsediments obtained by our microtremor study shows a clear relation betweenthe two parameters. During the 1938 earthquake, site effects played a prominent role due to the dimension of the source whose corner frequency wasabout 1 Hz.     

Determination of S-wave velocity structure using microtremors and spac method applied in Thessaloniki (Greece)

Array measurements of microtremors at 16 sites in the city of Thessaloniki were performed to estimate the V_s velocity of soil formations for site effect analysis. The spatial autocorrelation method was used to determine phase velocity dispersion curves in the frequency range from 0.8–1.5 to 6–7 Hz. A Rayleigh wave inversion technique (stochastic method) was subsequently applied to determine the Vs profiles at all the examined sites. The determination of V_s profiles reached a depth of 320 m. Comparisons with V_s values from cross-hole tests at the same sites proved the reliability of the SPAC method. The accuracy of the V_s profiles, the ability to reach large penetration depths in densely populated urban areas and its low cost compared to conventional geophysical prospecting, make Mictrotremor Exploration Method very attractive and useful for microzonation and site effects studies. An example of its application for the site characterization in Thessaloniki is presented herein.     

Tortoises and hares: dissolution,erosion and isostasy in landscape evolution

Denudation mechanisms differ fundamentally between limestone and silicate rock types, which are subject to very different rate thresholds and enhancers/inhibitors. Silicates are removed largely by erosion, the mechanical entrainment and transport of particles. This is a relatively high energy, and highly episodic, process which occurs only when a minimum threshold ?ow velocity is exceeded; it is inhibited by vegetation cover and favoured by strongly seasonal runoff. Limestone is removed largely by chemical dissolution at a rate directly proportional to runoff. Dissolution is a relatively low energy process that can occur at any ?ow velocity or in static water; in general it is enhanced by vegetation cover and non‐seasonality of runoff. These contrasting factors in the denudation of silicates versus limestone can produce strikingly uneven rates of surface lowering across a landscape, sometimes akin to the well known ‘tortoise and hare race’, where the slow and steady denudation of limestones may in the long term exceed the sometimes rapid, but often localized and episodic, erosion of silicates. Prolonged exposure of limestone to a humid temperate climate in a tectonically stable environment produces low‐relief corrosion plains in which limestone uplands are anomalous and, in most instances, due to recent unroo?ng from beneath a siliciclastic cover. In a highly seasonal or semi‐arid climate almost the exact inverse may develop, with ‘?ashy’ runoff and sparse vegetation favouring erosion rather than dissolution. Even under a constant humid climate progressive unroo?ng of a thick limestone unit within folded siliciclastics may lead to a topographic inversion over time, with the limestone outcrop always forming a topographic low ?anked by siliciclastic uplands. Valleys will be initiated on anticlinal crests, where the limestone is ?rst unroofed, but progressive lowering of the limestone causes these valleys to migrate to their ?nal position in the synclinal troughs. In humid climates isostatic compensation in response to slow, but continuous, denudation of extensive limestone outcrops may be a signi?cant factor in the development of relief on adjacent, more slowly eroding, silicate outcrops. Copyright © 2004 John Wiley & Sons, Ltd.     

鄂尔多斯盆地北部黄土沙漠区隐蔽薄气层预测

长庆油田公司针对黄土沙漠区低孔、低渗、薄气藏的地震目标处理,经过多年的探索与实践,已形成了一套适合该类气藏特点的储层横向预测技术。在榆林、乌审旗、苏里格庙和神木等地区的上古生界天然气勘探、开发阶段取得了显著效果,使长庆油田上古生界天然气储量连续八年保持稳定增长。本文详细剖析了该技术的应用效果及经验。     

Earthquake relocation and 3-dimensional crustal structure of P-wave velocity in cen-tral-western China

IntroductionUnderstandingthemechanismofcontinentalearthquakesisveryimportantforseismichaz-ardpreventionandearthquakeprediction.Themodernseismotectonictheoryandtheideaofearthquakepredictionaredevelopedmainlyfromthestudiesoninterplateearthquakes,whicharedifficulttoexplainthephenomenaofintraplateearthquakes,suchasthecontinentalearthquakesoccurredinChinesemainland.Whiletheinterplateearthquakesoccurredalongtheplatebounda-ries,theintraplateearthquakesdistributediffuselyintheinterioroftheplates.Thus…     

黄金带油田非均质储层特征及有利目标预测

为了加深对黄金带复杂断块油田进一步开发，针对该区断层复杂、储层变化大、油水关系复杂等特征，笔者综合运用岩心、分析化验及录井、测井资料，对其沉积特征、储层岩性、物性和非均质性进行详细研究。以JASON软件为平台，对含油砂体进行了追踪和预测，并在有利区提出井位部署，为该区的进一步开发提供了地质依据。     

非饱和土中镉离子传输模型参数反演

使用自行设计的的真空实验装置,采用一维实验室土柱实验方法在不同实验条件下进行了15组实验,并基于实验所得的穿透曲线,使用梯度正则化方法反演得到了镉离子在8种葡萄牙土样和1种比利时土样中相应的传输模型参数。同时通过数值模拟实验验证了梯度正则化方法的有效性和可靠性,进而证明所得结果是有效的。     

Detection of Methane Hydrate-Bearing Zones from Seismic Data

Abstract. The MITI Nankai Trough wells were drilled offshore Japan in the Tokai area in 1999 and 2000. The occurrence of methane hydrate was confirmed by various indicators in the borehole logs and from core data. These findings have a large impact on potential future Japanese energy resources and other related-scientific interests.
We first tried to find the methane hydrate-bearing zones using interval velocities derived from NMO velocity analysis. However, this analysis produced poor resolution. To achieve a more detailed delineation of the gas hydrate- and gas-bearing zones, we executed a seismic impedance inversion calibrated by the logs from two of the MITI Nankai Trough wells. Although these two wells are only about 90 m apart, we were able to produce an impedance section with fine detail by adopting a simple initial model and incorporating physical properties of the methane hydrate-bearing zones. The locations of the methane hydrate-bearing zones are readily apparent in the final section.     

Imaging the ramp–décollement geometry of the Chelungpu fault using coseismic GPS displacements from the 1999 Chi-Chi, Taiwan earthquake

We use coseismic GPS data from the 1999 Chi-Chi, Taiwan earthquake to estimate the subsurface shape of the Chelungpu fault that ruptured during the earthquake. Studies prior to the earthquake suggest a ramp–décollement geometry for the Chelungpu fault, yet many finite source inversions using GPS and seismic data assume slip occurred on the down-dip extension of the Chelungpu ramp, rather than on a sub-horizontal décollement. We test whether slip occurred on the décollement or the down-dip extension of the ramp using well-established methods of inverting GPS data for geometry and slip on faults represented as elastic dislocations. We find that a significant portion of the coseismic slip did indeed occur on a sub-horizontal décollement located at 8 km depth. The slip on the décollement contributes 21% of the total modeled moment release. We estimate the fault geometry assuming several different models for the distribution of elastic properties in the earth: homogeneous, layered, and layered with lateral material contrast across the fault. It is shown, however, that heterogeneity has little influence on our estimated fault geometry. We also investigate several competing interpretations of deformation within the E/W trending rupture zone at the northern end of the 1999 ground ruptures. We demonstrate that the GPS data require a 22- to 35-km-long lateral ramp at the northern end, contradicting other investigations that propose deformation is concentrated within 10 km of the Chelungpu fault. Lastly, we propose a simple tectonic model for the development of the lateral ramp.     

Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.