Three‐term inversion of prestack seismic data using a weighted l2, 1 mixed norm

We present a new inversion method to estimate, from prestack seismic data, blocky P‐ and S‐wave velocity and density images and the associated sparse reflectivity levels. The method uses the three‐term Aki and Richards approximation to linearise the seismic inversion problem. To this end, we adopt a weighted mixed l_{2, 1}‐norm that promotes structured forms of sparsity, thus leading to blocky solutions in time. In addition, our algorithm incorporates a covariance or scale matrix to simultaneously constrain P‐ and S‐wave velocities and density. This a priori information is obtained by nearby well‐log data. We also include a term containing a low‐frequency background model. The l_{2, 1} mixed norm leads to a convex objective function that can be minimised using proximal algorithms. In particular, we use the fast iterative shrinkage‐thresholding algorithm. A key advantage of this algorithm is that it only requires matrix–vector multiplications and no direct matrix inversion. The latter makes our algorithm numerically stable, easy to apply, and economical in terms of computational cost. Tests on synthetic and field data show that the proposed method, contrarily to conventional l₂‐ or l₁‐norm regularised solutions, is able to provide consistent blocky and/or sparse estimators of P‐ and S‐wave velocities and density from a noisy and limited number of observations.     

Moveout approximation for horizontal transversely isotropic and vertical transversely isotropic layered medium. Part I: 1D ray propagation‡

Anisotropy in subsurface geological models is primarily caused by two factors: sedimentation in shale/sand layers and fractures. The sedimentation factor is mainly modelled by vertical transverse isotropy (VTI), whereas the fractures are modelled by a horizontal transversely isotropic medium (HTI). In this paper we study hyperbolic and non‐hyperbolic normal reflection moveout for a package of HTI/VTI layers, considering arbitrary azimuthal orientation of the symmetry axis at each HTI layer. We consider a local 1D medium, whose properties change vertically, with flat interfaces between the layers. In this case, the horizontal slowness is preserved; thus, the azimuth of the phase velocity is the same for all layers of the package. In general, however, the azimuth of the ray velocity differs from the azimuth of the phase velocity. The ray azimuth depends on the layer properties and may be different for each layer. In this case, the use of the Dix equation requires projection of the moveout velocity of each layer on the phase plane. We derive an accurate equation for hyperbolic and high‐order terms of the normal moveout, relating the traveltime to the surface offset, or alternatively, to the subsurface reflection angle. We relate the azimuth of the surface offset to its magnitude (or to the reflection angle), considering short and long offsets. We compare the derived approximations with analytical ray tracing.     

Neuzeitlicher Schadstofftrend: Polycylische aromatische Kohlenwasserstoffe (PAK) in aquatischen und terrestrischen Sedimenten

A Recent Contamination Trend: Polycylic Aromatic Hydrocarbons (PAHs) in Aquatic and Terrestrial Sediments Concentrations of polycylic aromatic hydrocarbons (PAHs) measured in aquatic sediments of small, medium‐sized or very large inland lakes, of estuaries as well as in terrestrial sediments from Germany, the UK, Finland, and the USA were plotted in 10 trend curves over the period from 1800 to 1990. The segments of the bore cores were dated for their deposition age. For comparability of the results, the trend is preferentially shown by the guide parameter fluoranthene. Until the 19th century, a constant background value of C_Fluo = 10…50 μg/kg was observed in all regions. Then a steep increase in concentration followed, reaching in the aquatic sediments a maximum in the period 1940 to 1965, while the individual fluoranthene values ranged widely from 90 to 2400 μg/kg. In several regions, a more or less distinct decrease superseded this maximum, although this development did not appear in terrestrial sediments, in one estuary, and in remote forest lakes in Finland. The trend curves reflect the atmospheric deposition of PAHs in airborne dust and in some cases the import through runoff via rivers. Thus, the sediment profiles depict the development of air pollution by pyrogenic contaminants over more than a century. Sediment profiles from rivers can be used only with a high degree of reservation, because of the irregularity and low representativity of the sedimentation.     

前兆水位数据处理中的小波基及其参数选择方法研究

在分析不同小波基函数特征的基础上, 详细阐述了前兆水位数据处理中小波基及其参数的选择问题, 针对多个水位测项进行了实例分析和比较。 得到如下结论: ① 在检测前兆水位数据的突变点或不规则的突变部分, 选择db6、 bior2.6、 rbio2.6、 sym4和sym6小波基效果较好, 其中选用双正交小波bior2.6效果最佳; ② 在小波变换分解层数的选取方面: 第1～2层多适用于突变异常信号的获得, 第5～6层多适用于缓变异常信号的提取; 对于数据跳变异常显著的时间段, 第4～6层能较好地还原被高频信息掩盖的趋势变化过程。     

