北黄海重要历史地震参数复核

以文献记载资料和大连地震台的观测资料为主要依据,对1916年、1917年和1944年发生在北黄海鸭绿江口附近的3次重要地震参数进行复核确认,复核后的震中位于现代小地震活动相对密集的地段,且与大连台的观测结果更吻合.本文结果为合理识别发震构造和修订地震目录准备了条件.     

Inverse methods in hydrogeology: Evolution and recent trends

Parameter identification is an essential step in constructing a groundwater model. The process of recognizing model parameter values by conditioning on observed data of the state variable is referred to as the inverse problem. A series of inverse methods has been proposed to solve the inverse problem, ranging from trial-and-error manual calibration to the current complex automatic data assimilation algorithms. This paper does not attempt to be another overview paper on inverse models, but rather to analyze and track the evolution of the inverse methods over the last decades, mostly within the realm of hydrogeology, revealing their transformation, motivation and recent trends. Issues confronted by the inverse problem, such as dealing with multiGaussianity and whether or not to preserve the prior statistics are discussed.     

Misfit functionals in Laplace‐Fourier domain waveform inversion,with application to wide‐angle ocean bottom seismograph data

In seismic waveform inversion, non‐linearity and non‐uniqueness require appropriate strategies. We formulate four types of L2 normed misfit functionals for Laplace‐Fourier domain waveform inversion: i) subtraction of complex‐valued observed data from complex‐valued predicted data (the ‘conventional phase‐amplitude’ residual), ii) a ‘conventional phase‐only’ residual in which amplitude variations are normalized, iii) a ‘logarithmic phase‐amplitude’ residual and finally iv) a ‘logarithmic phase‐only’ residual in which the only imaginary part of the logarithmic residual is used. We evaluate these misfit functionals by using a wide‐angle field Ocean Bottom Seismograph (OBS) data set with a maximum offset of 55 km. The conventional phase‐amplitude approach is restricted in illumination and delineates only shallow velocity structures. In contrast, the other three misfit functionals retrieve detailed velocity structures with clear lithological boundaries down to the deeper part of the model. We also test the performance of additional phase‐amplitude inversions starting from the logarithmic phase‐only inversion result. The resulting velocity updates are prominent only in the high‐wavenumber components, sharpening the lithological boundaries. We argue that the discrepancies in the behaviours of the misfit functionals are primarily caused by the sensitivities of the model gradient to strong amplitude variations in the data. As the observed data amplitudes are dominated by the near‐offset traces, the conventional phase‐amplitude inversion primarily updates the shallow structures as a result. In contrast, the other three misfit functionals eliminate the strong dependence on amplitude variation naturally and enhance the depth of illumination. We further suggest that the phase‐only inversions are sufficient to obtain robust and reliable velocity structures and the amplitude information is of secondary importance in constraining subsurface velocity models.     

Transient process and optimal design of receiver coil for small‐loop transient electromagnetics

Distributed parameters of the receiver coils greatly affect transient electromagnetic signals over short time periods, causing a delay in the signal's effective sampling time and the loss of shallow exploration information. This paper investigates the influence of transient process on apparent resistivity calculation and analyses the relations between the error of apparent resistivity and receiver coil design. We find that, under the same effective area, different radii of the receiver coils lead to different levels of impact on the estimation of the apparent resistivity. An optimization model is proposed to determine the optimal receiver coil size that gives rise to the smallest estimation error of the apparent resistivity. The relationship between the optimal radius and the effective areas is developed, which serves as a guideline for the optimal receiver coil design. The results may provide a useful means for improving the accuracy of the small loop transient electromagnetic instrumentation for shallow‐depth mapping.     

随机天气模型参数化方案的研究及其模拟能力评估

