Post-stack impedance blocky inversion based on analytic solution of viscous acoustic wave equation

The conventional impedance inversion method ignores the attenuation effect, transmission loss and inter-layer multiple waves; the smooth-like regularization approach makes the corresponding impedance solution excessively smooth. Both fundamentally limit the resolution of impedance result and lead to the inadequate ability of boundary characterization. Therefore, a post-stack impedance blocky inversion method based on the analytic solution of viscous acoustic equation is proposed. Based on the derived recursive formula of reflections, the 1D viscous acoustic wave equation is solved analytically to obtain zero-offset full-wave field response. Applying chain rule, the analytical expression of the Fréchet derivative is derived for gradient-descent non-linear inversion. Combined with smooth constraints, the blocky constraints can be introduced into the Bayesian inference framework to obtain stable and well-defined inversion results. According to the above theory, we firstly use model data to analyse the influence of incompleteness of forward method on seismic response, and further verify the effectiveness of the proposed method. Then the Q-value sensitivity analysis of seismic trace is carried out to reduce the difficulty of Q-value estimation. Finally, the real data from Lower Congo Basin in West Africa indicate that the proposed approach provide the high-resolution and well-defined impedance result. As a supplement and development of linear impedance inversion method, the non-linear viscous inversion could recover more realistic and reliable impedance profiles.     

用遗传算法优化估计Van Genuchten方程参数

描述土壤水分特征曲线的Van Genuchten方程是一个非线性复杂方程,采用传统的方法对方程参数进行估计往往因为计算复杂和人为因素的影响,而使估计结果带有较大的误差。遗传算法是当前处理一般非线性数学模型优化的一种新的优秀方法。为此,采用实数编码多子种群遗传算法对方程的参数进行优化估计。算例结果表明用遗传算法优化估计Van Genuchten方程参数具有求解速度快、计算精度和自动化程度高、人为干扰因素小、通用性强等优点。因此,遗传算法可作为计算Van Genuchten方程参数的一种新方法。     

HDD扩孔钻具组合三维小挠度运动规律浅析

以HDD铺设管道直线段轨迹中的扩孔钻具组合为研究对象,将扩孔钻具简化为质量均匀分布的转子模型,对扩孔钻具进行了详细的受力分析。以转子动力学理论为基础,通过对扩孔钻具所受力及力矩表示并建立了力和力矩平衡方程,采用矢量分析方法推导出了在小挠度情况下扩孔钻具的三维控制微分方程。     

多模式热动力耦合诊断系统与暴雨预报

应用包含热动力耦合强迫的新型散度演化方程,研发了该方程的诊断计算方案和相应的软件系统,以多个数值预报模式产品为资料源,计算得到关于强降水激发、维持的热动力耦合诊断场,形成了一个多时效的诊断场序列,它包含了正压非平衡强迫和热动力耦合强迫对散度演化的作用,反映了暴雨激发和持续的数值特征,为探索暴雨激发维持机理和寻求暴雨预报指标提供了数据支持。多模式诊断场的综合分析是集合预报方法的一种实践,有利于提高预报准确率。我们对2009年的两次暴雨个例进行了诊断分析,结果表明：T213、T639模式的诊断场能够较为准确的预报出暴雨的开始和结束,有明显的预报指示意义。     

用于势场反演的特殊解法

本文以无旋场的积分路径为射线,以边界上已知的势值为输入、输出量,利用反拉东变换求一个梯度场的势函数,用反演的方法求解可以化为拉普拉斯方程的一类二阶偏微分方程。其思路新颖,原理独特,运算正确,边界条件处理方便,数值实现容易。在重力场和地电场的研究中可能得到应用。     

Equation of state of neutron star cores and spin down of isolated pulsars

We study possible impact of a softening of the equation of state by a phase transition, or appearance of hyperons, on the spin evolution of isolated pulsars. Numerical simulations are performed using exact 2-D simulations in general relativity. The equation of state of dense matter at supranuclear densities is poorly known. Therefore, the accent is put on the general correlations between evolution and equation of state, and mathematical strictness. General conjectures referring to the structure of the one-parameter families of stationary configurations are formulated. The interplay of the back bending phenomenon and stability with respect to axisymmetric perturbations is described. Changes of pulsar parameters in a corequake following instability are discussed, for a broad choice of phase transitions predicted by different theories of dense matter. The energy release in a corequake, at a given initial pressure, is shown to be independent of the angular momentum of collapsing configuration. This result holds for various types of phases transition, with and without metastability. We critically review observations of pulsars that could be relevant for the detection of the signatures of the phase transition in neutron star cores. This work was partially supported by the Polish MNiI Grant no. 1P03D-008-27.     

Development of a robust runoff-prediction model by fusing the Rational Equation and a modified SCS-CN method

Abstract