黏弹性双相介质参数反演的同伦方法

本文以介质的动态时域响应（位移、速度、加速度）作为反演的依据,对黏弹性双相介质的材料参数进行了反分析研究.根据附加条件,并考虑到由于测量误差和测量数据的不充分造成的反演问题的不适定性,将黏弹性双相介质的参数反演问题转化为非线性泛函的极小值问题,然后应用大范围收敛的同伦方法求出使泛函极小的根作为反问题的解.最后,以二维半空间黏弹性双相介质模型为例进行了数值反演分析,数值结果表明本文方法在参数的反演过程中目标函数稳定收敛,且具有一定的抗“噪声”干扰性能;与利用速度和加速度信息反演的结果相比,利用位移信息反演的结果具有较高的精度和抗噪声能力.     

The identification of conceptual hydrological models for surface water acidification

Ambiguity in parameter identification represents a potentially serious limitation to the application of models of surface water acidification. Previous work has concentrated on manipulation of two of the three factors affecting model identifiability, namely model structure and estimator properties. A new technique is proposed which uses different modes of response within the data to improve parameter identification. Preliminary results, obtained using the Birkenes model of surface water acidification, appear to show promise. The technique is robust in recovering model parameters from synthetic data, with and without error, and in assimilating problems of structural error.     

东海区带鱼年渔获量变动规律的探讨

本文根据东海区28a带鱼渔获量资料,采用参数模型法计算了东海区带鱼年渔获量变动规律,得出7a为一个变动周期的结论,井预测,若不采取有力措施保护资源,16a后资源将有枯竭的可能。     

高塔,高层建筑物振动与形变的超长摆综合观测

高大建筑物（如电视塔等）在外力作用和自身老化过程中出现的形变和振动频率特性变化的监测，是灾害预测和工程检验的重要内容。具有测站移动坐标系的超长摆观测系统，可以胜任这一超宽频（0～100HZ）振动和挠曲综合观测要求。文中导出超长摆观测方程xi＝－（R~2_iQ~2_s＋P~2_i－2R_iQ_sP_icosγ_s）~(1/2)Xmsin（ωt＋θ_i）并论证了在相应参数选择下，此系统具有足够的振幅放大系数（S_i≥0．8）和振动周期观测覆盖范围；对超低频一永久形变也具有真实的监测能力（S_i＝P_i）；具有监测频带宽、动态测量范围大、无漂移、抗干扰等优点。     

Assessment of radionuclide transport uncertainty in the unsaturated zone of Yucca Mountain

The present study assesses the uncertainty of flow and radionuclide transport in the unsaturated zone at Yucca Mountain using a Monte Carlo method. Matrix permeability, porosity, and sorption coefficient are considered random. Different from previous studies that assume distributions of the parameters, the distributions are determined in this study by applying comprehensive transformations and rigorous statistics to on-site measurements of the parameters. The distribution of permeability is further adjusted based on model calibration results. Correlation between matrix permeability and porosity is incorporated using the Latin Hypercube Sampling method. After conducting 200 Monte Carlo simulations of three-dimensional unsaturated flow and radionuclide transport for conservative and reactive tracers, the mean, variances, and 5th, 50th, and 95th percentiles for quantities of interest (e.g., matrix liquid saturation and water potential) are evaluated. The mean and 50th percentile are used as the mean predictions, and their associated predictive uncertainties are measured by the variances and the 5th and 95th percentiles (also known as uncertainty bounds). The mean predictions of matrix liquid saturation and water potential are in reasonable agreement with corresponding measurements. The uncertainty bounds include a large portion of the measurements, suggesting that the data variability can be partially explained by parameter uncertainty. The study illustrates propagation of predictive uncertainty of percolation flux, increasing downward from repository horizon to water table. Statistics from the breakthrough curves indicate that transport of the reactive tracer is delayed significantly by the sorption process, and prediction on the reactive tracer is of greater uncertainty than on the conservative tracer because randomness in the sorption coefficient increases the prediction uncertainty. Uncertainty in radionuclide transport is related to uncertainty in the percolation flux, suggesting that reducing the former entails reduction in the latter.     

1791年福建东山外海地震参数复核

本文依据校核后的历史资料,重新评定了1791年福建东山外海地震的影响场。再由影响场数据复核了该次地震的震中位置、震中烈度和震级值。得到的复核结果对修订已有的地震目录和合理确认该地区的发震构造将有所裨益。