文中介绍了随机天气模型 WGEN的基本结构及其模拟原理 ,并针对其中随机过程的统计结构特征和 GCMs输出要素的不同时空尺度特点 ,利用动态数据的参数化分析方法等统计学技术 ,确定了该模型参数的估计方法。同时基于蒙特卡罗数值计算原理 ,给出了 WGEN的随机试验方法 ,并通过模拟基准气候 ,从时间分布和空间场两方面对模型在中国东北地区的模拟效果及其能力进行了评估。结果表明 ,模型对于最高气温、最低气温、降水和辐射等要素均具有较好的模拟效果 ,模拟序列与观测序列的取值分布有较一致的概率特性。由此可以结合 GCMs大尺度网格上输出的月和年要素值 ,通过调控随机过程的参数 ,生成具有不同气候变率的 2× CO2 逐日气候变化情景 ,实现气候预测模式与气候影响模式的嵌套 ,进一步研究气候变率变化的可能影响。     

基于构造的地震活动参数扫描短期预测方法研究

应用山东省和江苏省地震局“九五”地质研究成果, 在华东地区共选择49条断裂和7个构造块体作为地震活动性研究的构造基础, 选用频度、 ΔT-T、 小震调制比3个参数, 对构造两侧各20 km宽度范围的地震资料进行参数时间扫描。 通过普查, 得到的构造优势映震距离为300 km; 得到了10个敏感构造及其映震特点。 通过对华东地区13个中强地震震例的构造异常的综合分析, 得到如下认识: ① 震前存在场的群体异常特征, 表明华东中强地震“前兆”是在构造应力作用下表现出的、 区域性的整体过程; ② 异常的构造常常分布在震中周围; ③ 地震常常发生在平静异常的构造附近; ④ 高频度异常常常发生在震中外围, 并常具有步进性; ⑤ 震中周围的构造在震前2～3年开始先后不断出现高频和平静异常, 平静后出现多处小震调制比高值是孕震中期向短期过渡的时间转折标志。     

Reflection-constrained 2D and 3D non-hyperbolic moveout analysis using particle swarm optimization

The objective of moveout parameter inversion is to derive sets of parameter models that can be used for moveout correction and stacking at each common midpoint location to increase the signal-to-noise ratio of the data and to provide insights into the kinematic characteristics of the data amongst other things. In this paper, we introduce a data-driven user-constrained optimization scheme that utilizes manual picks at a point on each reflector within a common midpoint gather to constrain the search space in which an optimization procedure can search for the optimal parameter sets at each reflection. The picks are used to create boundary curves which can be derived approximately via an optimization technique or analytically via the derivation of an analytical bounds function. In this paper, we derive analytical forms of bounds functions for four different moveout cases. These are normal moveout, non-hyperbolic moveout, azimuthally dependent normal moveout and azimuthally dependent non-hyperbolic moveout. The optimization procedure utilized here to search for the optimal moveout parameters is the particle swarm optimization technique. However, any metaheuristic optimization procedure could be modified to account for the constraints introduced in this paper. The technique is tested on two-layer synthetic models based on three of the four moveout cases discussed in this paper. It is also applied to an elastic forward modelled synthetic model called the HESS model, and finally to real 2D land data from Alaska. The resultant stacks show a marked improvement in the signal-to-noise ratio compared to the raw stacks. The results for the normal moveout, non-hyperbolic moveout and azimuthally dependent normal moveout tests suggest that the method is viable for said models. Results demonstrate that our method offers potential as an alternative to conventional parameter picking and inversion schemes, particularly for some cases where the number of parameters in the moveout approximation is 2 or greater.     

Multiwell study of seismic attenuation at the CO2CRC Otway project geosequestration site: Comparison of amplitude decay,centroid frequency shift and 1D waveform inversion methods