The objective of this study is to develop a Modified Rational Equation (MoRE) that combines the advantages of the Rational Equation (e.g. simplicity and global acceptance) and those of the standard US Department of Agriculture (USDA) Soil Conservation Service (SCS) curve number (CN) method (e.g. easy parameterization and extensive verification across the world). Herein, the hypothesis is that the MoRE is more accurate, consistent and robust than the SCS-CN method and its improved versions in predicting runoff in watersheds with limited data. The MoRE was designed to have a simple structure that is described by four intrinsic parameters: CN, permanent wilting point, field capacity and saturation soil moisture, and does not include initial abstraction as a variable. An evaluation of 77 USDA small agricultural watersheds indicated that CN of the MoRE has different physical meanings from CN of the SCS-CN method. The MoRE (mean Nash-Sutcliffe coefficient, E > 0.73) performed better than the SCS-CN (mean E < 0.32) and the four improved models (mean E < 0.56) in reproducing the runoff of the study watersheds. Performance of all six models varied greatly between watersheds, as well as between events, but was independent of watershed drainage area. However, the model performances tend to be better for watersheds and/or events with a runoff-to-rainfall ratio of between 0.1 and 0.3 than for those with a ratio outside this range. The MoRE has the most consistent and robust performance.

Editor D. Koutsoyiannis; Associate editor I. Nalbantis

Citation Wang, X., Liu, T., and Yang, W., 2012. Development of a robust runoff-prediction model by fusing the rational equation and a modified SCS-CN method. Hydrological Sciences Journal, 57 (6), 1118–1140.     

Hydrostatic compression of galenobismutite (PbBi<Subscript>2</Subscript>S<Subscript>4</Subscript>): elastic properties and high-pressure crystal chemistry

A single-crystal sample of galenobismutite was subjected to hydrostatic pressures in the range of 0.0001 and 9 GPa at room temperature using the diamond-anvil cell technique. A series of X-ray diffraction intensities were collected at ten distinct pressures using a CCD equipped 4-circle diffractometer. The crystal structure was refined to R1(|F₀| > 4σ) values of approximately 0.05 at all pressures. By fitting a third-order Birch-Murnaghan equation of state to the unit-cell volumes V ₀ = 700.6(2) ų, K ₀ = 43.9(7) GPa and dK/dP = 6.9(3) could be determined for the lattice compression. Both types of cations in galenobismutite have stereochemically active lone electron pairs, which distort the cation polyhedra at room pressure. The cation eccentricities decrease at higher pressure but are still pronounced at 9 GPa. Galenobismutite is isotypic with CaFe₂O₄ (CF) but moves away from the idealised CF-type structure during compression. Instead of the two octahedral cation sites and one bi-capped trigonal-prismatic site, PbBi₂S₄ attains a new high-pressure structure characterised by one octahedral site and two mono-capped trigonal-prismatic sites. Analyses of the crystal structure at high pressure confirm the preference of Bi for the octahedral site and the smaller one of the two trigonal-prismatic sites.     

THE INVERSE PROBLEM OF KLEIN-GORDON EQUATION BOUNDARY VALUE PROBLEM AND ITS APPLICATION IN DATA ASSIMILATION

Since the solution of elliptic partial differential equations continuously depends on the boundary condition, the Euler equation derived from variational method cannot be solved without boundary condition. It is often difficult to provide the exact boundary condition in the practical use of variational method. However, in some application problems such as the remote sensing data assimilation, the values can be easily obtained in the inner region of the domain. In this paper, the boundary condition is tried to be retrieved by using part solutions in the inner area. Firstly, the variational problem of remote sensing data assimilation within a circular area is established. The Klein-Gordon elliptic equation is derived from the Euler method of variational problems with assumed boundary condition. Secondly, a computer-friendly Green function is constructed for the Dirichlet problem of two-dimensional Klein-Gordon equation, with the formal solution according to Green formula. Thirdly, boundary values are retrieved by solving the optimal problem which is constructed according to the best approximation between formal solutions and high-accuracy measurements in the interior of the domain. Finally, the assimilation problem is solved on substituting the retrieved boundary values into the Klein-Gordon equation. It is a type of inverse problem in mathematics. The advantage of this method lies in that it overcomes the inherent instability of the inverse problem of Fredholm integral equation and alleviates the error introduced by artificial boundary condition in data fusion using variational method in the past.     

The Equation of State for Caspian Sea Waters

The density ρ of Caspian Sea waters was measured as a function of temperature (273.15–343.15) K at conductivity salinities of 7.8 and 11.3 using the Anton-Paar Densitometer. Measurements were also made on one of the samples (S = 11.38) diluted with water as a function of temperature (T = 273.15–338.15 K) and salinity (2.5–11.3). These latter results have been used to develop an equation of state for the Caspian Sea (σ = ±0.007 kg m⁻³)