At the CO2CRC Otway geosequestration site, the abundance of borehole seismic and logging data provides a unique opportunity to compare techniques of Q (measure of attenuation) estimation and validate their reliability. Specifically, we test conventional time-domain amplitude decay and spectral-domain centroid frequency shift methods versus the 1D waveform inversion constrained by well logs on a set of zero-offset vertical seismic profiles. The amplitude decay and centroid frequency shift methods of Q estimation assume that a seismic pulse propagates in a homogeneous medium and ignore the interference of the propagating wave with short-period multiples. The waveform inversion explicitly models multiple scattering and interference on a stack of thin layers using high-resolution data from sonic and density logs. This allows for stable Q estimation in small depth windows (in this study, 150 m), and separation of the frequency-dependent layer-induced scattering from intrinsic absorption. Besides, the inversion takes into account band-limited nature of seismic data, and thus, it is less dependent on the operating frequency bandwidth than on the other methods. However, all considered methods of Q estimation are unreliable in the intervals where subsurface significantly deviates from 1D geometry. At the Otway site, the attenuation estimates are distorted by sub-vertical faults close to the boreholes. Analysis of repeated vertical seismic profiles reveals that 15 kt injection of the CO₂-rich fluid into a thin saline aquifer at 1.5 km depth does not induce detectable absorption of P-waves at generated frequencies 5–150 Hz, most likely because the CO₂ plume in the monitoring well is thin, <15 m. At the Otway research site, strong attenuation Q ≈ 30–50 is observed only in shaly formations (Skull Creek Mudstone, Belfast Mudstone). Layer-induced scattering attenuation is negligible except for a few intervals, namely 500–650 m from the surface, and near the injection interval, at around 1400–1550 m, where Q_scat ≈ 50–65.     

On possible fluid detection in equivalent transversely isotropic media

We consider a transversely isotropic medium that is long-wave equivalent to a stack of thin, parallel, isotropic layers and is obtained using the Backus average. In such media, we analyse the relations among anisotropy parameters; Thomsen parameters, ε and δ, and a new parameter ϕ. We discuss the last parameter and show its essential properties; it equals 0 in the case of isotropy of equivalent medium and/or constant Lamé coefficient λ in layers. The second property occurs to make ϕ sensitive to variations of λ in thin-bedded sequences. According to Gassmann, in isotropic media the variation of fluid content affects only the Lamé coefficient λ, not μ; thus, the sensitivity to changes of λ is an essential property in the context of possible detection of fluids. We show algebraically and numerically that ϕ is more sensitive to these variations than ε or δ. Nevertheless, each of these parameters is dependent on the changes of μ; to understand this influence, we exhibit comprehensive tables that illustrate the behaviour of anisotropy parameters with respect to specific variations of λ and μ. The changes of μ in layers can be presented by the Thomsen parameter γ that depends on them solely. Hence, knowing the values of elasticity coefficients of equivalent transversely isotropic medium, we may compute ϕ and γ, and based on the aforementioned tables, we predict the expected variation of λ; in this way, we propose a new method of possible fluid detection. Also, we show that the prior approach of possible detection of fluids, proposed by Berryman et al., may be unreliable in specific cases. To establish our results, we use the Monte Carlo method; for the range and chosen variations of Lamé coefficients λ and μ – relevant to sandstones – we generate these coefficients in thin layers and, after the averaging process, we obtain an equivalent transversely isotropic medium. We repeat that process numerous times to get many equivalent transversely isotropic media, and – for each of them  – we compute their anisotropy parameters. We illustrate ϕ, ε and δ in the form of cross-plots that are relevant to the chosen variations of λ and μ. Additionally, we present a table with the computed ranges of anisotropy parameters that correspond to different variations of Lamé coefficients.     

钢板阻尼墙基于Bouc-Wen模型的模拟和工程应用研究

钢板阻尼墙抗震性能好,设计自由度大,在工程中使用越来越多。为了更高精度地模拟钢板阻尼墙的力学性能,提出基于Bouc-Wen模型的模拟分析方法。给出一种基于Bouc-Wen模型钢板阻尼墙参数辨识方法,根据试验数据和仿真结果的验证对比,文中提出的模拟方法可以较好地反映钢板阻尼墙的的力学性能曲线。该模拟方法可得到较好的仿真效果,对实际工程设计分析钢板阻尼墙对结构的抗震性能影响有高精度、数据化的帮助,可为进一步工程应用和设计分析提供参考。